Comparison of high throughput RNA sequences between Babesia bigemina and Babesia bovis revealed consistent differential gene expression that is required for the Babesia life cycle in the vertebrate and invertebrate hosts.

Bovine babesiosis caused by Babesia bigemina and Babesia bovis is an economically important disease that affects cattle worldwide. Both B. bigemina and B. bovis are transovarially transmitted by Rhipicephalus ticks. However, little is known regarding parasite gene expression during infection of the tick vector or mammalian host, which has limited the development of effective control strategies to alleviate the losses to the cattle industry. To understand Babesia gene regulation during tick and mammalian host infection, we performed high throughput RNA-sequencing using samples collected from calves and Rhipicephalus microplus ticks infected with B. bigemina. We evaluated gene expression between B. bigemina blood-stages and kinetes and compared them with previous B. bovis RNA-seq data. The results revealed similar patterns of gene regulation between these two tick-borne transovarially transmitted Babesia parasites. Like B. bovis, the transcription of several B. bigemina genes in kinetes exceeded a 1,000-fold change while a few of these genes had a >20,000-fold increase. To identify genes that may have important roles in B. bigemina and B. bovis transovarial transmission, we searched for genes upregulated in B. bigemina kinetes in the genomic datasets of B. bovis and non-transovarially transmitted parasites, Theileria spp. and Babesia microti. Using this approach, we identify genes that may be potential markers for transovarial transmission by B. bigemina and B. bovis. The findings presented herein demonstrate common Babesia genes linked to infection of the vector or mammalian host and may contribute to elucidating strategies used by the parasite to complete their life cycle.

Comparative analysis of gene expression between Babesia bovis blood stages and kinetes allowed by improved genome annotation.

Throughout their life cycle, Babesia parasites alternate between a mammalian host, where they cause babesiosis, and the tick vector. Transition between hosts results in distinct environmental signals that influence patterns of gene expression, consistent with the morphological and functional changes operating in the parasites during their life stages. In addition, comparing differential patterns of gene expression among mammalian and tick parasite stages can provide clues for developing improved methods of control. Hereby, we upgraded the genome assembly of Babesia bovis, a bovine hemoparasite, closing a 139 kbp gap, and used RNA-Seq datasets derived from mammalian blood and tick kinete stages to update the genome annotation. Of the originally annotated genes, 1,254 required structural changes, and 326 new genes were identified, leading to a different predicted proteome compared to the original annotation. Next, the RNA-Seq data was used to identify B. bovis genes that were differentially expressed in the vertebrate and arthropod hosts. In blood stages, 28% of the genes were upregulated up to 300 fold, whereas 26% of the genes in kinetes, a tick stage, were upregulated up to >19,000 fold. We thus discovered differentially expressed genes that may play key biological roles and serve as suitable targets for the development of vaccines to control bovine babesiosis.

Transcriptome dataset of Babesia bovis life stages within vertebrate and invertebrate hosts.

Babesia bovis is a hemoprotozoan parasite of cattle that has a complex life cycle within vertebrate and invertebrate hosts. In the mammalian host, B. bovis undergoes asexual reproduction while in the tick midgut, gametes are induced, fuse, and form zygotes. The zygote infects tick gut epithelial cells and transform into kinetes that are released into the hemolymph and invade other tick tissues such as the ovaries, resulting in transovarial transmission to tick offspring. To compare gene regulation between different B. bovis life stages, we collected parasites infecting bovine erythrocytes and tick hemolymph. Total RNA samples were isolated, and multiplexed libraries sequenced using paired-end 100 cycle reads of a HiSeq 2500. The data was normalized using the TMM method and analysed for significant differential expression using the generalized linear model likelihood ratio test (GLM LRT) in edgeR. To validate our datasets, ten genes were selected using NormFinder. Genes that had no significant fold change between the blood and tick stages in the RNA-Seq datasets were tested by quantitative PCR to determine their suitability as "housekeeping" genes. The normalized RNA-Seq data revealed genes upregulated during infection of the mammalian host or tick vector and six upregulated genes were validated by quantitative PCR. These datasets can help identify useful targets for controlling bovine babesiosis.