Evolutionary structure of Plasmodium falciparum major variant surface antigen genes in South America: Implications for epidemic transmission and surveillance

Strong founder effects resulting from human migration out of Africa have led to geographic variation in single nucleotide polymorphisms (SNPs) and microsatellites (MS) of the malaria parasite, Plasmodium falciparum. This is particularly striking in South America where two major founder populations of P. falciparum have been identified that are presumed to have arisen from the transatlantic slave trade. Given the importance of the major variant surface antigen of the blood stages of P. falciparum as both a virulence factor and target of immunity, we decided to investigate the population genetics of the genes encoding “Plasmodium falciparum Erythrocyte Membrane Protein 1” (PfEMP1) among several countries in South America, in order to evaluate the transmission patterns of malaria in this continent. Deep sequencing of the DBLα domain of var genes from 128 P. falciparum isolates from five locations in South America was completed using a 454 high throughput sequencing protocol. Striking geographic variation in var DBLα sequences, similar to that seen for SNPs and MS markers, was observed. Colombia and French Guiana had distinct var DBLα sequences, whereas Peru and Venezuela showed an admixture. The importance of such geographic variation to herd immunity and malaria vaccination is discussed.     

Age-specific patterns of DBLα var diversity can explain why residents of high malaria transmission areas remain susceptible to Plasmodium falciparum blood stage infection throughout life

Immunity to Plasmodium falciparum is non-sterilising, thus individuals residing in malaria-endemic areas are at risk of infection throughout their lifetime. Here we seek to find a genomic epidemiological explanation for why residents of all ages harbour blood stage infections despite lifelong exposure to P. falciparum in areas of high transmission. We do this by exploring, for the first known time, the age-specific patterns of diversity of variant antigen encoding (var) genes in the reservoir of infection. Microscopic and submicroscopic P. falciparum infections were analysed at the end of the wet and dry seasons in 2012–2013 for a cohort of 1541 residents aged from 1 to 91 years in an area characterised by high seasonal malaria transmission in Ghana. By sequencing the near ubiquitous Duffy-binding-like alpha domain (DBLα) that encodes immunogenic domains, we defined var gene diversity in an estimated 1096 genomes detected in sequential wet and dry season sampling of this cohort. Unprecedented var (DBLα) diversity was observed in all ages with 42,399 unique var types detected. There was a high degree of maintenance of types between seasons (>40% seen more than once), with many of the same types, especially upsA, appearing multiple times in isolates from different individuals. Children and adolescents were found to be significant reservoirs of var DBLα diversity compared with adults. Var repertoires within individuals were highly variable, with children having more related var repertoires compared to adolescents and adults. Individuals of all ages harboured multiple genomes with var repertoires unrelated to those infecting other hosts. High turnover of parasites with diverse isolate var repertoires was also observed in all ages. These age-specific patterns are best explained by variant-specific immune selection. The observed level of var diversity for the population was then used to simulate the development of variant-specific immunity to the diverse var types under conservative assumptions. Simulations showed that the extent of observed var diversity with limited repertoire relatedness was sufficient to explain why adolescents and adults in this community remain susceptible to blood stage infection, even with multiple genomes.     

Phylogeography of var gene repertoires reveals fine‐scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine‐scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high‐quality DBLα‐sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine‐scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.     

Multiple var2csa-Type PfEMP1 Genes Located at Different Chromosomal Loci Occur in Many Plasmodium falciparum Isolates

Background

The var2csa gene encodes a Plasmodium falciparum adhesion receptor which binds chondroitin sulfate A (CSA). This var gene is more conserved than other PfEMP1/var genes and is found in all P. falciparum isolates. In isolates 3D7, FCR3/It4 and HB3, var2csa is transcribed from a sub-telomeric position on the left arm of chromosome 12, but it is not known if this location is conserved in all parasites. Genome sequencing indicates that the var2csa gene is duplicated in HB3, but whether this is true in natural populations is uncertain.

Methodology/Principal Findings

To assess global variation in the VAR2CSA protein, sequence variation in the DBL2X region of var2csa genes in 54 P.falciparum samples was analyzed. Chromosome mapping of var2csa loci was carried out and a quantitative PCR assay was developed to estimate the number of var2csa genes in P.falciparum isolates from the placenta of pregnant women and from the peripheral circulation of other malaria patients. Sequence analysis, gene mapping and copy number quantitation in P.falciparum isolates indicate that there are at least two loci and that both var2csa-like genes can be transcribed. All VAR2CSA DBL2X domains fall into one of two distinct phylogenetic groups possessing one or the other variant of a large (∼26 amino acid) dimorphic motif, but whether either motif variant is linked to a specific locus is not known.

Conclusions/Significance

Two or more related but distinct var2csa-type PfEMP1/var genes exist in many P. falciparum isolates. One gene is on chromosome 12 but additional var2csa-type genes are on different chromosomes in different isolates. Multiplicity of var2csa genes appears more common in infected placentae than in samples from non-pregnant donors indicating a possible advantage of this genotype in pregnancy associated malaria.     

DNA secondary structures are associated with recombination in major Plasmodium falciparum variable surface antigen gene families

Many bacterial, viral and parasitic pathogens undergo antigenic variation to counter host immune defense mechanisms. In Plasmodium falciparum, the most lethal of human malaria parasites, switching of var gene expression results in alternating expression of the adhesion proteins of the Plasmodium falciparum-erythrocyte membrane protein 1 class on the infected erythrocyte surface. Recombination clearly generates var diversity, but the nature and control of the genetic exchanges involved remain unclear. By experimental and bioinformatic identification of recombination events and genome-wide recombination hotspots in var genes, we show that during the parasite’s sexual stages, ectopic recombination between isogenous var paralogs occurs near low folding free energy DNA 50-mers and that these sequences are heavily concentrated at the boundaries of regions encoding individual Plasmodium falciparum-erythrocyte membrane protein 1 structural domains. The recombinogenic potential of these 50-mers is not parasite-specific because these sequences also induce recombination when transferred to the yeast Saccharomyces cerevisiae. Genetic cross data suggest that DNA secondary structures (DSS) act as inducers of recombination during DNA replication in P. falciparum sexual stages, and that these DSS-regulated genetic exchanges generate functional and diverse P. falciparum adhesion antigens. DSS-induced recombination may represent a common mechanism for optimizing the evolvability of virulence gene families in pathogens.     

Ectopic Recombination of a Malaria var Gene during Mitosis Associated with an Altered var Switch Rate

The Plasmodium falciparum var multigene family encodes P. falciparum erythrocyte membrane protein 1, which is responsible for the pathogenic traits of antigenic variation and adhesion of infected erythrocytes to host receptors during malaria infection. Clonal antigenic variation of P. falciparum erythrocyte membrane protein 1 is controlled by the switching between exclusively transcribed var genes. The tremendous diversity of the var gene repertoire both within and between parasite strains is critical for the parasite's strategy of immune evasion. We show that ectopic recombination between var genes occurs during mitosis, providing P. falciparum with opportunities to diversify its var repertoire, even during the course of a single infection. We show that the regulation of the recombined var gene has been disrupted, resulting in its persistent activation although the regulation of most other var genes is unaffected. The var promoter and intron of the recombined var gene are not responsible for its atypically persistent activity, and we conclude that altered subtelomeric cis sequence is the most likely cause of the persistent activity of the recombined var gene.     

Genetic diversity and distribution patterns of PfEMP1 in Plasmodium falciparum isolates along the Thai-Myanmar border

Duffy binding–like domain (DBL) and cysteine-rich interdomain region (CIDR) domain genes of Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) encode malaria virulence proteins. The variants of these genes have been reported to be associated with severe/complicated malaria. The present study investigated the prevalence and distribution patterns of DBLα0.6/9, DBLα1.1, DBLα1 not var3 genes, DBLα2/α1.1/2/4/7, DBLβ12 & DBLβ3/5, DBLε8, CIDRα1.4, and CIDRα1.6 of P. falciparum isolates along the Thai-Myanmar border. The association between PfEMP1 variants and parasite density was also investigated. Two hundred and thirteen finger-prick dried blood spot (DBS) or whole blood samples were collected in 2007 and 2015, from patients with acute uncomplicated P. falciparum in Tak, Kanchanaburi, and Ranong provinces. Analysis of the variant genes was performed using polymerase chain reaction (PCR). The DBLs variant which was found at the highest and lowest frequencies in the three provinces were DBLα1 not var3 (72.77%), and DBLε8 (17.37%). The two CIDR domain variants were found at relatively lower frequencies compared with DBL domain variants (9.9% and 30.1%). P. falciparum isolates carrying the four PfEMP1 variants, i.e., DBLα0.6/9, DBLα1.1, DBLα2/α.1.1/2/4/7, and DBLε8 were found to be significantly associated with low parasitemia. Both DBLα0.6/9 and DBLα2/α1.1/2/4/7 variant genes which were present at high frequencies in this border area could be potential candidate markers for predicting P. falciparum hyperparasitemia and in this border area. Furthermore, the information could be exploited as candidate proteins for the development of an effective malaria vaccine in specific malaria-endemic areas.     

Generation of Antigenic Diversity in Plasmodium falciparum by Structured Rearrangement of Var Genes During Mitosis

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7–72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle.     

Plasmodium falciparum Erythrocyte Membrane Protein 1 Diversity in Seven Genomes – Divide and Conquer

The var gene encoded hyper-variable Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family mediates cytoadhesion of infected erythrocytes to human endothelium. Antibodies blocking cytoadhesion are important mediators of malaria immunity acquired by endemic populations. The development of a PfEMP1 based vaccine mimicking natural acquired immunity depends on a thorough understanding of the evolved PfEMP1 diversity, balancing antigenic variation against conserved receptor binding affinities. This study redefines and reclassifies the domains of PfEMP1 from seven genomes. Analysis of domains in 399 different PfEMP1 sequences allowed identification of several novel domain classes, and a high degree of PfEMP1 domain compositional order, including conserved domain cassettes not always associated with the established group A–E division of PfEMP1. A novel iterative homology block (HB) detection method was applied, allowing identification of 628 conserved minimal PfEMP1 building blocks, describing on average 83% of a PfEMP1 sequence. Using the HBs, similarities between domain classes were determined, and Duffy binding-like (DBL) domain subclasses were found in many cases to be hybrids of major domain classes. Related to this, a recombination hotspot was uncovered between DBL subdomains S2 and S3. The VarDom server is introduced, from which information on domain classes and homology blocks can be retrieved, and new sequences can be classified. Several conserved sequence elements were found, including: (1) residues conserved in all DBL domains predicted to interact and hold together the three DBL subdomains, (2) potential integrin binding sites in DBLα domains, (3) an acylation motif conserved in group A var genes suggesting N-terminal N-myristoylation, (4) PfEMP1 inter-domain regions proposed to be elastic disordered structures, and (5) several conserved predicted phosphorylation sites. Ideally, this comprehensive categorization of PfEMP1 will provide a platform for future studies on var/PfEMP1 expression and function.     

Defining the Antigenic Diversity of Plasmodium falciparum Apical Membrane Antigen 1 and the Requirements for a Multi-Allele Vaccine against Malaria

Apical Membrane Antigen 1 (AMA1) is a leading malaria vaccine candidate and a target of naturally-acquired human immunity. Plasmodium falciparum AMA1 is polymorphic and in vaccine trials it induces strain-specific protection. This antigenic diversity is a major roadblock to development of AMA1 as a malaria vaccine and understanding how to overcome it is essential. To assess how AMA1 antigenic diversity limits cross-strain growth inhibition, we assembled a panel of 18 different P. falciparum isolates which are broadly representative of global AMA1 sequence diversity. Antibodies raised against four well studied AMA1 alleles (W2Mef, 3D7, HB3 and FVO) were tested for growth inhibition of the 18 different P. falciparum isolates in growth inhibition assays (GIA). All antibodies demonstrated substantial cross-inhibitory activity against different isolates and a mixture of the four different AMA1 antibodies inhibited all 18 isolates tested, suggesting significant antigenic overlap between AMA1 alleles and limited antigenic diversity of AMA1. Cross-strain inhibition by antibodies was only moderately and inconsistently correlated with the level of sequence diversity between AMA1 alleles, suggesting that sequence differences are not a strong predictor of antigenic differences or the cross-inhibitory activity of anti-allele antibodies. The importance of the highly polymorphic C1-L region for inhibitory antibodies and potential vaccine escape was assessed by generating novel transgenic P. falciparum lines for testing in GIA. While the polymorphic C1-L epitope was identified as a significant target of some growth-inhibitory antibodies, these antibodies only constituted a minor proportion of the total inhibitory antibody repertoire, suggesting that the antigenic diversity of inhibitory epitopes is limited. Our findings support the concept that a multi-allele AMA1 vaccine would give broad coverage against the diversity of AMA1 alleles and establish new tools to define polymorphisms important for vaccine escape.     

A molecular epidemiological study of var gene diversity to characterize the reservoir of Plasmodium falciparum in humans in Africa

Background

The reservoir of Plasmodium infection in humans has traditionally been defined by blood slide positivity. This study was designed to characterize the local reservoir of infection in relation to the diverse var genes that encode the major surface antigen of Plasmodium falciparum blood stages and underlie the parasite''s ability to establish chronic infection and transmit from human to mosquito.

Methodology/Principal Findings

We investigated the molecular epidemiology of the var multigene family at local sites in Gabon, Senegal and Kenya which differ in parasite prevalence and transmission intensity. 1839 distinct var gene types were defined by sequencing DBLα domains in the three sites. Only 76 (4.1%) var types were found in more than one population indicating spatial heterogeneity in var types across the African continent. The majority of var types appeared only once in the population sample. Non-parametric statistical estimators predict in each population at minimum five to seven thousand distinct var types. Similar diversity of var types was seen in sites with different parasite prevalences.

Conclusions/Significance

Var population genomics provides new insights into the epidemiology of P. falciparum in Africa where malaria has never been conquered. In particular, we have described the extensive reservoir of infection in local African sites and discovered a unique var population structure that can facilitate superinfection through minimal overlap in var repertoires among parasite genomes. Our findings show that var typing as a molecular surveillance system defines the extent of genetic complexity in the reservoir of infection to complement measures of malaria prevalence. The observed small scale spatial diversity of var genes suggests that var genetics could greatly inform current malaria mapping approaches and predict complex malaria population dynamics due to the import of var types to areas where no widespread pre-existing immunity in the population exists.     

Plasmodium falciparum Variability and Immune Evasion Proceed from Antigenicity of Consensus Sequences from DBL6ε; Generalization to All DBL from VAR2CSA

We studied all consensus sequences within the four least ‘variable blocks’ (VB) present in the DBL6ε domain of VAR2CSA, the protein involved in the adhesion of infected red blood cells by Plasmodium falciparum that causes the Pregnancy-Associated Malaria (PAM). Characterising consensus sequences with respect to recognition of antibodies and percentage of responders among pregnant women living in areas where P. falciparum is endemic allows the identification of the most antigenic sequences within each VB. When combining these consensus sequences among four serotypes from VB1 or VB5, the most often recognized ones are expected to induce pan-reactive antibodies recognizing VAR2CSA from all plasmodial strains. These sequences are of main interest in the design of an immunogenic molecule. Using a similar approach than for DBL6ε, we studied the five other DBL and the CIDRpam from VAR2CSA, and again identified VB segments with highly conserved consensus sequences. In addition, we identified consensus sequences in other var genes expressed by non-PAM parasites. This finding paves the way for vaccine design against other pathologies caused by P. falciparum.     

A Network Approach to Analyzing Highly Recombinant Malaria Parasite Genes

The var genes of the human malaria parasite Plasmodium falciparum present a challenge to population geneticists due to their extreme diversity, which is generated by high rates of recombination. These genes encode a primary antigen protein called PfEMP1, which is expressed on the surface of infected red blood cells and elicits protective immune responses. Var gene sequences are characterized by pronounced mosaicism, precluding the use of traditional phylogenetic tools that require bifurcating tree-like evolutionary relationships. We present a new method that identifies highly variable regions (HVRs), and then maps each HVR to a complex network in which each sequence is a node and two nodes are linked if they share an exact match of significant length. Here, networks of var genes that recombine freely are expected to have a uniformly random structure, but constraints on recombination will produce network communities that we identify using a stochastic block model. We validate this method on synthetic data, showing that it correctly recovers populations of constrained recombination, before applying it to the Duffy Binding Like-α (DBLα) domain of var genes. We find nine HVRs whose network communities map in distinctive ways to known DBLα classifications and clinical phenotypes. We show that the recombinational constraints of some HVRs are correlated, while others are independent. These findings suggest that this micromodular structuring facilitates independent evolutionary trajectories of neighboring mosaic regions, allowing the parasite to retain protein function while generating enormous sequence diversity. Our approach therefore offers a rigorous method for analyzing evolutionary constraints in var genes, and is also flexible enough to be easily applied more generally to any highly recombinant sequences.     

Induction of Strain-Transcending Antibodies Against Group A PfEMP1 Surface Antigens from Virulent Malaria Parasites

Sequence diversity in pathogen antigens is an obstacle to the development of interventions against many infectious diseases. In malaria caused by Plasmodium falciparum, the PfEMP1 family of variant surface antigens encoded by var genes are adhesion molecules that play a pivotal role in malaria pathogenesis and clinical disease. PfEMP1 is a major target of protective immunity, however, development of drugs or vaccines based on PfEMP1 is problematic due to extensive sequence diversity within the PfEMP1 family. Here we identified the PfEMP1 variants transcribed by P. falciparum strains selected for a virulence-associated adhesion phenotype (IgM-positive rosetting). The parasites transcribed a subset of Group A PfEMP1 variants characterised by an unusual PfEMP1 architecture and a distinct N-terminal domain (either DBLα1.5 or DBLα1.8 type). Antibodies raised in rabbits against the N-terminal domains showed functional activity (surface reactivity with live infected erythrocytes (IEs), rosette inhibition and induction of phagocytosis of IEs) down to low concentrations (<10 µg/ml of total IgG) against homologous parasites. Furthermore, the antibodies showed broad cross-reactivity against heterologous parasite strains with the same rosetting phenotype, including clinical isolates from four sub-Saharan African countries that showed surface reactivity with either DBLα1.5 antibodies (variant HB3var6) or DBLα1.8 antibodies (variant TM284var1). These data show that parasites with a virulence-associated adhesion phenotype share IE surface epitopes that can be targeted by strain-transcending antibodies to PfEMP1. The existence of shared surface epitopes amongst functionally similar disease-associated P. falciparum parasite isolates suggests that development of therapeutic interventions to prevent severe malaria is a realistic goal.