Expression and strain variation of the novel "small open reading frame" (smorf) multigene family in Babesia bovis

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.

In Silico Survey and Characterization of Babesia microti Functional and Non-Functional Proteases

Human babesiosis caused by the intraerythrocytic apicomplexan Babesia microti is an expanding tick-borne zoonotic disease that may cause severe symptoms and death in elderly or immunocompromised individuals. In light of an increasing resistance of B. microti to drugs, there is a lack of therapeutic alternatives. Species-specific proteases are essential for parasite survival and possible chemotherapeutic targets. However, the repertoire of proteases in B. microti remains poorly investigated. Herein, we employed several combined bioinformatics tools and strategies to organize and identify genes encoding for the full repertoire of proteases in the B. microti genome. We identified 64 active proteases and 25 nonactive protease homologs. These proteases can be classified into cysteine (n = 28), serine (n = 21), threonine (n = 14), asparagine (n = 7), and metallopeptidases (n = 19), which, in turn, are assigned to a total of 38 peptidase families. Comparative studies between the repertoire of B. bovis and B. microti proteases revealed differences among sensu stricto and sensu lato Babesia parasites that reflect their distinct evolutionary history. Overall, this data may help direct future research towards our understanding of the biology and pathogenicity of Babesia parasites and to explore proteases as targets for developing novel therapeutic interventions.

Predicting Protein Therapeutic Candidates for Bovine Babesiosis Using Secondary Structure Properties and Machine Learning

Bovine babesiosis causes significant annual global economic loss in the beef and dairy cattle industry. It is a disease instigated from infection of red blood cells by haemoprotozoan parasites of the genus Babesia in the phylum Apicomplexa. Principal species are Babesia bovis, Babesia bigemina, and Babesia divergens. There is no subunit vaccine. Potential therapeutic targets against babesiosis include members of the exportome. This study investigates the novel use of protein secondary structure characteristics and machine learning algorithms to predict exportome membership probabilities. The premise of the approach is to detect characteristic differences that can help classify one protein type from another. Structural properties such as a protein’s local conformational classification states, backbone torsion angles ϕ (phi) and ψ (psi), solvent-accessible surface area, contact number, and half-sphere exposure are explored here as potential distinguishing protein characteristics. The presented methods that exploit these structural properties via machine learning are shown to have the capacity to detect exportome from non-exportome Babesia bovis proteins with an 86–92% accuracy (based on 10-fold cross validation and independent testing). These methods are encapsulated in freely available Linux pipelines setup for automated, high-throughput processing. Furthermore, proposed therapeutic candidates for laboratory investigation are provided for B. bovis, B. bigemina, and two other haemoprotozoan species, Babesia canis, and Plasmodium falciparum.

Applying Machine Learning to Predict the Exportome of Bovine and Canine Babesia Species That Cause Babesiosis

Babesia infection of red blood cells can cause a severe disease called babesiosis in susceptible hosts. Bovine babesiosis causes global economic loss to the beef and dairy cattle industries, and canine babesiosis is considered a clinically significant disease. Potential therapeutic targets against bovine and canine babesiosis include members of the exportome, i.e., those proteins exported from the parasite into the host red blood cell. We developed three machine learning-derived methods (two novel and one adapted) to predict for every known Babesia bovis, Babesia bigemina, and Babesia canis protein the probability of being an exportome member. Two well-studied apicomplexan-related species, Plasmodium falciparum and Toxoplasma gondii, with extensive experimental evidence on their exportome or excreted/secreted proteins were used as important benchmarks for the three methods. Based on 10-fold cross validation and multiple train-validation-test splits of training data, we expect that over 90% of the predicted probabilities accurately provide a secretory or non-secretory indicator. Only laboratory testing can verify that predicted high exportome membership probabilities are creditable exportome indicators. However, the presented methods at least provide those proteins most worthy of laboratory validation and will ultimately save time and money.

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.

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.

Babesia bovis Rad51 ortholog influences switching of ves genes but is not essential for segmental gene conversion in antigenic variation

The tick-borne apicomplexan parasite, Babesia bovis, a highly persistent bovine pathogen, expresses VESA1 proteins on the infected erythrocyte surface to mediate cytoadhesion. The cytoadhesion ligand, VESA1, which protects the parasite from splenic passage, is itself protected from a host immune response by rapid antigenic variation. B. bovis relies upon segmental gene conversion (SGC) as a major mechanism to vary VESA1 structure. Gene conversion has been considered a form of homologous recombination (HR), a process for which Rad51 proteins are considered pivotal components. This could make BbRad51 a choice target for development of inhibitors that both interfere with parasite genome integrity and disrupt HR-dependent antigenic variation. Previously, we knocked out the Bbrad51 gene from the B. bovis haploid genome, resulting in a phenotype of sensitivity to methylmethane sulfonate (MMS) and apparent loss of HR-dependent integration of exogenous DNA. In a further characterization of BbRad51, we demonstrate here that ΔBbrad51 parasites are not more sensitive than wild-type to DNA damage induced by γ-irradiation, and repair their genome with similar kinetics. To assess the need for BbRad51 in SGC, RT-PCR was used to observe alterations to a highly variant region of ves1α transcripts over time. Mapping of these amplicons to the genome revealed a significant reduction of in situ transcriptional switching (isTS) among ves loci, but not cessation. By combining existing pipelines for analysis of the amplicons, we demonstrate that SGC continues unabated in ΔBbrad51 parasites, albeit at an overall reduced rate, and a reduction in SGC tract lengths was observed. By contrast, no differences were observed in the lengths of homologous sequences at which recombination occurred. These results indicate that, whereas BbRad51 is not essential to babesial antigenic variation, it influences epigenetic control of ves loci, and its absence significantly reduces successful variation. These results necessitate a reconsideration of the likely enzymatic mechanism(s) underlying SGC and suggest the existence of additional targets for development of small molecule inhibitors.

B. bovis establishes highly persistent infections in cattle, in part by using cytoadhesion to avoid passage through the spleen. While protective, a host antibody response targeting the cytoadhesion ligand is quickly rendered ineffective by antigenic variation. In B. bovis, antigenic variation relies heavily upon segmental gene conversion (SGC), presumed to be a form of homologous recombination (HR), to generate variants. As Rad51 is generally considered essential to HR, we investigated its contribution to SGC. While diminishing the parasite’s capacity for HR-dependent integration of exogenous DNA, the loss of BbRad51 did not affect the parasite’s sensitivity to ionizing radiation, overall genome stability, or competence for SGC. Instead, loss of BbRad51 diminished the extent of in situ transcriptional switching (isTS) among ves gene loci, the accumulation of SGC recombinants, and the mean lengths of SGC sequence tracts. Given the overall reductions in VESA1 variability, compromise of the parasite’s capacity for in vivo persistence is predicted.

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.

Unravelling the cellular and molecular pathogenesis of bovine babesiosis: is the sky the limit?

The global impact of bovine babesiosis caused by the tick-borne apicomplexan parasites Babesia bovis, Babesia bigemina and Babesia divergens is vastly underappreciated. These parasites invade and multiply asexually in bovine red blood cells (RBCs), undergo sexual reproduction in their tick vectors (Rhipicephalus spp. for B. bovis and B. bigemina, and Ixodes ricinus for B. divergens) and have a trans-ovarial mode of transmission. Babesia parasites can cause acute and persistent infections to adult naïve cattle that can occur without evident clinical signs, but infections caused by B. bovis are associated with more severe disease and increased mortality, and are considered to be the most virulent agent of bovine babesiosis. In addition, babesiosis caused by B. divergens has an important zoonotic potential. The disease caused by B. bovis and B. bigemina can be controlled, at least in part, using therapeutic agents or vaccines comprising live-attenuated parasites, but these methods are limited in terms of their safety, ease of deployability and long-term efficacy, and improved control measures are urgently needed. In addition, expansion of tick habitats due to climate change and other rapidly changing environmental factors complicate efficient control of these parasites. While the ability to cause persistent infections facilitates transmission and persistence of the parasite in endemic regions, it also highlights their capacity to evade the host immune responses. Currently, the mechanisms of immune responses used by infected bovines to survive acute and chronic infections remain poorly understood, warranting further research. Similarly, molecular details on the processes leading to sexual reproduction and the development of tick-stage parasites are lacking, and such tick-specific molecules can be targets for control using alternative transmission blocking vaccines. In this review, we identify and examine key phases in the life-cycle of Babesia parasites, including dependence on a tick vector for transmission, sexual reproduction of the parasite in the midgut of the tick, parasite-dependent invasion and egression of bovine RBCs, the role of the spleen in the clearance of infected RBCs (IRBCs), and age-related disease resistance in cattle, as opportunities for developing improved control measures. The availability of integrated novel research approaches including "omics" (such as genomics, transcriptomics, and proteomics), gene modification, cytoadhesion assays, RBC invasion assays and methods for in vitro induction of sexual-stage parasites will accelerate our understanding of parasite vulnerabilities. Further, producing new knowledge on these vulnerabilities, as well as taking full advantage of existing knowledge, by filling important research gaps should result in the development of next-generation vaccines to control acute disease and parasite transmission. Creative and effective use of current and future technical and computational resources are needed, in the face of the numerous challenges imposed by these highly evolved parasites, for improving the control of this disease. Overall, bovine babesiosis is recognised as a global disease that imposes a serious burden on livestock production and human livelihood, but it largely remains a poorly controlled disease in many areas of the world. Recently, important progress has been made in our understanding of the basic biology and host-parasite interactions of Babesia parasites, yet a good deal of basic and translational research is still needed to achieve effective control of this important disease and to improve animal and human health.

Characterizing Cell Heterogeneity Using PCR Fingerprinting of Surface Multigene Families in Protozoan Parasites

Parasites counteract the action of the immune system and other environmental pressures by modulating and changing the composition of their cell surfaces. Surface multigene protein families are defined not only by highly variable regions in length and/or sequence exposed to the outer space but also by conserved sequences codifying for the signal peptide, hydrophobic C-terminal regions necessary for GPI modifications, as well as conserved UTR regions for mRNA regulation. The method here presented exploits these conserved signatures for characterizing variations in the mRNA expression of clonal cell populations of protozoan parasites using a combination of nested PCR amplification and capillary electrophoresis. With this workflow, in silico gels from isolated cell clones can be generated, thus providing an excellent tool for analyzing cellular heterogeneity in protozoan parasites.

Shared elements of host-targeting pathways among apicomplexan parasites of differing lifestyles

Apicomplexans are a diverse group of obligate parasites occupying different intracellular niches that require modification to meet the needs of the parasite. To efficiently manipulate their environment, apicomplexans translocate numerous parasite proteins into the host cell. Whereas some parasites remain contained within a parasitophorous vacuole membrane (PVM) throughout their developmental cycle, others do not, a difference that affects the machinery needed for protein export. A signal-mediated pathway for protein export into the host cell has been characterized in Plasmodium parasites, which maintain the PVM. Here, we functionally demonstrate an analogous host-targeting pathway involving organellar staging prior to secretion in the related bovine parasite, Babesia bovis, a parasite that destroys the PVM shortly after invasion. Taking into account recent identification of a similar signal-mediated pathway in the coccidian parasite Toxoplasma gondii, we suggest a model in which this conserved pathway has evolved in multiple steps from signal-mediated trafficking to specific secretory organelles for controlled secretion to a complex protein translocation process across the PVM.

Rhoptry-associated protein (rap-1) genes in the sheep pathogen Babesia sp. Xinjiang: Multiple transcribed copies differing by 3' end repeated sequences

Sheep babesiosis occurs mainly in tropical and subtropical areas. The sheep parasite Babesia sp. Xinjiang is widespread in China, and our goal is to characterize rap-1 (rhoptry-associated protein 1) gene diversity and expression as a first step of a long term goal aiming at developing a recombinant subunit vaccine. Seven different rap-1a genes were amplified in Babesia sp. Xinjiang, using degenerate primers designed from conserved motifs. Rap-1b and rap-1c gene types could not be identified. In all seven rap-1a genes, the 5' regions exhibited identical sequences over 936 nt, and the 3' regions differed at 28 positions over 147 nt, defining two types of genes designated α and β. The remaining 3' part varied from 72 to 360 nt in length, depending on the gene. This region consists of a succession of two to ten 36 nt repeats, which explains the size differences. Even if the nucleotide sequences varied, 6 repeats encoded the same stretch of amino acids. Transcription of at least four α and two β genes was demonstrated by standard RT-PCR.

The evolutionary dynamics of variant antigen genes in Babesia reveal a history of genomic innovation underlying host-parasite interaction

Babesia spp. are tick-borne, intraerythrocytic hemoparasites that use antigenic variation to resist host immunity, through sequential modification of the parasite-derived variant erythrocyte surface antigen (VESA) expressed on the infected red blood cell surface. We identified the genomic processes driving antigenic diversity in genes encoding VESA (ves1) through comparative analysis within and between three Babesia species, (B. bigemina, B. divergens and B. bovis). Ves1 structure diverges rapidly after speciation, notably through the evolution of shortened forms (ves2) from 5′ ends of canonical ves1 genes. Phylogenetic analyses show that ves1 genes are transposed between loci routinely, whereas ves2 genes are not. Similarly, analysis of sequence mosaicism shows that recombination drives variation in ves1 sequences, but less so for ves2, indicating the adoption of different mechanisms for variation of the two families. Proteomic analysis of the B. bigemina PR isolate shows that two dominant VESA1 proteins are expressed in the population, whereas numerous VESA2 proteins are co-expressed, consistent with differential transcriptional regulation of each family. Hence, VESA2 proteins are abundant and previously unrecognized elements of Babesia biology, with evolutionary dynamics consistently different to those of VESA1, suggesting that their functions are distinct.

Comparative transcriptome analysis of geographically distinct virulent and attenuated Babesia bovis strains reveals similar gene expression changes through attenuation

Background

Loss of virulence is a phenotypic adaptation commonly seen in prokaryotic and eukaryotic pathogens. This mechanism is not well studied, especially in organisms with multiple host and life cycle stages such as Babesia, a tick-transmitted hemoparasite of humans and animals. B. bovis, which infects cattle, has naturally occurring virulent strains that can be reliably attenuated in vivo. Previous studies suggest the virulence loss mechanism may involve post-genomic modification. We investigated the transcriptome profiles of two geographically distinct B. bovis virulent and attenuated strain pairs to better understand virulence loss and to gain insight into pathogen adaptation strategies.

Results

Expression microarray and RNA-sequencing approaches were employed to compare transcriptome profiles of two B. bovis strain pairs, with each pair consisting of a virulent parental and its attenuated derivative strain. Differentially regulated transcripts were identified within each strain pair. These included genes encoding for VESA1, SmORFs, undefined membrane and hypothetical proteins. The majority of individual specific gene transcripts differentially regulated within a strain were not shared between the two strains. There was a disproportionately greater number of ves genes upregulated in the virulent parental strains. When compared with their attenuated derivatives, divergently oriented ves genes were included among the upregulated ves genes in the virulent strains, while none of the upregulated ves genes in the attenuated derivatives were oriented head to head. One gene family whose specific members were consistently and significantly upregulated in expression in both attenuated strains was spherical body protein (SBP) 2 encoding gene where SBP2 truncated copies 7, 9 and 11 transcripts were all upregulated.

Conclusions

We conclude that ves heterodimer pair upregulation and overall higher frequency of ves gene expressions in the virulent strains is consistent with the involvement of this gene family in virulence. This is logical given the role of VESA1 proteins in cytoadherence of infected cells to endothelial cells. However, upregulation of some ves genes in the attenuated derivatives suggests that the consequence of upregulation is gene-specific. Furthermore, upregulation of the spherical body protein 2 gene family may play a role in the attenuated phenotype. Exactly how these two gene families may contribute to the loss or gain of virulence is discussed.

Differential expression of three members of the multidomain adhesion CCp family in Babesia bigemina, Babesia bovis and Theileria equi

Members of the CCp protein family have been previously described to be expressed on gametocytes of apicomplexan Plasmodium parasites. Knocking out Plasmodium CCp genes blocks the development of the parasite in the mosquito vector, making the CCp proteins potential targets for the development of a transmission-blocking vaccine. Apicomplexans Babesia bovis and Babesia bigemina are the causative agents of bovine babesiosis, and apicomplexan Theileria equi causes equine piroplasmosis. Bovine babesiosis and equine piroplasmosis are the most economically important parasite diseases that affect worldwide cattle and equine industries, respectively. The recent sequencing of the B. bovis and T. equi genomes has provided the opportunity to identify novel genes involved in parasite biology. Here we characterize three members of the CCp family, named CCp1, CCp2 and CCp3, in B. bigemina, B. bovis and T. equi. Using B. bigemina as an in vitro model, expression of all three CCp genes and proteins was demonstrated in temperature-induced sexual stages. Transcripts for all three CCp genes were found in vivo in blood stages of T. equi, and transcripts for CCp3 were detected in vivo in blood stages of B. bovis. However, no protein expression was detected in T. equi blood stages or B. bovis blood stages or B. bovis tick stages. Collectively, the data demonstrated a differential pattern of expression of three orthologous genes of the multidomain adhesion CCp family by B. bigemina, B. bovis and T. equi. The novel CCp members represent potential targets for innovative approaches to control bovine babesiosis and equine piroplasmosis.

Bioinformatic prediction of the exportome of Babesia bovis and identification of novel proteins in parasite-infected red blood cells

Babesia bovis is a pathogen of considerable economic significance to the livestock industry worldwide but the precise mechanisms by which this parasite causes disease in susceptible cattle remain poorly understood. It is clear, however, that alterations to the structure and function of red blood cells in which the parasites reside and replicate play an important role in pathogenesis and that these are secondary to the export of numerous, currently unknown and uncharacterised parasite-encoded proteins. Using a rational bioinformatic approach, we have identified a set of 362 proteins (117 of which are hypothetical) that we predict encompasses the B. bovis exportome. These exported proteins are likely to be trafficked to various cellular locations, with a subset destined for the red blood cell cytosol or the red blood cell cytoskeleton. These proteins are likely to play important roles in mediating the pathogenesis of babesiosis. We have selected three novel proteins and confirmed their predicted export and localisation within the host red blood cell by immunofluorescence using specific antibodies raised against these proteins. Complete characterisation of these novel exported parasite proteins will help elucidate their function within the host red blood cell and assist in identification of new therapeutic targets for babesiosis.

Bovine babesiosis in the 21st century: advances in biology and functional genomics

Bovine babesiosis caused by the protozoan parasite, Babesia bovis, remains a significant cause of avoidable economic losses to the livestock industry in many countries throughout the world. The molecular mechanisms underlying the pathophysiology of severe disease in susceptible cattle are not well understood and the tools available to study the biology of the parasite, including technologies for genetic manipulation, have only recently been developed. Recent availability of multiple parasite genomes and bioinformatic tools, in combination with the development of new biological reagents, will facilitate our better understanding of the parasite. This will ultimately assist in the identification of novel targets for the development of new therapeutics and vaccines. Here we describe some recent advances in Babesia research and highlight some important challenges for the future.