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

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.

Structural diversity, natural selection and intragenic recombination in the Plasmodium vivax merozoite surface protein 9 locus in Thailand

The merozoite surface protein 9 (MSP9) of malarial parasite forms co-ligand complex with the 19 kDa fragment of merozoite surface protein 1 (MSP1) prior to erythrocyte invasion. Interruption of this process could hamper subsequent asexual erythrocytic development of malaria parasites; therefore, these proteins are considered potential vaccine candidates. In Plasmodium vivax, MSP9 (PvMSP9) contains both conserved and polymorphic repetitive domains that were immunogenic upon natural malaria exposure and conferred protection in vaccination studies in animal models. To investigate the extent of sequence diversity at this locus, 104 P. vivax isolates from 4 major malaria endemic areas of Thailand were analyzed. Results revealed that pvmsp9 contained 3 repeat domains (R1-R3) flanked by conserved domains. Repeat domains exhibit extensive sequence and length variation, in which 14, 39 and 16 haplotypes for domains R1-R3, respectively, circulated in this country. Sequence diversity in pvmsp9 among P. vivax isolates from each endemic area displayed population structure. The extent of sequence diversity in pvmsp9 isolates from the provinces of Tak, Chanthaburi, Ubon Ratchathani and Prachuap Khiri Khan in northwestern, eastern, northeastern and southwestern areas, respectively, was almost comparable and was remarkably higher than that from Yala/Narathiwat population in southern Thailand. Evidence for intragenic recombination in this locus was observed within each P. vivax population except among isolates from Yala and Narathiwat. Synonymous nucleotide diversity significantly exceeded nonsynonymous nucleotide diversity in domains R2 and R3, indicating purifying selection. However, micro-scale signatures of positive and negative selections occurred in both conserved and repeat domains, implying two opposing forces, probably from functional or structural constraint and host immune pressure, could have influenced diversity at this locus. The immunodominant T and B cell epitopes so far identified were invariant or highly conserved across isolates. Further analysis of global isolates is warranted for vaccine design based on this protein.

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.

Polymorphism in merozoite surface protein-7E of Plasmodium vivax in Thailand: Natural selection related to protein secondary structure

Merozoite surface protein 7 (MSP-7) is a multigene family expressed during malaria blood-stage infection. MSP-7 forms complex with MSP-1 prior to merozoite egress from erythrocytes, and could affect merozoite invasion of erythrocytes. To characterize sequence variation in the orthologue in P. vivax (PvMSP-7), a gene member encoding PvMSP-7E was analyzed among 92 Thai isolates collected from 3 major endemic areas of Thailand (Northwest: Tak, Northeast: Ubon Ratchathani, and South: Yala and Narathiwat provinces). In total, 52 distinct haplotypes were found to circulate in these areas. Although population structure based on this locus was observed between each endemic area, no genetic differentiation occurred between populations collected from different periods in the same endemic area, suggesting spatial but not temporal genetic variation. Sequence microheterogeneity in both N- and C- terminal regions was predicted to display 4 and 6 α-helical domains, respectively. Signals of purifying selection were observed in α-helices II-X, suggesting structural or functional constraint in these domains. By contrast, α-helix-I spanning the putative signal peptide was under positive selection, in which amino acid substitutions could alter predicted CD4+ T helper cell epitopes. The central region of PvMSP-7E comprised the 5’-trimorphic and the 3’-dimorphic subregions. Positive selection was identified in the 3’ dimorphic subregion of the central domain. A consensus of intrinsically unstructured or disordered protein was predicted to encompass the entire central domain that contained a number of putative B cell epitopes and putative protein binding regions. Evidences of intragenic recombination were more common in the central region than the remainders of the gene. These results suggest that the extent of sequence variation, recombination events and selective pressures in the PvMSP-7E locus seem to be differentially affected by protein secondary structure.

Molecular approaches to determine the multiplicity of Plasmodium infections

Multiplicity of infection (MOI), also termed complexity of infection (COI), is defined as the number of genetically distinct parasite strains co-infecting a single host, which is an important indicator of malaria epidemiology. PCR-based genotyping often underestimates MOI. Next generation sequencing technologies provide much more accurate and genome-wide characterization of polyclonal infections. However, complete haplotype characterization of multiclonal infections remains a challenge due to PCR artifacts and sequencing errors, and requires efficient computational tools. In this review, the advantages and limitations of current molecular approaches to determine multiplicity of malaria parasite infection are discussed.

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.

Analysis of polymorphisms and selective pressures on ama1 gene in Plasmodium knowlesi isolates from Sabah, Malaysia

The apical membrane antigen-1 (AMA-1) of Plasmodium spp. is a merozoite surface antigen that is essential for the recognition and invasion of erythrocytes. Polymorphisms occurring in this surface antigen will cause major obstacles in developing effective malaria vaccines based on AMA-1. The objective of this study was to characterize ama1 gene in Plasmodium knowlesi isolates from Sabah. DNA was extracted from blood samples collected from Keningau, Kota Kinabalu and Kudat. The Pkama1 gene was amplified using nested PCR and subjected to bidirectional sequencing. Analysis of DNA sequence revealed that most of the nucleotide polymorphisms were synonymous and concentrated in domain I of PkAMA-1. Forteen haplotypes were identified based on amino acid variations and haplotype K5 was the most common haplotype. d_N/d_S ratios implied that purifying selection was prevalent in Pkama1 gene. Fu and Li's D and F values further provided evidence of negative selection acting on domain II of Pkama1. Lownucleotide diversitywas also detected for the Pkama1 sequences,which is similar to reports on Pkama1 from Peninsular Malaysia and Sarawak. The presence of purifying selection and low nucleotide diversity indicated that domain II of Pkama1 can be used as a target for vaccine development.

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.

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.

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.

The Plasmodium vivax merozoite surface protein 3β sequence reveals contrasting parasite populations in southern and northwestern Thailand

BACKGROUND

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

METHODOLOGY

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

PRINCIPAL FINDINGS

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77 ± 0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD=1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD=0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

CONCLUSION

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.

Invasion-inhibitory antibodies elicited by immunization with Plasmodium vivax apical membrane antigen-1 expressed in Pichia pastoris yeast

In a recent vaccine trial performed with African children, immunization with a recombinant protein based on Plasmodium falciparum apical membrane antigen 1 (AMA-1) conferred a significant degree of strain-specific resistance against malaria. To contribute to the efforts of generating a vaccine against Plasmodium vivax malaria, we expressed the ectodomain of P. vivax AMA-1 (PvAMA-1) as a secreted soluble protein in the methylotrophic yeast Pichia pastoris. Recognized by a high percentage of sera from individuals infected by P. vivax, this recombinant protein was found to have maintained its antigenicity. The immunogenicity of this protein was evaluated in mice using immunization protocols that included homologous and heterologous prime-boost strategies with plasmid DNA and recombinant protein. We used the following formulations containing different adjuvants: aluminum salts (Alum), Bordetella pertussis monophosphoryl lipid A (MPLA), flagellin FliC from Salmonella enterica serovar Typhimurium, saponin Quil A, or incomplete Freund's adjuvant (IFA). The formulations containing the adjuvants Quil A or IFA elicited the highest IgG antibody titers. Significant antibody titers were also obtained using a formulation developed for human use containing MPLA or Alum plus MPLA. Recombinant PvAMA-1 produced under "conditions of good laboratory practice" provided a good yield, high purity, low endotoxin levels, and no microbial contaminants and reproduced the experimental immunizations. Most relevant for vaccine development was the fact that immunization with PvAMA-1 elicited invasion-inhibitory antibodies against different Asian isolates of P. vivax. Our results show that AMA-1 expressed in P. pastoris is a promising antigen for use in future preclinical and clinical studies.

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.

Genetic Variation and Selection of Domain I of the Plasmodium vivax Apical Membrane Antigen-1(AMA-1) Gene in Clinical Isolates from Iran

BACKGROUND

Apical Membrane antigen 1 (AMA-1) is positioned on the surface of merozoite and it may play a role in attack to red blood cells. The main aim of present study was to determine the genetic variation, as well as, to detect of selection at domain I of AMA-1 gene Plasmodium vivax isolates in Iran.

METHODS

Blood samples were collected from 58 patients positive for P. vivax, mono infection and the domain I of AMA-1 gene was amplified by nested PCR and then sequenced.

RESULTS

A total 33 different haplotypes were identified among 58 Iranian sequences. The 23 new haplotypes were determined in this study that was not reported previously in other regions of the world. There were totally observed 36 point mutations at the nucleotide level in the analyzed sequences. Sequences analyses indicated 25 amino acid changes at 20 positions in which 5 sites demonstrated thrimorphic polymorphism and the others were dimorphic in the domain I of the Iranian PvAMA-1 isolates.

CONCLUSION

Our findings indicated relatively high level of allelic diversity at the domain I of PvAMA-1 among P.vivax isolates of Iran. Since, PvAMA-1 is considering as vaccine candidate antigen, these data provide valuable information for the development of a PvAMA-1 based malaria vaccine.

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.

Genomics, population genetics and evolutionary history of Plasmodium vivax

Plasmodium vivax is part of a highly diverse clade that includes several Plasmodium species found in nonhuman primates from Southeast Asia. The diversity of primate malarias in Asia is staggering; nevertheless, their origin was relatively recent in the evolution of Plasmodium. We discuss how humans acquired the lineage leading to P. vivax from a nonhuman primate determined by the complex geological processes that took place in Southeast Asia during the last few million years. We conclude that widespread population genomic investigations are needed in order to understand the demographic processes involved in the expansion of P. vivax in the human populations. India represents one of the few countries with widespread vivax malaria. Earlier studies have indicated high genetic polymorphism at antigenic loci and no evidence for geographic structuring. However, new studies using genetic markers in selectively neutral genetic regions indicate that Indian P. vivax presents complex evolutionary history but possesses features consistent with being part of the ancestral distribution range of this species. Such studies are possible due to the availability of the first P. vivax genome sequences. Next generation sequencing technologies are now paving the way for the sequencing of more P. vivax genomes that will dramatically increase our understanding of the unique biology of this species.

Natural acquisition of immunity to Plasmodium vivax: epidemiological observations and potential targets

Population studies show that individuals acquire immunity to Plasmodium vivax more quickly than Plasmodium falciparum irrespective of overall transmission intensity, resulting in the peak burden of P. vivax malaria in younger age groups. Similarly, actively induced P. vivax infections in malaria therapy patients resulted in faster and generally more strain-transcending acquisition of immunity than P. falciparum infections. The mechanisms behind the more rapid acquisition of immunity to P. vivax are poorly understood. Natural acquired immune responses to P. vivax target both pre-erythrocytic and blood-stage antigens and include humoral and cellular components. To date, only a few studies have investigated the association of these immune responses with protection, with most studies focussing on a few merozoite antigens (such as the Pv Duffy binding protein (PvDBP), the Pv reticulocyte binding proteins (PvRBPs), or the Pv merozoite surface proteins (PvMSP1, 3 & 9)) or the circumsporozoite protein (PvCSP). Naturally acquired transmission-blocking (TB) immunity (TBI) was also found in several populations. Although limited, these data support the premise that developing a multi-stage P. vivax vaccine may be feasible and is worth pursuing.

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.

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.

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.

CD8 T-cell activation in mice injected with a plasmid DNA vaccine encoding AMA-1 of the reemerging Korean Plasmodium vivax

Relatively little has been studied on the AMA-1 vaccine against Plasmodium vivax and on the plasmid DNA vaccine encoding P. vivax AMA-1 (PvAMA-1). In the present study, a plasmid DNA vaccine encoding AMA-1 of the reemerging Korean P. vivax has been constructed and a preliminary study was done on its cellular immunogenicity to recipient BALB/c mice. The PvAMA-1 gene was cloned and expressed in the plasmid vector UBpcAMA-1, and a protein band of approximately 56.8 kDa was obtained from the transfected COS7 cells. BALB/c mice were immunized intramuscularly or using a gene gun 4 times with the vaccine, and the proportions of splenic T-cell subsets were examined by fluorocytometry at week 2 after the last injection. The spleen cells from intramuscularly injected mice revealed no significant changes in the proportions of CD8(+) T-cells and CD4(+) T-cells. However, in mice immunized using a gene gun, significantly higher (P<0.05) proportions of CD8(+) cells were observed compared to UB vector-injected control mice. The results indicated that cellular immunogenicity of the plasmid DNA vaccine encoding AMA-1 of the reemerging Korean P. vivax was weak when it was injected intramuscularly; however, a promising effect was observed using the gene gun injection technique.

Genetic diversity of the Plasmodium vivax merozoite surface protein-5 locus from diverse geographic origins

Plasmodium vivax merozoite surface protein-5 (PvMsp-5), a potential vaccine candidate, is encoded by a two-exon single copy gene. We have conducted a comprehensive analysis of PvMsp-5 by sequencing the entire gene of four parasite populations from northwestern Thailand (n=73), southern Thailand (n=53), Indonesia (n=25) and Brazil (n=24), and five isolates from other endemic areas. Results reveal that exon I exhibits a significantly higher level of nucleotide diversity at both synonymous and nonsynonymous sites than exon II (p<0.01). Neutrality tests based on both intraspecific and interspecific nucleotide polymorphism have detected a signature of positive selection in exon I of all populations while substitutions in exon II mainly followed neutral expectation except that three residues in exon II of northwestern Thailand population appear to be positively selected using the Bayes Empirical Bayes method. Short imperfect repeats were identified in exon I at an equivalent region to its orthologue in P. knowlesi, supporting their close genetic relatedness. Significant levels of population subdivision were detected among most populations including those between northwestern and southern Thailand (p<10(-5)), implying absent or minimal gene flow between these populations. Importantly, evidences for intragenic recombination in PvMsp-5 were found in most populations except that from southern Thailand in which haplotype diversity and nucleotide diversity were exceptionally low. Results from Fu and Li's D*, F* and D and F tests suggested that PvMsp-5 of most P. vivax populations have been maintained by balancing selection whereas southern Thailand population could have gone through recent bottleneck events. These findings are concordant with a substantial reduction in the number of P. vivax cases in southern Thailand during the past decade, followed by a very recent population expansion. Therefore, spatio-temporal monitoring of parasite population genetics provides important implications for disease control.