Antigenic Switching Research Articles

Transcriptional Regulation of Telomeric Expression Sites and Antigenic Variation in Trypanosomes.

Introduction: Trypanosoma brucei uses antigenic variation to evade the host antibody clearance by periodically changing its Variant Surface Glycoprotein (VSGs) coat. T. brucei encode over 2,500 VSG genes and pseudogenes, however they transcribe only one VSG gene at time from one of the 20 telomeric Expression Sites (ESs). VSGs are transcribed in a monoallelic fashion by RNA polymerase I from an extranucleolar site named ES body. VSG antigenic switching occurs by transcriptional switching between telomeric ESs or by recombination of the VSG gene expressed. VSG expression is developmentally regulated and its transcription is controlled by epigenetic mechanisms and influenced by a telomere position effect.Conclusion: Here, we discuss 1) the molecular basis underlying transcription of telomeric ESs and VSG antigenic switching; 2) the current knowledge of VSG monoallelic expression; 3) the role of inositol phosphate pathway in the regulation of VSG expression and switching; and 4) the developmental regulation of Pol I transcription of procyclin and VSG genes.

Cis-acting DNA elements flanking the variable major protein expression site of Borrelia hermsii are required for murine persistence.

In Borrelia hermsii, antigenic variation occurs as a result of a nonreciprocal gene conversion event that places one of ~60 silent variable major protein genes downstream of a single, transcriptionally active promoter. The upstream homology sequence (UHS) and downstream homology sequence (DHS) are two putative cis‐acting DNA elements that have been predicted to serve as crossover points for homologous recombination. In this report, a targeted deletion/in cis complementation technique was used to directly evaluate the role for these elements in antigenic switching. The results demonstrate that deletion of the expression site results in an inability of the pathogen to relapse in immunocompetent mice, and that the utilized technique was successful in producing complemented mutants that are capable of antigenic switching. Additional complemented clones with mutations in the UHS and DHS of the expressed locus were then generated and evaluated for their ability to relapse in immunocompetent mice. Mutation of the UHS and inverted repeat sequence within the DHS rendered these mutants incapable of relapsing. Overall, the results establish the requirement of the inverted repeat of the DHS for antigenic switching, and support the importance of the UHS for B. hermsii persistence in the mammalian host.

Xer1-independent mechanisms of Vpma phase variation in Mycoplasma agalactiae are triggered by Vpma-specific antibodies

Despite their small genomes mycoplasmas maintain large multigene families devoted to surface antigenic variation. Although implicated as important factors for mycoplasma pathogenicity and persistence, the role of these antigenic switches in host immune evasion has never been unequivocally proven in these minimalist microbes. Mycoplasma agalactiae exhibits antigenic variation due to Xer1-mediated site-specific DNA inversions of vpma genes encoding abundant multiple surface lipoproteins. To evaluate the biological significance of Vpma oscillations the xer1 recombinase gene has been disrupted in earlier studies to abolish Vpma switching and to generate stable phase-locked mutants (PLMs) steadily expressing a single Vpma product. However, in previous animal infection studies, surprisingly these PLMs switched to new different Vpma phenotypes. The aim of the current study was to demonstrate the influence of anti-Vpma antibodies on change of Vpma expression in PLMs as well as on the wildtype strain. In in vitro assays it is shown that wild type M. agalactiae escapes the negative effects of Vpma-specific antibodies by high-frequency Xer1-mediated switching to alternative Vpma phenoytpes. Even for Xer1-disrupted PLMs that stably expressed the same Vpma for several in vitro generations, the presence of the corresponding Vpma-specific antibody caused repression of the target Vpma and induction of new Vpma phenotypes by novel complex vpma rearrangements like intragenic deletions and gene chimeras. These Xer1-independent vpma recombinations correlated very well with similar PLM switches observed in vivo in an earlier independent study, clearly demonstrating that Vpma phase variation is necessary to express ‘Vpma immune evasion proteins’ in order to escape the immune response and to survive in the immunocompetent host. The data clearly demonstrate that although the Xer1 recombinase is the sole factor responsible for Vpma switching of wild type M. agalactiae in vitro, other alternative molecular switches operate in its absence under the selective pressure of the immune response. Furthermore, this evasion from the immune attack of the host involves complex vpma rearrangements, a causal relationship that was so far never demonstrated for M. agalactiae, thereby illustrating novel features of its regulation under immune pressure. The results are anticipated to have a direct impact on understanding the in vivo role of surface antigenic variation systems and the immune evasion tactics of other pathogenic mycoplasma species.

Plasmodium knowlesi: a superb in vivo nonhuman primate model of antigenic variation in malaria.

Antigenic variation in malaria was discovered in Plasmodium knowlesi studies involving longitudinal infections of rhesus macaques (M. mulatta). The variant proteins, known as the P. knowlesi Schizont Infected Cell Agglutination (SICA) antigens and the P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) antigens, expressed by the SICAvar and var multigene families, respectively, have been studied for over 30 years. Expression of the SICA antigens in P. knowlesi requires a splenic component, and specific antibodies are necessary for variant antigen switch events in vivo. Outstanding questions revolve around the role of the spleen and the mechanisms by which the expression of these variant antigen families are regulated. Importantly, the longitudinal dynamics and molecular mechanisms that govern variant antigen expression can be studied with P. knowlesi infection of its mammalian and vector hosts. Synchronous infections can be initiated with established clones and studied at multi-omic levels, with the benefit of computational tools from systems biology that permit the integration of datasets and the design of explanatory, predictive mathematical models. Here we provide an historical account of this topic, while highlighting the potential for maximizing the use of P. knowlesi - macaque model systems and summarizing exciting new progress in this area of research.

A quorum sensing-independent path to stumpy development in Trypanosoma brucei.

For persistent infections of the mammalian host, African trypanosomes limit their population size by quorum sensing of the parasite-excreted stumpy induction factor (SIF), which induces development to the tsetse-infective stumpy stage. We found that besides this cell density-dependent mechanism, there exists a second path to the stumpy stage that is linked to antigenic variation, the main instrument of parasite virulence. The expression of a second variant surface glycoprotein (VSG) leads to transcriptional attenuation of the VSG expression site (ES) and immediate development to tsetse fly infective stumpy parasites. This path is independent of SIF and solely controlled by the transcriptional status of the ES. In pleomorphic trypanosomes varying degrees of ES-attenuation result in phenotypic plasticity. While full ES-attenuation causes irreversible stumpy development, milder attenuation may open a time window for rescuing an unsuccessful antigenic switch, a scenario that so far has not been considered as important for parasite survival.

Specific histone modifications play critical roles in the control of encystation and antigenic variation in the early-branching eukaryote Giardia lamblia

During evolution, parasitic microorganisms have faced the challenges of adapting to different environments to colonize a variety of hosts. Giardia lamblia, a common cause of intestinal disease, has developed fascinating strategies to adapt both outside and inside its host’s intestine, such as trophozoite differentiation into cyst and the switching of its major surface antigens. How gene expression is regulated during these adaptive processes remains undefined. Giardia lacks some typical eukaryotic features, like canonical transcription factors, linker histone H1, and complex promoter regions; suggesting that post-transcriptional and translational control of gene expression is essential for parasite survival. However, epigenetic factors may also play critical roles at the transcriptional level. Here, we describe the most common post-translational histone modifications; characterize enzymes involved in these reactions, and analyze their association with the Giardia’s differentiation processes. We present evidence that NAD+-dependent and NAD+-independent histone deacetylases regulate encystation; however, a unique NAD+-independent histone deacetylase modulate antigenic switching. The rates of acetylation of H4K8 and H4K16 are critical for encystation, whereas a decrease in acetylation of H4K8 and methylation of H3K9 occur preferentially during antigenic variation. These results show the complexity of the mechanisms regulating gene expression in this minimalistic protozoan parasite.

A transcription-independent epigenetic mechanism is associated with antigenic switching in Trypanosoma brucei.

Antigenic variation in Trypanosoma brucei relies on periodic switching of variant surface glycoproteins (VSGs), which are transcribed monoallelically by RNA polymerase I from one of about 15 bloodstream expression sites (BES). Chromatin of the actively transcribed BES is depleted of nucleosomes, but it is unclear if this open conformation is a mere consequence of a high rate of transcription, or whether it is maintained by a transcription-independent mechanism. Using an inducible BES-silencing reporter strain, we observed that chromatin of the active BES remains open for at least 24 hours after blocking transcription. This conformation is independent of the cell-cycle stage, but dependent upon TDP1, a high mobility group box protein. For two days after BES silencing, we detected a transient and reversible derepression of several silent BESs within the population, suggesting that cells probe other BESs before commitment to one, which is complete by 48 hours. FACS sorting and subsequent subcloning confirmed that probing cells are switching intermediates capable of returning to the original BES, switch to the probed BES or to a different BES. We propose that regulation of BES chromatin structure is an epigenetic mechanism important for successful antigenic switching.

The within-host dynamics of African trypanosome infections

African trypanosomes are single-celled protozoan parasites that are capable of long-term survival while living extracellularly in the bloodstream and tissues of mammalian hosts. Prolonged infections are possible because trypanosomes undergo antigenic variation—the expression of a large repertoire of antigenically distinct surface coats, which allows the parasite population to evade antibody-mediated elimination. The mechanisms by which antigen genes become activated influence their order of expression, most likely by influencing the frequency of productive antigen switching, which in turn is likely to contribute to infection chronicity. Superimposed upon antigen switching as a contributor to trypanosome infection dynamics is the density-dependent production of cell-cycle arrested parasite transmission stages, which limit the infection while ensuring parasite spread to new hosts via the bite of blood-feeding tsetse flies. Neither antigen switching nor developmental progression to transmission stages is driven by the host. However, the host can contribute to the infection dynamic through the selection of distinct antigen types, the influence of genetic susceptibility or trypanotolerance and the potential influence of host-dependent effects on parasite virulence, development of transmission stages and pathogenicity. In a zoonotic infection cycle where trypanosomes circulate within a range of host animal populations, and in some cases humans, there is considerable scope for a complex interplay between parasite immune evasion, transmission potential and host factors to govern the profile and outcome of infection.

Inositol phosphate pathway controls transcription of telomeric expression sites in trypanosomes

African trypanosomes evade clearance by host antibodies by periodically changing their variant surface glycoprotein (VSG) coat. They transcribe only one VSG gene at a time from 1 of about 20 telomeric expression sites (ESs). They undergo antigenic variation by switching transcription between telomeric ESs or by recombination of the VSG gene expressed. We show that the inositol phosphate (IP) pathway controls transcription of telomeric ESs and VSG antigenic switching in Trypanosoma brucei. Conditional knockdown of phosphatidylinositol 5-kinase (TbPIP5K) or phosphatidylinositol 5-phosphatase (TbPIP5Pase) or overexpression of phospholipase C (TbPLC) derepresses numerous silent ESs in T. brucei bloodstream forms. The derepression is specific to telomeric ESs, and it coincides with an increase in the number of colocalizing telomeric and RNA polymerase I foci in the nucleus. Monoallelic VSG transcription resumes after reexpression of TbPIP5K; however, most of the resultant cells switched the VSG gene expressed. TbPIP5K, TbPLC, their substrates, and products localize to the plasma membrane, whereas TbPIP5Pase localizes to the nucleus proximal to telomeres. TbPIP5Pase associates with repressor/activator protein 1 (TbRAP1), and their telomeric silencing function is altered by TbPIP5K knockdown. These results show that specific steps in the IP pathway control ES transcription and antigenic switching in T. brucei by epigenetic regulation of telomere silencing.

Quantitative proteomic analysis of Giardia duodenalis assemblage A: A baseline for host, assemblage, and isolate variation.

Giardia duodenalis is a gastrointestinal protozoan parasite of vertebrates and is a species complex comprised of eight assemblages, with the zoonotic assemblage A one of two subtypes infective for humans. With increasing genomic and transcriptomic data publicly available through the centralized giardiaDB.org, we have quantitatively analyzed the proteomes of eight G. duodenalis assemblage A strains (seven A1 and one A2) to provide a proteomic baseline to complement the available data. A nonredundant total of 1197 subassemblage A1 proteins and 719 subassemblage A2 proteins were identified with an average of 770 proteins in each strain. The eight strains were also searched against both assemblage A genome sequences (subassemblage A1 and A2 genomes) and demonstrated subassemblage specific differences in protein identifications, especially for variable gene families. Substantial differences were observed in the numbers and abundance in the variable surface protein family, and two different variable surface protein expression profiles that were independent of host origin, subassemblage, or geographic origin. We hypothesize that this variation in surface antigen switching events may be related to karotype and chromosomal variation, which would indicate an assemblage-independent mechanism of diversity generation in G. duodenalis. All MS data have been deposited in the ProteomeXchange with identifier PXD001272 (http://proteomecentral.proteomexchange.org/dataset/PXD001272).

DNA double-strand breaks and telomeres play important roles in trypanosoma brucei antigenic variation.

Human-infecting microbial pathogens all face a serious problem of elimination by the host immune response. Antigenic variation is an effective immune evasion mechanism where the pathogen regularly switches its major surface antigen. In many cases, the major surface antigen is encoded by genes from the same gene family, and its expression is strictly monoallelic. Among pathogens that undergo antigenic variation, Trypanosoma brucei (a kinetoplastid), which causes human African trypanosomiasis, Plasmodium falciparum (an apicomplexan), which causes malaria, Pneumocystis jirovecii (a fungus), which causes pneumonia, and Borrelia burgdorferi (a bacterium), which causes Lyme disease, also express their major surface antigens from loci next to the telomere. Except for Plasmodium, DNA recombination-mediated gene conversion is a major pathway for surface antigen switching in these pathogens. In the last decade, more sophisticated molecular and genetic tools have been developed in T. brucei, and our knowledge of functions of DNA recombination in antigenic variation has been greatly advanced. VSG is the major surface antigen in T. brucei. In subtelomeric VSG expression sites (ESs), VSG genes invariably are flanked by a long stretch of upstream 70-bp repeats. Recent studies have shown that DNA double-strand breaks (DSBs), particularly those in 70-bp repeats in the active ES, are a natural potent trigger for antigenic variation in T. brucei. In addition, telomere proteins can influence VSG switching by reducing the DSB amount at subtelomeric regions. These findings will be summarized and their implications will be discussed in this review.

Inositol phosphate regulation of variant surface glycoprotein expression in Trypanosoma brucei (232.3)

Trypanosoma brucei evades the host immune system by antigenic variation, expressing only one of >1,000 variant surface glycoprotein (VSG) genes through transcription at one of ~20 subtelomeric expression sites (ESs), thereby producing a homogenous surface coat. These coats vary periodically via VSG gene recombination or transcriptional ES switching, and the mechanisms regulating these processes are unknown. T. brucei has an inositol phosphate pathway encompassing 26 genes for inositol phosphate (IP) and phosphatidylinositol (PI) kinases, phosphatases, and metabolic enzymes. We generated T. brucei bloodstream form (BF) conditional null lines for genes in this pathway and found that PI(4) 5‐kinase and IP3 5‐phosphatase are essential for growth and infection of mice. Conditional repression of PI(4) 5‐kinase or IP3 5‐phosphatase, or overexpression of phospholipase C, resulted in expression of VSGs from multiple ESs, including BF and metacyclic ESs. Expression of procyclins or other genes encoded by RNA polymerase I, II or III were not affected, and repression of other IP pathway genes did not induce VSG expression. Repression of PI(4)5‐kinase or IP3(5)‐phosphatase also caused telomere repositioning in the nucleus; re‐expression of these enzymes generated cells that switched VSGs. Thus, T. brucei possesses an IP pathway which appears to be involved in regulation of VSG expression and antigenic switching.

Spleen-Dependent Regulation of Antigenic Variation in Malaria Parasites: Plasmodium knowlesi SICAvar Expression Profiles in Splenic and Asplenic Hosts

BackgroundAntigenic variation by malaria parasites was first described in Plasmodium knowlesi, which infects humans and macaque monkeys, and subsequently in P. falciparum, the most virulent human parasite. The schizont-infected cell agglutination (SICA) variant proteins encoded by the SICAvar multigene family in P. knowlesi, and Erythrocyte Membrane Protein-1 (EMP-1) antigens encoded by the var multigene family in P. falciparum, are expressed at the surface of infected erythrocytes, are associated with virulence, and serve as determinants of naturally acquired immunity. A parental P. knowlesi clone, Pk1(A+), and a related progeny clone, Pk1(B+)1+, derived by an in vivo induced variant antigen switch, were defined by the expression of distinct SICA variant protein doublets of 210/190 and 205/200 kDa, respectively. Passage of SICA[+] infected erythrocytes through splenectomized rhesus monkeys results in the SICA[-] phenotype, defined by the lack of surface expression and agglutination with variant specific antisera.Principal FindingsWe have investigated SICAvar RNA and protein expression in Pk1(A+), Pk1(B+)1+, and SICA[-] parasites. The Pk1(A+) and Pk1(B+)1+ parasites express different distinct SICAvar transcript and protein repertoires. By comparison, SICA[-] parasites are characterized by a vast reduction in SICAvar RNA expression, the lack of full-length SICAvar transcript signals on northern blots, and correspondingly, the absence of any SICA protein detected by mass spectrometry.SignificanceSICA protein expression may be under transcriptional as well as post-transcriptional control, and we show for the first time that the spleen, an organ central to blood-stage immunity in malaria, exerts an influence on these processes. Furthermore, proteomics has enabled the first in-depth characterization of SICA[+] protein phenotypes and we show that the in vivo switch from Pk1(A+) to Pk1(B+)1+ parasites resulted in a complete change in SICA profiles. These results emphasize the importance of studying antigenic variation in the context of the host environment.

The antigenic switching network of Plasmodium falciparum and its implications for the immuno-epidemiology of malaria

Antigenic variation in the human malaria parasite Plasmodium falciparum involves sequential and mutually exclusive expression of members of the var multi-gene family and appears to follow a non-random pattern. In this study, using a detailed in vitro gene transcription analysis of the culture-adapted HB3 strain of P. falciparum, we show that antigenic switching is governed by a global activation hierarchy favouring short and highly diverse genes in central chromosomal location. Longer and more conserved genes, which have previously been associated with severe infection in immunologically naive hosts, are rarely activated, however, implying an in vivo fitness advantage possibly through adhesion-dependent survival rates. We further show that a gene's activation rate is positively associated sequence diversity, which could offer important new insights into the evolution and maintenance of antigenic diversity in P. falciparum malaria. DOI:http://dx.doi.org/10.7554/eLife.01074.001.

Characterization of Free Exopolysaccharides Secreted by Mycoplasma mycoides Subsp. mycoides

Contagious bovine pleuropneumonia is a severe respiratory disease of cattle that is caused by a bacterium of the Mycoplasma genus, namely Mycoplasma mycoides subsp. mycoides (Mmm). In the absence of classical virulence determinants, the pathogenicity of Mmm is thought to rely on intrinsic metabolic functions and specific components of the outer cell surface. One of these latter, the capsular polysaccharide galactan has been notably demonstrated to play a role in Mmm persistence and dissemination. The free exopolysaccharides (EPS), also produced by Mmm and shown to circulate in the blood stream of infected cattle, have received little attention so far. Indeed, their characterization has been hindered by the presence of polysaccharide contaminants in the complex mycoplasma culture medium. In this study, we developed a method to produce large quantities of EPS by transfer of mycoplasma cells from their complex broth to a chemically defined medium and subsequent purification. NMR analyses revealed that the purified, free EPS had an identical β(1−>6)-galactofuranosyl structure to that of capsular galactan. We then analyzed intraclonal Mmm variants that produce opaque/translucent colonies on agar. First, we demonstrated that colony opacity was related to the production of a capsule, as observed by electron microscopy. We then compared the EPS extracts and showed that the non-capsulated, translucent colony variants produced higher amounts of free EPS than the capsulated, opaque colony variants. This phenotypic variation was associated with an antigenic variation of a specific glucose phosphotransferase permease. Finally, we conducted in silico analyses of candidate polysaccharide biosynthetic pathways in order to decipher the potential link between glucose phosphotransferase permease activity and attachment/release of galactan. The co-existence of variants producing alternative forms of galactan (capsular versus free extracellular galactan) and associated with an antigenic switch constitutes a finely tuned mechanism that may be involved in virulence.

DNA Break Site at Fragile Subtelomeres Determines Probability and Mechanism of Antigenic Variation in African Trypanosomes

Antigenic variation in African trypanosomes requires monoallelic transcription and switching of variant surface glycoprotein (VSG) genes. The transcribed VSG, always flanked by ‘70 bp’-repeats and telomeric-repeats, is either replaced through DNA double-strand break (DSB) repair or transcriptionally inactivated. However, little is known about the subtelomeric DSBs that naturally trigger antigenic variation in Trypanosoma brucei, the subsequent DNA damage responses, or how these responses determine the mechanism of VSG switching. We found that DSBs naturally accumulate close to both transcribed and non-transcribed telomeres. We then induced high-efficiency meganuclease-mediated DSBs and monitored DSB-responses and DSB-survivors. By inducing breaks at distinct sites within both transcribed and silent VSG transcription units and assessing local DNA resection, histone modification, G2/M-checkpoint activation, and both RAD51-dependent and independent repair, we reveal how breaks at different sites trigger distinct responses and, in ‘active-site’ survivors, different switching mechanisms. At the active site, we find that promoter-adjacent breaks typically failed to trigger switching, 70 bp-repeat-adjacent breaks almost always triggered switching through 70 bp-repeat recombination (∼60% RAD51-dependent), and telomere-repeat-adjacent breaks triggered switching through loss of the VSG expression site (25% of survivors). Expression site loss was associated with G2/M-checkpoint bypass, while 70 bp-repeat-recombination was associated with DNA-resection, γH2A-focus assembly and a G2/M-checkpoint. Thus, the probability and mechanism of antigenic switching are highly dependent upon the location of the break. We conclude that 70 bp-repeat-adjacent and telomere-repeat-adjacent breaks trigger distinct checkpoint responses and VSG switching pathways. Our results show how subtelomere fragility can generate the triggers for the major antigenic variation mechanisms in the African trypanosome.

Linking the antigen archive structure to pathogen fitness in African trypanosomes.

Systems that generate antigenic variation enable pathogens to evade host immune responses and are intricately interwoven with major pathogen traits, such as host choice, growth, virulence and transmission. Although much is understood about antigen switching at the molecular level, little is known about the cross-scale links between these molecular processes and the larger-scale within and between host population dynamics that they must ultimately drive. Inspired by the antigenic variation system of African trypanosomes, we apply modelling approaches to our expanding understanding of the organization and expression of antigen repertoires, and explore links across these scales. We predict how pathogen population processes are determined by underlying molecular genetics and infer resulting selective pressures on important emergent repertoire traits.

Antigenic Switching of Hepatitis B Virus by Alternative Dimerization of the Capsid Protein

Chronic hepatitis B virus (HBV) infection afflicts millions worldwide with cirrhosis and liver cancer. HBV e-antigen (HBeAg), a clinical marker for disease severity, is a nonparticulate variant of the protein (core antigen, HBcAg) that forms the building-blocks of capsids. HBeAg is not required for virion production, but is implicated in establishing immune tolerance and chronic infection. Here, we report the crystal structure of HBeAg, which clarifies how the short N-terminal propeptide of HBeAg induces a radically altered mode of dimerization relative to HBcAg (∼140° rotation), locked into place through formation of intramolecular disulfide bridges. This structural switch precludes capsid assembly and engenders a distinct antigenic repertoire, explaining why the two antigens are cross-reactive at the Tcell level (through sequence identity) but not at the B cell level (through conformation). The structure offers insight into how HBeAg may establish immune tolerance for HBcAg while evading its robust immunogenicity.

Insulator-like pairing elements regulate silencing and mutually exclusive expression in the malaria parasite Plasmodium falciparum

Plasmodium falciparum causes the deadliest form of human malaria. Its virulence is attributed to its ability to modify the infected RBC and to evade human immune attack through antigenic variation. Antigenic variation is achieved through tight regulation of antigenic switches between variable surface antigens named "P. falciparum erythrocyte membrane protein-1" encoded by the var multicopy gene family. Individual parasites express only a single var gene at a time, maintaining the remaining var genes in a transcriptionally silent state. Strict pairing between var gene promoters and a second promoter within an intron found in each var gene is required for silencing and counting of var genes by the mechanism that controls mutually exclusive expression. We have identified and characterized insulator-like DNA elements that are required for pairing var promoters and introns and thus are essential for regulating silencing and mutually exclusive expression. These elements, found in the regulatory regions of each var gene, are bound by distinct nuclear protein complexes. Any alteration in the specific, paired structure of these elements by either deletion or insertion of additional elements results in an unregulated var gene. We propose a model by which silencing and mutually exclusive expression of var genes is regulated by the precise arrangement of insulator-like DNA pairing elements.

A Statistically Rigorous Method for Determining Antigenic Switching Networks

Many vector-borne pathogens rely on antigenic variation to prolong infections and increase their likelihood of onward transmission. This immune evasion strategy often involves mutually exclusive switching between members of gene families that encode functionally similar but antigenically different variants during the course of a single infection. Studies of different pathogens have suggested that switching between variant genes is non-random and that genes have intrinsic probabilities of being activated or silenced. These factors could create a hierarchy of gene expression with important implications for both infection dynamics and the acquisition of protective immunity. Inferring complete switching networks from gene transcription data is problematic, however, because of the high dimensionality of the system and uncertainty in the data. Here we present a statistically rigorous method for analysing temporal gene transcription data to reconstruct an underlying switching network. Using artificially generated transcription profiles together with in vitro var gene transcript data from two Plasmodium falciparum laboratory strains, we show that instead of relying on data from long-term parasite cultures, accuracy can be greatly improved by using transcription time courses of several parasite populations from the same isolate, each starting with different variant distributions. The method further provides explicit indications about the reliability of the resulting networks and can thus be used to test competing hypotheses with regards to the underlying switching pathways. Our results demonstrate that antigenic switch pathways can be determined reliably from short gene transcription profiles assessing multiple time points, even when subject to moderate levels of experimental error. This should yield important new information about switching patterns in antigenically variable organisms and might help to shed light on the molecular basis of antigenic variation.