A Comparative Genomic Study of Attenuated and Virulent Strains of Babesia bigemina

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.

Identification of novel immune correlates of protection against acute bovine babesiosis by superinfecting cattle with in vitro culture attenuated and virulent Babesia bovis strains

The apicomplexan tickborne parasites Babesia bovis and B. bigemina are the major causative agents of bovine babesiosis, a disease that negatively affects the cattle industry and food safety around the world. The absence of correlates of protection represents one major impediment for the development of effective and sustainable vaccines against bovine babesiosis. Herein we superinfected cattle with attenuated and virulent strains of B. bovis to investigate immune correlates of protection against acute bovine babesiosis. Three 6-month-old Holstein calves were infected intravenously (IV) with the in vitro culture attenuated Att-S74-T3Bo B. bovis strain (106 infected bovine red blood cells (iRBC)/calf) while three age-matched Holstein calves were inoculated IV with normal RBC as controls (106 RBC/calf). All Att-S74-T3Bo-infected calves showed a significant increase in temperature early after inoculation but recovered without treatment. Att-S74-T3Bo-infected calves also developed: (a) monocytosis, neutropenia, and CD4+ lymphopenia in peripheral blood on days 3 to 7 post-inoculation; (b) significant levels of TNFα, CXCL10, IFNγ, IL-4, and IL-10 in sera at day 6 after infection; and (c) IgM and IgG against B. bovis antigens, starting at days 10 and 30 post-inoculation, respectively. At 46 days post-Att-S74-T3Bo inoculation, all experimental calves were infected IV with the homologous virulent B. bovis strain Vir-S74-T3Bo (107 iRBC/calf). All Att-S74-T3Bo-infected calves survived superinfection with Vir-S74-T3Bo without displaying signs of acute babesiosis. In contrast, control animals showed signs of acute disease, starting at day 10 post-Vir-S74-T3Bo infection, and two of them were humanely euthanized at days 13 and 14 after inoculation due to the severity of their symptoms. Also, control calves showed higher (P<0.05) parasite load in peripheral blood compared to animals previously exposed to Att-S74-T3Bo. No significant alterations in the profile of leukocytes and cytokines were observed in Att-S74-T3Bo-inoculated after Vir-S74-T3Bo infection. In conclusion, data demonstrate novel changes in the profile of blood immune cells and cytokine expression in peripheral blood that are associated with protection against acute bovine babesiosis. These identified immune correlates of protection may be useful for designing effective and sustainable vaccines against babesiosis in cattle.

A novel promising diagnostic candidate selected by screening the transcriptome of Babesia gibsoni (Wuhan isolate) asexual stages in infected beagles

BACKGROUND

Babesia gibsoni is one of the causative agents of canine babesiosis worldwide. Some dogs infected with B. gibsoni show severe clinical signs with progressive anemia, hemoglobinuria and splenomegaly. However, most infected dogs present a state of chronic infection and thereby may be a persistent pathogen carrier, increasing the risk of pathogen spreading. To date, little is known about this pathogen, with genomic and transcriptomic data in particular generally unavailable. This lack of knowledge extensively limits the development of effective diagnostic strategies and vaccines.

METHODS

High-throughput RNA sequencing of total RNA of B. gibsoni asexual stages collected from infected beagles was performed. The unigenes were annotated in seven databases. The genes were sorted according to their fragments per kilobase per million (FPKM) value, which was used as an indicator for expression level. The gene with the highest FPKM value was cloned from the genome of B. gibsoni and further tested for immunogenicity, cellular localization and efficacy as a potential diagnostic candidate for detecting B. gibsoni in sera collected from beagles.

RESULTS

A total of 62,580,653 clean reads were screened from the 64,336,475 raw reads, and the corresponding 70,134 transcripts and 36,587 unigenes were obtained. The gene with the highest FPKM value was screened from the unigenes; its full length was 1276 bp, and it was named BgP30. The BgP30 gene comprised three exons and two introns, with a 786-bp open reading frame, and encoded 261 amino acids with a predicted molecular weight of 30 kDa. The cellular localization assay confirmed the existence of P30 protein in B. gibsoni parasites. Moreover, P30 was detected in the serum of experimentally B. gibsoni-infected beagles, from 15 days up to 422 days post-infection, suggesting its usefulness as a diagnostic candidate for both acute and chronic infections.

CONCLUSIONS

We sequenced the transcriptome of B. gibsoni asexual stages for the first time. The BgP30 gene was highly expressed in the transcriptome screening experiments, with further studies demonstrating that it could induce immune response in B. gibsoni-infected dogs. These results lead us to suggest that bgP30 may be a good diagnostic candidate marker to detect both acute and chronic B. gibsoni infections.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

Mexican Strains of Anaplasma marginale: A First Comparative Genomics and Phylogeographic Analysis

The One Health approach looks after animal welfare and demands constant monitoring of the strains that circulate globally to prevent outbreaks. Anaplasma marginale is the etiologic agent of bovine anaplasmosis and is endemic worldwide. This study aimed to analyze, for the first time, the genetic diversity of seven Mexican strains of A. marginale and their relationship with other strains reported. The main features of A. marginale were obtained by characterizing all 24 genomes reported so far. Genetic diversity and phylogeography were analyzed by characterizing the msp1a gene and 5′-UTR microsatellite sequences and constructing a phylogenetic tree with 540 concatenated genes of the core genome. The Mexican strains show 15 different repeat sequences in six MSP1a structures and have phylogeographic relationships with strains from North America, South America, and Asia, which confirms they are highly variable. Based on our results, we encourage the performance of genome sequencing of A. marginale strains to obtain a high assembly level of molecular markers and the performance of extensive phylogeographic analysis. Undoubtedly, genomic surveillance helps build a picture of how a pathogen changes and evolves in geographical regions. However, we cannot discard the study of relationships pathogens establish with ticks and how they have co-evolved to establish themselves as a successful transmission system.