Universal primers suitable to assess population dynamics reveal apparent mutually exclusive transcription of the Babesia bovis ves1alpha gene

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.

Integration of DNA Repair, Antigenic Variation, Cytoadhesion, and Chance in Babesia Survival: A Perspective

Apicomplexan parasites live in hostile environments in which they are challenged chemically and their hosts attempt in many ways to kill them. In response, the parasites have evolved multiple mechanisms that take advantage of these challenges to enhance their survival. Perhaps the most impressive example is the evolutionary co-option of DNA repair mechanisms by the parasites as a means to rapidly manipulate the structure, antigenicity, and expression of the products of specific multigene families. The purpose of variant proteins that mediate cytoadhesion has long been thought to be primarily the avoidance of splenic clearance. Based upon known biology, I present an alternative perspective in which it is survival of the oxidative environment within which Babesia spp. parasites live that has driven integration of DNA repair, antigenic variation, and cytoadhesion, and speculate on how genome organization affects that integration. This perspective has ramifications for the development of parasite control strategies.

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.

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.

Unusual chromatin structure associated with monoparalogous transcription of the Babesia bovis ves multigene family

Rapid antigenic variation in Babesia bovis involves the variant erythrocyte surface antigen-1 (VESA1), a heterodimeric protein with subunits encoded by two branches of the ves multigene family. The ves1α and ves1β gene pair encoding VESA1a and 1b, respectively, are transcribed in a monoparalogous manner from a single locus of active ves transcription (LAT), just one of many quasi-palindromic ves loci. To determine whether this organization plays a role in transcriptional regulation, chromatin structure was first assessed. Limited treatment of isolated nuclei with micrococcal nuclease to assay nucleosomal patterning revealed a periodicity of 156-159 bp in both bulk chromatin and specific gene coding regions. This pattern also was maintained in the intergenic regions (IGr) of non-transcribed ves genes. In contrast, the LAT IGr adopts a unique pattern, yielding an apparent cluster of five closely-spaced hypersensitive sites flanked by regions of reduced nucleosomal occupancy. ves loci fall into three patterns of overall sensitivity to micrococcal nuclease or DNase I digestion, with only the LAT being consistently very sensitive. Non-transcribed ves genes are inconsistent in their sensitivity to the two enzymatic probes. Non-linear DNA structure in chromatin was investigated to determine whether unique structure arising as a result of the quasi-palindromic nature of the LAT may effect transcriptional control. The in vitro capacity of ves IGr sequences to adopt stable higher-order DNA structure is demonstrated here, but the presence of such structure in vivo was not supported. Based upon these results a working model is proposed for the chromatin structural remodeling responsible for the sequential expression of ves multigene family members from divergently-organized loci.

Characterization of the unusual bidirectional ves promoters driving VESA1 expression and associated with antigenic variation in Babesia bovis

Rapid clonal antigenic variation in Babesia bovis involves the variant erythrocyte surface antigen-1 (VESA1) protein expressed on the infected-erythrocyte surface. Because of the significance of this heterodimeric protein for demonstrated mechanisms of parasite survival and virulence, there is a need to understand how expression of the ves multigene family encoding this protein is controlled. As an initial step toward this goal, we present here initial characterization of the ves promoter driving transcription of VESA1a and -1b subunits. A series of transfection constructs containing various sequence elements from the in vivo locus of active ves transcription (LAT) were used to drive expression of the firefly luciferase gene in a dual luciferase-normalized assay. The results of this approach reveal the presence of two bidirectional promoter activities within the 434-bp intergenic region (IGr), influenced by putative regulatory sequences embedded within the flanking ves1α and ves1β genes. Repressor-like effects on the apposing gene were observed for intron 1 of both ves1α and ves1β. This effect is apparently not dependent upon intronic promoter activity and acts only in cis. The expression of genes within the ves family is likely modulated by local elements embedded within ves coding sequences outside the intergenic promoter region in concert with chromatin modifications. These results provide a framework to help us begin to understand gene regulation during antigenic variation in B. bovis.

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

Small open reading frame (smorf) genes comprise the second largest Babesia bovis multigene family. All known 44 variant smorf genes are located in close chromosomal proximity to ves1 genes, which encode proteins that mediate cytoadhesion and contribute to immune evasion. In this study, we characterised the general topology of smorf genes and investigated the gene repertoire, transcriptional profile and SMORF expression in two distinct strains, T2Bo and Mo7. Sequence analysis using degenerate primers identified additional smorf genes in each strain and demonstrated that the smorf gene repertoire varies between strains, with conserved and unique genes in both. Smorf genes have multiple semi-conserved and variable blocks, and a large hypervariable insertion in 20 of the 44 genes defines two major branches of the family, termed smorf A and smorf B. A total of 32 smorf genes are simultaneously transcribed in T2Bo strain B. bovis merozoites obtained from deep brain tissue of an acutely infected animal. SMORF peptide-specific antiserum bound in immunoblots to multiple proteins with a range of sizes predicted by smorf genes, confirming translation of smorf gene products from these transcripts. These results indicate that the smorf multigene family is larger than previously described and demonstrate that smorf genes are expressed and are undergoing variation, both within strains and in a lineage-specific pattern independent of strain specificity. The function of these novel proteins is unknown.

Recent insights into alteration of red blood cells by Babesia bovis: moovin' forward

Over the past decade or so, our understanding of the biology of apicomplexan parasites has increased dramatically, particularly in the case of malaria. Notable achievements are the availability of complete genome sequences, transcriptome and proteome profiles and the establishment of in vitro transfection techniques for asexual-stage malaria parasites. Interestingly, despite their major economic importance and striking similarities with malaria, Babesia parasites have been relatively ignored, but change is on the horizon. Here, we bring together recent work on Babesia bovis parasites which are beginning to unravel the molecular mechanisms that underlie the pathogenesis of babesiosis and highlight some opportunities and challenges that lie ahead.

The Babesia bovis VESA1 virulence factor subunit 1b is encoded by the 1beta branch of the ves multigene family

Babesia bovis, an intraerythrocytic parasite of cattle, establishes persistent infections of extreme duration. This is accomplished, at least in part, through rapid antigenic variation of a heterodimeric virulence factor, the variant erythrocyte surface antigen-1 (VESA1) protein. Previously, the VESA1a subunit was demonstrated to be encoded by a 1alpha member of the ves multigene family. Since its discovery the 1beta branch of this multigene family has been hypothesized to encode the VESA1b polypeptide, but formal evidence for this connection has been lacking. Here, we provide evidence that products of ves1beta genes are rapidly variant in antigenicity and size-polymorphic, matching known VESA1b polypeptides. Importantly, the ves1beta-encoded antigens are co-precipitated with VESA1a during immunoprecipitation with anti-VESA1a monoclonal antibodies, and antisera to ves1beta polypeptide co-precipitate VESA1a. Further, the ves1beta-encoded antigens significantly co-localize with VESA1a on the infected-erythrocyte membrane surface of live cells. These characteristics all match known properties of VESA1b, allowing us to conclude that the ves1beta gene divergently apposing the ves1beta gene within the locus of active ves transcription (LAT) encodes the 1b subunit of the VESA1 cytoadhesion ligand. However, the extent and stoichiometry of VESA1a and 1b co-localization on the surface of individual cells is quite variable, implicating competing effects on transcription, translation, or trafficking of the two subunits. These results provide essential information facilitating further investigation into this parasite virulence factor.