Bovine babesiosis, caused by the protozoan Babesia bigemina, is one of the most important hemoparasite diseases of cattle in Mexico and the world. An attenuated B. bigemina strain maintained under in vitro culture conditions has been used as a live attenuated vaccine; however, the biological mechanisms involved in attenuation are unknown. The objective of this study was to identify, through a comparative transcriptomics approach, the components of the B. bigemina virulent parasites that are differentially expressed in vivo, as opposed to those expressed by B. bigemina attenuated vaccine parasites when inoculated into naïve cattle. The biological material under study was obtained by inoculating spleen-intact cattle with infected erythrocytes containing either the attenuated strain or a virulent field strain. After RNA extraction, transcriptomic analysis (RNA-seq) was performed, followed by bioinformatic Differential Expression (DE) analysis and Gene Ontology (GO) term enrichment. The high-throughput sequencing results obtained by analyzing three biological replicates for each parasite strain ranged from 9,504,000 to 9,656,000, and 13,400,000 to 15,750,000 reads for the B. bigemina attenuated and virulent strains, respectively. At least 519 differentially expressed genes were identified in the analyzed strains. In addition, GO analysis revealed both similarities and differences across the three categories: cellular components, biological processes, and molecular functions. The attenuated strain of B. bigemina derived from in vitro culture presents global transcriptomic changes when compared to the virulent strain. Moreover, the obtained data provide insights into the potential molecular mechanisms associated with the attenuation or pathogenicity of each analyzed strain, offering molecular markers that might be associated with virulence or potential vaccine candidates.
Keywords: Babesia bigemina; differentially expressed genes; transcriptome profiling.