Orthologous gene sequences of merozoite surface protein 1 (MSP1) from Plasmodium reichenowi and P. gallinaceum confirm an ancient divergence of P. falciparum alleles

Plasmodium relictum MSP-1 capture antigen-based ELISA for detection of avian malaria antibodies in African penguins (Spheniscus demersus)

Avian malaria, caused by Plasmodium spp. and transmitted by mosquitos, is a leading cause of mortality of captive penguins. Antimalarial drugs are currently used to control infections in penguins. However, the effectiveness of treatment reduces significantly by the time the clinical signs appear, while early and unnecessary treatment interferes with development of protective immunity. Therefore, for suppressing parasitemia without affecting the development of immunity in captive penguins, antimalaria drugs need to be administered at the right time, which requires reliable diagnostic tools that can determine the levels of circulating antimalaria antibodies. In the present study, we have developed an enzyme-linked immunosorbent assay (ELISA) diagnostic assay based on the merozoite surface protein 1 (MSP-1) of P. relictum isolate SGS1 to specifically detect and relatively quantify antimalaria antibodies in penguins. We expressed and purified a truncated P. relictum isolate SGS1 MSP-1 and optimized its biotinylation and subsequent conjugation to streptavidin alkaline phosphatase for signal generation in ELISA. We tested the assay by analyzing sera obtained from penguins at the Baltimore Zoo, from Spring through Fall, and found that levels of detectable antibodies against MSP-1 varied seasonally for individual penguins, consistent with the expected seasonal variations in avian malaria prevalence. Corroboratively, we analyzed the sensitivity of the assay by titrating positive sera and found that the signal intensity generated was serum concentration-dependent, thus validating the ability of the assay to detect and relatively quantify the levels of antimalaria antibodies in penguin sera.

•

ELISA based on MSP1 for detection and quantification of antibodies against Plasmodium relictum in birds was developed.

•

Assay was validated to detect and quantify levels of antimalaria antibodies in infected penguins' sera.

•

Assay detected varied antibody levels against MSP-1 in penguin sera consistent with seasonal variations in malaria prevalence.

Classification for avian malaria parasite Plasmodium gallinaceum blood stages by using deep convolutional neural networks

The infection of an avian malaria parasite (Plasmodium gallinaceum) in domestic chickens presents a major threat to the poultry industry because it causes economic loss in both the quality and quantity of meat and egg production. Computer-aided diagnosis has been developed to automatically identify avian malaria infections and classify the blood infection stage development. In this study, four types of deep convolutional neural networks, namely Darknet, Darknet19, Darknet19-448 and Densenet201 are used to classify P. gallinaceum blood stages. We randomly collected a dataset of 12,761 single-cell images consisting of three parasite stages from ten-infected blood films stained by Giemsa. All images were confirmed by three well-trained examiners. The study mainly compared several image classification models and used both qualitative and quantitative data for the evaluation of the proposed models. In the model-wise comparison, the four neural network models gave us high values with a mean average accuracy of at least 97%. The Darknet can reproduce a superior performance in the classification of the P. gallinaceum development stages across any other model architectures. Furthermore, the Darknet has the best performance in multiple class-wise classification, with average values of greater than 99% in accuracy, specificity, and sensitivity. It also has a low misclassification rate (< 1%) than the other three models. Therefore, the model is more suitable in the classification of P. gallinaceum blood stages. The findings could help us create a fast-screening method to help non-experts in field studies where there is a lack of specialized instruments for avian malaria diagnostics.

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.

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.

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.

Plasmodium kinases as targets for new-generation antimalarials

There is an urgent need for the development of new antimalarial drugs with novel modes of actions. The malarial parasite, Plasmodium falciparum, has a relatively small kinome of <100 kinases, with many members exhibiting a high degree of structural divergence from their host counterparts. A number of Plasmodium kinases have recently been shown by reverse genetics to be essential for various parts of the complex parasitic life cycle, and are thus genetically validated as potential targets. Implementation of mass spectrometry-based phosphoproteomics approaches has informed on key phospho-signalling pathways in the parasite. In addition, global phenotypic screens have revealed a large number of putative protein kinase inhibitors with antimalarial potency. Taken together, these investigations point to the Plasmodium kinome as a rich source of potential new targets. In this review, we highlight recent progress made towards this goal.

Low level of sequence diversity at merozoite surface protein-1 locus of Plasmodium ovale curtisi and P. ovale wallikeri from Thai isolates

Background

The merozoite surface protein-1 (MSP-1) is a candidate target for the development of blood stage vaccines against malaria. Polymorphism in MSP-1 can be useful as a genetic marker for strain differentiation in malarial parasites. Although sequence diversity in the MSP-1 locus has been extensively analyzed in field isolates of Plasmodium falciparum and P. vivax, the extent of variation in its homologues in P. ovale curtisi and P. ovale wallikeri, remains unknown.

Methodology/Principal Findings

Analysis of the mitochondrial cytochrome b sequences of 10 P. ovale isolates from symptomatic malaria patients from diverse endemic areas of Thailand revealed co-existence of P. ovale curtisi (n = 5) and P. ovale wallikeri (n = 5). Direct sequencing of the PCR-amplified products encompassing the entire coding region of MSP-1 of P. ovale curtisi (PocMSP-1) and P. ovale wallikeri (PowMSP-1) has identified 3 imperfect repeated segments in the former and one in the latter. Most amino acid differences between these proteins were located in the interspecies variable domains of malarial MSP-1. Synonymous nucleotide diversity (π_S) exceeded nonsynonymous nucleotide diversity (π_N) for both PocMSP-1 and PowMSP-1, albeit at a non-significant level. However, when MSP-1 of both these species was considered together, π_S was significantly greater than π_N (p<0.0001), suggesting that purifying selection has shaped diversity at this locus prior to speciation. Phylogenetic analysis based on conserved domains has placed PocMSP-1 and PowMSP-1 in a distinct bifurcating branch that probably diverged from each other around 4.5 million years ago.

Conclusion/Significance

The MSP-1 sequences support that P. ovale curtisi and P. ovale wallikeri are distinct species. Both species are sympatric in Thailand. The low level of sequence diversity in PocMSP-1 and PowMSP-1 among Thai isolates could stem from persistent low prevalence of these species, limiting the chance of outcrossing at this locus.

Plasmodium falciparum merozoite surface protein 1 blocks the proinflammatory protein S100P

The malaria parasite, Plasmodium falciparum, and the human immune system have coevolved to ensure that the parasite is not eliminated and reinfection is not resisted. This relationship is likely mediated through a myriad of host-parasite interactions, although surprisingly few such interactions have been identified. Here we show that the 33-kDa fragment of P. falciparum merozoite surface protein 1 (MSP1(33)), an abundant protein that is shed during red blood cell invasion, binds to the proinflammatory protein, S100P. MSP1(33) blocks S100P-induced NFκB activation in monocytes and chemotaxis in neutrophils. Remarkably, S100P binds to both dimorphic alleles of MSP1, estimated to have diverged >27 Mya, suggesting an ancient, conserved relationship between these parasite and host proteins that may serve to attenuate potentially damaging inflammatory responses.

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.

Geographic differentiation of polymorphism in the Plasmodium falciparum malaria vaccine candidate gene SERA5

SERA5 is regarded as a promising malaria vaccine candidate of the most virulent human malaria parasite Plasmodium falciparum. SERA5 is a 120 kDa abundantly expressed blood-stage protein containing a papain-like protease. Since substantial polymorphism in blood-stage vaccine candidates may potentially limit their efficacy, it is imperative to fully investigate polymorphism of the SERA5 gene (sera5). In this study, we performed evolutionary and population genetic analysis of sera5. The level of inter-species divergence (kS=0.076) between P. falciparum and Plasmodium reichenowi, a closely related chimpanzee malaria parasite is comparable to that of housekeeping protein genes. A signature of purifying selection was detected in the proenzyme and enzyme domains. Analysis of 445 near full-length P. falciparum sera5 sequences from nine countries in Africa, Southeast Asia, Oceania and South America revealed extensive variations in the number of octamer repeat (OR) and serine repeat (SR) regions as well as substantial level of single nucleotide polymorphism (SNP) in non-repeat regions (2562 bp). Remarkably, a 14 amino acid sequence of SERA5 (amino acids 59-72) that is known to be the in vitro target of parasite growth inhibitory antibodies was found to be perfectly conserved in all 445 worldwide isolates of P. falciparum evaluated. Unlike other major vaccine target antigen genes such as merozoite surface protein-1, apical membrane antigen-1 or circumsporozoite protein, no strong evidence for positive selection was detected for SNPs in the non-repeat regions of sera5. A biased geographical distribution was observed in SNPs as well as in the haplotypes of the sera5 OR and SR regions. In Africa, OR- and SR-haplotypes with low frequency (<5%) and SNPs with minor allele frequency (<5%) were abundant and were mostly continent-specific. Consistently, significant genetic differentiation, assessed by the Wright's fixation index (Fst) of inter-population variance in allele frequencies, was detected for SNPs and both OR- and SR-haplotypes among almost all parasite populations. The exception was parasite populations between Tanzania and Ghana, suggesting frequent gene flow in Africa. The present study points to the importance of investigating whether biased geographical distribution for SNPs and repeat variants in the OR and SR regions affect the reactivity of human serum antibodies to variants.

Genomic diversity of the human intestinal parasite Entamoeba histolytica

Background

Results

Conclusions

In Silico Prediction of Evolutionarily Conserved GC-Rich Elements Associated with Antigenic Proteins of Plasmodium falciparum

The Plasmodium falciparum genome being AT-rich, the presence of GC-rich regions suggests functional significance. Evolution imposes selection pressure to retain functionally important coding and regulatory elements. Hence searching for evolutionarily conserved GC-rich, intergenic regions in an AT-rich genome will help in discovering new coding regions and regulatory elements. We have used elevated GC content in intergenic regions coupled with sequence conservation against P. reichenowi, which is evolutionarily closely related to P. falciparum to identify potential sequences of functional importance. Interestingly, ~30% of the GC-rich, conserved sequences were associated with antigenic proteins encoded by var and rifin genes. The majority of sequences identified in the 5′ UTR of var genes are represented by short expressed sequence tags (ESTs) in cDNA libraries signifying that they are transcribed in the parasite. Additionally, 19 sequences were located in the 3′ UTR of rifins and 4 also have overlapping ESTs. Further analysis showed that several sequences associated with var genes have the capacity to encode small peptides. A previous report has shown that upstream peptides can regulate the expression of var genes hence we propose that these conserved GC-rich sequences may play roles in regulation of gene expression.

The evolution of virulence in primate malaria parasites based on Bayesian reconstructions of ancestral states

Plasmodium parasites, the causative agents of malaria, are generally considered as harmful parasites, but many of them cause mild symptoms. Little is known about the evolutionary history and phylogenetic constraints that generate this interspecific variation in virulence due to uncertainties about the phylogenetic associations of parasites. Here, to account for such phylogenetic uncertainty, phylogenetic methods based on Bayesian statistics were followed in combination with sequence data from five genes to estimate the ancestral state of virulence in primate Plasmodium parasites. When recent parasites were categorised according to the damage caused to the host, Bayesian estimates of ancestral states indicated that the acquisition of a harmful host exploitation strategy is more likely to be a recent evolutionary event than a result of an ancient change in a character state altering virulence. On the contrary, there was more evidence for moderate host exploitation having a deep origin along the phylogenetic tree. Moreover, the evolution of host severity is determined by the phylogenetic relationships of parasites, as severity gains did not appear randomly on the evolutionary tree. Such phylogenetic constraints can be mediated by the acquisition of virulence genes. As the impact of a parasite on a host is the result of both the parasite's investment in reproduction and host sensitivity, virulence was also estimated by calculating peak parasitemia after eliminating host effects. A directional random-walk evolutionary model showed that the ancestral primate malarias reproduced at very low parasitemia in their hosts. Consequently, the extreme variation in the outcome of malaria infection in different host species can be better understood in light of the phylogeny of parasites.

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

Background

The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.

Results

We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1.

Conclusions

The present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.

Patterns of co-speciation and host switching in primate malaria parasites

Background

The evolutionary history of many parasites is dependent on the evolution of their hosts, leading to an association between host and parasite phylogenies. However, frequent host switches across broad phylogenetic distances may weaken this close evolutionary link, especially when vectors are involved in parasites transmission, as is the case for malaria pathogens. Several studies suggested that the evolution of the primate-infective malaria lineages may be constrained by the phylogenetic relationships of their hosts, and that lateral switches between distantly related hosts may have been occurred. However, no systematic analysis has been quantified the degree of phylogenetic association between primates and their malaria parasites.

Methods

Here phylogenetic approaches have been used to discriminate statistically between events due to co-divergence, duplication, extinction and host switches that can potentially cause historical association between Plasmodium parasites and their primate hosts. A Bayesian reconstruction of parasite phylogeny based on genetic information for six genes served as basis for the analyses, which could account for uncertainties about the evolutionary hypotheses of malaria parasites.

Results

Related lineages of primate-infective Plasmodium tend to infect hosts within the same taxonomic family. Different analyses testing for congruence between host and parasite phylogenies unanimously revealed a significant association between the corresponding evolutionary trees. The most important factor that resulted in this association was host switching, but depending on the parasite phylogeny considered, co-speciation and duplication may have also played some additional role. Sorting seemed to be a relatively infrequent event, and can occur only under extreme co-evolutionary scenarios. The concordance between host and parasite phylogenies is heterogeneous: while the evolution of some malaria pathogens is strongly dependent on the phylogenetic history of their primate hosts, the congruent evolution is less emphasized for other parasite lineages (e.g. for human malaria parasites). Estimation of ancestral states of host use along the phylogenetic tree of parasites revealed that lateral transfers across distantly related hosts were likely to occur in several cases. Parasites cannot infect all available hosts, and they should preferentially infect hosts that provide a similar environment for reproduction. Marginally significant evidence suggested that there might be a consistent variation within host ranges in terms of physiology.

Conclusion

The evolution of primate malarias is constrained by the phylogenetic associations of their hosts. Some parasites can preserve a great flexibility to infect hosts across a large phylogenetic distance, thus host switching can be an important factor in mediating host ranges observed in nature. Due to this inherent flexibility and the potential exposure to various vectors, the emergence of new malaria disease in primates including humans cannot be predicted from the phylogeny of parasites.

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.

Strain-specific immunity may drive adaptive polymorphism in the merozoite surface protein 1 of the rodent malaria parasite Plasmodium chabaudi

Clinical immunity against malaria is slow to develop, poorly understood and strongly strain-specific. Understanding how strain-specific immunity develops and identifying the parasite antigens involved is crucial to developing effective vaccines against the disease. In previous experiments we have shown that strain-specific protective immunity (SSPI) exists between genetically distinct strains (cloned lines) of the rodent malaria parasite Plasmodium chabaudi chabaudi in mice [Cheesman, S., Raza, A., Carter, R., 2006. Mixed strain infections and strain-specific protective immunity in the rodent malaria parasite P. chabaudi chabaudi in mice. Infect. Immun. 74, 2996-3001]. In two subsequent studies, we identified the highly polymorphic Merozoite Surface Protein 1 (MSP-1) as being the principal candidate molecule for the control of SSPI against P. c. chabaudi malaria [Martinelli et al., 2005; Pattaradilokrat, S., Cheesman, S.J., Carter R., 2007. Linkage group selection: towards identifying genes controlling strain-specific protective immunity in malaria. PLoS ONE 2(9):e857]. In the present study, we sequenced the whole msp1 gene of several genetically distinct strains of P. chabaudi and found high levels of genetic diversity. Protein sequence alignments reveal extensive allelic polymorphism between the P. chabaudi strains, concentrated primarily within five regions of the protein. The 3'-end sequence region, encoding the C-terminal 21 kDa region (MSP-1(21)), which is analogous and homologous to MSP-1(19) of Plasmodium falciparum, appears to have been subject to balancing selection. We have found that the strains with the lowest sequence identity at MSP-1(21) (i.e. AS/CB and AJ/CB) induce robust and reciprocal SSPI in experimental mice. In contrast, two strains that do not induce reciprocal SSPI are identical at the 21 kDa region. Final identification of the region(s) controlling SSPI will provide important information to help guide decisions about MSP-1 based vaccines.

Patterns of polymorphism in genomic regions flanking three highly polymorphic surface antigens in Plasmodium falciparum

Many surface antigens of the human malaria parasite Plasmodium falciparum show extraordinary diversity, with different alleles being so divergent as to be unalignable in some coding regions. To better understand the population history and modes of selection on such loci, we sequenced genomic regions flanking the highly polymorphic genes merozoite surface protein-1, merozoite surface protein-2, and circumsporozoite protein, from reference isolates of P. falciparum. Diversity was much lower in genomic flanking regions than in the coding sequences. Average pairwise nucleotide diversity for these regions was 0.00088, similar to other genomic regions not thought to be evolving under balancing selection, suggesting against balancing selection acting on promoter regions of these genes. Most observed polymorphisms were singletons. A higher ratio of SNPs to indels than previously reported for P. falciparum was observed. An 11 bp repeat upstream of msp2 showed an intriguing pattern of polymorphism possibly suggestive of purifying selection on total allele length.

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.

Recent independent evolution of msp1 polymorphism in Plasmodium vivax and related simian malaria parasites

The Plasmodium MSP-1 is a promising malaria vaccine candidate. However, the highly polymorphic nature of the MSP-1 gene (msp1) presents a potential obstacle for effective vaccine development. To investigate the evolutionary history of msp1 polymorphism in P. vivax, we construct phylogenetic trees of msp1 from P. vivax and related monkey malaria parasite species. All P. vivax msp1 alleles cluster in the P. vivax lineage and are not distributed among other species. Similarly, all P. cynomolgi msp1 alleles cluster in the P. cynomolgi lineage. This suggests that, in contrast to presumed ancient origin of P. falciparum msp1 polymorphism, the origin of P. vivax msp1 polymorphism is relatively recent. We observed positive selection in the P. vivax lineage but not in P. cynomolgi. Also, positive selection acts on different regions of msp1 in P. vivax and P. falciparum. This study shows that the evolutionary history of msp1 differs greatly among parasite lineages.

Functional characterization of both MAP kinases of the human malaria parasite Plasmodium falciparum by reverse genetics

The kinome of the human malaria parasite Plasmodium falciparum includes two genes encoding mitogen-activated protein kinase (MAPK) homologues, pfmap-1 and pfmap-2, but no clear orthologue of the MAPK kinase (MAPKK) family, raising the question of the mode of activation and function of the plasmodial MAPKs. Functional studies in the rodent malaria model Plasmodium berghei recently showed the map-2 gene to be dispensable for asexual growth and gametocytogenesis, but essential for male gametogenesis in the mosquito vector. Here, we demonstrate by using a reverse genetics approach that the map-2 gene is essential for completion of the asexual cycle of P. falciparum, an unexpected result in view of the non-essentiality of the orthologous gene for P. berghei erythrocytic schizogony. This validates Pfmap-2 as a potential target for chemotherapeutic intervention. In contrast, the other P. falciparum MAPK, Pfmap-1, is required neither for in vitro schizogony and gametocytogenesis in erythrocytes, nor for gametogenesis and sporogony in the mosquito vector. However, Pfmap-2 protein levels are elevated in pfmap-1(-) parasites, suggesting that Pfmap-1 fulfils an important function in asexual parasites that necessitates compensatory adaptation in parasites lacking this enzyme.

Allelic dimorphism-associated restriction of recombination in Plasmodium falciparum msp1

Allelic dimorphism is a characteristic feature of the Plasmodium falciparum msp1 gene encoding the merozoite surface protein 1, a strong malaria vaccine candidate. Meiotic recombination is a major mechanism for the generation of msp1 allelic diversity. Potential recombination sites have previously been mapped to specific regions within msp1 (a 5' 1-kb region and a 3' 0.4-kb region) with no evidence for recombination events in a central 3.5-kb region. However, evidence for the lack of recombination events is circumstantial and inconclusive because the number of msp1 sequences analysed is limited, and the frequency of recombination events has not been addressed previously in a high transmission area, where the frequency of meiotic recombination is expected to be high. In the present study, we have mapped potential allelic recombination sites in 34 full-length msp1 sequences, including 24 new sequences, from various geographic origins. We also investigated recombination events in blocks 6 to 16 by population genetic analysis of P. falciparum populations in Tanzania, where malaria transmission is intense. The results clearly provide no evidence of recombination events occurring between the two major msp1 allelic types, K1-type and Mad20-type, in the central region, but do show recombination events occurring throughout the entire gene within sequences of the Mad20-type. Thus, the present study indicates that allelic dimorphism of msp1 greatly affects inter-allelic recombination events, highlighting a unique feature of allelic diversity of P. falciparum msp1.

Evolution of allelic dimorphism in malarial surface antigens

The extensive sequence variation in most surface antigens of Plasmodium falciparum is one of the major factors why clinical immunity to malaria develops only after repeated infections with the same species over several years. For some P. falciparum surface antigens, all observed alleles clearly fall into two allelic classes, with divergence between classes dwarfing divergence within classes. We discuss the ways in which such allelic dimorphism deviates from the expected shape of the genealogy of genes under either neutral evolution or standard balancing selection, and present a simple test, based on coalescent theory, to detect this deviation in samples of DNA sequences. We review previous hypotheses for the origin and evolution of allelic dimorphism in malarial antigens and discuss the difficulties of explaining the available data under these proposals. We conclude by offering several possible classes of explanations for allelic dimorphism, which are worthy of further theoretical and empirical exploration.

AA. A comparative study of the genetic diversity of the 42kDa fragment of the merozoite surface protein 1 in Plasmodium falciparum and P. vivax

We investigated the genetic diversity of the 42kDa fragment of the merozoite surface protein 1 (MSP-1) antigen in Plasmodium falciparum and P. vivax, as well as in non-human primate malarial parasites. This fragment undergoes a proteolytic cleavage generating two fragments of 19kDa (MSP-1(19)) and 33kDa (MSP-1(33)) that are critical in erythrocyte invasion. We found that overall the MSP-1(33) fragment exhibits greater genetic diversity than the MSP-1(19) regardless of the species. We have found evidence for positive natural selection only in the human malaria parasites by comparing the rate of non-synonymous versus synonymous substitutions. In addition, we found clear differences between the two major human malaria parasites. In the case of P. falciparum, positive natural selection is acting on the MSP-1(19) region while the MSP-1(33) is neutral or under purifying selection. The opposite pattern was observed in P. vivax. Our results suggest different roles of this antigen in the host-parasite immune interaction in each of the major human malarial parasites.

Plasmodium falciparum: worldwide sequence diversity and evolution of the malaria vaccine candidate merozoite surface protein-2 (MSP-2)

We examined patterns and putative mechanisms of sequence diversification in the merozoite surface protein-2 (MSP-2) of Plasmodium falciparum, a major dimorphic malaria vaccine candidate antigen, by analyzing 448 msp-2 alleles from all continents. We describe several nucleotide replacements, insertion and deletion events, frameshift mutations, and proliferations of repeat units that generate the extraordinary diversity found in msp-2 alleles. We discuss the role of positive selection exerted by naturally acquired type- and variant-specific immunity in maintaining the observed levels of polymorphism and suggest that this is the most likely explanation for the significant excess of nonsynonymous nucleotide replacements found in dimorphic msp-2 domains. Hybrid sequences created by meiotic recombination between alleles of different dimorphic types were observed in few (3.1%) isolates, mostly from Africa. We found no evidence for an extremely ancient origin of allelic dimorphism at the msp-2 locus, predating P. falciparum speciation, in contrast with recent findings for other surface malarial antigens.