- Open Access
Analysis of Pvama1 genes from China-Myanmar border reveals little regional genetic differentiation of Plasmodium vivax populations
Parasites & Vectorsvolume 9, Article number: 614 (2016)
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.
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.
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.
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.
In South and Southeast Asia, Latin America and Oceania, Plasmodium vivax is the major malaria species; around 2.5 billion people are living in areas of P. vivax transmission . In the Greater Mekong Subregion (GMS) where countries are pursuing regional malaria elimination (www.apmen.org), P. vivax is often sympatric with P. falciparum, P. ovale and P. malariae , although P. vivax has become the predominant species in recent years . With the ability to form hypnozoites that are responsible for relapses, P. vivax presents a major challenge for malaria elimination. Within the GMS, malaria is distributed very unevenly; malaria transmission is concentrated along international borders, whereas central plains are mostly free from malaria [3, 4]. The more intensified control efforts in this region have led to a further reduction of malaria incidence, creating isolated areas or pockets of high malaria prevalence separated by areas with extremely low endemicity or malaria-free zones. In China, autochthonous malaria incidence is mostly located in counties bordering with Myanmar, where malaria burden is the highest in the GMS [5, 6]. In these border regions, cross-border human migration as a major source of malaria introduction presents a significant challenge to the malaria elimination course [7, 8]. Since control efforts are expected to have great impacts on the genetic diversity of the parasite populations [9, 10], tracking their spatial and temporal dynamics may provide timely measurement of the progress of regional malaria elimination.
The genetic diversity of antigens in malaria parasites has been extensively studied not only because of their importance as malaria vaccine candidates , but also due to their usefulness as molecular markers for differentiating parasite populations. Several P. vivax proteins, including Duffy-binding protein (DBP), apical membrane antigen 1 (AMA1), and merozoite surface proteins (MSPs), have been selected as vaccine candidates for their essential functions during erythrocytes invasion and their antigenicity in natural host immune response [12–15]. Among them, AMA1 has been identified as an essential target of the host immune system, and considered an attractive malaria vaccine candidate [15–17]. The ama1 gene has been extensively studied in a number of Plasmodium species . As a type I transmembrane protein, AMA1 is secreted by microneme organelles. Together with RON proteins, AMA1 is involved in merozoite reorientation and tight junction formation during the invasion process [19–24]. Antibodies raised against the AMA1 ectodomain have been shown to inhibit erythrocyte invasion, and AMA1 immunization protects against malaria infection [15, 25–27]. The ectodomain of AMA1 was divided into three subdomains referred to as Domain I (DI), Domain II (DII) and Domain III (DIII) based on the conserved cysteine residues . DI harbors higher levels of genetic variation compared to DII and DIII, suggesting this domain is a target of the host immune system . Within DI of PfAMA1, eight polymorphic amino acids located in the cluster 1 loop (c1L) were identified as the targets of allele-specific, protective immune response . Evidence of diversifying selection was observed in DII of AMA1 in some studies such as in the Sri Lankan parasites, suggesting that this region may also be targeted by host immunity [31, 32]. In addition, serological studies showed that DII is the most immunogenic of the three domains . Due to the highly polymorphic feature of the ama1 gene, it has been used as a molecular marker for population genetic studies [34, 35].
Although the genetic diversity of P. vivax ama1 (Pvama1) has been extensively studied with samples from Asia, Oceania and South America [36–39], there is a relative lack of data from Southeast Asia, where parasite populations are now being restricted to the international borders as the regional malaria elimination program is unfolding. Thus, in this study, we aimed to investigate the genetic diversity of Pvama1 gene from the China-Myanmar border area. We wanted to determine whether this potentially isolated parasite population was genetically different from those of other endemic regions. Furthermore, we wanted to learn whether the intensified malaria control efforts resulted in evident reduction of genetic diversity of the parasite population as observed elsewhere in Southeast Asia . Our study revealed a considerably high level of genetic diversity of the Pvama1 gene in the China-Myanmar border parasite population.
Study sites and isolates
Parasite sample collection was carried out from April 2011 to October 2012 around the township of Laiza, which is located in the northeast Kachin State, along the China-Myanmar border . This region has a subtropical climate with most of the precipitation occurring in June-August (the rainy season). The majority of the local inhabitants are ethnic Kachin (called Jingpo in China) with farming being their major occupation. In recent years, there has been a considerable decline of annual malaria incidence and P. vivax has become the predominant parasite species [40, 41]. Like other international borders within the GMS, the China-Myanmar border also has experienced an increase in cross-border migration with large migrant populations engaging in international trade, logging, quarry, plantation and construction activities, which may play an important role for introducing malaria [8, 42]. Filter paper blood spots were collected as our passive malaria case detection efforts in the township hospital and clinics servicing surrounding villages and two camps for internally displaced people. P. vivax infections were based on microscopy of Giemsa-stained thick smears. The study protocol was approved by the IRB of Pennsylvania State University and by the Health Department of Kachin.
DNA extraction, amplification and sequencing of the PvAMA1 gene
Genomic DNA from dried filter paper blood spots of slide-positive samples was extracted by using a QIAamp DNA Blood Mini kit (QIAGEN, Hilden, Germany). A DNA fragment of the Pvama1 gene encoding the extracellular region (1,290 bp, nucleotide region (nt) 148–1,437 of the Sal-I sequence; GenBank accession no. AF063138) was amplified and sequenced with six oligonucleotide primers listed in Additional file 1: Table S1. The amplifications were performed in a 20 μl volume reaction containing 1× KOD-Plus-Neo buffer, 200 μM dNTPs, 1 mM MgSO4, 250 nM of each primer, 0.4 units of KOD Plus neo polymerase (Toyobo, Osaka, Japan), and 1 μl genomic DNA as template. The reaction was run at 94 °C for 5 min, followed by 45 cycles of 94 °C for 15 s, 56 °C for 15 s, and 68 °C for 90 s, and extension at 68 °C for 5 min. PCR fragments were analyzed by electrophoresis on a 1.2% agarose gels. Direct sequencing of purified PCR fragments was carried out on both strands with primers listed in Additional file 1: Table S1 by using the ABI BigDye™ Terminator Reaction Ready kit (Applied Biosystems, CA, USA).
Sequence assembling and polymorphism analysis
Of the 78 P. vivax infections amplified by PCR, five samples showing dual peaks suggestive of a mixed infection were excluded from further analysis. Pvama1 was successfully sequenced from 73 samples (sequences submitted to GenBank under accession nos. KX495505–KX495577). A single contiguous 1,290 bp of Pvama1 (nt 148–1,437 and codons 50–479) was derived for each of the 73 Pvama1 sequences, which include DI (codons 94–247), DII (codons 265–363) and DIII (codons 388–451). The consensus of the sequences were aligned to the P. vivax Sal-I strain by using the CLUSTAL W program in MEGA6.0 . Seven primate-adapted P. vivax isolates (Belem, Chesson I, India VII, Indonesia XIX, North Korea, Palo Alto and Simium; GenBank accession nos EU395595–EU395601) were also analyzed as additional reference sequences . Aligned sequences were exported as FASTA alignment for statistical analysis using the DnaSP v5.10.01 software . Sequencing was repeated to confirm singletons. The numbers of segregating sites (S), the total number of mutations (η), the average number of pairwise nucleotide differences (k), nucleotide diversity (π), the number of haplotypes (H) and haplotype diversity (Hd), and their corresponding standard deviation were computed using the options available in DnaSP v5.10.01 software . The distribution of nucleotide diversity (π) across the DI-III of Pvama1 gene was analyzed using the sliding window approach.
Haplotype network construction, linkage disequilibrium (LD) and F ST analysis
A haplotype network based on the Pvama1 sequence was constructed by using the NETWORK software Version 188.8.131.52 with the Median-Joining method . A total of 534 publically available Pvama1 sequences representing seven P. vivax populations were retrieved from the GenBank: South Korea , India , Sri Lanka , Thailand , Iran , Venezuela , and Papua New Guinea (PNG) . LD was performed to calculate the minimum number of recombination events (R M)  and to give an estimation of the recombination parameter C . The relationship between LD and distance between nucleotide sites was plotted by using indices D’  and R 2 . To assess the proportion of genetic variance due to population subdivision, Wright’s F ST among populations were calculated by using DnaSP v5.10.01 .
Statistical and phylogenetic analysis
To determine departure from neutrality, Tajima’s D, Fu and Li’s D* and F*, and McDonald-Kreitman indices were calculated using DnaSP v5.10.01 and via a sliding window method . In Tajima’s D test, departure from neutrality is determined by comparing the values of θ (estimated nucleotide diversity) derived from π (observed average pairwise nucleotide diversity) and the total number of segregating sites (S) . The Fu and Li’s D* and F* tests determine departure from neutrality by differences between estimation of θ derived from the number of singletons and that based on either the total number of mutations (D*) or the average pairwise diversity (F*) . The McDonald-Kreitman test  was applied by using a single P. cynomolgi ama1 sequence (GenBank accession no. X86099; ) as the outgroup for comparison with the P. vivax populations. It allows for determination of the ratio of synonymous substitutions to nonsynonymous substitutions between and within species. A two-tailed Fisher’s exact test was computed to determine the statistical significance (P < 0.05). A phylogenetic analysis of Pvama1 sequences (nt 322–737 relative to Sal-I strain) with Myanmar samples  was performed by the neighbor-joining (NJ) method in the MEGA 6.0 software . The bootstrap method with 1,000 replications was used to assess the reliability of the gene tree.
Structural modelling of PvAMA1
The potential B cell epitopes of DI-III in PvAMA1 were predicted by using the ABCpred server . In order to consider a given region as a valid linear B-cell epitope, a threshold of 0.7 was used to predict the 16 residues. The prediction of intrinsically unstructured/disordered regions (IURs) in proteins, which is highly desirable in the design of vaccines and diagnostic tests, was carried out using the RONN server . Residues with an output score above 0.5 were regarded as disordered. The overlapped regions between B-cell epitopes and IURs, as well as mutation sites across DI-III were mapped on a three dimensional structure of PvAMA1 (Protein Data Bank ID: 1W8K; ). ViewerLite 4.2 was used to visualize the surface distribution of mapped IURs, the polymorphism, and B-cell epitopes on the crystal structure of PvAMA1.
Genetic diversity of Pvama1 among P. vivax isolates from the China-Myanmar border
The 1,290 bp sequence corresponding to the ectodomain (nt 148–1,437 encoding amino acids 50–479 relative to the Sal-I strain) of Pvama1 was amplified from 73 P. vivax samples acquired from the China-Myanmar border. Of the 73 Pvama1 sequences, there were 41 haplotypes, giving overall haplotype diversity (Hd) of 0.958 (Table 1). The average number of pairwise nucleotide differences (k) for the entire 1,290 bp sequenced region, DI, DII, and DIII were 10.895, 7.433, 1.385, and 0.419, respectively (Table 1). A total of 46 SNPs was detected, including 13 synonymous and 33 non-synonymous (Table 2). Nucleotide diversity (π) at DI, DII and DIII of Pvama1 among 73 examined P. vivax isolates was 0.016, 0.005, and 0.002, respectively, giving an overall π value of the sequenced region of 0.008 (Table 1). A sliding window plot of π with a window of 90 bp and step size of 3 bp revealed values ranging from 0.000 to 0.023 (Fig. 1a). The ectodomain contained 43 polymorphic sites, of which four (nt 162, 957, 1377, 1412) were novel compared to previous reports [32, 37, 38, 47–50]. The majority of polymorphic sites were identified in DI (26 sites) and DII (5 sites) as compared to two sites in DIII (Fig. 1b). Of these polymorphic sites, five were trimorphic, while the rest were dimorphic, displaying only two alternative amino acids (Fig. 1b). Within DI, 10 of the 26 Pvama1 mutations were mapped to the same positions as the Pfama1 c1, c1L, c2, and c3 clusters. Among them, four mapped to the c1 cluster, and one was located within the Pvama1 c1L cluster (Additional file 2: Figure S1).
Genetic diversity of Pvama1 in the eight worldwide parasite populations
The Pvama1 sequences obtained from the China-Myanmar border (n = 73) were then compared to previously published sequences from seven other geographic regions and primate-adapted P. vivax isolates (nt 148–1,437 encoding amino acids 50–479). A total of 117 mutations were identified among the 615 global and reference Pvama1 sequences, including 38 synonymous SNPs (SP SNPs) and 79 non-synonymous SNPs (NS SNPs) (Table 2). The majority of the NS SNPs (53.2%) were clustered within DI, resulting in a peak nucleotide diversity (π = 0.020) for this region (Additional file 3: Figure S2). A total of 308 haplotypes were identified among these global isolates, demonstrating an extremely high level of haplotype diversity (Hd = 0.988). Nucleotide diversity was the highest in the Iranian and Thai populations and lowest in the South Korean population (Table 2). Despite that the majority of haplotypes was unique to the China-Myanmar border and Myanmar populations, Pvama1 sequences from the China-Myanmar border population did not form distinct clades with the Myanmar isolates  by phylogenetic analyses (Additional file 4: Figure S3).
Recombination and linkage disequilibrium
Analysis of the ectodomain of Pvama1 from the China-Myanmar border samples provided estimates of the minimum number of recombination events of six, while values of the recombination parameter C between adjacent sites and per gene were 0.043 and 25.5, respectively (Table 3). These values within the China-Myanmar border population were higher compared to estimates from the South Korea, India, Sri Lanka and Venezuela populations, but lower than those of the Thailand, PNG and Iran populations, respectively (Table 3). Figure 2 shows the relationship between R 2, D’ with distance between sites. This LD analysis showed a decline of the LD index (R 2) with increasing nucleotide sites distance within the Pvama1 gene of the China-Myanmar border isolates, indicating a high meiotic recombination rate (Fig. 2).
F ST analysis
In order to understand the distribution of diversity across global populations, F ST values of the China-Myanmar border population and seven worldwide populations with full-length ectodomain sequences were evaluated . A high level of genetic differentiation (F ST = 0.47) was detected between the South Korean population and the China-Myanmar border population, whereas a low to moderate level of genetic differentiation was detected among other worldwide populations (ranging from 0.03 to 0.24). A moderate range of F ST values (0.13–0.16) were detected when comparing the China-Myanmar border populations with the India, Iran, Venezuela, and PNG populations, but a much lower genetic difference was identified when comparing the China-Myanmar border population with the Thai population (F ST = 0.03) (Table 4). Low F ST values (0.01–0.10) were also observed when comparing part of Pvama1 DI (nt 322–737) among Myanmar, Thai and China-Myanmar border populations (Additional file 5: Table S2).
Evidence of positive diversifying selection on Pvama1 gene
To examine whether natural selection contributed to the generation of the diversity in Pvama1 of the China-Myanmar border P. vivax population, Tajima’s D and Fu and Li’s D* and F* were performed. The values of Tajima’s D and Fu and Li’s D* and F* indices for the entire ectodomain and each domain separately are displayed in Table 1. Even though these values were not statistically significant, the positive value indicates deviation from neutral evolution and the tendency for positive diversifying selection. Meanwhile, a sliding window plot depicted significant positive values in DI, suggesting positive diversifying selection in this region (Fig. 3). Examination of the two sequence sets using a McDonald-Kreitman test showed that the entire sequenced region, particularly within DI, had significantly more non-synonymous substitutions than expected from comparison with P. cynomolgi (P < 0.0001 and P = 0.0037). This result suggests that polymorphisms found for Pvama1, especially in the DI, are maintained by diversifying selection, presumably due to host immune pressure.
Haplotype network reconstruction
On the basis of the 34 polymorphisms in the Pvama1 DI global isolates, 162 unique haplotypes were identified in 615 sequences, of which 54.9% was singleton. Individual haplotype prevalence ranged from 0.6 to 45.1%. Of the reference strain Pvama1 sequences, only the Belem haplotype was detected, but was restricted to the Thai population. No haplotype matched that of the Sal-I strain. Twelve of the 27 China-Myanmar border haplotypes (44.4%) were shared with other populations, of which 91.7% (11/12) were identical to some of the Pvama1 haplotypes observed in the Thai population. Haplotype 21 is the predominant haplotype both in the China-Myanmar border and Thai populations, with a frequency of 24.7% and 23.8%, respectively. Intriguingly, an independent haplotype analysis performed among the China-Myanmar border, Myanmar, and Thai populations indicates that Haplotype 21 is also the only haplotype shared among these three populations (Additional file 6: Figure S4). Haplotype 10 is shared among the China-Myanmar border, Indian, Sri Lankan, Thai and Iranian populations, with a frequency of 2.7, 10.0, 30.4, 2.6, and 10.3%, respectively. Haplotype 30 is the only haplotype shared between the China-Myanmar border population and the outside Asian population. The haplotype network, drawn by excluding the 89 singletons and ≤ 2% frequency haplotypes from the analysis, showed that clusters from the South American, Oceania, and Asian populations overlapped (Fig. 4).
Analysis of PvAMA1 structure and antigenicity
To determine the distribution of the B-cell epitopes, IURs, and the polymorphic residues found within the China-Myanmar border population, the modeled structure of PvAMA1 was used. The potential B-cell epitopes were identified across the ectodomain of PvAMA1, which also overlap with the locations of most of the SNPs (Additional file 7: Figure S5). Also, the putative IURs are distributed in five regions (amino acids - aa: 121–141, 171–183, 291–309, 392–406, and 427–441) within all three domains of the PvAMA1 protein with no polymorphic residues detected in aa 171–183 and 291–309. In Fig. 5, the overlap between the putative B-cell epitopes and IURs have been demonstrated in the 3-D model of PvAMA1 in yellow (with polymorphic residues, aa: 121–141, 392–406, and 427–441) and red colors (without polymorphic residues, aa: 171–183 and 291–309), respectively. The polymorphic residues clustered on one side of the PvAMA1 protein. Interestingly, the two overlapping non-polymorphic residues of the B-cell epitopes and IURs are located on the opposing faces of the PvAMA1 protein (Fig. 5).
Consistent with worldwide reports, we detected high genetic diversity in the ectodomain of Pvama1 sequences obtained in 73 P. vivax isolates from the China-Myanmar border. Globally, substantial geographic variability has been observed, with lower diversity in East Asian (Korean, π = 0.006 ) and South American (Venezuelan, π = 0.007 ) populations, compared to higher diversity in Southeast Asian (Thai, π = 0.009 ) and West Asian (Iranian, π = 0.010 ) populations. As genetic diversity of P. vivax populations varies significantly due to transmission intensity [62–64], it is likely that the higher genetic diversity observed in Thailand and Iran was due to higher P. vivax transmission intensity than in Korea and Venezuela. The detected high level of genetic diversity in Pvama1 in parasite population from the China-Myanmar border was not in line with the rapid decline of malaria incidence in recent years. While it is expected that genetic diversity diminishes following the reduction in malaria incidence, the relationship between the two aspects is complex and could vary according to many epidemiological parameters. In Thailand, for example, P. vivax population in the western border was highly diverse as compared to being almost clonal in the southern border [9, 10]. It remains to be determined whether these differences were due to different histories of malaria control efforts and levels of parasite introduction from neighboring highly endemic areas. Nevertheless, with the overall increased control efforts in the GMS and more isolated parasite populations on the international borders, future longitudinal tracking of parasite genetic diversity in the local parasite populations as well as in potential sources of parasite introduction will provide answers to these questions. Meanwhile, the use of more than one genetic marker such as a molecular barcode may provide a better resolution of the parasite populations .
The LD indices decline rapidly with the increasing distance between sites. Considering the high values of recombination parameter, C (C = 4Nr), these results indicate that recombination is very frequent in the Pvama1 gene of diverse geographic populations [32, 37, 38, 47–50]. An estimated moderate level of the recombination parameter C (both between adjacent site, Ra and per gene, Rb) was detected for Pvama1 (0.043 and 25.5) in the China-Myanmar border population. This C value observed for Pvama1 is much higher than that of Pfama1 (0.003 and 3.8) from a sympatric population of P. falciparum , for which the effective population size (N e ) is only 1.3 × 103 when calculated by using the average recombination rate of P. falciparum (~6 × 10-7 Morgans/base) . We conclude from these data that recombination events in Pvama1 of the China-Myanmar border population are common, and that this area appears to have a considerably larger effective population of P. vivax than P. falciparum. In P. falciparum, strong LD and low genetic diversity often are observed in low-transmission areas . However, such a spectrum of population structures may not fit for P. vivax populations, as we observed LD and high-level haplotype diversity (0.958) of Pvama1 in the China-Myanmar border, where malaria transmission is low . We speculate that two factors might contribute to the observed high genetic diversity in Pvama1 in the border regions. The ability to relapse in P. vivax increases the prevalence of mixed infections and chances for genetic recombination in the mosquito vectors. In addition, heavy human population flow across the porous international borders further favors mixing of parasite populations and increasing parasite genetic diversity.
Pairwise F ST comparisons only detected a moderate level of genetic differentiation (F ST = 0.23) among global populations. This may reflect the functional constraints of this protein and similar balancing selection from human immunity on this protein. Despite this, major geographical genetic differentiation does exist. For instance, the Korean parasite population was drastically different from other parasite populations from the world. The Thai P. vivax populations collected from three regions (northwestern, eastern and southern) from 1996 and 2007 also showed significant population subdivision , which may be due to temporal and spatial separations among these populations. Later studies confirmed that the southern Thai population had experienced a genetic bottleneck with significantly reduced genetic diversity [9, 10]. Yet, very little genetic differentiation was observed between the combined Thai parasite population and the China-Myanmar border population, despite that our parasite samples were acquired more recently. This demonstrates the persistence of highly diverse P. vivax populations with only slight genetic differentiation in different regions of the GMS, suggesting the lack of major gene flow barriers in this region. Together with the fact that similar allelic forms of Pvama1 are circulating in the Myanmar, China-Myanmar border, and Thai P. vivax populations [48, 59], the low genetic differentiation between these regions support the idea that a general malaria elimination or control strategy may be applicable for these endemic regions.
The pattern of polymorphisms observed in the ectodomain as well as DI of the China-Myanmar border Pvama1 showed that both domains are under diversifying selection, suggesting a significant departure from neutrality as confirmed by the McDonald and Kreitman test. These results are consistent with previous reports of Pvama1 DI from various global isolates [32, 37, 38, 47–50]. Although DII of Pvama1 has been reported as highly immunogenic [33, 69], and there is evidence of positive selection in DII of Sri Lankan P. vivax populations , no evidence was found for diversifying selection on DII and DIII of the China-Myanmar border Pvama1 as confirmed by neutrality tests, which is in agreement with several other reports [37, 50]. These results suggest that DI is the dominant target of the host immune response.
A total of 308 haplotypes, with a haplotype diversity of 0.988, was identified in global isolates when analyzing the ectodomain of Pvama1. This extreme polymorphism indicates that each haplotype might be recognized as immunologically distinct by the immune system. Further observation of Pvama1 DI showed that only haplotype 30 in the China-Myanmar border is shared with populations outside of Asia, and haplotype 10 is the only one shared among the majority of analyzed Asian populations. Together with the complex network analysis results, covering diversity of P. vivax will be a great challenge for a PvAMA1-based vaccine design. Although highly diverse, the observation that the majority of the China-Myanmar border Pvama1 haplotypes (44.4%) are shared with Thai populations and that the dominant haplotype (H21) is identical in both populations suggested little differentiation of parasite populations within GMS, and a PvAMA1-based malaria vaccine may be effective in this entire region. Notably, the vaccine strain Sal-I haplotype was not identified in any of the analyzed population in the current study. Thus, the dominant Pvama1 alleles from field parasite populations should be considered in the future studies.
Compared to the Sal-I sequence, 46 SNPs resulting in 43 amino acid substitutions were identified in the China-Myanmar border Pvama1. The majority of the diversity clustered within Pvama1 DI, as has been reported previously for P. falciparum and P. vivax [29, 37, 50]. The 90.7% (39/43) polymorphic sites within the ectodomain identified in the China-Myanmar border Pvama1 are comparable to those previously reported [32, 37, 38, 47–50], indicating that these sites may be under strong natural selection. Four clusters within DI (c1, c1L, c2 and c3) were mapped to the surface of the PfAMA1 protein and are potentially associated with antigenic escape [70, 71]. In P. falciparum, the c1 (positions 187 to 231) and c1L (positions 196–207) clusters contain the most polymorphic residues [70, 71]. However, only four Pvama1 polymorphic sites of the China-Myanmar border isolates were located in the Pfama1 c1 cluster. These mutation sites are comparable with those observed in the PNG isolates , supporting the hypothesis that the presentation of the PvAMA1’s DI loops is different from that of PfAMA1. Structural modeling of PvAMA1 from the China-Myanmar border isolates revealed that all of the polymorphic residues are mapped to one surface of the PvAMA1 protein, suggesting this is the side face exposed to the host immune system . Previous studies by using invasion-inhibitory monoclonal antibodies specific for PfAMA1 revealed that DII of PvAMA1 is a target of the protective immune response. Furthermore, IUR region (aa 290–307) in the DII of PvAMA1 was found to be an important antigenic region during natural human infections . Comparable to these previous reports, amino acid sequence conservation has also been identified in China-Myanmar border isolates, with no polymorphic sites in aa 171–183 (in DI) and 291–309 (in DII) of B-cell epitope and IUR overlapped region have been identified in all 73 sequenced samples. The current results support the inclusion of amino acid sequences aa 291–309 in DII in PvAMA1-based vaccine.
Our study demonstrated that the genetic diversity of the Pvama1 in 73 P. vivax isolates from the China-Myanmar border area, as well as global isolates, was exceptionally high, and that Pvama1 DI is the dominant target of positive diversifying selection. Furthermore, a low level of genetic differentiation among the Myanmar, China-Myanmar border and Thai populations suggests high levels of gene flow within the GMS.
Cluster 1 loop
- F ST :
Intrinsically unstructured/disordered regions
- MK :
McDonald Kreitman test
Merozoite surface protein 1
- NS SNPs:
Papua New Guinea
Plasmodium vivax apical membrane antigen 1
- SP SNPs:
Gething PW, Elyazar IR, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl Trop Dis. 2012;6:e1814.
Zhou X, Huang JL, Njuabe MT, Li SG, Chen JH, Zhou XN. A molecular survey of febrile cases in malaria-endemic areas along China-Myanmar border in Yunnan province, People’s Republic of China. Parasite. 2014;21:27.
Hewitt S, Delacollette C, Chavez I. Malaria situation in the Greater Mekong Subregion. Southeast Asian J Trop Med Public Health. 2013;44 Suppl 1:46–72. discussion 306–307.
Cui L, Yan G, Sattabongkot J, Cao Y, Chen B, Chen X, et al. Malaria in the Greater Mekong Subregion: heterogeneity and complexity. Acta Trop. 2012;121:227–39.
Feng J, Xiao H, Xia Z, Zhang L, Xiao N. Analysis of malaria epidemiological characteristics in the People’s Republic of China, 2004–2013. Am J Trop Med Hyg. 2015;93:293–9.
Feng J, Xiao H, Zhang L, Yan H, Feng X, Fang W, Xia Z. The Plasmodium vivax in China: decreased in local cases but increased imported cases from Southeast Asia and Africa. Sci Rep. 2015;5:8847.
Wang D, Li S, Cheng Z, Xiao N, Cotter C, Hwang J, et al. Transmission risk from imported Plasmodium vivax malaria in the China-Myanmar border region. Emerg Infect Dis. 2015;21:1861–4.
Zhou G, Sun L, Xia R, Duan Y, Xu J, Yang H, et al. Clinical malaria along the China-Myanmar border, Yunnan Province, China, January 2011–August 2012. Emerg Infect Dis. 2014;20:675–8.
Putaporntip C, Miao J, Kuamsab N, Sattabongkot J, Sirichaisinthop J, Jongwutiwes S, Cui L. The Plasmodium vivax merozoite surface protein 3beta sequence reveals contrasting parasite populations in southern and northwestern Thailand. PLoS Negl Trop Dis. 2014;8:e3336.
Jongwutiwes S, Putaporntip C, Hughes AL. Bottleneck effects on vaccine-candidate antigen diversity of malaria parasites in Thailand. Vaccine. 2010;28:3112–7.
Rich SM, Ferreira MU, Ayala FJ. The origin of antigenic diversity in Plasmodium falciparum. Parasitol Today. 2000;16:390–6.
Tsuboi T, Kappe SH, Al-Yaman F, Prickett MD, Alpers M, Adams JH. Natural variation within the principal adhesion domain of the Plasmodium vivax duffy binding protein. Infect Immun. 1994;62:5581–6.
de Sousa TN, Kano FS, de Brito CF, Carvalho LH. The Duffy binding protein as a key target for a Plasmodium vivax vaccine: lessons from the Brazilian Amazon. Mem Inst Oswaldo Cruz. 2014;109:608–17.
Versiani FG, Almeida ME, Mariuba LA, Orlandi PP, Nogueira PA. N-terminal Plasmodium vivax merozoite surface protein-1, a potential subunit for malaria vivax vaccine. Clin Dev Immunol. 2013;2013:965841.
Remarque EJ, Faber BW, Kocken CH, Thomas AW. Apical membrane antigen 1: a malaria vaccine candidate in review. Trends Parasitol. 2008;24:74–84.
Mitchell GH, Thomas AW, Margos G, Dluzewski AR, Bannister LH. Apical membrane antigen 1, a major malaria vaccine candidate, mediates the close attachment of invasive merozoites to host red blood cells. Infect Immun. 2004;72:154–8.
Dutta S, Haynes JD, Moch JK, Barbosa A, Lanar DE. Invasion-inhibitory antibodies inhibit proteolytic processing of apical membrane antigen 1 of Plasmodium falciparum merozoites. Proc Natl Acad Sci U S A. 2003;100:12295–300.
Preiser P, Kaviratne M, Khan S, Bannister L, Jarra W. The apical organelles of malaria merozoites: host cell selection, invasion, host immunity and immune evasion. Microbes Infect. 2000;2:1461–77.
Gaur D, Mayer DC, Miller LH. Parasite ligand-host receptor interactions during invasion of erythrocytes by Plasmodium merozoites. Int J Parasitol. 2004;34:1413–29.
Richard D, MacRaild CA, Riglar DT, Chan JA, Foley M, Baum J, et al. Interaction between Plasmodium falciparum apical membrane antigen 1 and the rhoptry neck protein complex defines a key step in the erythrocyte invasion process of malaria parasites. J Biol Chem. 2010;285:14815–22.
Lamarque M, Besteiro S, Papoin J, Roques M, Vulliez-Le Normand B, Morlon-Guyot J, et al. The RON2-AMA1 interaction is a critical step in moving junction-dependent invasion by apicomplexan parasites. PLoS Pathog. 2011;7:e1001276.
Galinski MR, Barnwell JW. Plasmodium vivax: merozoites, invasion of reticulocytes and considerations for malaria vaccine development. Parasitol Today. 1996;12:20–9.
Narum DL, Thomas AW. Differential localization of full-length and processed forms of PF83/AMA-1 an apical membrane antigen of Plasmodium falciparum merozoites. Mol Biochem Parasitol. 1994;67:59–68.
Bannister LH, Hopkins JM, Dluzewski AR, Margos G, Williams IT, Blackman MJ, et al. Plasmodium falciparum apical membrane antigen 1 (PfAMA-1) is translocated within micronemes along subpellicular microtubules during merozoite development. J Cell Sci. 2003;116:3825–34.
Kusi KA, Faber BW, Thomas AW, Remarque EJ. Humoral immune response to mixed PfAMA1 alleles; multivalent PfAMA1 vaccines induce broad specificity. PLoS One. 2009;4:e8110.
Stowers AW, Kennedy MC, Keegan BP, Saul A, Long CA, Miller LH. Vaccination of monkeys with recombinant Plasmodium falciparum apical membrane antigen 1 confers protection against blood-stage malaria. Infect Immun. 2002;70:6961–7.
Anders RF, Crewther PE, Edwards S, Margetts M, Matthew ML, Pollock B, Pye D. Immunisation with recombinant AMA-1 protects mice against infection with Plasmodium chabaudi. Vaccine. 1998;16:240–7.
Pizarro JC, Vulliez-Le Normand B, Chesne-Seck ML, Collins CR, Withers-Martinez C, et al. Crystal structure of the malaria vaccine candidate apical membrane antigen 1. Science. 2005;308:408–11.
Polley SD, Conway DJ. Strong diversifying selection on domains of the Plasmodium falciparum apical membrane antigen 1 gene. Genetics. 2001;158:1505–12.
Ouattara A, Takala-Harrison S, Thera MA, Coulibaly D, Niangaly A, Saye R, et al. Molecular basis of allele-specific efficacy of a blood-stage malaria vaccine: vaccine development implications. J Infect Dis. 2013;207:511–9.
Dias S, Somarathna M, Manamperi A, Escalante AA, Gunasekera AM, Udagama PV. 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. Vaccine. 2011;29:7491–504.
Gunasekera AM, Wickramarachchi T, Neafsey DE, Ganguli I, Perera L, Premaratne PH, et al. Genetic diversity and selection at the Plasmodium vivax apical membrane antigen-1 (PvAMA-1) locus in a Sri Lankan population. Mol Biol Evol. 2007;24:939–47.
Mufalo BC, Gentil F, Bargieri DY, Costa FT, Rodrigues MM, Soares IS. Plasmodium vivax apical membrane antigen-1: comparative recognition of different domains by antibodies induced during natural human infection. Microbes Infect. 2008;10:1266–73.
Joshi H. Markers for population genetic analysis of human plasmodia species, P. falciparum and P. vivax. J Vector Borne Dis. 2003;40:78–83.
Cheng Q, Saul A. Sequence analysis of the apical membrane antigen I (AMA-1) of Plasmodium vivax. Mol Biochem Parasitol. 1994;65:183–7.
Rodrigues MH, Rodrigues KM, Oliveira TR, Comodo AN, Rodrigues MM, Kocken CH, et al. Antibody response of naturally infected individuals to recombinant Plasmodium vivax apical membrane antigen-1. Int J Parasitol. 2005;35:185–92.
Arnott A, Mueller I, Ramsland PA, Siba PM, Reeder JC, Barry AE. Global population structure of the genes encoding the malaria vaccine candidate, Plasmodium vivax apical membrane antigen 1 (PvAMA1). PLoS Negl Trop Dis. 2013;7:e2506.
Rajesh V, Elamaran M, Vidya S, Gowrishankar M, Kochar D, Das A. Plasmodium vivax: genetic diversity of the apical membrane antigen-1 (AMA-1) in isolates from India. Exp Parasitol. 2007;116:252–6.
Chung JY, Chun EH, Chun JH, Kho WG. Analysis of the Plasmodium vivax apical membrane antigen-1 gene from re-emerging Korean isolates. Parasitol Res. 2003;90:325–9.
Li N, Parker DM, Yang Z, Fan Q, Zhou G, Ai G, et al. Risk factors associated with slide positivity among febrile patients in a conflict zone of north-eastern Myanmar along the China-Myanmar border. Malar J. 2013;12:361.
Zhou G, Lo E, Zhong D, Wang X, Wang Y, Malla S, et al. Impact of interventions on malaria in internally displaced persons along the China-Myanmar border: 2011–2014. Malar J. 2016;15:471.
Jitthai N. Migration and malaria. Southeast Asian J Trop Med Public Health. 2013;44 Suppl 1:166–200. discussion 306–167.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Ntumngia FB, McHenry AM, Barnwell JW, Cole-Tobian J, King CL, Adams JH. Genetic variation among Plasmodium vivax isolates adapted to non-human primates and the implication for vaccine development. Am J Trop Med Hyg. 2009;80:218–27.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.
Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Kang JM, Lee J, Cho PY, Moon SU, Ju HL, Ahn SK, et al. Population genetic structure and natural selection of apical membrane antigen-1 in Plasmodium vivax Korean isolates. Malar J. 2015;14:455.
Putaporntip C, Jongwutiwes S, Grynberg P, Cui L, Hughes AL. Nucleotide sequence polymorphism at the apical membrane antigen-1 locus reveals population history of Plasmodium vivax in Thailand. Infect Genet Evol. 2009;9:1295–300.
Zakeri S, Sadeghi H, Mehrizi AA, Djadid ND. Population genetic structure and polymorphism analysis of gene encoding apical membrane antigen-1 (AMA-1) of Iranian Plasmodium vivax wild isolates. Acta Trop. 2013;126:269–79.
Ord RL, Tami A, Sutherland CJ. ama1 genes of sympatric Plasmodium vivax and P. falciparum from Venezuela differ significantly in genetic diversity and recombination frequency. PLoS One. 2008;3:e3366.
Hudson RR, Kaplan NL. Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985;111:147–64.
Hudson RR. Estimating the recombination parameter of a finite population model without selection. Genet Res. 1987;50:245–50.
Lewontin RC. The interaction of selection and linkage. I. General considerations; heterotic models. Genetics. 1964;49:49–67.
Hill WG, Robertson A. Linkage disequilibrium in finite populations. Theor Appl Genet. 1968;38:226–31.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.
McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991;351:652–4.
Dutta S, Malhotra P, Chauhan VS. Sequence analysis of apical membrane antigen 1 (AMA-1) of Plasmodium cynomolgi bastianelli. Mol Biochem Parasitol. 1995;73:267–70.
Moon SU, Na BK, Kang JM, Kim JY, Cho SH, Park YK, et al. Genetic polymorphism and effect of natural selection at domain I of apical membrane antigen-1 (AMA-1) in Plasmodium vivax isolates from Myanmar. Acta Trop. 2010;114:71–5.
Saha S, Raghava GP. Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins. 2006;65:40–8.
Yang ZR, Thomson R, McNeil P, Esnouf RM. RONN: the bio-basis function neural network technique applied to the detection of natively disordered regions in proteins. Bioinformatics. 2005;21:3369–76.
Osier FH, Weedall GD, Verra F, Murungi L, Tetteh KK, Bull P, et al. Allelic diversity and naturally acquired allele-specific antibody responses to Plasmodium falciparum apical membrane antigen 1 in Kenya. Infect Immun. 2010;78:4625–33.
Karunaweera ND, Ferreira MU, Munasinghe A, Barnwell JW, Collins WE, King CL, et al. Extensive microsatellite diversity in the human malaria parasite Plasmodium vivax. Gene. 2008;410:105–12.
Imwong M, Nair S, Pukrittayakamee S, Sudimack D, Williams JT, Mayxay M, et al. Contrasting genetic structure in Plasmodium vivax populations from Asia and South America. Int J Parasitol. 2007;37:1013–22.
Baniecki ML, Faust AL, Schaffner SF, Park DJ, Galinsky K, Daniels RF, et al. Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections. PLoS Negl Trop Dis. 2015;9:e0003539.
Zhu X, Zhao Z, Feng Y, Li P, Liu F, Liu J, et al. Genetic diversity of the Plasmodium falciparum apical membrane antigen I gene in parasite population from the China-Myanmar border area. Infect Genet Evol. 2016;39:155–62.
Su X, Ferdig MT, Huang Y, Huynh CQ, Liu A, You J, et al. A genetic map and recombination parameters of the human malaria parasite Plasmodium falciparum. Science. 1999;286:1351–3.
Anderson TJ, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Mollinedo R, et al. Microsatellite markers reveal a spectrum of population structures in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000;17:1467–82.
Gentil F, Bargieri DY, Leite JA, Francoso KS, Patricio MB, Espindola NM, et al. A recombinant vaccine based on domain II of Plasmodium vivax Apical Membrane Antigen 1 induces high antibody titres in mice. Vaccine. 2010;28:6183–90.
Dutta S, Lee SY, Batchelor AH, Lanar DE. Structural basis of antigenic escape of a malaria vaccine candidate. Proc Natl Acad Sci U S A. 2007;104:12488–93.
Takala SL, Coulibaly D, Thera MA, Batchelor AH, Cummings MP, Escalante AA, et al. Extreme polymorphism in a vaccine antigen and risk of clinical malaria: implications for vaccine development. Sci Transl Med. 2009;1:2ra5.
We would like to thank patients from the surrounding villages for participating in this study and providing the blood samples.
This work was supported by grants from NIAID, National Institutes of Health (U19AI089672) and grants from National Natural Science Foundation of China (Grant no. 81301455).
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
XZ, YC and LC designed the studies and drafted the manuscript. PZ and SW took part in the genomic DNA extraction, PCR amplification and writing of the manuscript. JL, FL, and JW participated in sequence analysis and writing of the manuscript. ZY, QF and GY performed epidemiology study, sample collection and writing of the manuscript. YC and LC revised the final manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study protocol was approved by the IRB of Pennsylvania State University and by the Health Department of Kachin. Written informed consent/assent was obtained from all participants.
Oligonucleotides used for amplification and sequencing of Pvama1 gene sequences. (DOCX 14 kb)
Alignment of P. falciparum, P. vivax, and P. cynomolgi AMA1 protein sequences. P. falciparum 3D7 (PlasmoDB ID: PF3D7_1133400), P. vivax Sal-I (GenBank access no. AF063138) and P. cynomolgi (GenBank access no. X86099) were aligned using MUltiple Sequence Comparison by Log-Expectation (MUSCLE, http://www.ebi.ac.uk/Tools/mas/muscle/). Gaps are indicating by dashes. Conserved amino acids among the three analyzed strains are marked by asterisks. Red bold types indicate conserved cysteine residues that divide the ectodomain of AMA1. Domains I-III is shown in grey, green and yellow color, respectively. Boxes indicate c1, c2, and c3 clusters. The c1L cluster region is marked with underlines. Light blue bold type indicates residues that are polymorphic in P. vivax of China-Myanmar border isolates. (TIF 2305 kb)
Patterns of nucleotide diversity and amino acid polymorphisms of Pvama1. Sliding window plot of nucleotide diversity (π) and amino acid polymorphism of Pvama1 ectodomain in 607 global isolates and 8 reference strains were shown. The π value was calculated using DnaSP v5.10.01 with window length of 90 bp and step size of 3 bp Domain I is highlighted in red color. Nucleotide and amino acid positions are after the Sal-I sequence. (TIF 183 kb)
Phylogenetic analysis of China-Myanmar border and Myanmar Pvama1 sequences. The tree was constructed by analyzing DI (nt, 322–737) with 73 Pvama1 sequences obtained from China-Myanmar border and 24 Pvama1 sequences obtained from Myanmar P. vivax isolates (GenBank accession nos. FJ157247, FJ157249, FJ157251, FJ157253–FJ157269 and FJ157282–FJ157285) using a neighbor-joining method. The bootstrap method with 1,000 replications was used to construct the gene tree. Red circle indicate sequences obtained from Myanmar. (TIF 1803 kb)
Genetic differentiation (F ST ) of the Pvama1 gene among Thai, Myanmar and China-Myanmar border populations across DI. (DOCX 14 kb)
The proportion of Pvama1 haplotypes variation observed in China-Myanmar border, Myanmar, and Thai populations. The size of the pies reflects the frequency of a particular haplotype. The lengths of the lines connecting the pies, measured from their centers, are in proportion to the number of base pair substitutions separating the haplotypes. Color of each pie represents different country. Abbreviation: H21, haplotype 21. (TIF 831 kb)
Location of current study detected SNPs in the predicted B-cell epitopes and IURs in PvAMA1. Predictions of B-cell epitopes and IURs were performed by using ABCpred and RONN server, respectively. The cysteine residues have been shown in red bold. DI = 94–247 aa; DII = 265–363 aa; DIII = 388–451 aa; SNPs in the current study were marked in grey shadow; underlined shows B-cell epitopes; Bold overbars show IURs. (TIF 323 kb)