Genetic Diversity and Selection at the Plasmodium vivax Apical Membrane Antigen-1 (PvAMA-1) Locus in a Sri Lankan Population

Genetic diversity and natural selection of Plasmodium vivax reticulocyte invasion genes in Ecuador

BACKGROUND

Knowledge of the diversity of invasion ligands in malaria parasites in endemic regions is essential to understand how natural selection influences genetic diversity of these ligands and their feasibility as possible targets for future vaccine development. In this study the diversity of four genes for merozoite invasion ligands was studied in Ecuadorian isolates of Plasmodium vivax.

METHODS

Eighty-eight samples from P. vivax infected individuals from the Coast and Amazon region of Ecuador were obtained between 2012 and 2015. The merozoite invasion genes pvmsp-1-19, pvdbpII, pvrbp1a-2 and pvama1 were amplified, sequenced, and compared to the Sal-1 strain. Polymorphisms were mapped and genetic relationships between haplotypes were determined.

RESULTS

Only one nonsynonymous polymorphism was detected in pvmsp-1-19, while 44 nonsynonymous polymorphisms were detected in pvdbpII, 56 in pvrbp1a-2 and 33 in pvama1. While haplotypes appeared to be more related within each area of study and there was less relationship between parasites of the coastal and Amazon regions of the country, diversification processes were observed in the two Amazon regions. The highest haplotypic diversity for most genes occurred in the East Amazon of the country. The high diversity observed in Ecuadorian samples is closer to Brazilian and Venezuelan isolates, but lower than reported in other endemic regions. In addition, departure from neutrality was observed in Ecuadorian pvama1. Polymorphisms for pvdbpII and pvama1 were associated to B-cell epitopes.

CONCLUSIONS

pvdbpII and pvama1 genetic diversity found in Ecuadorian P. vivax was very similar to that encountered in other malaria endemic countries with varying transmission levels and segregated by geographic region. The highest diversity of P. vivax invasion genes in Ecuador was found in the Amazonian region. Although selection appeared to have small effect on pvdbpII and pvrbp1a-2, pvama1 was influenced by significant balancing selection.

Population genetic analyses inferred a limited genetic diversity across the pvama-1 DI domain among Plasmodium vivax isolates from Khyber Pakhtunkhwa regions of Pakistan

Background

Plasmodium vivax apical membrane antigen-1 (pvama-1) is an important vaccine candidate against Malaria. The genetic composition assessment of pvama-1 from wide-range geography is vital to plan the antigen based vaccine designing against Malaria.

Methods

The blood samples were collected from 84 P. vivax positive malaria patients from different districts of Khyber Pakhtunkhwa (KP) province of Pakistan. The highly polymorphic and immunogenic domain-I (DI) region of pvama-1 was PCR amplified and DNA sequenced. The QC based sequences raw data filtration was done using DNASTAR package. The downstream population genetic analyses were performed using MEGA4, DnaSP, Arlequin v3.5 and Network.5 resources.

Results

The analyses unveiled total 57 haplotypes of pvama-1 (DI) in KP samples with majorly prevalent H-14 and H-5 haplotypes. Pairwise comparative population genetics analyses identified limited to moderate genetic distinctions among the samples collected from different districts of KP, Pakistan. In context of worldwide available data, the KP samples depicted major genetic differentiation against the Korean samples with Fst = 0.40915 (P-value = 0.0001), while least distinction was observed against Indian and Iranian samples. The statistically significant negative values of Fu and Li’s D* and F* tests indicate the evidence of population expansion and directional positive selection signature. The slow LD decay across the nucleotide distance in KP isolates indicates low nucleotide diversity. In context of reference pvama-1 sequence, the KP samples were identified to have 09 novel non-synonymous single nucleotide polymorphisms (nsSNPs), including several trimorphic and tetramorphic substitutions. Few of these nsSNPs are mapped within the B-cell predicted epitopic motifs of the pvama-1, and possibly modulate the immune response mechanism.

Conclusion

Low genetic differentiation was observed across the pvama-1 DI among the P. vivax isolates acquired from widespread regions of KP province of Pakistan. The information may implicate in future vaccine designing strategies based on antigenic features of pvama-1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07798-1.

Genetic Diversity and Selection of Plasmodium vivax Apical Membrane Antigen-1 in China-Myanmar Border of Yunnan Province, China, 2009-2016

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is an important vaccine candidate for vivax malaria. However, antigenic variation within PvAMA-1 is a major obstacle to the design of a global protective malaria vaccine. In this study, we analyzed the genetic polymorphism and selection of the PvAMA-1 gene from 152 P. vivax isolates from imported cases to China, collected in the China–Myanmar border (CMB) area in Yunnan Province (YP) during 2009–2011 (n = 71) and 2014–2016 (n = 81), in comparison with PvAMA-1 gene information from Myanmar (n = 73), collected from public data. The overall nucleotide diversity of the PvAMA-1 gene from the 152 YP isolates was 0.007 with 76 haplotypes identified (Hd = 0.958). Results from the population structure suggested three groups among the YP and Myanmar isolates with optimized clusters value of K = 7. In addition, YP (2014–2016) isolates generally lacked some K components that were commonly found in YP (2009–2011) and Myanmar. Meanwhile, PvAMA-1 domain I is found to be the dominant target of positive diversifying selection and most mutation loci were found in this domain. The mutation frequencies of D107N/A, R112K/T, K120R, E145A, E277K, and R438H in PvAMA-1 were more than 70% in the YP isolates. In conclusion, high genetic diversity and positive selection were found in the PvAMA-1 gene from YP isolates, which are significant findings for the design and development of PvAMA-1-based malaria vaccine.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Antibodies Against the Plasmodium vivax Apical Membrane Antigen 1 From the Belem Strain Share Common Epitopes Among Other Worldwide Variants

Malaria is a human parasitic disease distributed in many tropical countries and caused by various Plasmodium species. Plasmodium vivax has the largest geographical distribution of the Plasmodium species and is predominant in the Americas, including Brazil. Only a small number of P. vivax vaccine formulations have successfully reached clinical trials relative to their P. falciparum counterparts. One of the candidate antigens for a blood-stage P. vivax vaccine is apical membrane antigen 1 (PvAMA-1). Due to the worldwide distribution of Plasmodium parasites, a high degree of variability has been detected in this antigen sequence, representing a considerable challenge to the development of a universal vaccine against malaria. In this study, we evaluated how PvAMA-1 polymorphisms influence vaccine-derived immune responses in P. vivax malaria. To this end, we expressed 9 recombinant protein representatives of different PvAMA-1 allelic variants in the yeast Pichia pastoris: Belem, Chesson I, Sal-1, Indonesia XIX, SK0814, TC103, PNG_05_ESP, PNG_62_MU, and PNG_68_MAS. After protein expression and purification, we evaluated the breadth of the immune responses derived from malaria-exposed individuals from the Amazon region. From 611 serum samples of malaria-exposed individuals, 53.68% of them reacted against the PvAMA-1 Belem through ELISA. Positive samples were further tested against recombinant proteins representing the other PvAMA-1 allelic variants. Whereas Sal-1, Chesson I and SK0814 variants were highly recognized by tested serum samples, Indonesia XIX, TC103, PNG_05_ESP, PNG_62_MU, and PNG_68_MAS were only slightly recognized. Moreover, polyclonal sera derived from C57BL/6 mice immunized with the PvAMA-1 Belem protein predominantly recognized Belem, Sal-1, Chesson I, SK0814, and Indonesia XIX through ELISA. Last, ELISA-based competition assays demonstrated that a previous interaction between anti-Belem polyclonal serum and Sal-1, Chesson I, SK0814, or Indonesia XIX proteins could further inhibit antibody binding to the Belem variant. Our human and mouse data suggest the presence of common epitopes or cross-reactivity between Belem, Sal-1, Chesson I, and SK0814 variants. Although the PvAMA-1 Belem variant induces strain-transcendent antibodies, PvAMA-1 variants from Thailand and Papua New Guinea may need to be included in a universal vaccine formulation to achieve protection against P. vivax malaria.

Plasmodium vivax AMA1: Implications of distinct haplotypes for immune response

In Brazil, Plasmodium vivax infection accounts for around 80% of malaria cases. This infection has a substantial impact on the productivity of the local population as the course of the disease is usually prolonged and the development of acquired immunity in endemic areas takes several years. The recent emergence of drug-resistant strains has intensified research on alternative control methods such as vaccines. There is currently no effective available vaccine against malaria; however, numerous candidates have been studied in the past several years. One of the leading candidates is apical membrane antigen 1 (AMA1). This protein is involved in the invasion of Apicomplexa parasites into host cells, participating in the formation of a moving junction. Understanding how the genetic diversity of an antigen influences the immune response is highly important for vaccine development. In this study, we analyzed the diversity of AMA1 from Brazilian P. vivax isolates and 19 haplotypes of P. vivax were found. Among those sequences, 33 nonsynonymous PvAMA1 amino acid sites were identified, whereas 20 of these sites were determined to be located in predicted B-cell epitopes. Nonsynonymous mutations were evaluated for their influence on the immune recognition of these antigens. Two distinct haplotypes, 5 and 16, were expressed and evaluated for reactivity in individuals from northern Brazil. Both PvAMA1 variants were reactive. Moreover, the IgG antibody response to these two PvAMA1 variants was analyzed in an exposed but noninfected population from a P. vivax endemic area. Interestingly, over 40% of this population had antibodies recognizing both variants. These results have implications for the design of a vaccine based on a polymorphic antigen.

Plasmodium vivax is the most abundant Plasmodium species in Brazil. While this species has been neglected for many years, the recent emergence of drug-resistant strains and the absence of a vaccine intensified the efforts for a better control method. Naturally acquired immune response analysis is a useful tool for understanding the antigenicity of Plasmodium proteins and evaluating the potential of a vaccine candidate. In this study, the genetic variability of one of the leading P. vivax vaccine candidates (PvAMA1) was analyzed. Two distinct variants were expressed and the antibody response was evaluated in infected and noninfected individuals in the Brazilian Amazon. This improved understanding of the magnitude and dynamics of the antibody response will contribute to the knowledge of a vaccine candidate and open new perspectives in vivax malaria vaccine development.

Population structure of avian malaria parasites

The geographic distribution of genetic diversity in malaria parasite populations (Apicomplexa: Haemosporida) presumably influences local patterns of virulence and the evolution of host-resistance, but little is known about population genetic structure in these parasites. We assess the distribution of genetic diversity in the partial Domain I of apical membrane antigen 1 (AMA1) in three mtDNA-defined lineages of avian Plasmodium to determine spatial population structure and host-parasite genetic relationships. We find that one parasite lineage is genetically differentiated in association with a single host genus and among some locations, but not with respect to other hosts. Two other parasite lineages are undifferentiated with respect to host species but exhibit geographic differentiation that is inconsistent with shared geographic barriers or with isolation-by-distance. Additional differentiation within two other lineages is unassociated with host species or location; in one case, we tentatively interpret this differentiation as the result of mitochondrial introgression from one of the lineages into a second lineage. More sampling of nuclear genetic diversity within populations of avian Plasmodium is needed to rule out coinfection as a possible confounding factor. If coinfections are not responsible for these findings, further assessment is needed to determine the frequency of mitonuclear discordance and its implications for defining parasite lineages based on mitochondrial genetic variation.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at Genbank https://www.ncbi.nlm.nih.gov/genbank/, accession numbers MK965548-MK965653 and MK929797-MK930264.

Identification of conserved peptides containing B-cell epitopes of Babesia bovis AMA-1 and their potential as diagnostics candidates

The apical membrane antigen 1 (AMA-1) is a protein of the micronemes that is present in all organisms of the phylum Apicomplexa; it has been shown that AMA-1 plays an essential role for parasite invasion to target cells. It has been reported that AMA-1 is conserved among different isolates of Babesia; however, it is unknown whether the protein contains conserved B-cell epitopes and whether these epitopes are recognized by antibodies from cattle in endemic areas. In this research, using an in silico analysis, four peptides were designed containing exposed and conserved linear B-cell epitopes from the extracellular region of Babesia bovis AMA-1. The selected peptides were chemically synthesized, and then each peptide was emulsified and used to immunize two bovines per peptide. The antibodies produced against these peptides were able to recognize intra-erythrocytic parasites in an IFAT, except peptide 4, which was insoluble. The synthetic peptides were covalently fixed to the wells of an ELISA plate and incubated with sera from B. bovis naturally infected cattle. Peptides P2AMA and P3AMA were recognized by the sera of naturally infected cattle from different regions of Mexico. Statistical analysis showed that the ELISA test for peptides P2AMA and P3AMA had a concordance of 91.2% and 61.1% compared to the IFAT, a sensitivity of 94.56% and 71.74%, and a specificity of 76.19% and 14.2%, respectively. The presence of antibodies in bovine sera from endemic areas that bind to the identified peptides indicates that AMA-1 from B. bovis has conserved B-cell epitopes involved in the immune response under natural conditions. However, to propose their use as vaccine or diagnostics candidates, a further characterization of the humoral immune response elicited in cattle by these peptides is needed.

Limited genetic variability of Cytauxzoon felis apical membrane antigen-1 (ama1) from domestic cats and bobcats

Background

Cytauxzoon felis is a tick-transmitted apicomplexan that causes cytauxzoonosis in domestic cats (Felis catus). Even with intensive care, the mortality rate of acute cytauxzoonosis approaches 40% in domestic cats, while bobcats (Lynx rufus), the natural intermediate host of C. felis, remain clinically asymptomatic. However, multiple reports of domestic cats surviving acute disease without any treatment exist. One hypothesis for survival of these cats is infection with unique C. felis genotypes of lower pathogenicity. Prior studies have identified genetically distinct C. felis isolates containing polymorphisms within internal transcribed spacer regions (ITS) of the rRNA operon. However, these polymorphisms do not correlate with the clinical outcome of cytauxzoonosis, and so additional genetic markers are needed to test this hypothesis. We selected C. felis apical membrane antigen-1 (ama1) as a potential genetic marker of differential pathogenicity. AMA1 is a vaccine candidate for relatives of C. felis within Plasmodium spp.; however its historically high level of genetic polymorphism has resulted in escape from vaccine-induced immunity. While such diversity has hindered vaccine development, the expected polymorphism within the ama1 gene may be useful to evaluate population genetics.

Results

A 677 bp sequence of the C. felis ama1 gene was PCR-amplified from 84 domestic cats and 9 bobcats and demonstrated 99.9% sequence identity across all samples. A single nucleotide polymorphism (SNP) was identified in domestic cats and bobcats with evidence for co-infection with both genotypes identified in two domestic cats. The prevalence of the two genotypes varied with geographical distribution in domestic cats. Nucleotide diversity (π) and haplotype diversity (H) were calculated for C. felis ama1 and ama1 of related apicomplexans to assess genetic diversity. Based on these values (π = 0.00067 and H = 0.457), the diversity of the C. felis ama1 gene region analyzed is considerably lower than what is documented in related apicomplexans.

Conclusions

In surprising contrast to related apicomplexans, our results support that the sequence of the C. felis ama1 gene is highly conserved. While lack of genetic diversity limits utility of C. felis AMA1 as a genetic marker for clinical outcome, it supports further investigation as a vaccine candidate for cytauxzoonosis.

Genetic sequence characterization and naturally acquired immune response to Plasmodium vivax Rhoptry Neck Protein 2 (PvRON2)

Background

The genetic diversity of malaria antigens often results in allele variant-specific immunity, imposing a great challenge to vaccine development. Rhoptry Neck Protein 2 (PvRON2) is a blood-stage antigen that plays a key role during the erythrocyte invasion of Plasmodium vivax. This study investigates the genetic diversity of PvRON2 and the naturally acquired immune response to P. vivax isolates.

Results

Here, the genetic diversity of PvRON2_1828–2080 and the naturally acquired humoral immune response against PvRON2_1828–2080 in infected and non-infected individuals from a vivax malaria endemic area in Brazil was reported. The diversity analysis of PvRON2_1828–2080 revealed that the protein is conserved in isolates in Brazil and worldwide. A total of 18 (19%) patients had IgG antibodies to PvRON2_1828–2080. Additionally, the analysis of the antibody response in individuals who were not acutely infected with malaria, but had been infected with malaria in the past indicated that 32 patients (33%) exhibited an IgG immune response against PvRON2.

Conclusions

PvRON2 was conserved among the studied isolates. The presence of naturally acquired antibodies to this protein in the absence of the disease suggests that PvRON2 induces a long-term antibody response. These results indicate that PvRON2 is a potential malaria vaccine candidate.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2543-7) contains supplementary material, which is available to authorized users.

On the Evolution and Function of Plasmodium vivax Reticulocyte Binding Surface Antigen (pvrbsa)

The RBSA protein is encoded by a gene described in Plasmodium species having tropism for reticulocytes. Since this protein is antigenic in natural infections and can bind to target cells, it has been proposed as a potential candidate for an anti-Plasmodium vivax vaccine. However, genetic diversity (a challenge which must be overcome for ensuring fully effective vaccine design) has not been described at this locus. Likewise, the minimum regions mediating specific parasite-host interaction have not been determined. This is why the rbsa gene’s evolutionary history is being here described, as well as the P. vivax rbsa (pvrbsa) genetic diversity and the specific regions mediating parasite adhesion to reticulocytes. Unlike what has previously been reported, rbsa was also present in several parasite species belonging to the monkey-malaria clade; paralogs were also found in Plasmodium parasites invading reticulocytes. The pvrbsa locus had less diversity than other merozoite surface proteins where natural selection and recombination were the main evolutionary forces involved in causing the observed polymorphism. The N-terminal end (PvRBSA-A) was conserved and under functional constraint; consequently, it was expressed as recombinant protein for binding assays. This protein fragment bound to reticulocytes whilst the C-terminus, included in recombinant PvRBSA-B (which was not under functional constraint), did not. Interestingly, two PvRBSA-A-derived peptides were able to inhibit protein binding to reticulocytes. Specific conserved and functionally important peptides within PvRBSA-A could thus be considered when designing a fully-effective vaccine against P. vivax.

Analysis of Plasmodium vivax Apical Membrane Antigen-1 (PvAMA-1) Haplotypes among Iranian Isolates

Plasmodium vivax apical membrane antigen-1(PvAMA-1) is a surface protein with polymorphic sites. This study was aimed to analyze the polymorphic amino acid residues at PvAMA-1 in different infected age groups. 92 blood samples were collected from the south and southeast of Iran. The DNA coding for the domain I (DI), DII, and partial DIII of this antigen was amplified by Nested-PCR, and sequenced. Nucleotide mutations were found in 49 sites and based on the amino acid sequence, 30 variable sites were detected. Age distribution of malaria cases showed that the majority of the patients were between 10 to 30 years old. The scattering plot haplotypes by age showed an increasing incidence rate with age during childhood, whereas, incidence was the lowest in patients under five years old. Comparison of the polymorphic sites of PvAMA-1 in Iranian isolates with those found in other geographic regions of the world indicated nine common variable positions. In addition, a significant dependence was found between some particular substitutions and age categories. Dependence between particular substitutions and age groups suggests that certain residues in AMA-1 are responsible for clinical attacks in different ages, likely as a result of host immune pressure. The crystal structure of the PvAMA-1 showed that the amino acid substitutions that changed the protein charge were exclusively located in loops and turns where, the interactions with antibodies could occur. These data provide the necessary information for an AMA-1 based malaria vaccine design to be effective across all ages.

Low levels of genetic diversity associated with evidence of negative selection on the Babesia bovis apical membrane antigen 1 from parasite populations in Thailand

Babesia bovis, a parasite infecting cattle and buffalo, continues to spread throughout the developing world. The babesial vaccine was developed to be a sustainable alternative treatment to control the parasite. However, genetic diversity is a major obstacle for designing and developing a safe and effective vaccine. The apical membrane antigen 1 (AMA-1) is considered to be a potential vaccine candidate antigen among immunogenic genes of B. bovis. To gain a more comprehensive understanding of B. bovis AMA-1 (BbAMA-1), three B. bovis DNA samples were randomly selected to characterize in order to explore genetic diversity and natural selection and to predict the antigen epitopes. The sequence analysis revealed that BbAMA-1 has a low level of polymorphism and is highly conserved (95.46-99.94%) among Thai and global isolates. The majority of the polymorphic sites were observed in domains I and III. Conversely, domain II contained no polymorphic sites. We report the first evidence of strong negative or purifying selection across the full length of the gene, especially in domain I, by demonstrating a significant excess of the average number of synonymous (dS) over the non-synonymous (dN) substitutions. Finally, we also predict the linear and conformational B-cell epitope. The predicted B-cell epitopes appeared to be involved with the amino acid changes. Collectively, the results suggest that the conserved BbAMA-1 may be used to detect regional differences in the B. bovis parasite. Importantly, the limitation of BbAMA-1 diversity under strong negative selection indicates strong functional constraints on this gene. Thus, the gene could be a valuable target vaccine candidate antigen.

Temporal genetic changes in Plasmodium vivax apical membrane antigen 1 over 19 years of transmission in southern Mexico

Background

Mexico advanced to the pre-elimination phase in 2009 due to a significant reduction in malaria cases, and since 2000, Plasmodium vivax is the only species transmitted. During the last two decades, malaria transmission has been mostly local and isolated to a few regions. It is important to gain further insights into the impact of control measures on the parasite population structure. Hence, the aim of the current study was to determine detailed changes in P. vivax genetic diversity and population structure based on analysing the gene that encodes the apical membrane antigen 1 (pvama1). This analysis covered from control to pre-elimination (1993–2011) in a hypo-endemic region in southern Mexico.

Results

The 213 pvama1_I-II sequences presently analysed were grouped into six periods of three years each. They showed low genetic diversity, with 15 haplotypes resolved. Among the DNA sequences, there was a gradual decrease in genetic diversity, the number of mixed genotype infections and the intensity of positive selection, in agreement with the parallel decline in malaria cases. At the same time, linkage disequilibrium (R²) increased. The three-dimensional haplotype network revealed that pvama1_I-II haplotypes were separated by 1–11 mutational steps, and between one another by 0–3 unsampled haplotypes. In the temporal network, seven haplotypes were detected in at least two of the six-time layers, and only four distinct haplotypes were evidenced in the pre-elimination phase. Structure analysis indicated that three subpopulations fluctuated over time. Only 8.5% of the samples had mixed ancestry. In the pre-elimination phase, subpopulation P1 was drastically reduced, and the admixture was absent.

Conclusions

The results suggest that P. vivax in southern Mexico evolved based on local adaptation into three “pseudoclonal” subpopulations that diversified at the regional level and persisted over time, although with varying frequency. Control measures and climate events influenced the number of malaria cases and the genetic structure. The sharp decrease in parasite diversity and other related genetic parameters during the pre-elimination phase suggests that malaria elimination is possible in the near future. These results are useful for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2156-y) contains supplementary material, which is available to authorized users.

Analysis of Pvama1 genes from China-Myanmar border reveals little regional genetic differentiation of Plasmodium vivax populations

Background

With the premise of diminishing parasite genetic diversity following the reduction of malaria incidence, the analysis of polymorphic antigenic markers may provide important information about the impact of malaria control on local parasite populations. Here we evaluated the genetic diversity of Plasmodium vivax apical membrane antigen 1 (Pvama1) gene in a parasite population from the China-Myanmar border and compared it with global P. vivax populations.

Methods

We performed evolutionary analysis to examine the genetic diversity, natural selection, and population differentiation of 73 Pvama1 sequences acquired from the China-Myanmar border as well as 615 publically available Pvama1 sequences from seven global P. vivax populations.

Results

A total of 308 Pvama1 haplotypes were identified among the global P. vivax isolates. The overall nucleotide diversity of Pvama1 gene among the 73 China-Myanmar border parasite isolates was 0.008 with 41 haplotypes being identified (Hd = 0.958). Domain I (DI) harbored the majority (26/33) of the polymorphic sites. The McDonald Kreitman test showed a significant positive selection across the ectodomain and the DI of Pvama1. The fixation index (F_ST) estimation between the China-Myanmar border, Thailand (0.01) and Myanmar (0.10) showed only slight geographical genetic differentiation. Notably, the Sal-I haplotype was not detected in any of the analyzed global isolates, whereas the Belem strain was restricted to the Thai population. The detected mutations are mapped outside the overlapped region of the predicted B-cell epitopes and intrinsically unstructured/disordered regions.

Conclusions

This study revealed high levels of genetic diversity of Pvama1 in the P. vivax parasite population from the China-Myanmar border with DI displaying stronger diversifying selection than other domains. There were low levels of population subdivision among parasite populations from the Greater Mekong Subregion.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1899-1) contains supplementary material, which is available to authorized users.

Characterizing the impact of sustained sulfadoxine/pyrimethamine use upon the Plasmodium falciparum population in Malawi

BACKGROUND

Malawi experienced prolonged use of sulfadoxine/pyrimethamine (SP) as the front-line anti-malarial drug, with early replacement of chloroquine and delayed introduction of artemisinin-based combination therapy. Extended use of SP, and its continued application in pregnancy is impacting the genomic variation of the Plasmodium falciparum population.

METHODS

Whole genome sequence data of P. falciparum isolates covering 2 years of transmission within Malawi, alongside global datasets, were used. More than 745,000 SNPs were identified, and differences in allele frequencies between countries assessed, as well as genetic regions under positive selection determined.

RESULTS

Positive selection signals were identified within dhps, dhfr and gch1, all components of the parasite folate pathway associated with SP resistance. Sitting predominantly on a dhfr triple mutation background, a novel copy number increase of ~twofold was identified in the gch1 promoter. This copy number was almost fixed (96.8% frequency) in Malawi samples, but found at less than 45% frequency in other African populations, and distinct from a whole gene duplication previously reported in Southeast Asian parasites.

CONCLUSIONS

SP resistance selection pressures have been retained in the Malawian population, with known resistance dhfr mutations at fixation, complemented by a novel gch1 promoter duplication. The effects of the duplication on the fitness costs of SP variants and resistance need to be elucidated.

Diversity and evolutionary genetics of the three major Plasmodium vivax merozoite genes participating in reticulocyte invasion in southern Mexico

Background

Reported malaria cases in the Americas had been reduced to about one-half million by 2012. To advance towards elimination of this disease, it is necessary to gain insights into how the malaria parasite is evolving, including the emergence, spread and persistence of new haplotypes in affected regions. In here, the genetic diversity of the three major P. vivax merozoite genes was analyzed.

Methods

From P. vivax-infected blood samples obtained in southern Mexico (SMX) during 2006–2007, nucleotide sequences were achieved for: the 42 kDa carboxyl fragment of the merozoite surface protein-1 (msp1₄₂), domains I-II of the apical membrane antigen-1 (ama1_I-II), and domain II of the Duffy binding protein (dbp_II). Gene polymorphism was examined and haplotype networks were developed to depict parasite relationships in SMX. Then genetic diversity, recombination and natural selection were analyzed and the degree of differentiation was determined as F_ST values.

Results

The diversity of P. vivax merozoite genes in SMX was less than that of parasites from other geographic origins, with dbp_II < ama1_I-II < msp1₄₂. Ama1_I-II and msp1₄₂ exposed the more numerous haplotypes exclusive to SMX. While, all dbp_II haplotypes from SMX were separated from one to three mutational steps, the networks of ama1_I-II and msp1₄₂ were more complex; loops and numerous mutational steps were evidenced, likely due to recombination. Sings of local diversification were more evident for msp1₄₂. Sixteen combined haplotypes were determined; one of these haplotypes not detected in 2006 was highly frequent in 2007. The Rm value was higher for msp1₄₂than for ama1_I-II, being insignificant for dbp_II. The dN-dS value was highly significant for ama1_I-II and lesser so for dbp_II. The F_ST values were higher for dbp_II than msp1₄₂, and very low for ama1_I-II.

Conclusions

In SMX, P. vivax ama1_I-II, dbp_II and msp1₄₂ demonstrated limited diversity, and exhibited a differentiated parasite population. The results suggest that differential intensities of selective forces are operating on these gene fragments, and probably related to their timing, length of exposure and function during reticulocyte adhesion and invasion. Therefore, these finding are essential for mono and multivalent vaccine development and for epidemiological surveillance.

Electronic supplementary material

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

Genetic polymorphism in domain I of the apical membrane antigen-1 among Plasmodium knowlesi clinical isolates from Peninsular Malaysia

The simian malaria parasite Plasmodium knowlesi is now recognized as a species that can cause human malaria. The first report of large scale human knowlesi malaria was in 2004 in Malaysia Borneo. Since then, hundreds of human knowlesi malaria cases have been reported in Southeast Asia. The present study investigates the genetic polymorphism of P. knowlesi DI domain of the apical membrane antigen-1 (AMA-1), a protein considered as a promising vaccine candidate for malaria. The DI domain of AMA-1 gene of P. knowlesi clinical isolates from Peninsular Malaysia was amplified by PCR, cloned into Escherichia coli, then sequenced and analysed. Ninety-seven DI domain sequences were obtained. Comparison at the nucleotide level against P. knowlesi strain H as reference sequence showed 21 synonymous and 25 nonsynonymous mutations. Nonetheless, nucleotide sequence analysis revealed low genetic diversity of the DI domain, and it was under purifying (negative) selection. At the amino acid level, 26 different haplotypes were identified and 2 were predominant haplotypes (H1, H2) with high frequencies. Phylogenetic analysis revealed that the 26 haplotypes could be clustered into 2 distinct groups (I and II). Members of the groups were basically derived from haplotypes H1 and H2, respectively.

Population genetic structure and natural selection of apical membrane antigen-1 in Plasmodium vivax Korean isolates

Background

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is a leading candidate antigen for blood stage malaria vaccine. However, antigenic variation is a major obstacle in the development of an effective vaccine based on this antigen. In this study, the genetic structure and the effect of natural selection of PvAMA-1 among Korean P. vivax isolates were analysed.

Methods

Blood samples were collected from 66 Korean patients with vivax malaria. The entire PvAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. The PvAMA-1 sequence of each isolate was sequenced and the polymorphic characteristics and effect of natural selection were analysed using the DNASTAR, MEGA4, and DnaSP programs.

Results

Thirty haplotypes of PvAMA-1, which were further classified into seven different clusters, were identified in the 66 Korean P. vivax isolates. Domain II was highly conserved among the sequences, but substantial nucleotide diversity was observed in domains I and III. The difference between the rates of non-synonymous and synonymous mutations suggested that the gene has evolved under natural selection. No strong evidence indicating balancing or positive selection on PvAMA-1 was identified. Recombination may also play a role in the resulting genetic diversity of PvAMA-1.

Conclusions

This study is the first comprehensive analysis of nucleotide diversity across the entire PvAMA-1 gene using a single population sample from Korea. Korean PvAMA-1 had limited genetic diversity compared to PvAMA-1 in global isolates. The overall pattern of genetic polymorphism of Korean PvAMA-1 differed from other global isolates and novel amino acid changes were also identified in Korean PvAMA-1. Evidences for natural selection and recombination event were observed, which is likely to play an important role in generating genetic diversity across the PvAMA-1. These results provide useful information for the understanding the population structure of P. vivax circulating in Korea and have important implications for the design of a vaccine incorporating PvAMA-1.

Molecular Evolution of PvMSP3α Block II in Plasmodium vivax from Diverse Geographic Origins

Block II of Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is conserved and has been proposed as a potential candidate for a malaria vaccine. The present study aimed to compare sequence diversity in PvMSP3a block II at a local microgeographic scale in a village as well as from larger geographic regions (countries and worldwide). Blood samples were collected from asymptomatic carriers of P. vivax in a village at the western border of Thailand and PvMSP3α was amplified and sequenced. For population genetic analysis, 237 PvMSP3α block II sequences from eleven P. vivax endemic countries were analyzed. PvMSP3α sequences from 20 village-level samples revealed two length variant types with one type containing a large deletion in block I. In contrast, block II was relatively conserved; especially, some non-synonymous mutations were extensively shared among 11 parasite populations. However, the majority of the low-frequency synonymous variations were population specific. The conserved pattern of nucleotide diversity in block II sequences was probably due to functional/structural constraints, which were further supported by the tests of neutrality. Notably, a small region in block II that encodes a predicted B cell epitope was highly polymorphic and showed signs of balancing selection, signifying that this region might be influenced by the immune selection and may serve as a starting point for designing multi-antigen/stage epitope based vaccines against this parasite.

Babesia divergens apical membrane antigen-1 (BdAMA-1): A poorly polymorphic protein that induces a weak and late immune response

Babesiosis is an important veterinary and zoonotic tick borne disease caused by the hemoprotozoan Babesia spp. which infects red blood cell of its vertebrate host. In order to control the infection, vaccination that targets molecules involved in the invasion process of red blood cells could provide a good alternative to chemotherapy. Among these molecules, Apical Membrane Antigen-1 (AMA-1) has been described as an excellent vaccine candidate in Plasmodium spp. In this paper, we have investigated AMA-1 of Babesia divergens (BdAMA-1) as vaccine candidate by evaluating its polymorphism and by studying the humoral response against BdAMA-1 of sheep experimentally infected with B. divergens. Polymorphism of BdAMA-1 was investigated by sequencing the corresponding gene of 9 B. divergens isolates from different geographical areas in France. Two Bdama-1 haplotypes (A and B) could be defined based on 2 non-synonymous point mutations. In silico prediction of linear epitopes revealed that the antigenicity of the 2 haplotypes is very similar. Antibody production against the extracellular domain of BdAMA-1 is weak and late, between 1 and 5 months after the inoculation of parasites. Both IgG1 and IgG2 are components of the anti-BdAMA-1 response. These results indicate that while BdAMA-1 may not be an immuno-dominant antigen, it could induce a mixed type 1 and type 2 immune response. In light of these results, the potential of BdAMA-1 as vaccine candidate is discussed.

Inferring natural selection signals in Plasmodium vivax-encoded proteins having a potential role in merozoite invasion

Detecting natural selection signals in Plasmodium parasites antigens might be used for identifying potential new vaccine candidates. Fifty-nine Plasmodium vivax-Sal-I genes encoding proteins having a potential role in invasion were used as query for identifying them in recent P. vivax strain genome sequences and two closely-related Plasmodium species. Several measures of DNA sequence variation were then calculated and selection signatures were detected by using different approaches. Our results may be used for determining which genes expressed during P. vivax merozoite stage could be prioritised for further population genetics or functional studies for designing a P. vivax vaccine which would avoid allele-specific immune responses.

Structure and dynamics of apical membrane antigen 1 from Plasmodium falciparum FVO

Apical membrane antigen 1 (AMA1) interacts with RON2 to form a protein complex that plays a key role in the invasion of host cells by malaria parasites. Blocking this protein-protein interaction represents a potential route to controlling malaria and related parasitic diseases, but the polymorphic nature of AMA1 has proven to be a major challenge to vaccine-induced antibodies and peptide inhibitors exerting strain-transcending inhibitory effects. Here we present the X-ray crystal structure of AMA1 domains I and II from Plasmodium falciparum strain FVO. We compare our new structure to those of AMA1 from P. falciparum 3D7 and Plasmodium vivax. A combination of normalized B factor analysis and computational methods has been used to investigate the flexibility of the domain I loops and how this correlates with their roles in determining the strain specificity of human antibody responses and inhibitory peptides. We also investigated the domain II loop, a key region involved in inhibitor binding, by comparison of multiple AMA1 crystal structures. Collectively, these results provide valuable insights that should contribute to the design of strain-transcending agents targeting P. falciparum AMA1.

On the road to eliminate malaria in Sri Lanka: lessons from history, challenges, gaps in knowledge and research needs

Malaria is one of the most important tropical diseases that has caused devastation throughout the history of mankind. Malaria eradication programmes in the past have had many positive effects but failed to wipe out malaria from most tropical countries, including Sri Lanka. Encouraged by the impressive levels of reduction in malaria case numbers during the past decade, Sri Lanka has launched a programme to eliminate malaria by year 2014. This article reviews the historical milestones associated with the malaria eradication programme that failed subsequently and the events that led to the launch of the ongoing malaria elimination plans at national-level and its strategies that are operational across the entire country. The existing gaps in knowledge are also discussed together with the priority areas for research to fill in these gaps that are posing as challenges to the envisaged goal of wiping out malaria from this island nation.

Global Population Structure of the Genes Encoding the Malaria Vaccine Candidate, Plasmodium vivax Apical Membrane Antigen 1 (PvAMA1)

Background

The Plasmodium vivax Apical Membrane Antigen 1 (PvAMA1) is a promising malaria vaccine candidate, however it remains unclear which regions are naturally targeted by host immunity and whether its high genetic diversity will preclude coverage by a monovalent vaccine. To assess its feasibility as a vaccine candidate, we investigated the global population structure of PvAMA1.

Methodology and Principal Findings

New sequences from Papua New Guinea (PNG, n = 102) were analysed together with published sequences from Thailand (n = 158), India (n = 8), Sri Lanka (n = 23), Venezuela (n = 74) and a collection of isolates from disparate geographic locations (n = 8). A total of 92 single nucleotide polymorphisms (SNPs) were identified including 22 synonymous SNPs and 70 non-synonymous (NS) SNPs. Polymorphisms and signatures of balancing (positive Tajima's D and low F_ST values) selection were predominantly clustered in domain I, suggesting it is a dominant target of protective immune responses. To estimate global antigenic diversity, haplotypes comprised of (i) non-singleton (n = 40) and (ii) common (≥10% minor allele frequency, n = 23) polymorphic amino acid sites were then analysed revealing a total of 219 and 210 distinct haplotypes, respectively. Although highly diverse, the 210 haplotypes comprised of only common polymorphisms were grouped into eleven clusters, however substantial geographic differentiation was observed, and this may have implications for the efficacy of PvAMA1 vaccines in different malaria-endemic areas. The PNG haplotypes form a distinct group of clusters not found in any other geographic region. Vaccine haplotypes were rare and geographically restricted, suggesting potentially poor efficacy of candidate PvAMA1 vaccines.

Conclusions

It may be possible to cover the existing global PvAMA1 diversity by selection of diverse alleles based on these analyses however it will be important to first define the relationships between the genetic and antigenic diversity of this molecule.

Traditionally misclassified as benign and neglected as a research priority, it is now understood that P. vivax is an increasingly important cause of human malaria. This important human pathogen poses an enormous obstacle to malaria control and elimination efforts due its broad geographic distribution, ability to cause recurring episodes of malaria after long periods of inactivity and extreme biodiversity. Vaccines are an essential component of global malaria control and elimination campaigns but the diversity of malaria antigens is thought to be a major cause of vaccine failure. Furthermore, at present the majority of current vaccine research is directed toward P. falciparum. The aims of this study were to investigate the global diversity of the P. vivax vaccine candidate, Apical Membrane Antigen 1 (PvAMA1), to determine the feasibility of designing a globally effective PvAMA1 vaccine and to determine which region of PvAMA1 is targeted by host immune responses, in order to identify the most promising vaccine candidates. We report that PvAMA1 diversity is extremely high, and that PvAMA1 domain I is a dominant target of host immune responses. These analyses of PvAMA1 diversity from several geographic regions provide a framework to guide development of a broadly efficacious P. vivax vaccine.

A Large Plasmodium vivax Reservoir and Little Population Structure in the South Pacific

Introduction

The importance of Plasmodium vivax in malaria elimination is increasingly being recognized, yet little is known about its population size and population genetic structure in the South Pacific, an area that is the focus of intensified malaria control.

Methods

We have genotyped 13 microsatellite markers in 295 P. vivax isolates from four geographically distinct sites in Papua New Guinea (PNG) and one site from Solomon Islands, representing different transmission intensities.

Results

Diversity was very high with expected heterozygosity values ranging from 0.62 to 0.98 for the different markers. Effective population size was high (12′872 to 19′533 per site). In PNG population structuring was limited with moderate levels of genetic differentiation. F_ST values (adjusted for high diversity of markers) were 0.14–0.15. Slightly higher levels were observed between PNG populations and Solomon Islands (F_ST = 0.16).

Conclusions

Low levels of population structure despite geographical barriers to transmission are in sharp contrast to results from regions of low P. vivax endemicity. Prior to intensification of malaria control programs in the study area, parasite diversity and effective population size remained high.

Population genetic structure and polymorphism analysis of gene encoding apical membrane antigen-1 (AMA-1) of Iranian Plasmodium vivax wild isolates

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is a major candidate antigen for human malaria vaccine. In the present study, polymorphism of pvama-1 among Iranian isolates was investigated to generate useful information on this vaccine candidate antigen, which is required for the rational design of a vaccine against P. vivax. Blood samples were collected from P. vivax-infected Iranian patients during 2009-2010. Of 99 collected isolates, 37 were analyzed for almost the entire pvama-1 gene using sequencing. The overall nucleotide diversity (π) was 0.00826 ± 0.0004 and the majority of polymorphic sites were identified in domain I (DI) of the pvama-1 gene. Neutrality analysis using Tajima's D, Fu and Li's D* and F* and McDonald Kreitman tests showed a significant positive departure from neutral substitution patterns, indicating a possible balancing selection across the entire ectodomain and DI sequences of pvama-1 gene. However, no evidence was found for the balancing selection in DII and DIII regions of Iranian PvAMA-1. Also, 29 haplotypes with different frequencies were identified and the overall haplotype diversity was 0.982 ± 0.012. Epitope mapping prediction of PvAMA-1 showed the potential B-cell epitopes across DI-DIII overlap with E145K, P210S, R249H, G253E, K352E, R438H and N445D mutations; however, no mutation has been found in intrinsically unstructured/disordered regions. The fixation index (Fst) estimation between Iran and the closest geographical sites such as India (0.0707) showed a slight geographical genetic differentiation; however, the Fst estimation between Iran and Thailand (0.1253) suggested a moderate geographical isolation. In summary, genetic investigation in pvama-1 among Iranian P. vivax isolates indicates that this antigen showed limited antigenic diversity and most of the detected mutations are located outside B-cell epitopes. Therefore, the present results have significant implications in understanding the nature of P. vivax population circulating in Iran as well as in providing useful information for malaria vaccine development based on this antigen.

Population genetics, sequence diversity and selection in the gene encoding the Plasmodium falciparum apical membrane antigen 1 in clinical isolates from the south-east of Iran

The Plasmodium falciparum apical membrane antigen1 (AMA1) is a leading malaria vaccine candidate antigen. In the present investigation, for the first time, the almost full length of the ama1 gene covering domain I (DI), DII and DIII was PCR amplified and sequenced in 21 P. falciparum isolates collected from the southeastern parts of Iran. The result showed the low genetic diversity of Iranian PfAMA1 with 11 PfAMA1 haplotypes in which nine out of 11 haplotypes are novel and have been reported for the first time. The Iranian P. falciparum population indicated a moderate level of genetic differentiation. The difference among the rates of non-synonymous and synonymous mutations, Tajima's D and McDonald-Kreitman tests suggested that the diversity at DI is due to positive natural selection. In addition, recombination contributes to the diversity of Iranian PfAMA1 and this is supported by the decline of the linkage disequilibrium index R(2) with increasing the nucleotide distance. The highly polymorphic residues (positions: 187, 197, 200, 230 and 243) were polymorphic; however, most of the SNPs in non-polymorphic residues were conserved except the residue at position 395. Nevertheless, no mutation was found in the DII loop of the Iranian PfAMA1, indicating that it is subjected to purifying selection. In conclusion, the low genetic diversity in PfAMA1 among Iranian isolates supports and provides valuable information for the development of a PfAMA1-based malaria vaccine.

Genetic diversity of Plasmodium vivax and Plasmodium falciparum in Honduras

BACKGROUND

Understanding the population structure of Plasmodium species through genetic diversity studies can assist in the design of more effective malaria control strategies, particularly in vaccine development. Central America is an area where malaria is a public health problem, but little is known about the genetic diversity of the parasite's circulating species. This study aimed to investigate the allelic frequency and molecular diversity of five surface antigens in field isolates from Honduras.

METHODS

Five molecular markers were analysed to determine the genotypes of Plasmodium vivax and Plasmodium falciparum from endemic areas in Honduras. Genetic diversity of ama-1, msp-1 and csp was investigated for P. vivax, and msp-1 and msp-2 for P. falciparum. Allelic frequencies were calculated and sequence analysis performed.

RESULTS AND CONCLUSION

A high genetic diversity was observed within Plasmodium isolates from Honduras. A different number of genotypes were elucidated: 41 (n = 77) for pvama-1; 23 (n = 84) for pvcsp; and 23 (n = 35) for pfmsp-1. Pvcsp sequences showed VK210 as the only subtype present in Honduran isolates. Pvmsp-1 (F2) was the most polymorphic marker for P. vivax isolates while pvama-1 was least variable. All three allelic families described for pfmsp-1 (n = 30) block 2 (K1, MAD20, and RO33), and both allelic families described for the central domain of pfmsp-2 (n = 11) (3D7 and FC27) were detected. However, K1 and 3D7 allelic families were predominant. All markers were randomly distributed across the country and no geographic correlation was found. To date, this is the most complete report on molecular characterization of P. vivax and P. falciparum field isolates in Honduras with regards to genetic diversity. These results indicate that P. vivax and P. falciparum parasite populations are highly diverse in Honduras despite the low level of transmission.

The evolutionary consequences of blood-stage vaccination on the rodent malaria Plasmodium chabaudi

A candidate malaria vaccine promoted the evolution of more virulent malaria parasites in mice.

Malaria vaccine developers are concerned that antigenic escape will erode vaccine efficacy. Evolutionary theorists have raised the possibility that some types of vaccine could also create conditions favoring the evolution of more virulent pathogens. Such evolution would put unvaccinated people at greater risk of severe disease. Here we test the impact of vaccination with a single highly purified antigen on the malaria parasite Plasmodium chabaudi evolving in laboratory mice. The antigen we used, AMA-1, is a component of several candidate malaria vaccines currently in various stages of trials in humans. We first found that a more virulent clone was less readily controlled by AMA-1-induced immunity than its less virulent progenitor. Replicated parasites were then serially passaged through control or AMA-1 vaccinated mice and evaluated after 10 and 21 rounds of selection. We found no evidence of evolution at the ama-1 locus. Instead, virulence evolved; AMA-1-selected parasites induced greater anemia in naïve mice than both control and ancestral parasites. Our data suggest that recombinant blood stage malaria vaccines can drive the evolution of more virulent malaria parasites.

Vaccination can drive the evolution of pathogens. Most obviously, molecules targeted by vaccine-induced immunity can change. Such evolution makes vaccines less effective. A different possibility is that more virulent pathogens are favored in vaccinated hosts. In that case, vaccination would create pathogens that cause more harm to unvaccinated individuals. To test this idea, we studied a rodent malaria parasite in laboratory mice immunized with a component of malaria vaccines currently in human trials. We found that a more virulent parasite clone was less well controlled by vaccine-induced immunity than was its less virulent ancestor. We then passaged parasites through sham- or vaccinated mice to study how the parasites might evolve after multiple rounds of infection of mouse hosts. The parasite molecule targeted by the vaccine did not change during this process. Instead, the parasites became more virulent if they evolved in vaccinated hosts. Our data suggest that some vaccines can drive the evolution of more virulent parasites.

SNP genotyping identifies new signatures of selection in a deep sample of West African Plasmodium falciparum malaria parasites

We used a high-density single-nucleotide polymorphism array to genotype 75 Plasmodium falciparum isolates recently collected from Senegal and The Gambia to search for signals of selection in this malaria endemic region. We found little geographic or temporal stratification of the genetic diversity among the sampled parasites. Through application of the iHS and REHH haplotype-based tests for positive selection, we found evidence of recent selective sweeps at a known drug resistance locus, at several known antigenic loci, and at several genomic regions not previously identified as sites of recent selection. We discuss the value of deep population-specific genomic analyses for identifying selection signals within sampled endemic populations of parasites, which may correspond to local selection pressures such as distinctive therapeutic regimes or mosquito vectors.

Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite outside of Africa, and its range extends well into the temperate zones. Previous studies provided evidence for vivax population differentiation, but temperate vivax parasites were not well represented in these analyses. Here we address this deficit by using complete mitochondrial (mt) genome sequences to elucidate the broad genetic diversity and population structure of P. vivax from temperate regions in East and Southeast Asia.

RESULTS

From the complete mtDNA sequences of 99 clinical samples collected in China, Myanmar and Korea, a total of 30 different haplotypes were identified from 26 polymorphic sites. Significant differentiation between different East and Southeast Asian parasite populations was observed except for the comparison between populations from Korea and southern China. Haplotype patterns and structure diversity analysis showed coexistence of two different groups in East Asia, which were genetically related to the Southeast Asian population and Myanmar population, respectively. The demographic history of P. vivax, examined using neutrality tests and mismatch distribution analyses, revealed population expansion events across the entire P. vivax range and the Myanmar population. Bayesian skyline analysis further supported the occurrence of ancient P. vivax population expansion.

CONCLUSIONS

This study provided further resolution of the population structure and evolution of P. vivax, especially in temperate/warm-temperate endemic areas of Asia. The results revealed divergence of the P. vivax populations in temperate regions of China and Korea from other populations. Multiple analyses confirmed ancient population expansion of this parasite. The extensive genetic diversity of the P. vivax populations is consistent with phenotypic plasticity of the parasites, which has implications for malaria control.

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.

Understanding the population genetics of Plasmodium vivax is essential for malaria control and elimination

Traditionally, infection with Plasmodium vivax was thought to be benign and self-limiting, however, recent evidence has demonstrated that infection with P. vivax can also result in severe illness and death. Research into P. vivax has been relatively neglected and much remains unknown regarding the biology, pathogenesis and epidemiology of this parasite. One of the fundamental factors governing transmission and immunity is parasite diversity. An understanding of parasite population genetic structure is necessary to understand the epidemiology, diversity, distribution and dynamics of natural P. vivax populations. In addition, studying the population structure of genes under immune selection also enables investigation of the dynamic interplay between transmission and immunity, which is crucial for vaccine development. A lack of knowledge regarding the transmission and spread of P. vivax has been particularly highlighted in areas where malaria control and elimination programmes have made progress in reducing the burden of Plasmodium falciparum, yet P. vivax remains as a substantial obstacle. With malaria elimination back on the global agenda, mapping of global and local P. vivax population structure is essential prior to establishing goals for elimination and the roll-out of interventions. A detailed knowledge of the spatial distribution, transmission and clinical burden of P. vivax is required to act as a benchmark against which control targets can be set and measured. This paper presents an overview of what is known and what is yet to be fully understood regarding P. vivax population genetics, as well as the importance and application of P. vivax population genetics studies.

Evaluation of the genetic diversity of domain II of Plasmodium vivax Apical Membrane Antigen 1 (PvAMA-1) and the ensuing strain-specific immune responses in patients from Sri Lanka

Antigenic polymorphism displayed by malaria parasites is a skewed schema to escape the host immune system. The prevailing genetic diversity at domain II of the Plasmodium vivax Apical Membrane Antigen-1 (Pvama-1DII) was characterized in 64 single clone P. vivax isolates from Sri Lanka, where unstable malaria prevails with low intensity. In Sri Lanka, the Pvama-1DII gene showed meager meiotic recombination with the enclosure of single nucleotide polymorphisms (SNPs). Eleven amino acid (a.a.) variant positions defined 21 a.a. haplotypes with 9 unique to the island, where the predominant haplotype, H1, was identical to the reference Salvador I strain. A further 376 globally dispersed isolates defined 38 a.a. haplotypes (H22-H59), with 4 and 26 haplotypes exclusive to India and Thailand, respectively. The phylogenetic tree revealed no clustering, where most isolates had a very recent common origin. The polymorphism detected in PvAMA-1DII B and T cell epitopes evidenced an immune evasion mechanism exploited by the parasite. Majority of Sri Lankan patients developed antibody responses to both conformational and linear B cell epitopes. The ensuing strain-specific immunity due to extensive antigenic polymorphism was evaluated by aligning a.a. sequences of PvAMA-1DII with the homologous total (IgM+IgG) antibody responses assayed by in-house established indirect ELISAs against 7 PvAMA-1DII overlapping synthetic peptides, P01-P07. While the antibody responses to P01-P03, P06, P07 harbouring P. vivax clinical isolates with polymorphic a.a. haplotype to Sal I was clearly strain-transcending (cross-reactive), individuals with isolates identical to the Sal I strain observed varying antibody prevalence against the seven PvAMA-1DII Sal-I synthetic peptides, with the highest prevalence detected against P04. Synthetic peptide P04, spanning a.a. positions 302-324 of the PvAMA-1DII of the Sal I strain that included the epitope recognized by the invasion inhibitory 4G2 monoclonal antibody of PfAMA-1, was highly conserved in all 440 local and global P. vivax isolates examined. A functional role for this region is reinforced by the highly immunogenic nature of P04, and could point towards a presumably "protective" anti-P04 antibody response that elicited an isotype switch from IgM to IgG, with increasing exposure to malaria exclusively in endemic residents. Thus the conserved and seemingly "protective" nature of the domain II loop of PvAMA-1 makes it a putative contender to be included in a cocktail vaccine against P. vivax asexual erythrocytic stages in Sri Lanka.

Identification of a highly antigenic linear B cell epitope within Plasmodium vivax apical membrane antigen 1 (AMA-1)

Apical membrane antigen 1 (AMA-1) is considered to be a major candidate antigen for a malaria vaccine. Previous immunoepidemiological studies of naturally acquired immunity to Plasmodium vivax AMA-1 (PvAMA-1) have shown a higher prevalence of specific antibodies to domain II (DII) of AMA-1. In the present study, we confirmed that specific antibody responses from naturally infected individuals were highly reactive to both full-length AMA-1 and DII. Also, we demonstrated a strong association between AMA-1 and DII IgG and IgG subclass responses. We analyzed the primary sequence of PvAMA-1 for B cell linear epitopes co-occurring with intrinsically unstructured/disordered regions (IURs). The B cell epitope comprising the amino acid sequence 290–307 of PvAMA-1 (SASDQPTQYEEEMTDYQK), with the highest prediction scores, was identified in domain II and further selected for chemical synthesis and immunological testing. The antigenicity of the synthetic peptide was identified by serological analysis using sera from P. vivax-infected individuals who were knowingly reactive to the PvAMA-1 ectodomain only, domain II only, or reactive to both antigens. Although the synthetic peptide was recognized by all serum samples specific to domain II, serum with reactivity only to the full-length protein presented 58.3% positivity. Moreover, IgG reactivity against PvAMA-1 and domain II after depletion of specific synthetic peptide antibodies was reduced by 18% and 33% (P = 0.0001 for both), respectively. These results suggest that the linear epitope SASDQPTQYEEEMTDYQK is highly antigenic during natural human infections and is an important antigenic region of the domain II of PvAMA-1, suggesting its possible future use in pre-clinical studies.

Genetic diversity of Plasmodium vivax Duffy Binding Protein II (PvDBPII) under unstable transmission and low intensity malaria in Sri Lanka

Elucidating the genetic diversity of the Duffy Binding Protein II (PvDBPII), a leading vaccine candidate for vivax malaria, in different geographical settings is vital. In Sri Lanka malaria transmission is unstable with low intensity. A relatively high level of allelic diversity, with 27 polymorphic nucleotides and 33 (aa) haplotypes was detected among the PvdbpII gene in 100 local Plasmodium vivax isolates collected from two hypoendemic areas, and from a non endemic area of the country. Mutations, recombination and balancing selection seem to maintain the observed local allelic diversity of PvdbpII. Lack of gene flow was evidenced by high Fst values between the two endemic study sites. Some of the aa polymorphisms may alter the binding and expression capacity of predicted T cell epitopes in PvDBPII. Of the 8 binding inhibitory linear B cell epitopes, 2 (H2 and M1) in the vicinity of the exact binding region of PvDBPII appeared to be highly conserved in Sri Lankan, Iran and Colombian isolates, while H3, M2, M3 and L3 neutralizing epitopes seem to be polymorphic globally, with H1 and L2 conserved in Colombian, South Korean and Iran isolates. In comparison to the reference Sal-1 strain, among 402 world-wide isolates (302 global and 100 local), 121 aa polymorphisms and 138 haplotypes were recorded of which 3 aa polymorphisms and 21 haplotypes seem to be unique to Sri Lanka. PvdbpII phylogeny suggests that local P. vivax parasites represent a sample of the global population. The ubiquitous presence of some PvDBPII aa haplotypes among both local and global P. vivax isolates may aid future vaccination strategies based on PvDBPII.

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in Brazilian Plasmodium vivax populations

Background

Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP_II), which is the most variable segment of the protein.

Methods

To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP_II in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP_II, and T- and B-cell epitopes were localized on the 3-D structure.

Results

The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP_II, and (ii) PvDBP_II appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium.

Conclusions

This study shows that some polymorphisms of PvDBP_II are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

Genetic polymorphism and effect of natural selection at domain I of apical membrane antigen-1 (AMA-1) in Plasmodium vivax isolates from Myanmar

Malaria is endemic or hypoendemic in Myanmar and the country still contributes to the high level of malaria deaths in South-East Asia. Although information on the nature and extent of population diversity within malaria parasites in the country is essential not only for understanding the epidemic situation but also to establish a proper control strategy, very little data is currently available on the extent of genetic polymorphisms of the malaria parasites in Myanmar. In this study, we analyzed the genetic polymorphism and natural selection at domain I of the apical membrane antigen-1 (AMA-1) among Plasmodium vivax Myanmar isolates. A total of 34 distinguishable haplotypes were identified among the 76 isolates sequenced. Comparison with the previously available PvAMA-1 sequences in the GenBank database revealed that 21 of them were new haplotypes that have never been reported till date. The difference between the rate of nonsynonymous (dN) and synonymous (dS) mutations was positive (dN-dS, 0.013+/-0.005), suggesting the domain I is under positive natural selection. The Tajima's D statistics was found to be -0.74652, suggesting that the gene has evolved under population size expansion and/or positive selection. The minimum recombination events were also high, indicating that recombination may occur within the domain I resulting in allelic diversity of PvAMA-1. Our results collectively suggest that PvAMA-1 displays high genetic polymorphism among Myanmar P. vivax isolates with highly diversifying selection at domain I. These results have significant implications in understanding the nature of P. vivax population circulating in Myanmar as well as providing useful information for malaria vaccine development based on this antigen.

Nucleotide sequence polymorphism at the apical membrane antigen-1 locus reveals population history of Plasmodium vivax in Thailand

Apical membrane antigen-1 is a candidate for inclusion in a vaccine for the human malaria parasite Plasmodium vivax. We collected 231 complete sequences of the gene encoding this antigen (pvama-1) from three regions of Thailand, the most extensive collection to date of sequences at this locus. The domain II loop (previously mentioned as a potential vaccine component) was almost completely conserved, with a single amino acid variant (I313R) observed in a single sequence. The 3' portion of the gene (domain II through the stop codon) showed significantly lower nucleotide diversity than the 5' portion (start codon through domain I); and a given domain I sequence might be found in a haplotype with more than one domain II sequence. These results imply a hotspot of recombination between domains I and II. We found significant geographic subdivision among the three regions of Thailand (NW, East, and South) in which collections were made in 2007. Numbers of P. vivax infections have experienced overall declines since 1990 in all three regions; but the decline has been most recent in the NW, and there has been a rebound in numbers of infections in the South since 2000. Consistent with population history, amino acid sequence diversity was greatest in the NW. The South, which had by far the lowest sequence diversity of the three regions, showed signs of a population that has expanded from a small number of founders after a bottleneck.

Limited global diversity of the Plasmodium vivax merozoite surface protein 4 gene

Merozoite surface proteins (MSPs) of the malaria parasites are major candidates for vaccine development targeting asexual blood stages. However, the diverse antigenic repertoire of these antigens that induce strain-specific protective immunity in human is a major challenge for vaccine design and often determines the efficacy of a vaccine. Here we further assessed the genetic diversity of Plasmodium vivax MSP4 (PvMSP4) protein using 195 parasite samples collected mostly from Thailand, Indonesia and Brazil. Overall, PvMSP4 is highly conserved with only eight amino acid substitutions. The majority of the haplotype diversity was restricted to the two short tetrapeptide repeat arrays in exon 1 and 2, respectively. Selection and neutrality tests indicated that exon 1 and the entire coding region of PvMSP4 were under purifying selection. Despite the limited nucleotide polymorphism of PvMSP4, significant genetic differentiation among the three major parasite populations was detected. Moreover, microgeographical heterogeneity was also evident in the parasite populations from different endemic areas of Thailand.

Polymorphism at the apical membrane antigen 1 gene (AMA1) of the malaria parasite Plasmodium falciparum in a Vietnamese population

The patterns of molecular evolution of the most diverse region of the apical membrane antigen 1 (AMA1) gene in Plasmodium falciparum from a Vietnamese subpopulation (Bao Loc) were investigated. Within the Bao Loc population, the sequenced gene region showed relatively high allelic and nucleotide diversity (0.985 and 0.02694, respectively). Further, the level of population recombination was substantial, resulting in a significant decay of linkage disequilibrium along the gene region. The results suggest that AMA1 is a useful genetic marker for studying the relationships between adaptation of parasite populations (to the human host immune system) and malaria epidemiology.

ama1 genes of sympatric Plasmodium vivax and P. falciparum from Venezuela differ significantly in genetic diversity and recombination frequency

Background

We present the first population genetic analysis of homologous loci from two sympatric human malaria parasite populations sharing the same human hosts, using full-length sequences of ama1 genes from Plasmodium vivax and P. falciparum collected in the Venezuelan Amazon.

Methodology/Principal Findings

Significant differences between the two species were found in genetic diversity at the ama1 locus, with 18 distinct haplotypes identified among the 73 Pvama1 sequences obtained, compared to 6 unique haplotypes from 30 Pfama1 sequences, giving overall diversity estimates of h = 0.9091, and h = 0.538 respectively. Levels of recombination were also found to differ between the species, with P. falciparum exhibiting very little recombination across the 1.77kb sequence. In contrast, analysis of patterns of nucleotide substitutions provided evidence that polymorphisms in the ama1 gene of both species are maintained by balancing selection, particularly in domain I. The two distinct population structures observed are unlikely to result from different selective forces acting upon the two species, which share both human and mosquito hosts in this setting. Rather, the highly structured P. falciparum population appears to be the result of a population bottleneck, while the much less structured P. vivax population is likely to be derived from an ancient pool of diversity, as reflected in a larger estimate of effective population size for this species. Greatly reduced mosquito transmission in 1997, due to low rainfall prior to the second survey, was associated with far fewer P. falciparum infections, but an increase in P. vivax infections, probably due to hypnozoite activation.

Conclusions/Significance

The relevance of these findings to putative competitive interactions between these two important human pathogen species is discussed. These results highlight the need for future control interventions to employ strategies targeting each of the parasite species present in endemic areas.

Plasmodium vivax apical membrane antigen-1: comparative recognition of different domains by antibodies induced during natural human infection

The Apical Membrane Antigen-1 (AMA-1) of Plasmodium sp. has been suggested as a vaccine candidate against malaria. This protein seems to be involved in merozoite invasion and its extra-cellular portion contains three distinct domains: DI, DII, and DIII. Previously, we described that Plasmodium vivax AMA-1 (PvAMA-1) ectodomain is highly immunogenic in natural human infections. Here, we expressed each domain, separately or in combination (DI-II or DII-III), as bacterial recombinant proteins to map immunodominant epitopes within the PvAMA-1 ectodomain. IgG recognition was assessed by ELISA using sera of P. vivax-infected individuals collected from endemic regions of Brazil or antibodies raised in immunized mice. The frequencies of responders to recombinant proteins containing the DII were higher than the others and similar to the ones observed against the PvAMA-1 ectodomain. Moreover, ELISA inhibition assays using the PvAMA-1 ectodomain as substrate revealed the presence of many common epitopes within DI-II that are recognized by human immune antibodies. Finally, immunization of mice with the PvAMA-1 ectodomain induced high levels of antibodies predominantly to DI-II. Together, our results indicate that DII is particularly immunogenic during natural human infections, thus indicating that this region could be used as part of an experimental sub-unit vaccine to prevent vivax malaria.

Polymorphism at the apical membrane antigen 1 locus reflects the world population history of Plasmodium vivax

Background

In malaria parasites (genus Plasmodium), ama-1 is a highly polymorphic locus encoding the Apical Membrane Protein-1, and there is evidence that the polymorphism at this locus is selectively maintained. We tested the hypothesis that polymorphism at the ama-1 locus reflects population history in Plasmodium vivax, which is believed to have originated in Southeast Asia and is widely geographically distributed. In particular, we tested for a signature of the introduction of P. vivax into the New World at the time of the European conquest and African slave trade and subsequent population expansion.

Results

One hundred and five ama-1 sequences were generated and analyzed from samples from six different Brazilian states and compared with database sequences from the Old World. Old World populations of P. vivax showed substantial evidence of population substructure, with high sequence divergence among localities at both synonymous and nonsynonymous sites, while Brazilian isolates showed reduced diversity and little population substructure.

Conclusion

These results show that genetic diversity in P. vivax AMA-1 reflects population history, with population substructure characterizing long-established Old World populations, whereas Brazilian populations show evidence of loss of diversity and recent population expansion.

Note

Nucleotide sequence data reported is this paper are available in the GenBank™ database under the accession numbers EF031154 – EF031216 and EF057446 – EF057487