Up-regulated expression of spherical body protein 2 truncated copy 11 in Babesia bovis is associated with reduced cytoadhesion to vascular endothelial cells

Pushing the frontiers of babesiosis research: in vitro culture and gene editing

Babesiosis is a tick-borne parasitic disease that poses a significant risk to both animal and human health. A comprehensive understanding of Babesia biology necessitates the application of advanced laboratory techniques. This review explores recent advancements in gene editing technologies of Babesia, emphasizing the foundational importance of in vitro culture systems. We highlight the historical challenges encountered in establishing effective in vitro culture and discuss the need for optimizing these methods to enhance gene editing efficiency. Here, we describe recent progress in Babesia transfection, different gene manipulation systems, and the applications of gene editing. This review aims to provide essential insights and technical guidance for future studies in Babesia genetics, highlighting the transformative potential of gene manipulation in combating this important parasitic disease.

ves1α genes expression is the major determinant of Babesia bovis-infected erythrocytes cytoadhesion to endothelial cells

Babesia bovis causes the most pathogenic form of babesiosis in cattle, resulting in high mortality in naive adults. This parasite invades red blood cells (RBCs) within the bovine hosts where they multiply and produce clinical disease. Babesia bovis exports numerous proteins into invaded RBCs changing its properties. Thus, the infected RBCs (iRBCs) are capable to cytoadhere in the microvasculature of internal organs and brain, leading to respiratory distress, neurologic signs, and mortality. Variant Erythrocyte Surface Antigen 1 (VESA1) is one of those exported proteins by B. bovis which represents a major virulence factor due to its central role in immune evasion by antigenic variation and intravascular parasite sequestration. VESA1 is a heterodimer protein encoded by ves1α and ves1β multigene family and localized on the ridges, the focal point for cytoadhesion. To gain further insights into the molecular mechanisms of cytoadhesion of B. bovis, we panned the parasites with bovine brain microvasculature endothelial cells, which resulted in obtaining several clones with different cytoadherence abilities. The transcriptome analysis of 2 high and 2 low cytoadherent clones revealed that ves1α sequences were diversified, likely resulting from genomic recombination. On the other hand, ves1β sequences were almost identical among these 4 clones. Insertion and expression of ves1α of a clone with high binding into ef-1α locus of a low binding clone increased cytoadherence confirming the role of ves1α suggested by our transcriptome data. Whole genome sequencing of cytoadherent clones revealed active locus of ves1 on chromosome 2. These results suggest that VESA1a proteins encoded by ves1α genes determine the cytoadherence strength of B. bovis and they are in the active site for recombination.

Critical role of Babesia bovis spherical body protein 3 in ridge formation on infected red blood cells

Babesia bovis, an apicomplexan intraerythrocytic protozoan parasite, causes serious economic loss to cattle industries around the world. Infection with this parasite leads to accumulation of infected red blood cells (iRBCs) in the brain microvasculature that results in severe clinical complications known as cerebral babesiosis. Throughout its growth within iRBCs, the parasite exports various proteins to the iRBCs that lead to the formation of protrusions known as "ridges" on the surface of iRBCs, which serve as sites for cytoadhesion to endothelial cells. Spherical body proteins (SBPs; proteins secreted from spherical bodies, which are organelles specific to Piroplasmida) are exported into iRBCs, and four proteins (SBP1-4) have been reported to date. In this study, we elucidated the function of SBP3 using an inducible gene knockdown (KD) system. Localization of SBP3 was assessed by immunofluorescence assay, and only partial colocalization was detected between SBP3 and SBP4 inside the iRBCs. In contrast, colocalization was observed with VESA-1, which is a major parasite ligand responsible for the cytoadhesion. Immunoelectron microscopy confirmed localization of SBP3 at the ridges. SBP3 KD was performed using the glmS system, and effective KD was confirmed by Western blotting, immunofluorescence assay, and RNA-seq analysis. The SBP3 KD parasites showed severe growth defect suggesting its importance for parasite survival in the iRBCs. VESA-1 on the surface of iRBCs was scarcely detected in SBP3 KD parasites, whereas SBP4 was still detected in the iRBCs. Moreover, abolition of ridges on the iRBCs and reduction of iRBCs cytoadhesion to the bovine brain endothelial cells were observed in SBP3 KD parasites. Immunoprecipitation followed by mass spectrometry analysis detected the host Band 3 multiprotein complex, suggesting an association of SBP3 with iRBC cytoskeleton proteins. Taken together, this study revealed the vital role of SBP3 in ridge formation and its significance in the pathogenesis of cerebral babesiosis.

Babesia gibsoni Whole-Genome Sequencing, Assembling, Annotation, and Comparative Analysis

The intracellular protozoan parasite Babesia gibsoni infects canine erythrocytes and causes babesiosis. The hazards to animal health have increased due to the rise of B. gibsoni infections and medication resistance. However, the lack of high-quality full-genome sequencing sets has expanded the obstacles to the development of pathogeneses, drugs, and vaccines. In this study, the whole genome of B. gibsoni was sequenced, assembled, and annotated. The genomic size of B. gibsoni was 7.94 Mbp in total. Four chromosomes with the size of 0.69 Mb, 2.10 Mb, 2.77 Mb, and 2.38 Mb, respectively, 1 apicoplast (28.4 Kb), and 1 mitochondrion (5.9 Kb) were confirmed. KEGG analysis revealed 2,641 putative proteins enriched on 316 pathways, and GO analysis showed 7,571 annotations of the nuclear genome in total. Synteny analysis showed a high correlation between B. gibsoni and B. bovis. A new divergent point of B. gibsoni occurred around 297.7 million years ago, which was earlier than that of B. bovis, B. ovata, and B. bigemina. Orthology analysis revealed 22 and 32 unique genes compared to several Babesia spp. and apicomplexan species. The metabolic pathways of B.gibsoni were characterized, pointing to a minimal size of the genome. A species-specific secretory protein SA1 and 19 homologous genes were identified. Selected specific proteins, including apetala 2 (AP2) factor, invasion-related proteins BgAMA-1 and BgRON2, and rhoptry function proteins BgWH_04g00700 were predicted, visualized, and modeled. Overall, whole-genome sequencing provided molecular-level support for the diagnosis, prevention, clinical treatment, and further research of B. gibsoni. IMPORTANCE The whole genome of B. gibsoni was first sequenced, annotated, and disclosed. The key part of genome composition, four chromosomes, was comparatively analyzed for the first time. A full-scale phylogeny evolution analysis based on the whole-genome-wide data of B. gibsoni was performed, and a new divergent point on the evolutionary path was revealed. In previous reports, molecular studies were often limited by incomplete genomic data, especially in key areas like life cycle regulation, metabolism, and host-pathogen interaction. With the whole-genome sequencing of B. gibsoni, we provide useful genetic data to encourage the exploration of new terrain and make it feasible to resolve the theoretical and practical problems of babesiosis.

Comparative genomics reveals unique features of two Babesia motasi subspecies: Babesia motasi lintanensis and Babesia motasi hebeiensis

Parasites of the Babesia genus are prevalent worldwide and infect a wide diversity of domestic animals and humans. Herein, using Oxford Nanopore Technology and Illumina sequencing technologies, we sequenced two Babesia sub-species, Babesia motasi lintanensis and Babesia motasi hebeiensis. We identified 3,815 one-to-one ortholog genes that are specific to ovine Babesia spp. Phylogenetic analysis reveals that the two B. motasi subspecies form a distinct clade from other Piroplasma spp. Consistent with their phylogenetic position, comparative genomic analysis reveals that these two ovine Babesia spp. share higher colinearity with Babesia bovis than with Babesia microti. Concerning the speciation date, B. m. lintanensis split from B. m. hebeiensis approximately 17 million years ago. Genes correlated to transcription, translation, protein modification and degradation, as well as differential/specialized gene family expansions in these two subspecies may favor adaptation to vertebrate and tick hosts. The close relationship between B. m. lintanensis and B. m. hebeiensis is underlined by a high degree of genomic synteny. Compositions of most invasion, virulence, development, and gene transcript regulation-related multigene families, including spherical body protein, variant erythrocyte surface antigen, glycosylphosphatidylinositol anchored proteins, and transcription factor Apetala 2 genes, is largely conserved, but in contrast to this conserved situation, we observe major differences in species-specific genes that may be involved in multiple functions in parasite biology. For the first time in Babesia spp., we find abundant fragments of long terminal repeat-retrotransposons in these two species. We provide fundamental information to characterize the genomes of B. m. lintanensis and B. m. hebeiensis, providing insights into the evolution of B. motasi group parasites.

Molecular characterization and genetic diversity of Babesia bovis and Babesia bigemina of cattle in Thailand

Babesia bovis and B. bigemina are the most common tick-borne parasites that cause bovine babesiosis which effects livestock production, leading to economic losses in tropical and subtropical areas of the world. The aims of this study were to determine the molecular detection, genetic diversity and antigenicity prediction of B. bovis based on spherical body protein 2 (sbp-2) gene and B. bigemina based on rhoptry-associated protein 1a (rap-1a) gene in cattle in Thailand. By PCR assay, the molecular detection of B. bovis and B. bigemina infection revealed levels of 2.58% (4/155) and 5.80% (9/155), respectively. The phylograms showed that B. bovis sbp-2 and B. bigemina rap-1a sequences displayed 5 and 3 clades with similarity ranging between 85.53 to 100% and 98.28 to 100%, respectively, when compared within Thailand strain. Diversity analysis of sbp-2 and rap-1a sequences showed 18 and 4 haplotypes, respectively. The entropy analysis illustrated 104 and 7 polymorphic sites of sbp-2 and rap-1a nucleic acid sequences, respectively, while those of sbp-2 and rap-1a amino acid sequences showed 46 and 4 high entropy peaks, respectively. Motifs analysis exhibited the distribution and conservation among sbp-2 and rap-1a sequences. The continuous and discontinuous B-cell epitopes have also been evaluated in this work. Therefore, our findings may be used to ameliorate the understanding inputs of molecular phylogeny, genetic diversity and antigenicity of B. bovis and B. bigemina Thailand stains.

Advances in understanding red blood cell modifications by Babesia

Babesia are tick-borne protozoan parasites that can infect livestock, pets, wildlife animals, and humans. In the mammalian host, they invade and multiply within red blood cells (RBCs). To support their development as obligate intracellular parasites, Babesia export numerous proteins to modify the RBC during invasion and development. Such exported proteins are likely important for parasite survival and pathogenicity and thus represent candidate drug or vaccine targets. The availability of complete genome sequences and the establishment of transfection systems for several Babesia species have aided the identification and functional characterization of exported proteins. Here, we review exported Babesia proteins; discuss their functions in the context of immune evasion, cytoadhesion, and nutrient uptake; and highlight possible future topics for research and application in this field.

Establishment of a Transient and Stable Transfection System for Babesia duncani Using a Homologous Recombination Strategy

Genetic modification provides an invaluable molecular tool to dissect the biology and pathogenesis of pathogens. However, no report is available about the genetic modification of Babesia duncani, a pathogen responsible for human babesiosis that is widespread in North America, suggesting the necessity to develop a genetic manipulation method to improve the strategies for studying and understanding the biology of protozoan pathogens. The establishment of a genetic modification method requires promoters, selectable markers, and reporter genes. Here, the double-copy gene elongation factor-1α (ef-1α) and its promoters were amplified by conventional PCR and confirmed by sequencing. We established a transient transfection system by using the ef-1αB promoter and the reporter gene mCherry and achieved stable transfection through homologous recombination to integrate the selection marker hDHFR-eGFP into the parasite genome. The potential of this genetic modification method was tested by knocking out the thioredoxin peroxidase-1 (TPX-1) gene, and under the drug pressure of 5 nM WR99210, 96.3% of the parasites were observed to express green fluorescence protein (eGFP) by flow cytometry at day 7 post-transfection. Additionally, the clone line of the TPX-1 knockout parasite was successfully obtained by the limiting dilution method. This study provided a transfection method for B. duncani, which may facilitate gene function research and vaccine development of B. duncani.

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.

Identification of a novel variant erythrocyte surface antigen-1 (VESA1) in Babesia orientalis

Babesia orientalis, belonging to the phylum Apicomplexa, is mainly accountable for water buffalo babesiosis, which adversely affected the livestock industry in China. Variant erythrocyte surface antigen-1 (VESA1), an antigen that helps infected erythrocytes to escape from host immune responses, was first reported in Babesia bovis. Various VESA1 proteins have also been characterized in other Babesia species. Nevertheless, there is no research on the identification and characterization of VESA1 proteins in Babesia orientalis. In this study, the BoVESA1 gene was amplified from both gDNA and cDNA. The results revealed that it is an intronless gene with a full length of 753 bp, encoding a protein of 250 amino acids with a predicted molecular weight of 28 kDa. The coding sequence (CDS) was cloned into the pGEX-6p-1 vector using a homologous recombination kit and expressed as a glutathione-S-transferase (GST)-fusion protein with a molecular weight of 53 kDa. The tertiary structure of BoVESA1 was predicted using the I-TASSER software. The recombinant protein was subjected to western blotting; the immunogenicity of recombinant BoVESA1 (rBoVESA1) was identified by incubating it with B. orientalis-positive serum. The native BoVESA1 was identified using the lysates of B. orientalis-infected water buffalo erythrocytes incubated with the anti-rBoVESA1 mouse serum. The results showed a band of ~ 28 kDa, which is similar to the predicted size. Immunofluorescence assay (IFA) using anti-rBoVESA1 serum probed indicated a strong signal in the infected RBCs, while the negative control showed no signal. In conclusion, the VESA1 protein was first identified in B. orientalis. This study facilitated further investigation of B. orientalis, and the results indicated that BoVESA1 may serve as a potential candidate antigen for diagnosis and detection of B. orientalis infection.

A framework for signaling throughout the life cycle of Babesia species

Babesia species are tick-borne intracellular parasites that infect the red blood cells of their mammalian host, leading to severe or fatal disease. Babesia spp. infect a wide range of mammalian species and cause a significant economic burden globally, predominantly through disease in cattle. Several Babesia spp. are increasingly being recognized as zoonotic pathogens of humans. Babesia spp. have complex life cycles involving multiple stages in the tick and the mammalian host. The parasite utilizes complex signaling pathways during replication, egress, and invasion in each of these stages. They must also rapidly respond to their environment when switching between the mammalian and tick stages. This review will focus on the signaling pathways and environmental stimuli that Babesia spp. utilize in the bloodstream and for transmission to the tick, with an emphasis on the role of phosphorylation- and calcium-based signaling during egress and invasion. The expanding availability of in vitro and in vivo culture systems, genomes, transcriptomes, and transgenic systems available for a range of Babesia spp. should encourage further biological and translational studies of these ubiquitous parasites.

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.

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.

Establishment of a stable transfection method in Babesia microti and identification of a novel bidirectional promoter of Babesia microti

Babesia microti, an emerging human pathogen, is primarily transmitted through a bite of an infected tick and blood transfusions in human. Stable transfection technique has been reported in many protozoan parasites over the past few years. However, in vivo transient and stable transfection method has not been established for Babesia microti. Here, for the first time, we present a method of transient as well as stable transfection of the Babesia microti (B. microti) in the in vivo conditions. We have identified a novel promoter of B. microti. We also demonstrated that Plasmodium berghei DHFR promoter is recognized and functional in B. microti. We show that BM-CTQ41297 promoter control the expression of two genes, which are present on either side and thus represents a bi-functional promoter in B. microti. The predicted promoter activity values using Promoter 2.0 program is higher for BM- CTQ41297 promoter than strong promoters such as β-actin, ef-1β, and many other promoters. Furthermore, we discovered a non-essential locus for the genetic manipulation of the parasite, allowing us to stably integrate foreign genes; GFP, mCherry, into the B. microti. The transfection using an electroporation method and genetic manipulation of B. microti is now achievable and it is possible to obtain transfected viable parasites under in vivo growing conditions. The growth curve analysis of transfected and WT B. microti are similar indicating no defects in the transgenic parasites. This study will enable other researchers in understanding the B. microti biology, host modulation and diverse parasite developmental stages using reverse genetics and holds great potential to identify novel drug targets and vaccine development.

Identification and antigenicity of the Babesia caballi spherical body protein 4 (SBP4)

BACKGROUND

The tick-borne intra-erythrocytic apicomplexan Babesia caballi is one of the etiological agents of equine babesiosis, an economically important disease of equids in most tropical and subtropical areas of the world. Discovering candidate antigens for improved diagnostic tools and vaccines remains needed for controlling equine babesiosis. This study describes the B. caballi sbp4 (Bcsbp4) gene and protein (BcSBP4) and analyzes its antigenicity in infected equids.

METHODS

BLAST searches of an uncurated B. caballi assembly genome using the B. bovis SBP4 as a query were carried out, followed by PCR amplification and sequencing of a newly identified BcSBP4. Characterization of this novel gene and protein was performed by bioinformatics analysis, western blots, immunofluorescence (IFA) and an in vitro neutralization test using anti SBP4 peptide antibodies. Antigenicity of recombinant BcSBP4 (rBcSBP4) was tested with sera from field animals (n = 18) using an indirect ELISA (iELISA).

RESULTS

Babesia caballi genome searches using B. bovis SBP4 as a query allowed identification of a novel gene termed Bcsbp4. The Bcsbp4 gene encodes for a protein of 30.58 kDa, which is fully conserved among B. caballi isolates from USA and Egypt. Bioinformatics analysis indicates that BcSBP4 contains a signal peptide and lacks additional transmembrane domains. Expression of BcSBP4 in blood stages of B. caballi was confirmed by western blot and IFA using antibodies against synthetic peptides representing putative B-cell epitopes of BcSBP4 predicted by in silico analysis. In vitro neutralization tests using anti-BcSBP4 peptide antibodies showed a marginal, but statistically significant inhibitory effect on the infectivity of B. caballi merozoites in horse red blood cells. Sera from eight B. caballi-infected equids, but none out of ten negative equid control sera, gave a positive signal in an rBcSBP4 based iELISA.

CONCLUSIONS

The Bcsbp4 gene is expressed in B. caballi blood stages. The BcSBP4 protein is a potential candidate for developing a novel serological test that could detect B. caballi infection in equids in tropical and subtropical countries worldwide.

Transient Transfection of the Zoonotic Parasite Babesia microti

The development of genetic manipulation techniques has been reported in many protozoan parasites over the past few years. However, these techniques have not been established for Babesia microti. Here, we report the first successful transient transfection of B. microti. The plasmids containing the firefly luciferase reporter gene were transfected into B. microti by an AMAXA 4D Nucleofection system. Twenty-four-hour synchronization, the 5'-actin promoter, program FA100, and 50 μg of plasmid DNA constituted the best conditions for the transient transfection of B. microti. This finding is the first step towards a stable transfection method for B. microti, which may contribute to a better understanding of the biology of the parasite.

Interplay between Attenuation- and Virulence-Factors of Babesia bovis and Their Contribution to the Establishment of Persistent Infections in Cattle

Bovine babesiosis is an acute and persistent tick-borne global disease caused mainly by the intraerythrocytic apicomplexan parasites Babesia bovis and B. bigemina. B. bovis infected erythrocytes sequester in blood capillaries of the host (cytoadhesion), causing malaria-like neurological signs. Cytoadhesion and antigenic variation in B. bovis are linked to the expression of members of the Variant Erythrocyte Surface Antigen (VESA) gene family. Animals that survive acute B. bovis infection and those vaccinated with attenuated strains remain persistently infected, suggesting that B. bovis parasites use immune escape mechanisms. However, attenuated B. bovis parasites do not cause neurological signs in vaccinated animals, indicating that virulence or attenuation factors play roles in modulating parasite virulence phenotypes. Artificial overexpression of the SBP2t11 protein, a defined attenuation factor, was associated with reduced cytoadhesion, suggesting a role for this protein as a key modulator of virulence in the parasite. Hereby, we propose a model that might be functional in the modulation of B. bovis virulence and persistence that relies on the interplay among SBP2t, VESA proteins, cytoadhesion, and the immune responses of the host. Elucidation of mechanisms used by the parasite to establish persistent infection will likely contribute to the design of new methods for the control of bovine babesiosis.

Variable and Variant Protein Multigene Families in Babesia bovis Persistence

Cattle infected with Babesia bovis face a bifurcated fate: Either die of the severe acute infection, or survive and carry for many years a highly persistent but generally asymptomatic infection. In this review, the author describes known and potential contributions of three variable or highly variant multigene-encoded families of proteins to persistence in the bovine host, and the mechanisms by which variability arises among these families. Ramifications arising from this variability are discussed.