Genomic resources for a unique, low-virulence Babesia taxon from China

Chromosome-level genome assembly of Babesia caballi reveals diversity of multigene families among Babesia species

BACKGROUND

Babesia caballi is an intraerythrocytic parasite from the phylum Apicomplexa, capable of infecting equids and causing equine piroplasmosis. However, since there is limited genome information available on B. caballi, molecular mechanisms involved in host specificity and pathogenicity of this species have not been fully elucidated yet.

RESULTS

Genomic DNA from a B. caballi subclone was purified and sequenced using both Illumina and Nanopore technologies. The resulting assembled sequence consisted of nine contigs with a size of 12.9 Mbp, rendering a total of 5,910 protein-coding genes. The phylogenetic tree of Apicomplexan species was reconstructed using 263 orthologous genes. We identified 481 ves1-like genes and named "ves1c". In contrast, expansion of the major facilitator superfamily (mfs) observed in closely related B. bigemina and B. ovata species was not found in B. caballi. A set of repetitive units containing an open reading frame with a size of 297 bp was also identified.

CONCLUSIONS

We present a chromosome-level genome assembly of B. caballi. Our genomic data may contribute to estimating gene expansion events involving multigene families and exploring the evolution of species from this genus.

Whole genome sequence and diversity in multigene families of Babesia ovis

Ovine babesiosis, caused by Babesia ovis, is an acute, lethal, and endemic disease worldwide and causes a huge economic loss to animal industry. Pathogen genome sequences can be utilized for selecting diagnostic markers, drug targets, and antigens for vaccine development; however, those for B. ovis have not been available so far. In this study, we obtained a draft genome sequence for B. ovis isolated from an infected sheep in Turkey. The genome size was 7.81 Mbp with 3,419 protein-coding genes. It consisted of 41 contigs, and the N50 was 526 Kbp. There were 259 orthologs identified among eight Babesia spp., Plasmodium falciparum, and Toxoplasma gondii. A phylogeny was estimated on the basis of the orthologs, which showed B. ovis to be closest to B. bovis. There were 43 ves genes identified using hmm model as well. They formed a discriminating cluster to other ves multigene family of Babesia spp. but showed certain similarities to those of B. bovis, B. caballi, and Babesia sp. Xinjiang, which is consistent with the phylogeny. Comparative genomics among B. ovis and B. bovis elucidated uniquely evolved genes in these species, which may account for the adaptation.

Small Ruminant Piroplasmosis: High Prevalence of Babesia aktasi n. sp. in Goats in Türkiye

Small ruminant piroplasmosis is the hemoparasitic infection of sheep and goats caused by Babesia and Theileria species responsible for clinical infections with high mortality outcomes. The disease is transmitted by ixodid ticks and prevalent in the tropical and subtropical regions of the world, including Türkiye. A prevalence survey, using molecular methods, is conducted in this study to determine the frequency of newly defined Babesia aktasi n. sp. and other tick-borne piroplasm species in small ruminants in Turkiye. A total of 640 blood samples from sheep (n = 137) and goats (n = 503) were analyzed by nested PCR-based reverse line blot (RLB) hybridization. The results show that 32.3% (207/640) of apparently healthy, small ruminants are infected with three Theileria and two Babesia species. Babesia aktasi n. sp. was the most prevalent species in goats, with 22.5% of samples being positive, followed by B. ovis (4%), T. ovis (2.8%), T. annulata (2.6%), and Theileria sp. (0.6%). None of the sheep samples were positive for Babesia aktasi n. sp.; however, 51.8% were infected with T. ovis. In conclusion, the findings reveal that B. aktasi n. sp. is highly prevalent in goats, but absent in sheep. In future studies, experimental infections will determine whether B. aktasi n. sp. is infectious to sheep, as well as its pathogenicity in small ruminants.

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.

Phylogenetic analyses of the mitochondrial, plastid, and nuclear genes of Babesia sp. Mymensingh and its naming as Babesia naoakii n. sp

Background

The recently discovered Babesia sp. Mymensingh, which causes clinical bovine babesiosis, has a wide geographical distribution. We investigated the phylogenetic position of Babesia sp. Mymensingh using its mitochondrial, plastid, and nuclear genes. Based on morphological and molecular data, Babesia sp. Mymensingh is a unique species and we named it as Babesia naoakii n. sp.

Methods

A blood DNA sample from a Babesia sp. Mymensingh-infected cow was subjected to genome sequencing to obtain the sequences of mitochondrial, plastid, and nuclear genes. Six phylogenetic trees were then constructed with (1) concatenated amino acid sequences of cytochrome oxidase subunit I, cytochrome oxidase subunit III, and cytochrome b genes of the mitochondrial genome; (2) 16S rRNA of the plastid genome; (3) nucleotide sequences of the elongation factor Tu gene of the plastid genome; (4) ITS1-5.8S rRNA-ITS2; (5) concatenated nucleotide sequences of 89 nuclear genes; and (6) concatenated amino acid sequences translated from the 89 nuclear genes.

Results

In all six phylogenetic trees, B. naoakii n. sp. formed a sister clade to the common ancestor of Babesia bigemina and B. ovata. The concatenated nuclear genes of B. naoakii n. sp. and their translated amino acid sequences shared lower identity scores with the sequences from B. bigemina (82.7% and 84.7%, respectively) and B. ovata (83.5% and 85.5%, respectively) compared with the identity scores shared between the B. bigemina and B. ovata sequences (86.3% and 87.9%, respectively).

Conclusions

Our study showed that B. naoakii n. sp. occupies a unique phylogenetic position distinct from existing Babesia species. Our findings, together with morphological differences, identify B. naoakii n. sp. as a distinct parasite species.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05374-9.

Comparative genomic analysis of Babesia duncani responsible for human babesiosis

Background

Human babesiosis, caused by parasites of the genus Babesia, is an emerging and re-emerging tick-borne disease that is mainly transmitted by tick bites and infected blood transfusion. Babesia duncani has caused majority of human babesiosis in Canada; however, limited data are available to correlate its genomic information and biological features.

Results

We generated a B. duncani reference genome using Oxford Nanopore Technology (ONT) and Illumina sequencing technology and uncovered its biological features and phylogenetic relationship with other Apicomplexa parasites. Phylogenetic analyses revealed that B. duncani form a clade distinct from B. microti, Babesia spp. infective to bovine and ovine species, and Theileria spp. infective to bovines. We identified the largest species-specific gene family that could be applied as diagnostic markers for this pathogen. In addition, two gene families show signals of significant expansion and several genes that present signatures of positive selection in B. duncani, suggesting their possible roles in the capability of this parasite to infect humans or tick vectors.

Conclusions

Using ONT sequencing and Illumina sequencing technologies, we provide the first B. duncani reference genome and confirm that B. duncani forms a phylogenetically distinct clade from other Piroplasm parasites. Comparative genomic analyses show that two gene families are significantly expanded in B. duncani and may play important roles in host cell invasion and virulence of B. duncani. Our study provides basic information for further exploring B. duncani features, such as host-parasite and tick-parasite interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01361-9.

Recent Advances in Molecular Genetic Tools for Babesia

Development of in vitro culture and completion of genome sequencing of several Babesia parasites promoted the efforts to establish transfection systems for these parasites to dissect the gene functions. It has been more than a decade since the establishment of first transfection for Babesia bovis, the causative agent of bovine babesiosis. However, the number of genes that were targeted by genetic tools in Babesia parasites is limited. This is partially due to the low efficiencies of these methods. The recent adaptation of CRISPR/Cas9 for genome editing of Babesia bovis can accelerate the efforts for dissecting this parasite’s genome and extend the knowledge on biological aspects of erythrocytic and tick stages of Babesia. Additionally, glmS ribozyme as a conditional knockdown system is available that could be used for the characterization of essential genes. The development of high throughput genetic tools is needed to dissect the function of multigene families, targeting several genes in a specific pathway, and finally genome-wide identification of essential genes to find novel drug targets. In this review, we summarized the current tools that are available for Babesia and the genes that are being targeted by these tools. This may draw a perspective for the future development of genetic tools and pave the way for the identification of novel drugs or vaccine targets.

Stable transfection system for Babesia sp. Xinjiang

BACKGROUND

Stable transfection systems have been described in many protozoan parasites, including Plasmodium falciparum, Cryptosporidium parvum, Babesia bovis, Babesia ovata, and Babesia gibsoni. For Babesia sp. Xinjiang (Bxj), which is the causative pathogen of ovine babesiosis and mainly prevails across China, the platform of those techniques remains absent. Genetic manipulation techniques are powerful tools to enhance our knowledge on parasite biology, which may provide potential drug targets and diagnostic markers.

METHODS

We evaluated the inhibition efficiency of blasticidin (BSD) and WR99210 to Bxj. Then, a plasmid was constructed bearing selectable marker BSD, green fluorescent protein (GFP) gene, and rhoptry-associated protein-1 3' terminator region (rap 3' TR). The plasmid was integrated into the elongation factor-1 alpha (ef-1α) site of Bxj genome by cross-over homologous recombination technique. Twenty μg of plasmid was transfected into Bxj merozoites. Subsequently, drug selection was performed 24 h after transfection to generate transfected parasites.

RESULTS

Transfected parasite lines, Bxj-c1, Bxj-c2, and Bxj-c3, were successfully obtained after transfection, drug selection, and colonization. Exogenous genes were integrated into the Bxj genome, which were confirmed by PCR amplification and sequencing. In addition, results of western blot (WB) and indirect immunofluorescence assay (IFA) revealed that GFP-BSD had expressed for 11 months.

CONCLUSIONS

In our present study, stable transfection system for Bxj was successfully developed. We anticipate that this platform will greatly facilitate basic research of Bxj.

Systematic Comparison of the Performances of De Novo Genome Assemblers for Oxford Nanopore Technology Reads From Piroplasm

Background

Emerging long reads sequencing technology has greatly changed the landscape of whole-genome sequencing, enabling scientists to contribute to decoding the genetic information of non-model species. The sequences generated by PacBio or Oxford Nanopore Technology (ONT) be assembled de novo before further analyses. Some genome de novo assemblers have been developed to assemble long reads generated by ONT. The performance of these assemblers has not been completely investigated. However, genome assembly is still a challenging task.

Methods and Results

We systematically evaluated the performance of nine de novo assemblers for ONT on different coverage depth datasets. Several metrics were measured to determine the performance of these tools, including N50 length, sequence coverage, runtime, easy operation, accuracy of genome and genomic completeness in varying depths of coverage. Based on the results of our assessments, the performances of these tools are summarized as follows: 1) Coverage depth has a significant effect on genome quality; 2) The level of contiguity of the assembled genome varies dramatically among different de novo tools; 3) The correctness of an assembled genome is closely related to the completeness of the genome. More than 30× nanopore data can be assembled into a relatively complete genome, the quality of which is highly dependent on the polishing using next generation sequencing data.

Conclusion

Considering the results of our investigation, the advantage and disadvantage of each tool are summarized and guidelines of selecting assembly tools are provided under specific conditions.

N-Glycosylation in Piroplasmids: Diversity within Simplicity

N-glycosylation has remained mostly unexplored in Piroplasmida, an order of tick-transmitted pathogens of veterinary and medical relevance. Analysis of 11 piroplasmid genomes revealed three distinct scenarios regarding N-glycosylation: Babesia sensu stricto (s.s.) species add one or two N-acetylglucosamine (NAcGlc) molecules to proteins; Theileria equi and Cytauxzoon felis add (NAcGlc)2-mannose, while B. microti and Theileria s.s. synthesize dolichol-P-P-NAcGlc and dolichol-P-P-(NAcGlc)2 without subsequent transfer to proteins. All piroplasmids possess the gene complement needed for the synthesis of the N-glycosylation substrates, dolichol-P and sugar nucleotides. The oligosaccharyl transferase of Babesia species, T. equi and C. felis, is predicted to be composed of only two subunits, STT3 and Ost1. Occurrence of short N-glycans in B. bovis merozoites was experimentally demonstrated by fluorescence microscopy using a NAcGlc-specific lectin. In vitro growth of B. bovis was significantly impaired by tunicamycin, an inhibitor of N-glycosylation, indicating a relevant role for N-glycosylation in this pathogen. Finally, genes coding for N-glycosylation enzymes and substrate biosynthesis are transcribed in B. bovis blood and tick stages, suggesting that this pathway is biologically relevant throughout the parasite life cycle. Elucidation of the role/s exerted by N-glycans will increase our understanding of these successful parasites, for which improved control measures are needed.

Whole genome sequencing of Theileria parva using target capture

Protozoan parasite isolation and purification are laborious and time-consuming processes required for high quality genomic DNA used in whole genome sequencing. The objective of this study was to capture whole Theileria parva genomes directly from cell cultures and blood samples using RNA baits. Cell culture material was bait captured or sequenced directly, while blood samples were all captured. Baits had variable success in capturing T. parva genomes from blood samples but were successful in cell cultures. Genome mapping uncovered extensive host contamination in blood samples compared to cell cultures. Captured cell cultures had over 81 fold coverage for the reference genome compared to 0-33 fold for blood samples. Results indicate that baits are specific to T. parva, are a good alternative to conventional methods and thus ideal for genomic studies. This study also reports the first whole genome sequencing of South African T. parva.

Novel Babesia bovis exported proteins that modify properties of infected red blood cells

Babesia bovis causes a pathogenic form of babesiosis in cattle. Following invasion of red blood cells (RBCs) the parasite extensively modifies host cell structural and mechanical properties via the export of numerous proteins. Despite their crucial role in virulence and pathogenesis, such proteins have not been comprehensively characterized in B. bovis. Here we describe the surface biotinylation of infected RBCs (iRBCs), followed by proteomic analysis. We describe a multigene family (mtm) that encodes predicted multi-transmembrane integral membrane proteins which are exported and expressed on the surface of iRBCs. One mtm gene was downregulated in blasticidin-S (BS) resistant parasites, suggesting an association with BS uptake. Induced knockdown of a novel exported protein encoded by BBOV_III004280, named VESA export-associated protein (BbVEAP), resulted in a decreased growth rate, reduced RBC surface ridge numbers, mis-localized VESA1, and abrogated cytoadhesion to endothelial cells, suggesting that BbVEAP is a novel virulence factor for B. bovis.

A Historical Overview of Research on Babesia orientalis, a Protozoan Parasite Infecting Water Buffalo

Babesiosis is a globally important zoonotic disease caused by tick-borne intraerythrocytic protozoan of the genus Babesia (phylum apicomplexa). In China, there are five species that infect cattle buffalo and cause great economic loss, which include Babesia bigemina, B. bovis, B. major, B. ovata, and B. orientalis. Among them, B. orientalis is the most recently identified new Babesia species epidemic in China. This review summarized the work done in the past 33 years to give an overview of what learned about this parasite. This parasitic protozoan was found in 1984 in Central and South China and then named as B. orientalis in 1997 based on its differences in transmitting host, morphology, pathogenicity and characteristics of in vitro cultivation when compared with B. bigemina and B. bovis. It was found that Rhipicephalus haemaphysaloides is the transmitting vector and water buffalo is the only reported host. Phylogenetic analysis based on the 18S rRNA gene also confirmed that B. orientalis is a new species. After species verification, four diagnostic methods including semi-nest PCR, loop-mediated isothermal amplification assay, reverse line blot hybridization assay, and real-time PCR were established for lab and field use purposes. Genomic sequencing was conducted and the complete genomes of mitochondria and apicoplast were annotated. Future work will be focused on developing effective vaccines, identifying drug targets and screening useful drugs for controlling B. orientalis in water buffalo.