Evolution of allelic dimorphism in malarial surface antigens

Role for gene conversion in the evolution of cell-surface antigens of the malaria parasite Plasmodium falciparum

While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.

Most Plasmodiophora brassicae Populations in Single Canola Root Galls from Alberta Fields are Mixtures of Multiple Strains

Clubroot, caused by Plasmodiophora brassicae, is an important disease on canola in Alberta, Canada. The pathogen is grouped into pathotypes according to their virulence reaction on differential hosts. Multiple pathotypes or strains are known exist in one field, one plant, or even one gall. This study was conducted with the objective of testing the prevalence of the coexistence of multiple strains in a single gall. In all, 79 canola clubroot galls were collected from 22 fields across Northern Alberta in 2018. Genomic DNA extracted from these single galls was analyzed using RNase H-dependent PCR (rhPCR). The rhPCR primers were designed to amplify a partial sequence of a dimorphic gene, with one primer pair specific to one sequence and the other primer pair specific to the alternative sequence. The amplification of both sequences from DNA obtained from a single gall would indicate that it contains two different P. brassicae strains. The rhPCR analyses indicated that the P. brassicae populations in 50 of the 79 galls consisted of more than one strain. This result emphasizes the need for cautious interpretation of results when a single-gall population is subject to pathotyping or being used as inoculum in plant pathology research. It also confirms that the maintenance of pathotype diversity within single root galls is a common occurrence which has implications for the durability, and stewardship, of single-gene host resistance.

Genetic Diversity of Plasmodium vivax Causing Epidemic Malaria in the Republic of Korea

Plasmodium vivax is more challenging to control and eliminate than P. falciparum due to its more asymptomatic infections with low parasite densities making diagnosis more difficult, in addition to its unique biological characteristics. The potential re-introduction of incidence cases, either through borders or via human migrations, is another major hurdle to sustained control and elimination. The Republic of Korea has experienced re-emergence of vivax malaria in 1993 but is one of the 32 malaria-eliminating countries to-date. Despite achieving successful nationwide control and elimination of vivax malaria, the evolutionary characteristics of vivax malaria isolates in the Republic of Korea have not been fully understood. In this review, we present an overview of the genetic variability of such isolates to increase understanding of the epidemiology, diversity, and dynamics of vivax populations in the Republic of Korea.

Molecular epidemiology and evolution of drug-resistant genes in the malaria parasite Plasmodium falciparum in southwestern Nigeria

Malaria is an age-old disease of human kind living in the tropical and sub-tropical regions of the globe, with Africa contributing the highest incidence of morbidity and mortality. Among many hurdles, evolution and spread of drug-resistant Plasmodium falciparum parasites constitute major challenges to malaria control and elimination. Information on molecular epidemiology and pattern of evolution of genes conferring resistance to different antimalarials are needed to track the route of the spread of resistant parasites and also to inform if the drug-resistant genes are adapted in the population following the Darwinian model of evolution. In the present study, we have followed molecular methods to detect both the known and emerging mutations in three genes (Pfcrt, Pfdhfr and Pfdhps) of P. falciparum conferring resistance to chloroquine and sulfadoxine-pyrimethamine from two different states (Edo: meso-endemic and Lagos: hypo-endemic) in southwestern Nigeria. High diversities in haplotypes and nucleotides in genes responsible for chloroquine (Pfcrt) and sulfadoxine (Pfdhps) resistance are recorded. About 96% of Pfdhfr and Pfdhps gene in both the meso- and hypo- endemic areas were mutant type, followed by 61% in Pfcrt gene. Many unique haplotypes of Pfdhps and Pfcrt were found to be segregated in these two populations. One particular mutant haplotype of Pfdhfr (AIRNI) was found to be in very high frequency in both Lagos and Edo. While the net haplotype diversity was highest in Pfdhps (0.81 in Lagos, 0.87 in Edo), followed by Pfcrt (0.69 in Lagos, 0.65 in Edo); highest number of haplotype was found in Pfdhps with 13 distinct haplotypes, followed by seven in Pfcrt and four in Pfdhfr gene. Moreover, detection of strong linkage among mutations of Pfcrt and Pfdhfr and feeble evidence for balancing selection in Pfdhps are indicative of evolutionary potential of mutation in genes responsible for drug resistance in Nigerian populations of P. falciparum.

Evolution of Fseg/Cseg dimorphism in region III of the Plasmodium falciparum eba-175 gene

The 175-kDa erythrocyte binding antigen (EBA-175) of the malaria parasite Plasmodium falciparum is important for its invasion into human erythrocytes. The primary structure of eba-175 is divided into seven regions, namely I to VII. Region III contains highly divergent dimorphic segments, termed Fseg and Cseg. The allele frequencies of segmental dimorphism within a P. falciparum population have been extensively examined; however, the molecular evolution of segmental dimorphism is not well understood. A comprehensive comparison of nucleotide sequences among 32 P. falciparum eba-175 alleles identified in our previous study, two Plasmodium reichenowi, and one P. gaboni orthologous alleles obtained from the GenBank database was conducted to uncover the origin and evolutionary processes of segmental dimorphism in P. falciparum eba-175. In the eba-175 nucleotide sequence derived from a P. reichenowi CDC strain, both Fseg and Cseg were found in region III, which implies that the original eba-175 gene had both segments, and deletions of F- and C-segments generated Cseg and Fseg alleles, respectively. We also confirmed the presence of allele with Fseg and Cseg in another P. reichenowi strain (SY57), by re-mapping short reads obtained from the GenBank database. On the other hand, the segmental sequence of eba-175 ortholog in P. gaboni was quite diverged from those of the other species, suggesting that the original eba-175 dimorphism of P. falciparum can be traced back to the stem linage of P. falciparum and P. reichenowi. Our findings suggest that Fseg and Cseg alleles are derived from a single eba-175 allele containing both segments in the ancestral population of P. falciparum and P. reichenowi, and that the allelic dimorphism of eba-175 was shaped by the independent emergence of similar dimorphic lineage in different species that has never been observed in any evolutionary mode of allelic dimorphism at other loci in malaria genomes.

Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum

The malaria parasite Plasmodium falciparum has a great capacity for evolutionary adaptation to evade host immunity and develop drug resistance. Current understanding of parasite evolution is impeded by the fact that a large fraction of the genome is either highly repetitive or highly variable and thus difficult to analyze using short-read sequencing technologies. Here, we describe a resource of deep sequencing data on parents and progeny from genetic crosses, which has enabled us to perform the first genome-wide, integrated analysis of SNP, indel and complex polymorphisms, using Mendelian error rates as an indicator of genotypic accuracy. These data reveal that indels are exceptionally abundant, being more common than SNPs and thus the dominant mode of polymorphism within the core genome. We use the high density of SNP and indel markers to analyze patterns of meiotic recombination, confirming a high rate of crossover events and providing the first estimates for the rate of non-crossover events and the length of conversion tracts. We observe several instances of meiotic recombination within copy number variants associated with drug resistance, demonstrating a mechanism whereby fitness costs associated with resistance mutations could be compensated and greater phenotypic plasticity could be acquired.

Out of Africa: origins and evolution of the human malaria parasites Plasmodium falciparum and Plasmodium vivax

Plasmodium falciparum and Plasmodium vivax account for more than 95% of all human malaria infections, and thus pose a serious public health challenge. To control and potentially eliminate these pathogens, it is important to understand their origins and evolutionary history. Until recently, it was widely believed that P. falciparum had co-evolved with humans (and our ancestors) over millions of years, whilst P. vivax was assumed to have emerged in southeastern Asia following the cross-species transmission of a parasite from a macaque. However, the discovery of a multitude of Plasmodium spp. in chimpanzees and gorillas has refuted these theories and instead revealed that both P. falciparum and P. vivax evolved from parasites infecting wild-living African apes. It is now clear that P. falciparum resulted from a recent cross-species transmission of a parasite from a gorilla, whilst P. vivax emerged from an ancestral stock of parasites that infected chimpanzees, gorillas and humans in Africa, until the spread of the protective Duffy-negative mutation eliminated P. vivax from human populations there. Although many questions remain concerning the biology and zoonotic potential of the P. falciparum- and P. vivax-like parasites infecting apes, comparative genomics, coupled with functional parasite and vector studies, are likely to yield new insights into ape Plasmodium transmission and pathogenesis that are relevant to the treatment and prevention of human malaria.

Merozoite surface protein-1 genetic diversity in Plasmodium malariae and Plasmodium brasilianum from Brazil

Background

The merozoite surface protein 1 (MSP1) gene encodes the major surface antigen of invasive forms of the Plasmodium erythrocytic stages and is considered a candidate vaccine antigen against malaria. Due to its polymorphisms, MSP1 is also useful for strain discrimination and consists of a good genetic marker. Sequence diversity in MSP1 has been analyzed in field isolates of three human parasites: P. falciparum, P. vivax, and P. ovale. However, the extent of variation in another human parasite, P. malariae, remains unknown. This parasite shows widespread, uneven distribution in tropical and subtropical regions throughout South America, Asia, and Africa. Interestingly, it is genetically indistinguishable from P. brasilianum, a parasite known to infect New World monkeys in Central and South America.

Methods

Specific fragments (1 to 5) covering 60 % of the MSP1 gene (mainly the putatively polymorphic regions), were amplified by PCR in isolates of P. malariae and P. brasilianum from different geographic origin and hosts. Sequencing of the PCR-amplified products or cloned PCR fragments was performed and the sequences were used to construct a phylogenetic tree by the maximum likelihood method. Data were computed to give insights into the evolutionary and phylogenetic relationships of these parasites.

Results

Except for fragment 4, sequences from all other fragments consisted of unpublished sequences. The most polymorphic gene region was fragment 2, and in samples where this region lacks polymorphism, all other regions are also identical. The low variability of the P. malariae msp1 sequences of these isolates and the identification of the same haplotype in those collected many years apart at different locations is compatible with a low transmission rate. We also found greater diversity among P. brasilianum isolates compared with P. malariae ones. Lastly, the sequences were segregated according to their geographic origins and hosts, showing a strong genetic and geographic structure.

Conclusions

Our data show that there is a low level of sequence diversity and a possible absence of allelic dimorphism of MSP1 in these parasites as opposed to other Plasmodium species. P. brasilianum strains apparently show greater divergence in comparison to P. malariae, thus P. malariae could derive from P. brasilianum, as it has been proposed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-1238-8) contains supplementary material, which is available to authorized users.

The Plasmodium gaboni genome illuminates allelic dimorphism of immunologically important surface antigens in P. falciparum

In the deadly human malaria parasite Plasmodium falciparum, several major merozoite surface proteins (MSPs) show a striking pattern of allelic diversity called allelic dimorphism (AD). In AD, the vast majority of observed alleles fall into two highly divergent allelic classes, with recombinant alleles being rare or not observed, presumably due to repression by natural selection (recombination suppression, or RS). The three AD loci, merozoite surface proteins (MSPs) 1, 2, and 6, along with MSP3, which also exhibits RS among four allelic classes, can be collectively called AD/RS. The causes of AD/RS and the evolutionary history of allelic diversity at these loci remain mysterious. The few available sequences from a single closely related chimpanzee parasite, P. reichenowi, have suggested that for 3/4 loci, AD/RS is an ancient state that has been retained in P. falciparum since well before the P. falciparum-P. reichenowi ancestor. On the other hand, based on comparative sequence analysis, we recently suggested that (i) AD/RS P. falciparum loci have undergone interallelic recombination over longer evolutionary times (on the timescale of recent speciation events), and thus (ii) AD/RS may be a recent phenomenon. The recent publication of genomic sequencing efforts for P. gaboni, an outgroup to P. falciparum and P. reichenowi, allows for improved reconstruction of the evolutionary history of these loci. In this work, I report genic sequence for P. gaboni for all four AD/RS P. falciparum loci (MSP1, 2, 3, and 6). Comparison of these sequences with available P. falciparum and P. reichenowi data strengthens the evidence for interallelic recombination over the evolutionary history of these species and also strengthens the case that AD/RS at these loci is ancient. Combined with previous results, these data provide evidence that AD/RS at different loci has evolved at several different times in the evolutionary history of P. falciparum: (i) before the P. gaboni-P. falciparum divergence, for much of MSP1 and MSP3; (ii) between the P. gaboni-P. falciparum and P. reichenowi-P. falciparum divergences, for the 5' end of the AD region of MSP6 and block 3 of MSP1; (iii) near the P. reichenowi-P. falciparum divergence, for the 3' end of the AD region of MSP6; and (iv) after the P. reichenowi-P. falciparum divergence, for MSP2. Based on these results, I suggest a new hypothesis for long-term evolutionary maintenance of AD/RS by recombination within allelic groups.

Plasmodium knowlesi genome sequences from clinical isolates reveal extensive genomic dimorphism

Plasmodium knowlesi is a newly described zoonosis that causes malaria in the human population that can be severe and fatal. The study of P. knowlesi parasites from human clinical isolates is relatively new and, in order to obtain maximum information from patient sample collections, we explored the possibility of generating P. knowlesi genome sequences from archived clinical isolates. Our patient sample collection consisted of frozen whole blood samples that contained excessive human DNA contamination and, in that form, were not suitable for parasite genome sequencing. We developed a method to reduce the amount of human DNA in the thawed blood samples in preparation for high throughput parasite genome sequencing using Illumina HiSeq and MiSeq sequencing platforms. Seven of fifteen samples processed had sufficiently pure P. knowlesi DNA for whole genome sequencing. The reads were mapped to the P. knowlesi H strain reference genome and an average mapping of 90% was obtained. Genes with low coverage were removed leaving 4623 genes for subsequent analyses. Previously we identified a DNA sequence dimorphism on a small fragment of the P. knowlesi normocyte binding protein xa gene on chromosome 14. We used the genome data to assemble full-length Pknbpxa sequences and discovered that the dimorphism extended along the gene. An in-house algorithm was developed to detect SNP sites co-associating with the dimorphism. More than half of the P. knowlesi genome was dimorphic, involving genes on all chromosomes and suggesting that two distinct types of P. knowlesi infect the human population in Sarawak, Malaysian Borneo. We use P. knowlesi clinical samples to demonstrate that Plasmodium DNA from archived patient samples can produce high quality genome data. We show that analyses, of even small numbers of difficult clinical malaria isolates, can generate comprehensive genomic information that will improve our understanding of malaria parasite diversity and pathobiology.

A new model for the origins of allelic dimorphism in Plasmodium falciparum

In his landmark 1987 study of the merozoite surface protein-1 locus in Plasmodium falciparum, Kazuyuki Tanabe and coauthors introduced the phenomenon of allelic dimorphism, in which antigenic diversity is arranged into two maximally diverged haplotypes. Further work has extended this finding to other loci in P. falciparum. Each of the loci at which allelic dimorphism is observed encodes major surface antigens of blood-stage malaria parasites, and is consequently a major vaccine target, thus understanding the origins and implications of allelic dimorphism is of crucial importance. Here we examine the essential features of allelic dimorphism in dimorphic malarial surface antigens. From sequence analysis, we conclude that the ancestral population may have been recombining/multimorphic rather than dimorphic. We hypothesize a pathway to allelic dimorphism in which an ancestral allele-rich recombining population could have undergone a severe population bottleneck, putatively caused by the lateral transfer of P. falciparum from apes to humans. This bottleneck produced a reduction in allelic diversity, favoring the survival of the most divergent alleles, which in turn led to recombination suppression by strong natural selection against recombinants.

Hybridization between previously isolated ancestors may explain the persistence of exactly two ancient lineages in the genome of the oyster parasite Perkinsus marinus

Theory predicts that neutral genetic variation accumulates within populations to a level determined by gains through mutation and losses by genetic drift. This balance results in a characteristic distribution of allelic variation with the maximum allelic difference determined by effective population size. Here, we report a striking departure from these expectations in the form of allelic dimorphism, observed at the majority of seven loci examined in Perkinsus marinus, an important oyster parasite that causes Dermo disease. DNA sequences were collected from five loci flanking microsatellite repeats and two loci coding for superoxide dismutase enzymes that may mediate the parasite's interaction with its host. Based on 474 sequences, sampled across 5000 km of the eastern United States coastline, no more than two alleles were observed at each locus (discounting singletons). Depending on the locus, the common allele ranged in overall frequency from 72% to 92%. At each locus the two alleles differed substantially (3.8% sequence difference, on average), and the among-locus variance in divergences was not sufficient to reject a simultaneous origin for all dimorphisms using approximate Bayesian methods. Dimorphic alleles were estimated to have diverged from a common ancestral allele at least 0.9 million years ago. Across these seven loci, only five other alleles were ever observed, always as singletons and differing from the dimorphic alleles by no more than two nucleotides. Free recombination could potentially have shuffled these dimorphisms into as many as 243 multilocus combinations, but the existence of only ten combinations among all samples strongly supports low recombination frequencies and is consistent with the observed absence of intragenic recombination. We consider several demographic and evolutionary hypotheses to explain these patterns. Few can be conclusively rejected with the present data, but we advance a recent hybridization of ancient divergent lineages scenario as the most parsimonious.

Hypervariable antigen genes in malaria have ancient roots

Background

The var genes of the human malaria parasite Plasmodium falciparum are highly polymorphic loci coding for the erythrocyte membrane proteins 1 (PfEMP1), which are responsible for the cytoaherence of P. falciparum infected red blood cells to the human vasculature. Cytoadhesion, coupled with differential expression of var genes, contributes to virulence and allows the parasite to establish chronic infections by evading detection from the host’s immune system. Although studying genetic diversity is a major focus of recent work on the var genes, little is known about the gene family's origin and evolutionary history.

Results

Using a novel hidden Markov model-based approach and var sequences assembled from additional isolates and species, we are able to reveal elements of both the early evolution of the var genes as well as recent diversifying events. We compare sequences of the var gene DBLα domains from divergent isolates of P. falciparum (3D7 and HB3), and a closely-related species, Plasmodium reichenowi. We find that the gene family is equally large in P. reichenowi and P. falciparum -- with a minimum of 51 var genes in the P. reichenowi genome (compared to 61 in 3D7 and a minimum of 48 in HB3). In addition, we are able to define large, continuous blocks of homologous sequence among P. falciparum and P. reichenowi var gene DBLα domains. These results reveal that the contemporary structure of the var gene family was present before the divergence of P. falciparum and P. reichenowi, estimated to be between 2.5 to 6 million years ago. We also reveal that recombination has played an important and traceable role in both the establishment, and the maintenance, of diversity in the sequences.

Conclusions

Despite the remarkable diversity and rapid evolution found in these loci within and among P. falciparum populations, the basic structure of these domains and the gene family is surprisingly old and stable. Revealing a common structure as well as conserved sequence among two species also has implications for developing new primate-parasite models for studying the pathology and immunology of falciparum malaria, and for studying the population genetics of var genes and associated virulence phenotypes.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

Differential transcription of two highly divergent gut-expressed Bm86 antigen gene homologues in the tick Rhipicephalus appendiculatus (Acari: Ixodida)

The transcriptional control of gene expression is not well documented in the Arthropoda. We describe transcriptional analysis of two exceptionally divergent homologues (Ra86) of the Bm86 gut antigen from Rhipicephalus appendiculatus. Bm86 forms the basis of a commercial vaccine for the control of Rhipicephalus (Boophilus) microplus. The R. appendiculatus Ra86 proteins contain 654 and 693 amino acids, with only 80% amino acid sequence identity. Reverse-transcription PCR of gut cDNA showed transcription of only one genotype in individual female ticks. PCR amplification of 3' untranslated sequences from genomic DNA indicated that both variants could be encoded within a single genome. When both variants were present, one of the two Ra86 genotypes was transcriptionally dominant.

Limited variation in vaccine candidate Plasmodium falciparum Merozoite Surface Protein-6 over multiple transmission seasons

BACKGROUND

Plasmodium falciparum Merozoite Surface Protein-6 (PfMSP6) is a component of the complex proteinacious coat that surrounds P. falciparum merozoites. This location, and the presence of anti-PfMSP6 antibodies in P. falciparum-exposed individuals, makes PfMSP6 a potential blood stage vaccine target. However, genetic diversity has proven to be a major hurdle for vaccines targeting other blood stage P. falciparum antigens, and few endemic field studies assessing PfMSP6 gene diversity have been conducted. This study follows PfMSP6 diversity in the Peruvian Amazon from 2003 to 2006 and is the first longitudinal assessment of PfMSP6 sequence dynamics.

METHODS

Parasite DNA was extracted from 506 distinct P. falciparum infections spanning the transmission seasons from 2003 to 2006 as part of the Malaria Immunology and Genetics in the Amazon (MIGIA) cohort study near Iquitos, Peru. PfMSP6 was amplified from each sample using a nested PCR protocol, genotyped for allele class by agarose gel electrophoresis, and sequenced to detect diversity. Allele frequencies were analysed using JMP v.8.0.1.0 and correlated with clinical and epidemiological data collected as part of the MIGIA project.

RESULTS

Both PfMSP6 allele classes, K1-like and 3D7-like, were detected at the study site, confirming that both are globally distributed. Allele frequencies varied significantly between transmission seasons, with 3D7-class alleles dominating and K1-class alleles nearly disappearing in 2005 and 2006. There was a significant association between allele class and village location (p-value = 0.0008), but no statistically significant association between allele class and age, sex, or symptom status. No intra-allele class sequence diversity was detected.

CONCLUSIONS

Both PfMSP6 allele classes are globally distributed, and this study shows that allele frequencies can fluctuate significantly between communities separated by only a few kilometres, and over time in the same community. By contrast, PfMSP6 was highly stable at the sequence level, with no SNPs detected in the 506 samples analysed. This limited diversity supports further investigation of PfMSP6 as a blood stage vaccine candidate, with the clear caveat that any such vaccine must either contain both alleles or generate cross-protective responses that react against both allele classes. Detailed immunoepidemiology studies are needed to establish the viability of these approaches before PfMSP6 advances further down the vaccine development pipeline.

Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates.

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.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7. Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5' and 3' segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition, we uncovered two surprising patterns of evolution in repetitive sequence. First, in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi.

Differential selective pressures on the merozoite surface protein 2 locus of Plasmodium falciparum in a low endemic area

195 Plasmodium falciparum merozoite surface protein-2 alleles collected in Tak Province, Thailand, in 1996 and 2006 revealed extremely limited sequence polymorphism except in the variable (V) region, which defines the two allelic families 3D7 and FC27. This pattern is most easily explained by repeated inter-allelic gene conversion events homogenizing alleles outside the V region. Comparison of synonymous and nonsynonymous differences in V regions within allelic families supported the hypothesis that amino acid sequence polymorphism in this region is selectively favored. The pattern of sequence differentiation supported the hypothesis that repeats in the V region have evolved by concerted evolution in the 3D7 family but not in the FC27 family. In the FC27 family two alleles of relatively high frequency were the most common V-region alleles in both 1996 and 2006, while 3D7 alleles constituted a significantly greater proportion of the sequences collected in 2006 (56.1%) than of those collected in 1996 (28.9%). These changes in the frequencies of 3D7 alleles may reflect increased intensity of selection on the P. falciparum population in Thailand as a result of effective control measures that have sharply reduced the incidence of malaria infection.

Diversity and evolution of the rhoph1/clag multigene family of Plasmodium falciparum

A complex of high-molecular-mass proteins (PfRhopH) of the human malaria parasite Plasmodium falciparum induces host protective immunity and therefore is a candidate for vaccine development. Understanding the level of polymorphism and the evolutionary processes is important for advancements in both vaccine design and knowledge of the evolution of cell invasion in this parasite. In the present study, we sequenced the entire open reading frames of seven genes encoding the proteins of the PfRhopH complex (rhoph2, rhoph3, and five rhoph1/clag gene paralogs). We found that four rhoph1/clag genes (clag2, 3.1, 3.2, and 8) were highly polymorphic. Amino acid substitutions and indels are predominantly clustered around amino acid positions 1000-1200 of these four rhoph1/clag genes. An excess of nonsynonymous substitutions over synonymous substitutions was detected for clag8 and 9, indicating positive selection. The McDonald-Kreitman test with a Plasmodium reichenowi orthologous sequence also supports positive selection on clag8. Based on the ratio of interspecific genetic distance to intraspecific distance, the time to the most recent common ancestor of the clag2 and 8 polymorphisms was estimated to be 1.89 and 0.87 million years ago, respectively, assuming divergence of P. falciparum and P. reichenowi 6 million years ago. In addition to a copy number polymorphism, gene conversion events were detected for the rhoph1/clag genes on chromosome 3, which likely play a role in increasing the diversity of each locus. Our results indicate that a high diversity of the PfRhopH1/Clag multigene family is maintained by diversifying selection forces over a considerably long period.