Shared themes of antigenic variation and virulence in bacterial, protozoal, and fungal infections

An immune memory-structured SIS epidemiological model for hyperdiverse pathogens

A hyperdiverse class of pathogens of humans and wildlife, including the malaria parasite Plasmodium falciparum, relies on multigene families to encode antigenic variation. As a result, high (asymptomatic) prevalence is observed despite high immunity in local populations under high-transmission settings. The vast diversity of "strains" and genes encoding this variation challenges the application of established models for the population dynamics of such infectious diseases. Agent-based models have been formulated to address theory on strain coexistence and structure, but their complexity can limit application to gain insights into population dynamics. Motivated by P. falciparum malaria, we develop an alternative formulation in the form of a structured susceptible-infected-susceptible population model in continuous time, where individuals are classified not only by age, as is standard, but also by the diversity of parasites they have been exposed to and retain in their specific immune memory. We analyze the population dynamics and bifurcation structure of this system of partial-differential equations, showing the existence of alternative steady states and an associated tipping point with transmission intensity. We attribute the critical transition to the positive feedback between parasite genetic diversity and force of infection. Basins of attraction show that intervention must drastically reduce diversity to prevent a rebound to high infection levels. Results emphasize the importance of explicitly considering pathogen diversity and associated specific immune memory in the population dynamics of hyperdiverse epidemiological systems. This statement is discussed in a more general context for ecological competition systems with hyperdiverse trait spaces.

Giardia VSPAS7 protein attenuates Giardia intestinalis-induced host macrophage pyroptosis

BACKGROUND

The unicellular protozoan parasite Giardia intestinalis, which primarily infects humans and animals such as cattle and sheep, is having a major negative impact on public health. Giardia is able to evade the recognition and elimination of the host immune system because of the trophozoite surface and extracellular vesicles (EVs) covered by variant-specific surface proteins (VSPs). As key proteins for immune evasion, whether VSPs can regulate Giardia-induced pyroptosis and promote Giardia evasion of host immune responses has not been reported.

METHODS

To examine the role of Giardia VSPAS7 on Giardia-induced activation of the signaling pathway, secretion of pro-inflammatory cytokines, pyroptosis and the mechanism involved, we constructed the pcDNA3.1-vspas7 expression plasmid and transfected this plasmid into mouse macrophages. Key proteins for pyroptosis, IL-1β secretion and LDH release were detected in pcDNA3.1-vspas7-transfected wild-type (WT) cells and NLRP3-deficient cells by western blot, ELISA and LDH assays, respectively. The interactions of Giardia VSPAS7 and mouse NLRP3 were examined using immunofluorescence assays (IFA), co-immunoprecipitation (Co-IP) and bimolecular fluorescence complementation (BiFC) assays.

RESULTS

VSPAS7 could decrease the levels of phosphorylated-p65 (P-p65), P-IκBα and P-ERK caused by Giardia and reduce the production levels of Giardia-induced pro-inflammatory cytokine IL-6, IL-12 p40 and TNF-α. The results showed that VSPAS7 inhibited Giardia-mediated activation of NF-κB, ERK/MAPK signaling and secretion of pro-inflammatory cytokines. Furthermore, VSPAS7 suppressed Giardia-induced macrophage pyroptosis by reducing GSDMD cleavage, caspase-1 activation, IL-1β secretion and LDH release. We further found that VSPAS7 could interact with mouse NLRP3 directly, and in NLRP3-deficient cells the suppression of Giardia-induced macrophage pyroptosis by VSPAS7 was significantly attenuated.

CONCLUSIONS

Overall, VSPAS7 could inhibit Giardia-induced activation of signaling pathways and pyroptosis in host macrophages, allowing Giardia evasion of host immune responses. Studies on Giardia VSP-mediated immune evasion provide an important theoretical basis for in-depth studies on Giardia pathogenicity.

Identification of G-quadruplex DNA sequences in SARS-CoV2

G-quadruplex structure or Putative Quadruplex Sequences (PQSs) are abundant in human, microbial, DNA, or RNA viral genomes. These sequences in RNA viral genome play critical roles in integration into human genome as LTR (Long Terminal Repeat), genome replication, chromatin rearrangements, gene regulation, antigen variation (Av), and virulence. Here, we investigated whether the genome of SARS-CoV2, an RNA virus, contained such potential G-quadruplex structures. Using bioinformatic tools, we searched for such sequences and found thirty-seven (forward strand (twenty-five) + reverse strand (Twelve)) QGRSs (Quadruplex forming G-Rich Sequences)/PQSs in SARS-CoV2 genome. These sequences are dispersed mainly in the upstream of SARS-CoV2 genes. We discuss whether existing PQS/QGRS ligands could inhibit the SARS-CoV2 replication and gene transcription as has been observed in other RNA viruses. Further experimental validation would determine the role of these G-quadruplex sequences in SARS-CoV2 genome function to survive in the host cells and identify therapeutic agents to destabilize these PQSs/QGRSs.

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.

An antigenic diversification threshold for falciparum malaria transmission at high endemicity

In malaria and several other important infectious diseases, high prevalence occurs concomitantly with incomplete immunity. This apparent paradox poses major challenges to malaria elimination in highly endemic regions, where asymptomatic Plasmodium falciparum infections are present across all age classes creating a large reservoir that maintains transmission. This reservoir is in turn enabled by extreme antigenic diversity of the parasite and turnover of new variants. We present here the concept of a threshold in local pathogen diversification that defines a sharp transition in transmission intensity below which new antigen-encoding genes generated by either recombination or migration cannot establish. Transmission still occurs below this threshold, but diversity of these genes can neither accumulate nor recover from interventions that further reduce it. An analytical expectation for this threshold is derived and compared to numerical results from a stochastic individual-based model of malaria transmission that incorporates the major antigen-encoding multigene family known as var. This threshold corresponds to an “innovation” number we call R_div; it is different from, and complementary to, the one defined by the classic basic reproductive number of infectious diseases, R₀, which does not readily is better apply under large and dynamic strain diversity. This new threshold concept can be exploited for effective malaria control and applied more broadly to other pathogens with large multilocus antigenic diversity.

The vast diversity of the falciparum malaria parasite, as seen by the immune system of hosts in high transmission regions, underlies both high prevalence of asymptomatic infections and partial protection to re-infection despite previous exposure. This large antigenic diversity of the parasite challenges control and elimination efforts. We propose a threshold quantity for antigenic innovation, we call R_div, measuring the potential of transmission to accumulate new antigenic variants over time. When R_div is pushed below one by reduced transmission intensity, new genes encoding this variation can no longer accumulate, resulting in a lower number of strains and facilitating further intervention. This innovation number can be applied to other infectious diseases with fast turnover of antigens, where large standing diversity similarly opposes successful intervention.

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.

Insights into the genetic variation profile of tprK in Treponema pallidum during the development of natural human syphilis infection

Background

Although the tprK gene of Treponema pallidum are thought to play a critical role in the pathogenesis of syphilis, the profile of variations in tprK during the development of human syphilis infection have remained unclear.

Methods/Principal findings

Through next-generation sequencing, we compared the tprK gene of 14 secondary syphilis patients with that of 14 primary syphilis patients, and the results showed an increased number of variants within the seven V regions of the tprK gene in the secondary syphilis samples. The length of the sequences within each V region also presented a 3-bp changing pattern. Interestingly, the frequencies of predominant sequences within the V regions in the secondary syphilis samples were generally decreased compared with those found in the primary syphilis samples, particularly in the V7 region, where a frequency below 60% was found in up to 57% (8/14) of all secondary samples compared with 7% (1/14) of all primary samples. Moreover, the number of minor variants distributed between frequencies of 10 and 49.9% was increased. The alignment of all amino acid sequences within each V region of the primary and secondary syphilis samples revealed that some amino acid sequences, particularly the amino acid sequences IASDGGAIKH and IASEDGSAGNLKH in V1, were highly stable. Additionally, the amino acid sequences in V6 also exhibited notable intrastrain heterogeneity and were likely to form a strain-specific pattern at the interstrain level.

Conclusions

The identification of different profiles of the tprK gene in primary and secondary syphilis patients indicated that the tprK gene of T. pallidum undergoes constant variation to result in the best adaptation to the host. The highly stable peptides found in V1 are likely promising potential vaccine components. The highly heterogenetic regions (e.g., V6) could help to understand the role of tprK in immune evasion.

Antigenic variation of the TprK antigen has been acknowledged to explain the persistence of Treponema pallidum in the host, however, the profile of variations in tprK during the development of human syphilis infection has not been well characterized. Here, we performed next-generation sequencing to compare the variations in tprK between primary and secondary syphilis samples. The profiles of tprK in the samples at different stages showed differences. A higher amount of pool variants within seven V regions was found in the secondary syphilis samples, and the frequencies of their predominant sequences generally decreased with increases in the number of minor variants with frequencies in the range of 10 to 49.9%. However, the length of variable sequences within the V regions of tprK in the secondary syphilis samples also presented a 3-bp changing pattern. Notably, the amino acid sequences IASDGGAIKH and IASEDGSAGNLKH in V1 not only presented a high proportion of interstrain sharing but also were found at a relatively high frequency (above 80%) in the populations. The sequences in V6 of the samples demonstrated substantial variability at the intra- and interstrain levels. These findings could provide insights into the potential syphilis vaccine components and the role of TprK in immune evasion.

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.

G-Quadruplexes: More Than Just a Kink in Microbial Genomes

G-quadruplexes (G4s) are noncanonical nucleic acid secondary structures formed by guanine-rich DNA and RNA sequences. In this review we aim to provide an overview of the biological roles of G4s in microbial genomes with emphasis on recent discoveries. G4s are enriched and conserved in the regulatory regions of microbes, including bacteria, fungi, and viruses. Importantly, G4s in hepatitis B virus (HBV) and hepatitis C virus (HCV) genomes modulate genes crucial for virus replication. Recent studies on Epstein-Barr virus (EBV) shed light on the role of G4s within the microbial transcripts as cis-acting regulatory signals that modulate translation and facilitate immune evasion. Furthermore, G4s in microbial genomes have been linked to radioresistance, antigenic variation, recombination, and latency. G4s in microbial genomes represent novel therapeutic targets for antimicrobial therapy.

Interactions between host immune response and antigenic variation that control Borrelia burgdorferi population dynamics

The population dynamics of pathogens within hosts result from interactions between host immune responses and mechanisms of the pathogen to evade or resist immune responses. Vertebrate hosts have evolved adaptive immune responses to eliminate the infection, while many pathogens evade immune clearance through altering surface antigens. Such interactions can result in a characteristic pattern of pathogen population dynamics within hosts consisting of population growth after infection, rapid population decline following specific immune responses, followed by persistence at low densities during a chronic infection stage. Despite the medical importance of chronic infections, little is known about the conditions of the interactions between variable antigens and the adaptive immune system that cause the characteristic pathogen population dynamics. Using the vls antigenic variation system of the Lyme disease pathogen, Borrelia burgdorferi, as a model system, we investigated conditions of the interaction between the antigenic variation system and the adaptive immune response that can explain the within-host population dynamics of B. burgdorferi using mathematical modelling. This characteristic population dynamic pattern can be explained by models that assume a variable immune removal rate of antibody-bound B. burgdorferi. However, models with a constant immune removal rate could reproduce the rapid population decline of B. burgdorferi populations but not their long-term persistence within hosts using parameter values determined by fitting empirical data. The model predictions, along with the assumptions about the interactions between B. burgdorferi and the immune response, can be tested experimentally to estimate the likelihood that each mechanism affects B. burgdorferi population dynamics in real infections.

Correlation between antigenicity and variability in the vls antigenic variation system of Borrelia burgdorferi

Many parasites have evolved antigenic variation systems that alter surface proteins in order to evade recognition by presently expressed antibodies and subsequent death. Although the amino acid positions in antigens to which antibodies most commonly target are expected to be the most variable, this assumption has not been investigated. Using the vls antigenic variation system of Borrelia burgdorferi as a model, we first investigated this assumption computationally and then developed a sensitive immunoassay to experimentally validate the computational results. There was a strong correlation between variability at an amino acid position and each of the computational metrics associated with antibody reactivity. However, empirical measures of antibody reactivity were not consistently greater at the variable amino acid positions than at the invariant amino acid positions. The inconsistent experimental support for this hypothesis suggests that the biological effect of variability at an amino acid position is obfuscated by other factors.

DNA Recombination Strategies During Antigenic Variation in the African Trypanosome

Survival of the African trypanosome in its mammalian hosts has led to the evolution of antigenic variation, a process for evasion of adaptive immunity that has independently evolved in many other viral, bacterial and eukaryotic pathogens. The essential features of trypanosome antigenic variation have been understood for many years and comprise a dense, protective Variant Surface Glycoprotein (VSG) coat, which can be changed by recombination-based and transcription-based processes that focus on telomeric VSG gene transcription sites. However, it is only recently that the scale of this process has been truly appreciated. Genome sequencing of Trypanosoma brucei has revealed a massive archive of >1000 VSG genes, the huge majority of which are functionally impaired but are used to generate far greater numbers of VSG coats through segmental gene conversion. This chapter will discuss the implications of such VSG diversity for immune evasion by antigenic variation, and will consider how this expressed diversity can arise, drawing on a growing body of work that has begun to examine the proteins and sequences through which VSG switching is catalyzed. Most studies of trypanosome antigenic variation have focused on T. brucei, the causative agent of human sleeping sickness. Other work has begun to look at antigenic variation in animal-infective trypanosomes, and we will compare the findings that are emerging, as well as consider how antigenic variation relates to the dynamics of host-trypanosome interaction.

The Global Transcription Factor Lrp Controls Virulence Modulation in Xenorhabdus nematophila

The bacterium Xenorhabdus nematophila engages in phenotypic variation with respect to pathogenicity against insect larvae, yielding both virulent and attenuated subpopulations of cells from an isogenic culture. The global regulatory protein Lrp is necessary for X. nematophila virulence and immunosuppression in insects, as well as colonization of the mutualistic host nematode Steinernema carpocapsae, and mediates expression of numerous genes implicated in each of these phenotypes. Given the central role of Lrp in X. nematophila host associations, as well as its involvement in regulating phenotypic variation pathways in other bacteria, we assessed its function in virulence modulation. We discovered that expression of lrp varies within an isogenic population, in a manner that correlates with modulation of virulence. Unexpectedly, although Lrp is necessary for optimal virulence and immunosuppression, cells expressing high levels of lrp were attenuated in these processes relative to those with low to intermediate lrp expression. Furthermore, fixed expression of lrp at high and low levels resulted in attenuated and normal virulence and immunosuppression, respectively, and eliminated population variability of these phenotypes. These data suggest that fluctuating lrp expression levels are sufficient to drive phenotypic variation in X. nematophila.

IMPORTANCE

Many bacteria use cell-to-cell phenotypic variation, characterized by distinct phenotypic subpopulations within an isogenic population, to cope with environmental change. Pathogenic bacteria utilize this strategy to vary antigen or virulence factor expression. Our work establishes that the global transcription factor Lrp regulates phenotypic variation in the insect pathogen Xenorhabdus nematophila, leading to attenuation of virulence and immunosuppression in insect hosts. Unexpectedly, we found an inverse correlation between Lrp expression levels and virulence: high levels of expression of Lrp-dependent putative virulence genes are detrimental for virulence but may have an adaptive advantage in other aspects of the life cycle. Investigation of X. nematophila phenotypic variation facilitates dissection of this phenomenon in the context of a naturally occurring symbiosis.

Transcriptome sequencing and analysis of Plasmodium gallinaceum reveals polymorphisms and selection on the apical membrane antigen-1

Background

Plasmodium erythrocyte invasion genes play a key role in malaria parasite transmission, host-specificity and immuno-evasion. However, the evolution of the genes responsible remains understudied. Investigating these genes in avian malaria parasites, where diversity is particularly high, offers new insights into the processes that confer malaria pathogenesis. These parasites can pose a significant threat to birds and since birds play crucial ecological roles they serve as important models for disease dynamics. Comprehensive knowledge of the genetic factors involved in avian malaria parasite invasion is lacking and has been hampered by difficulties in obtaining nuclear data from avian malaria parasites. Thus the first Illumina-based de novo transcriptome sequencing and analysis of the chicken parasite Plasmodium gallinaceum was performed to assess the evolution of essential Plasmodium genes.

Methods

White leghorn chickens were inoculated intravenously with erythrocytes containing P. gallinaceum. cDNA libraries were prepared from RNA extracts collected from infected chick blood and sequencing was run on the HiSeq2000 platform. Orthologues identified by transcriptome sequencing were characterized using phylogenetic, ab initio protein modelling and comparative and population-based methods.

Results

Analysis of the transcriptome identified several orthologues required for intra-erythrocytic survival and erythrocyte invasion, including the rhoptry neck protein 2 (RON2) and the apical membrane antigen-1 (AMA-1). Ama-1 of avian malaria parasites exhibits high levels of genetic diversity and evolves under positive diversifying selection, ostensibly due to protective host immune responses.

Conclusion

Erythrocyte invasion by Plasmodium parasites require AMA-1 and RON2 interactions. AMA-1 and RON2 of P. gallinaceum are evolutionarily and structurally conserved, suggesting that these proteins may play essential roles for avian malaria parasites to invade host erythrocytes. In addition, host-driven selection presumably results in the high levels of genetic variation found in ama-1 of avian Plasmodium species. These findings have implications for investigating avian malaria epidemiology and population dynamics. Moreover, this work highlights the P. gallinaceum transcriptome as an important public resource for investigating the diversity and evolution of essential Plasmodium genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-382) contains supplementary material, which is available to authorized users.

Hidden diversity in honey bee gut symbionts detected by single-cell genomics

Microbial communities in animal guts are composed of diverse, specialized bacterial species, but little is known about how gut bacteria diversify to produce genetically and ecologically distinct entities. The gut microbiota of the honey bee, Apis mellifera, presents a useful model, because it consists of a small number of characteristic bacterial species, each showing signs of diversification. Here, we used single-cell genomics to study the variation within two species of the bee gut microbiota: Gilliamella apicola and Snodgrassella alvi. For both species, our analyses revealed extensive variation in intraspecific divergence of protein-coding genes but uniformly high levels of 16S rRNA similarity. In both species, the divergence of 16S rRNA loci appears to have been curtailed by frequent recombination within populations, while other genomic regions have continuously diverged. Furthermore, gene repertoires differ markedly among strains in both species, implying distinct metabolic capabilities. Our results show that, despite minimal divergence at 16S rRNA genes, in situ diversification occurs within gut communities and generates bacterial lineages with distinct ecological niches. Therefore, important dimensions of microbial diversity are not evident from analyses of 16S rRNA, and single cell genomics has potential to elucidate processes of bacterial diversification.

Gut microbial communities are often complex, consisting of bacteria from divergent phyla as well as multiple strains of each of the constituent species. But because the composition of these communities is typically assessed using 16S rRNA analyses, little is known about genomic changes associated with in situ diversification of bacterial lineages in animal guts. We undertook a single-cell genomic approach to investigate the diversification within two species of the gut microbiota of honey bees. Each species exhibited a surprisingly high level of genomic variation, despite uniformity in the 16S rRNA sequences. Our data indicate that genetically and ecologically distinct lineages can evolve in the gut of the same host species in the presence of frequent recombination at 16S rRNA genes. These findings parallel observations from mammals, suggesting that in situ diversification of a few bacterial lineages is a common pattern in the evolution of gut communities.

Genome evolution in trypanosomatid parasites

A decade of genome sequencing has transformed our understanding of how trypanosomatid parasites have evolved and provided fresh impetus to explaining the origins of parasitism in the Kinetoplastida. In this review, I will consider the many ways in which genome sequences have influenced our view of genomic reduction in trypanosomatids; how species-specific genes, and the genomic domains they occupy, have illuminated the innovations in trypanosomatid genomes; and how comparative genomics has exposed the molecular mechanisms responsible for innovation and adaptation to a parasitic lifestyle.

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

BACKGROUND

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

RESULTS

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409-1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

CONCLUSIONS

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

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.

Natural selection promotes antigenic evolvability

The hypothesis that evolvability - the capacity to evolve by natural selection - is itself the object of natural selection is highly intriguing but remains controversial due in large part to a paucity of direct experimental evidence. The antigenic variation mechanisms of microbial pathogens provide an experimentally tractable system to test whether natural selection has favored mechanisms that increase evolvability. Many antigenic variation systems consist of paralogous unexpressed 'cassettes' that recombine into an expression site to rapidly alter the expressed protein. Importantly, the magnitude of antigenic change is a function of the genetic diversity among the unexpressed cassettes. Thus, evidence that selection favors among-cassette diversity is direct evidence that natural selection promotes antigenic evolvability. We used the Lyme disease bacterium, Borrelia burgdorferi, as a model to test the prediction that natural selection favors amino acid diversity among unexpressed vls cassettes and thereby promotes evolvability in a primary surface antigen, VlsE. The hypothesis that diversity among vls cassettes is favored by natural selection was supported in each B. burgdorferi strain analyzed using both classical (dN/dS ratios) and Bayesian population genetic analyses of genetic sequence data. This hypothesis was also supported by the conservation of highly mutable tandem-repeat structures across B. burgdorferi strains despite a near complete absence of sequence conservation. Diversification among vls cassettes due to natural selection and mutable repeat structures promotes long-term antigenic evolvability of VlsE. These findings provide a direct demonstration that molecular mechanisms that enhance evolvability of surface antigens are an evolutionary adaptation. The molecular evolutionary processes identified here can serve as a model for the evolution of antigenic evolvability in many pathogens which utilize similar strategies to establish chronic infections.

Cryptosporidium pathogenicity and virulence

Cryptosporidium is a protozoan parasite of medical and veterinary importance that causes gastroenteritis in a variety of vertebrate hosts. Several studies have reported different degrees of pathogenicity and virulence among Cryptosporidium species and isolates of the same species as well as evidence of variation in host susceptibility to infection. The identification and validation of Cryptosporidium virulence factors have been hindered by the renowned difficulties pertaining to the in vitro culture and genetic manipulation of this parasite. Nevertheless, substantial progress has been made in identifying putative virulence factors for Cryptosporidium. This progress has been accelerated since the publication of the Cryptosporidium parvum and C. hominis genomes, with the characterization of over 25 putative virulence factors identified by using a variety of immunological and molecular techniques and which are proposed to be involved in aspects of host-pathogen interactions from adhesion and locomotion to invasion and proliferation. Progress has also been made in the contribution of host factors that are associated with variations in both the severity and risk of infection. Here we provide a review comprised of the current state of knowledge on Cryptosporidium infectivity, pathogenesis, and transmissibility in light of our contemporary understanding of microbial virulence.

Duplication and concerted evolution in a master sex determiner under balancing selection

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

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.

Telomere length affects the frequency and mechanism of antigenic variation in Trypanosoma brucei

Trypanosoma brucei is a master of antigenic variation and immune response evasion. Utilizing a genomic repertoire of more than 1000 Variant Surface Glycoprotein-encoding genes (VSGs), T. brucei can change its protein coat by “switching” from the expression of one VSG to another. Each active VSG is monoallelically expressed from only one of approximately 15 subtelomeric sites. Switching VSG expression occurs by three predominant mechanisms, arguably the most significant of which is the non-reciprocal exchange of VSG containing DNA by duplicative gene conversion (GC). How T. brucei orchestrates its complex switching mechanisms remains to be elucidated. Recent work has demonstrated that an exogenous DNA break in the active site could initiate a GC based switch, yet the source of the switch-initiating DNA lesion under natural conditions is still unknown. Here we investigated the hypothesis that telomere length directly affects VSG switching. We demonstrate that telomerase deficient strains with short telomeres switch more frequently than genetically identical strains with long telomeres and that, when the telomere is short, switching preferentially occurs by GC. Our data supports the hypothesis that a short telomere at the active VSG expression site results in an increase in subtelomeric DNA breaks, which can initiate GC based switching. In addition to their significance for T. brucei and telomere biology, the findings presented here have implications for the many diverse pathogens that organize their antigenic genes in subtelomeric regions.

A broad array of human pathogens (including bacteria, fungi and parasites) vary the proteins on their cell surface to escape the immune response of their hosts. This process, called antigenic variation, relies on a repertoire of variant protein encoding genes in the genome and the organism's ability to accurately switch from the expression of one variant gene to another. A common theme in both the diversification of these variant genes and the mechanisms required for their expression is that they are often located near the ends of chromosomes. The ends of chromosomes are protected by structures called telomeres. Regions near the telomere are referred to as subtelomeric and are commonly thought to be comparatively unstable DNA sites. It is therefore intriguing that organisms that rely on antigenic variation for survival would organize their critical survival genes in these sites. Trypanosoma brucei is a model organism for the study of antigenic variation. The causative agent of African sleeping sickness, this unicellular parasite possesses an antigenic repertoire of unparalleled diversity, which can only be expressed from specific subtelomeric sites. Here we use the power of the T. brucei model to investigate the effect of telomere length on antigenic variation.

Impact of mosquito bites on asexual parasite density and gametocyte prevalence in asymptomatic chronic Plasmodium falciparum infections and correlation with IgE and IgG titers

An immunomodulatory role of arthropod saliva has been well documented, but evidence for an effect on Plasmodium sp. infectiousness remains controversial. Mosquito saliva may orient the immune response toward a Th2 profile, thereby priming a Th2 response against subsequent antigens, including Plasmodium. Orientation toward a Th1 versus a Th2 profile promotes IgG and IgE proliferation, respectively, where the former is crucial for the development of an efficient antiparasite immune response. Here we assessed the direct effect of mosquito bites on the density of Plasmodium falciparum asexual parasites and the prevalence of gametocytes in chronic, asymptomatic infections in a longitudinal cohort study of seasonal transmission. We additionally correlated these parasitological measures with IgE and IgG antiparasite and anti-salivary gland extract titers. The mosquito biting density was positively correlated with the asexual parasite density but not asexual parasite prevalence and was negatively correlated with gametocyte prevalence. Individual anti-salivary gland IgE titers were also negatively correlated with gametocyte carriage and were strongly positively correlated with antiparasite IgE titers, consistent with the hypothesis that mosquito bites predispose individuals to develop an IgE antiparasite response. We provide evidence that mosquito bites have an impact on asymptomatic infections and differentially so for the production of asexual and sexual parasites. An increased research focus on the immunological impact of mosquito bites during asymptomatic infections is warranted, to establish whether strategies targeting the immune response to saliva can reduce the duration of infection and the onward transmission of the parasite.

After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selection

Many of the most virulent bacterial pathogens show low genetic diversity and sexual isolation. Accordingly, Mycobacterium tuberculosis, the deadliest human pathogen, is thought to be clonal and evolve by genetic drift. Yet, its genome shows few of the concomitant signs of genome degradation. We analyzed 24 genomes and found an excess of genetic diversity in regions encoding key adaptive functions including the type VII secretion system and the ancient horizontally transferred virulence-related regions. Four different approaches showed evident signs of recombination in M. tuberculosis. Recombination tracts add a high density of polymorphisms, and many are thus predicted to arise from outside the clade. Some of these tracts match Mycobacterium canettii sequences. Recombination introduced an excess of non-synonymous diversity in general and even more in genes expected to be under positive or diversifying selection, e.g., cell wall component genes. Mutations leading to non-synonymous SNPs are effectively purged in MTBC, which shows dominance of purifying selection. MTBC mutation bias toward AT nucleotides is not compensated by biased gene conversion, suggesting the action of natural selection also on synonymous changes. Together, all of these observations point to a strong imprint of recombination and selection in the genome affecting both non-synonymous and synonymous positions. Hence, contrary to some other pathogens and previous proposals concerning M. tuberculosis, this lineage may have come out of its ancestral bottleneck as a very successful pathogen that is rapidly diversifying by the action of mutation, recombination, and natural selection.

Characterization and tissue-specific expression patterns of the Plasmodium chabaudi cir multigene family

BACKGROUND

Variant antigens expressed on the surface of parasitized red blood cells (pRBCs) are important virulence factors of malaria parasites. Whereas Plasmodium falciparum erythrocyte membrane proteins 1 (PfEMP1) are responsible for sequestration of mature parasites, little is known about putative ligands mediating cytoadherence to host receptors in other Plasmodium species. Candidates include members of the pir superfamily found in the human parasite Plasmodium vivax (vir), in the simian pathogen Plasmodium knowlesi (kir) and in the rodent malarias Plasmodium yoelii (yir), Plasmodium berghei (bir) and Plasmodium chabaudi (cir). The aim of this study was to reveal a potential involvement of cir genes in P. chabaudi sequestration.

METHODS

Subfamilies of cir genes were identified by bioinformatic analyses of annotated sequence data in the Plasmodium Genome Database. In order to examine tissue-specific differences in the expression of cir mRNAs, RT-PCR with subfamily-specific primers was used. In total, 432 cDNA clones derived from six different tissues were sequenced to characterize the transcribed cir gene repertoire. To confirm differences in transcription profiles of cir genes, restriction fragment length polymorphism (RFLP) analyses were performed to compare different host tissues and to identify changes during the course of P. chabaudi infections in immunocompetent mice.

RESULTS

The phylogenetic analysis of annotated P. chabaudi putative CIR proteins identified two major subfamilies. Comparison of transcribed cir genes from six different tissues revealed significant differences in the frequency clones belonging to individual cir gene subgroups were obtained from different tissues. Further hints of difference in the transcription of cir genes in individual tissues were obtained by RFLP. Whereas only minimal changes in the transcription pattern of cir genes could be detected during the developmental cycle of the parasites, switching to expression of other cir genes during the course of an infection was observed around or after peak parasitemia.

CONCLUSIONS

The tissue-specific expression of cir mRNAs found in this study indicates correlation between expression of CIR antigens and distribution of parasites in inner organs. Together with comparable results for other members of the pir superfamily this suggests a role of cir and other pir genes in antigenic variation and sequestration of malaria parasites.

Antigenic variation in Plasmodium falciparum malaria involves a highly structured switching pattern

Many pathogenic bacteria, fungi, and protozoa achieve chronic infection through an immune evasion strategy known as antigenic variation. In the human malaria parasite Plasmodium falciparum, this involves transcriptional switching among members of the var gene family, causing parasites with different antigenic and phenotypic characteristics to appear at different times within a population. Here we use a genome-wide approach to explore this process in vitro within a set of cloned parasite populations. Our analyses reveal a non-random, highly structured switch pathway where an initially dominant transcript switches via a set of switch-intermediates either to a new dominant transcript, or back to the original. We show that this specific pathway can arise through an evolutionary conflict in which the pathogen has to optimise between safeguarding its limited antigenic repertoire and remaining capable of establishing infections in non-naïve individuals. Our results thus demonstrate a crucial role for structured switching during the early phases of infections and provide a unifying theory of antigenic variation in P. falciparum malaria as a balanced process of parasite-intrinsic switching and immune-mediated selection.

The malaria parasite Plasmodium falciparum avoids recognition and clearance by the immune system by sequentially switching between members of the var multi-gene family which encode the immunodominant surface proteins PfEMP1. However, some mechanism must exist to prevent rapid exposure of the pathogen's entire antigenic repertoire as this would quickly terminate the infection. It has previously been shown that the immune system can play an important role in orchestrating the sequential display of variants once an infection is established; however this does not explain how repertoire exhaustion is avoided in the initial phases of infection before an immune response has been established. Here we show that P. falciparum has evolved a highly structured switching pattern to prevent repertoire exhaustion in the early stages of infection without compromising the ability to establish new infections among partially immune individuals.

Towards delineating functions within the fasciola secreted cathepsin l protease family by integrating in vivo based sub-proteomics and phylogenetics

BACKGROUND

fasciola hepatica, along with Fasciola gigantica, is the causative agent of fasciolosis, a foodborne zoonotic disease affecting grazing animals and humans worldwide. Pathology is directly related to the release of parasite proteins that facilitate establishment within the host. The dominant components of these excretory-secretory (ES) products are also the most promising vaccine candidates, the cathepsin L (Cat L) protease family.

METHODOLOGY/PRINCIPAL FINDINGS

the sub-proteome of Cat L proteases from adult F. hepatica ES products derived from in vitro culture and in vivo from ovine host bile were compared by 2-DE. The individual Cat L proteases were identified by tandem mass spectrometry with the support of an in-house translated liver fluke EST database. The study reveals plasticity within the CL1 clade of Cat L proteases; highlighted by the identification of a novel isoform and CL1 sub-clade, resulting in a new Cat L phylogenetic analysis including representatives from other adult Cat L phylogenetic clades. Additionally, for the first time, mass spectrometry was shown to be sufficiently sensitive to reveal single amino acid polymorphisms in a resolved 2-DE protein spot derived from pooled population samples.

CONCLUSIONS/SIGNIFICANCE

we have investigated the sub-proteome at the population level of a vaccine target family using the Cat L proteases from F. hepatica as a case study. We have confirmed that F. hepatica exhibits more plasticity in the expression of the secreted CL1 clade of Cat L proteases at the protein level than previously realised. We recommend that superfamily based vaccine discovery programmes should screen parasite populations from different host populations and, if required, different host species via sub-proteomic assay in order to confirm the relative expression at the protein level prior to the vaccine development phase.

Molecular and evolutionary bases of within-patient genotypic and phenotypic diversity in Escherichia coli extraintestinal infections

Although polymicrobial infections, caused by combinations of viruses, bacteria, fungi and parasites, are being recognised with increasing frequency, little is known about the occurrence of within-species diversity in bacterial infections and the molecular and evolutionary bases of this diversity. We used multiple approaches to study the genomic and phenotypic diversity among 226 Escherichia coli isolates from deep and closed visceral infections occurring in 19 patients. We observed genomic variability among isolates from the same site within 11 patients. This diversity was of two types, as patients were infected either by several distinct E. coli clones (4 patients) or by members of a single clone that exhibit micro-heterogeneity (11 patients); both types of diversity were present in 4 patients. A surprisingly wide continuum of antibiotic resistance, outer membrane permeability, growth rate, stress resistance, red dry and rough morphotype characteristics and virulence properties were present within the isolates of single clones in 8 of the 11 patients showing genomic micro-heterogeneity. Many of the observed phenotypic differences within clones affected the trade-off between self-preservation and nutritional competence (SPANC). We showed in 3 patients that this phenotypic variability was associated with distinct levels of RpoS in co-existing isolates. Genome mutational analysis and global proteomic comparisons in isolates from a patient revealed a star-like relationship of changes amongst clonally diverging isolates. A mathematical model demonstrated that multiple genotypes with distinct RpoS levels can co-exist as a result of the SPANC trade-off. In the cases involving infection by a single clone, we present several lines of evidence to suggest diversification during the infectious process rather than an infection by multiple isolates exhibiting a micro-heterogeneity. Our results suggest that bacteria are subject to trade-offs during an infectious process and that the observed diversity resembled results obtained in experimental evolution studies. Whatever the mechanisms leading to diversity, our results have strong medical implications in terms of the need for more extensive isolate testing before deciding on antibiotic therapies.

We investigated whether an infection is a site of pathogen within-species diversity. Our results indicate that there is indeed extensive diversity during human extraintestinal infections by Escherichia coli. This diversity was of two types, not mutually exclusive, as we found that patients were infected either by several distinct E. coli clones or by members of a single clone that exhibit micro-heterogeneity. The high degree of phenotypic diversity, including antibiotic resistance, suggests that there is no uniform selection pressure leading to a single fitter clone during an infection. We discuss a possible mechanism and a mathematical model that explains these unexpected results. Our data suggest that the evolution of diversity in the course of an infection and in in vitro experimental evolution in the absence of host immune selective pressure may have many parallels. Whatever the mechanisms leading to diversity, our results have strong medical implications in terms of the need for more extensive isolate testing before deciding on antibiotic therapies.

Helicobacter pylori HopE and HopV porins present scarce expression among clinical isolates

AIM

To evaluate how widely Helicobacter pylori (H. pylori) HopE and HopV porins are expressed among Chilean isolates and how seroprevalent they are among infected patients in Chile.

METHODS

H. pylori hopE and hopV genes derived from strain CHCTX-1 were cloned by polymerase chain reaction (PCR), sequenced and expressed in Escherichia coli AD494 (DE3). Gel-purified porins were used to prepare polyclonal antibodies. The presence of both genes was tested by PCR in a collection of H. pylori clinical isolates and their expression was detected in lysates by immunoblotting. Immune responses against HopE, HopV and other H. pylori antigens in sera from infected and non-infected patients were tested by Western blotting using these sera as first antibody on recombinant H. pylori antigens.

RESULTS

PCR and Western blotting assays revealed that 60 and 82 out of 130 Chilean isolates carried hopE and hopV genes, respectively, but only 16 and 9, respectively, expressed these porins. IgG serum immunoreactivity evaluation of 69 H. pylori-infected patients revealed that HopE and HopV were infrequently recognized (8.7% and 10.1% respectively) compared to H. pylori VacA (68.1%) and CagA (59.5%) antigens. Similar values were detected for IgA serum immunoreactivity against HopE (11.6%) and HopV (10.5%) although lower values for VacA (42%) and CagA (17.4%) were obtained when compared to the IgG response.

CONCLUSION

A scarce expression of HopE and HopV among Chilean isolates was found, in agreement with the infrequent seroconversion against these antigens when tested in infected Chilean patients.

Structure-function investigation of vsp serotypes of the spirochete Borrelia hermsii

Background

Relapsing fever (RF) spirochetes are notable for multiphasic antigenic variation of polymorphic outer membrane lipoproteins, a phenomenon responsible for immune evasion. An additional role in tissue localization is suggested by the finding that isogenic serotypes 1 (Bt1) and 2 (Bt2) of the RF spirochete Borrelia turicatae, which differ only in the Vsp they express, exhibit marked differences in clinical disease severity and tissue localization during infection.

Methodology/Principal Findings

Here we used known vsp DNA sequences encoding for B. turicatae and Borrelia hermsii Vsp proteins with variable regions and then studied whether there are differences in disease expression and tissue localization of their corresponding serotypes during mouse infection. For sequence and structural comparisons we focused exclusively on amino acid residues predicted to project away from the spirochetes surface, referred to as the Vsp dome. Disease severity and tissue localization were studied during persistent infection with individual or mixed serotypes in SCID mice. The results showed that all Vsp domes clustered into 3 main trunks, with the domes for B. turicatae Vsp1 (BtVsp1) and BtVsp2 clustering into separate ones. B. hermsii serotypes whose Vsp domes clustered with the BtVsp1 dome were less virulent but localized to the brain more. The BtVsp2 dome was the oddball among all and Bt2 was the only serotype that caused severe arthritis.

Conclusion/Significance

These findings indicate that there is significant variability in Vsp dome structure, disease severity, and tissue localization among serotypes of B. hermsii.

High frequency of hotspot mutations in core genes of Escherichia coli due to short-term positive selection

Core genes comprising the ubiquitous backbone of bacterial genomes are not subject to frequent horizontal transfer and generally are not thought to contribute to the adaptive evolution of bacterial pathogens. We determined, however, that at least one-third and possibly more than one-half of the core genes in Escherichia coli genomes are targeted by repeated replacement substitutions in the same amino acid positions-hotspot mutations. Occurrence of hotspot mutations is driven by positive selection, as their rate is significantly higher than expected by random chance alone, and neither intragenic recombination nor increased mutability can explain the observed patterns. Also, commensal E. coli strains have a significantly lower frequency of mutated genes and mutations per genome than pathogenic strains. E. coli strains causing extra-intestinal infections accumulate hotspot mutations at the highest rate, whereas the highest total number of mutated genes has been found among Shigella isolates, suggesting the pathoadaptive nature of such mutations. The vast majority of hotspot mutations are of recent evolutionary origin, implying short-term positive selection, where adaptive mutations emerge repeatedly but are not sustained in natural circulation for long. Such pattern of dynamics is consistent with source-sink model of virulence evolution.

Sir2 paralogues cooperate to regulate virulence genes and antigenic variation in Plasmodium falciparum

Cytoadherance of Plasmodium falciparum-infected erythrocytes in the brain, organs and peripheral microvasculature is linked to morbidity and mortality associated with severe malaria. Parasite-derived P. falciparum Erythrocyte Membrane Protein 1 (PfEMP1) molecules displayed on the erythrocyte surface are responsible for cytoadherance and undergo antigenic variation in the course of an infection. Antigenic variation of PfEMP1 is achieved by in situ switching and mutually exclusive transcription of the var gene family, a process that is controlled by epigenetic mechanisms. Here we report characterisation of the P. falciparum silent information regulator's A and B (PfSir2A and PfSir2B) and their involvement in mutual exclusion and silencing of the var gene repertoire. Analysis of P. falciparum parasites lacking either PfSir2A or PfSir2B shows that these NAD⁺-dependent histone deacetylases are required for silencing of different var gene subsets classified by their conserved promoter type. We also demonstrate that in the absence of either of these molecules mutually exclusive expression of var genes breaks down. We show that var gene silencing originates within the promoter and PfSir2 paralogues are involved in cis spreading of silenced chromatin into adjacent regions. Furthermore, parasites lacking PfSir2A but not PfSir2B have considerably longer telomeric repeats, demonstrating a role for this molecule in telomeric end protection. This work highlights the pivotal but distinct role for both PfSir2 paralogues in epigenetic silencing of P. falciparum virulence genes and the control of pathogenicity of malaria infection.

The unicellular parasite Plasmodium falciparum is the cause of the most severe form of malaria and is responsible for 300 million infections and ∼2 million deaths a year. Infected erythrocytes clump and block capillaries in the peripheral circulation, the brain, and placenta and are a major contributor to the pathology of malaria. A parasite-derived protein displayed on the surface of the infected erythrocyte is responsible for erythrocyte clumping in capillaries. Although 60 subtelomeric var genes can encode different versions of this “sticky” capillary-binding protein, only one protein is expressed at a time, and switches in expression between these genes causes variation of this pathogenic molecule enabling the parasite to evade the immune system. Here we identify two chromatin-modifying proteins that cooperate to mediate silencing and mutual exclusive expression of var genes. These proteins are thus important virulence factors of the malaria-causing parasite.

Investigation into two Sir2 histone deacetylases in the malaria-causing parasite revealstrans-acting epigenetic factors control mutually exclusive expression of a major subtelomeric virulence gene family.

Genetics of antigenic variation in Plasmodium falciparum

Malaria caused by the protozoan parasite Plasmodium falciparum is characterized by long-term, persistent infections that can last for many months. The ability of this parasite to avoid clearance by the human immune system is dependent on its capacity to continuously alter the surface exposed antigenic proteins that that are vulnerable to antibody recognition and attack, a process called antigenic variation. Significant work in recent years has contributed to our understanding of the mechanisms underlying this process, including the genes encoding the antigenic proteins and the DNA sequence elements that control their expression. In addition, the epigenetic "marks" that are associated with activation and silencing of individual genes have been extensively characterized. These studies have led to a model that includes multiple layers of regulation that ultimately lead to the tight coordination of expression of the genes responsible for antigenic variation by malaria parasites. Here we review some more recent data that adds additional complexity to our understanding of these regulatory layers.

Protective role of the lung collectins surfactant protein A and surfactant protein D in airway inflammation

The acute inflammatory airway response is characterized by a time-dependent onset followed by active resolution. Emerging evidence suggests that epithelial cells of the proximal and distal air spaces release host defense mediators that can facilitate both the initiation and the resolution part of inflammatory airway changes. These molecules, also known as the hydrophilic surfactant proteins (surfactant protein [SP]-A and SP-D) belong to the class of collagenous lectins (collectins). The collectins are a small family of soluble pattern recognition receptors containing collagenous regions and C-type lectin domains. SP-A and SP-D are most abundant in the lung. Because of their structural uniqueness, specific localization, and functional versatility, lung collectins are important players of the pulmonary immune responses. Recent studies in our laboratory and others indicated significant associations of lung collectin levels with acute and chronic airway inflammation in both animal models and patients, suggesting the usefulness of these molecules as disease biomarkers. Research on wild-type and mutant recombinant molecules in vivo and in vitro showed that SP-A and SP-D bind carbohydrates, lipids, and nucleic acids with a broad-spectrum specificity and initiate phagocytosis of inhaled pathogens as well as apoptotic cells. Investigations on gene-deficient and conditional overexpresser mice indicated that lung collectins also directly modulate innate immune cell function and T-cell-dependent inflammatory events. Thus, these molecules have a unique, dual-function capacity to induce pathogen elimination and control proinflammatory mechanisms, suggesting a potential suitability for therapeutic prevention and treatment of chronic airway inflammation. This article reviews evidence supporting that the lung collectins play an immune-protective role and are essential for maintenance of the immunologic homeostasis in the lung.

Loss of allergen 1 confers a hypervirulent phenotype that resembles mucoid switch variants of Cryptococcus neoformans

Microbial survival in a host is usually dependent on the ability of a pathogen to undergo changes that promote escape from host defense mechanisms. The human-pathogenic fungus Cryptococcus neoformans undergoes phenotypic switching in vivo that promotes persistence in tissue. By microarray and real-time PCR analyses, the allergen 1 gene (ALL1) was found to be downregulated in the hypervirulent mucoid switch variant, both during logarithmic growth and during intracellular growth in macrophages. The ALL1 gene encodes a small cytoplasmic protein that is involved in capsule formation. Growth of an all1Delta gene deletion mutant was normal. Similar to cells of the mucoid switch variant, all1Delta cells produced a larger polysaccharide capsule than cells of the smooth parent and the complemented strain produced, and the enlarged capsule inhibited macrophage phagocytosis. The mutant exhibited a modest defect in capsule induction compared to all of the other variants. In animal models the phenotype of the all1Delta mutant mimicked the hypervirulent phenotype of the mucoid switch variant, which is characterized by decreased host survival and elevated intracranial pressure. Decreased survival is likely the result of both an ineffective cell-mediated immune response and impaired phagocytosis by macrophages. Consequently, we concluded that, unlike loss of most virulence-associated genes, where loss of gene function results in attenuated virulence, loss of the ALL1 gene enhances virulence by altering the host-pathogen interaction and thereby impairing clearance. Our data identified the first cryptococcal gene associated with elevated intracranial pressure and support the hypothesis that an environmental opportunistic pathogen has modified its virulence in vivo by epigenetic downregulation of gene function.

varDB: a pathogen-specific sequence database of protein families involved in antigenic variation

Infectious diseases are a major threat to global public health and prosperity. The causative agents consist of a suite of pathogens, ranging from bacteria to viruses, including fungi, helminthes and protozoa. Although these organisms are extremely varied in their biological structure and interactions with the host, they share similar methods of evading the host immune system. Antigenic variation and drift are mechanisms by which pathogens change their exposed epitopes while maintaining protein function. Accordingly, these traits enable pathogens to establish chronic infections in the host. The varDB database was developed to serve as a central repository of protein and nucleotide sequences as well as associated features (e.g. field isolate data, clinical parameters, etc.) involved in antigenic variation. The data currently contained in varDB were mined from GenBank as well as multiple specialized data repositories (e.g. PlasmoDB, GiardiaDB). Family members and ortholog groups were identified using a hierarchical search strategy, including literature/author-based searches and HMM profiles. Included in the current release are>29,00 sequences from 39 gene families from 25 different pathogens. This resource will enable researchers to compare antigenic variation within and across taxa with the goal of identifying common mechanisms of pathogenicity to assist in the fight against a range of devastating diseases.

AVAILABILITY

varDB is freely accessible at http://www.vardb.org/

Trypanosoma brucei BRCA2 acts in antigenic variation and has undergone a recent expansion in BRC repeat number that is important during homologous recombination

Antigenic variation in Trypanosoma brucei has selected for the evolution of a massive archive of silent Variant Surface Glycoprotein (VSG) genes, which are activated by recombination into specialized expression sites. Such VSG switching can occur at rates substantially higher than background mutation and is dependent on homologous recombination, a core DNA repair reaction. A key regulator of homologous recombination is BRCA2, a protein that binds RAD51, the enzyme responsible for DNA strand exchange. Here, we show that T. brucei BRCA2 has undergone a recent, striking expansion in the number of BRC repeats, a sequence element that mediates interaction with RAD51. T. brucei BRCA2 mutants are shown to be significantly impaired in antigenic variation and display genome instability. By generating BRCA2 variants with reduced BRC repeat numbers, we show that the BRC expansion is crucial in determining the efficiency of T. brucei homologous recombination and RAD51 localization. Remarkably, however, this appears not to be a major determinant of the activation of at least some VSG genes.

Plasmodium falciparum antigenic variation. Mapping mosaic var gene sequences onto a network of shared, highly polymorphic sequence blocks

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) is a potentially important family of immune targets, encoded by an extremely diverse gene family called var. Understanding of the genetic organization of var genes is hampered by sequence mosaicism that results from a long history of non-homologous recombination. Here we have used software designed to analyse social networks to visualize the relationships between large collections of short var sequences tags sampled from clinical parasite isolates. In this approach, two sequences are connected if they share one or more highly polymorphic sequence blocks. The results show that the majority of analysed sequences including several var-like sequences from the chimpanzee parasite Plasmodium reichenowi can be either directly or indirectly linked together in a single unbroken network. However, the network is highly structured and contains putative subgroups of recombining sequences. The major subgroup contains the previously described group A var genes, previously proposed to be genetically distinct. Another subgroup contains sequences found to be associated with rosetting, a parasite virulence phenotype. The mosaic structure of the sequences and their division into subgroups may reflect the conflicting problems of maximizing antigenic diversity and minimizing epitope sharing between variants while maintaining their host cell binding functions.

Virulence factor rtx in Legionella pneumophila, evidence suggesting it is a modular multifunctional protein

BACKGROUND

The repeats in toxin (Rtx) are an important pathogenicity factor involved in host cells invasion of Legionella pneumophila and other pathogenic bacteria. Its role in escaping the host immune system and cytotoxic activity is well known. Its repeated motives and modularity make Rtx a multifunctional factor in pathogenicity.

RESULTS

The comparative analysis of rtx gene among 6 strains of L. pneumophila showed modularity in their structures. Among compared genomes, the N-terminal region of the protein presents highly dissimilar repeats with functionally similar domains. On the contrary, the C-terminal region is maintained with a fashionable modular configuration, which gives support to its proposed role in adhesion and pore formation. Despite the variability of rtx among the considered strains, the flanking genes are maintained in synteny and similarity.

CONCLUSION

In contrast to the extracellular bacteria Vibrio cholerae, in which the rtx gene is highly conserved and flanking genes have lost synteny and similarity, the gene region coding for the Rtx toxin in the intracellular pathogen L. pneumophila shows a rapid evolution. Changes in the rtx could play a role in pathogenicity. The interplay of the Rtx toxin with host membranes might lead to the evolution of new variants that are able to escape host cell defences.

Alterations in local chromatin environment are involved in silencing and activation of subtelomeric var genes in Plasmodium falciparum

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), encoded by the var gene family, undergoes antigenic variation and plays an important role in chronic infection and severe malaria. Only a single var gene is transcribed per parasite, and epigenetic control mechanisms are fundamental in this strategy of mutually exclusive transcription. We show that subtelomeric upsB var gene promoters carried on episomes are silenced by default, and that promoter activation is sufficient to silence all other family members. However, they are active by default when placed downstream of a second active var promoter, underscoring the significance of local chromatin environment and nuclear compartmentalization in var promoter regulation. Native chromatin covering the SPE2-repeat array in upsB promoters is resistant to nuclease digestion, and insertion of these regulatory elements into a heterologous promoter causes local alterations in nucleosomal organization and promoter repression. Our findings suggest a common logic underlying the transcriptional control of all var genes, and have important implications for our understanding of the epigenetic processes involved in the regulation of this major virulence gene family.

Trypanosoma brucei homologous recombination is dependent on substrate length and homology, though displays a differential dependence on mismatch repair as substrate length decreases

Homologous recombination functions universally in the maintenance of genome stability through the repair of DNA breaks and in ensuring the completion of replication. In some organisms, homologous recombination can perform more specific functions. One example of this is in antigenic variation, a widely conserved mechanism for the evasion of host immunity. Trypanosoma brucei, the causative agent of sleeping sickness in Africa, undergoes antigenic variation by periodic changes in its variant surface glycoprotein (VSG) coat. VSG switches involve the activation of VSG genes, from an enormous silent archive, by recombination into specialized expression sites. These reactions involve homologous recombination, though they are characterized by an unusually high rate of switching and by atypical substrate requirements. Here, we have examined the substrate parameters of T. brucei homologous recombination. We show, first, that the reaction is strictly dependent on substrate length and that it is impeded by base mismatches, features shared by homologous recombination in all organisms characterized. Second, we identify a pathway of homologous recombination that acts preferentially on short substrates and is impeded to a lesser extent by base mismatches and the mismatch repair machinery. Finally, we show that mismatches during T. brucei recombination may be repaired by short-patch mismatch repair.

Temporal expression analysis of the Borrelia burgdorferi paralogous gene family 54 genes BBA64, BBA65, and BBA66 during persistent infection in mice

Members of the Borrelia burgdorferi paralogous gene family 54 (pgf 54) are regulated by conditions simulating mammalian infection and are thought to be instrumental in borrelial host survival and pathogenesis. To explore the activities of these genes in vivo, a comprehensive analysis of pgf 54 genes BBA64, BBA65, and BBA66 was performed to assess the genetic stability, host antibody responses, and kinetics of gene expression in the murine model of persistent infection. DNA sequencing of pgf 54 genes obtained from re-isolates at 1 year postinfection demonstrated that all genes of this family are stable and do not undergo recombination to generate variant antigens during persistent infection. Antibodies against BBA64 and BBA66 appeared soon after infection and were detectable throughout the infection, suggesting that there was gene expression during infection. However, quantitative reverse transcription-PCR revealed that BBA64 gene expression was considerably decreased in Borrelia residing in the mouse ear tissue compared to the expression in cultured spirochetes by 20 days postinfection and that the levels of expression remained low throughout the infection. Conversely, transcription of the BBA65 and BBA66 genes was increased, and both of these genes were continuously expressed until 100 days postinfection; this was followed by periods of differential expression late in infection. The expression profile of the BBA64 gene suggests that this gene has an important role during tick-to-host transmission and early infection, whereas the expression profile of the BBA65 and BBA66 genes suggests that these genes have a role in persistent infection. The differential regulation of pgf 54 genes observed during infection may help confer a survival advantage during persistent infection, influencing mechanisms for B. burgdorferi dissemination, tissue tropism, or evasion of the adaptive immune response.

Phenotypic switching and its implications for the pathogenesis of Cryptococcus neoformans

Phenotypic switching has been described in several strains of Cryptococcus neoformans. It occurs in vivo during chronic infection and is associated with differential gene expression and changes in virulence. The switch involves changes in the polysaccharide capsule and cell wall that affect the yeast's ability to resist phagocytosis. In addition, the phenotypic switch variants elicit qualitatively different inflammatory responses in the host. The host's immune response ultimately affects selection of the switch variants in animal models of chronic cryptococcosis. The biological relevance of phenotypic switching is demonstrated in several murine infection models and further underlines the importance of phenotypic switching in the setting of human disease. This includes the association of switching and poor outcome in chronic infection, the ability of the mucoid variant of strain RC-2 (RC-2 MC) but not the smooth variant (RC-2 SM) to promote increased intracranial pressure in a rat model, and lastly the observation that antifungal interventions can promote the selection of more virulent switch variants during chronic murine infection.

Extreme sequence divergence but conserved ligand-binding specificity in Streptococcus pyogenes M protein

Many pathogenic microorganisms evade host immunity through extensive sequence variability in a protein region targeted by protective antibodies. In spite of the sequence variability, a variable region commonly retains an important ligand-binding function, reflected in the presence of a highly conserved sequence motif. Here, we analyze the limits of sequence divergence in a ligand-binding region by characterizing the hypervariable region (HVR) of Streptococcus pyogenes M protein. Our studies were focused on HVRs that bind the human complement regulator C4b-binding protein (C4BP), a ligand that confers phagocytosis resistance. A previous comparison of C4BP-binding HVRs identified residue identities that could be part of a binding motif, but the extended analysis reported here shows that no residue identities remain when additional C4BP-binding HVRs are included. Characterization of the HVR in the M22 protein indicated that two relatively conserved Leu residues are essential for C4BP binding, but these residues are probably core residues in a coiled-coil, implying that they do not directly contribute to binding. In contrast, substitution of either of two relatively conserved Glu residues, predicted to be solvent-exposed, had no effect on C4BP binding, although each of these changes had a major effect on the antigenic properties of the HVR. Together, these findings show that HVRs of M proteins have an extraordinary capacity for sequence divergence and antigenic variability while retaining a specific ligand-binding function.

Many pathogens have evolved mechanisms to evade host immunity. In one such mechanism, the sequence of a surface protein varies among different strains of a pathogen. This sequence variability represents an apparent paradox, because the variable protein must retain an important function. The authors studied this problem in Streptococcus pyogenes, a major human pathogen. The surface-localized M protein of this bacterium varies extensively in sequence between bacterial strains, allowing immune escape. Nevertheless, the most variable part of the M protein commonly binds a human plasma protein. By hijacking this human protein the bacteria evade attack by complement an important part of the innate immune system. Comparison of the ligand-binding region in different M proteins showed that these regions lack a shared amino acid sequence motif. Thus, a variable protein can retain a ligand-binding function in the absence of a conserved binding motif. Evidence is also presented that a single amino acid change in the variable region may cause a major antigenic change, providing a selective advantage for the bacteria. Together, these data bear witness to the extraordinary ability of pathogens to escape host immunity, without losing ability to cause disease.

Phenotypic switching in Cryptococcus neoformans

Phenotypic switching has been described in serotype A and D strains of Cryptococcus neoformans. It occurs in vivo during chronic infection and is associated with differential gene expression and changes in virulence. The switch involves changes in the polysaccharide capsule and cell wall that affect the yeast's ability to resist phagocytosis. In addition, the phenotypic switch variants elicit qualitatively different inflammatory responses in the host. In animal models of chronic cryptococosis, the immune response of the host ultimately determines which of the switch variants are selected and maintained. The importance of phenotypic switching is further underscored by several findings that are relevant in the setting of human disease. These include the ability of the mucoid colony variant of RC-2 (RC-2 MC) but not the smooth variant (RC-2 SM) to promote increased intracerebral pressure in a rat model of cryptococcal meningitis. Furthermore, chemotherapeutic and immunological antifungal interventions can promote the selection of the RC-2 MC variant during chronic murine infection.

Distinct roles for two RAD51-related genes in Trypanosoma brucei antigenic variation

In Trypanosoma brucei, DNA recombination is crucial in antigenic variation, a strategy for evading the mammalian host immune system found in a wide variety of pathogens. T.brucei has the capacity to encode >1000 antigenically distinct variant surface glycoproteins (VSGs). By ensuring that only one VSG is expressed on the cell surface at one time, and by periodically switching the VSG gene that is expressed, T.brucei can evade immune killing for prolonged periods. Much of VSG switching appears to rely on a widely conserved DNA repair pathway called homologous recombination, driven by RAD51. Here, we demonstrate that T.brucei encodes a further five RAD51-related proteins, more than has been identified in other single-celled eukaryotes to date. We have investigated the roles of two of the RAD51-related proteins in T.brucei, and show that they contribute to DNA repair, homologous recombination and RAD51 function in the cell. Surprisingly, however, only one of the two proteins contributes to VSG switching, suggesting that the family of diverged RAD51 proteins present in T.brucei have assumed specialized functions in homologous recombination, analogous to related proteins in metazoan eukaryotes.

Replication dynamics of Mycobacterium tuberculosis in chronically infected mice

The dynamics of host-pathogen interactions have important implications for the design of new antimicrobial agents to treat chronic infections such as tuberculosis (TB), which is notoriously refractory to conventional drug therapy. In the mouse model of TB, an acute phase of exponential bacterial growth in the lungs is followed by a chronic phase characterized by relatively stable numbers of bacteria. This equilibrium could be static, with little ongoing replication, or dynamic, with continuous bacterial multiplication balanced by bacterial killing. A static model predicts a close correspondence between "viable counts" (live bacteria) and "total counts" (live plus dead bacteria) in the lungs over time. A dynamic model predicts the divergence of total counts and viable counts over time due to the accumulation of dead bacteria. Here, viable counts are defined as bacterial CFU enumerated by plating lung homogenates; total counts are defined as bacterial chromosome equivalents (CEQ) enumerated by using quantitative real-time PCR. We show that the viable and total bacterial counts in the lungs of chronically infected mice do not diverge over time. Rapid degradation of dead bacteria is unlikely to account for the stability of bacterial CEQ numbers in the lungs over time, because treatment of mice with isoniazid for 8 weeks led to a marked reduction in the number of CFU without reducing the number of CEQ. These observations support the hypothesis that the stable number of bacterial CFU in the lungs during chronic infection represents a static equilibrium between host and pathogen.

Genomic analysis of Bacteroides fragilis reveals extensive DNA inversions regulating cell surface adaptation

Bacteroides are predominant human colonic commensals, but the principal pathogenic species, Bacteroides fragilis (BF), lives closely associated with the mucosal surface, whereas a second major species, Bacteroides thetaiotaomicron (BT), concentrates within the colon. We find corresponding differences in their genomes, based on determination of the genome sequence of BF and comparative analysis with BT. Both species have acquired two mechanisms that contribute to their dominance among the colonic microbiota: an exceptional capability to use a wide range of dietary polysaccharides by gene amplification and the capacity to create variable surface antigenicities by multiple DNA inversion systems. However, the gene amplification for polysaccharide assimilation is more developed in BT, in keeping with its internal localization. In contrast, external antigenic structures can be changed more systematically in BF. Thereby, at the mucosal surface, where microbes encounter continuous attack by host defenses, BF evasion of the immune system is favored, and its colonization and infectious potential are increased.