Skip to main content

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

Abstract

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 (R2) 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.

Background

Plasmodium vivax is the most prevalent malaria species in Latin America, the Middle East, South and Southeast Asia, Oceania and The Horn of Africa [1]. In these regions, more than 2.5 billion people are at risk [1], and approximately 13.8 million cases were reported in 2014 [2]. During the 2000–2014 period in Latin America, malaria cases gradually diminished from 1.2 million to 390,000 cases, representing a 74% decline. Argentina, Paraguay and Costa Rica are in the elimination phase, while Belize, Dominican Republic, Ecuador, El Salvador and Mexico are in the pre-elimination phase [2].

In Mexico, the annual number of malaria cases fluctuated between 20,000 and 130,000 in the 1980s. Since the late 90’s, the improvement of strategies and efforts by the malaria national control program has driven a continuous decrease in malaria incidence, with only P. vivax causing cases of autochthonous transmission. The Pacific coast is the location of almost all transmission, mainly near the boundaries between Mexico and Guatemala as well as between the states of Chihuahua, Sinaloa and Sonora, as well as Durango and Nayarit. In the former case, along the Mexican side of the border is the State of Chiapas, and transmission occurs at the pacific side and in the Lacandon rainforest [3, 4].

Plasmodium vivax is haploid in the human host [5], which facilitates the haplotype assembly when analysing human blood stages. It has been found that this species is indeed genetically diverse worldwide [6, 7] and that the haplotype frequency can be modified by evolutionary forces such as mutation, natural selection, migration and genetic drift [8,9,10]. The understanding of how these forces induce P. vivax populations to differentiate and expand in time and space [11] are relevant for epidemiological surveillance and essential for the design and evaluation of control and elimination strategies. To study the genetic diversity of P. vivax, different molecular markers have been used, including microsatellites [7, 12,13,14,15], mitochondrial genes [16,17,18] and genes encoding antigenic molecules [19,20,21,22]. Many blood stage antigens especially those involved in immune evasion are highly polymorphic; they evolve more rapidly than neutral markers and allow high resolution of parasite haplotypes [23].

The apical membrane antigen 1 (ama1) gene is located on chromosome 9 and encodes an integral membrane protein type 1, considered as a strong vaccine candidate (AMA1). This antigen has an ectodomain that plays a key role in the reorientation of the merozoite during the invasion of red blood cells [24, 25]. The ectodomain comprises three domains: domain I, domain II and domain III [26, 27]. In domain DI of P. vivax, most studied, higher nonsynonymous than synonymous nucleotide changes have been detected in isolates from distinct regions such as Asia [28,29,30,31], Oceania [22], the Middle East [32] and Latin America [33].

The P. vivax population is genetically unique in southern Mexico. By using microsatellite markers, three subpopulations related to mosquito specificity and geographical distribution were identified [12]. Accordingly, the ookinete Pvs25-28 polymorphism was associated with vector compatibility [34]. The icb5-6 blocks of pvmsp1 revealed new hybrid lineages that diversified, possibly because of recombination between divergent haplotypes from South America and Asia [35]. In additional, the analysis of Pvama1 I-II in a limited number of samples from 2006–2007, exposed that this gene fragment had moderate genetic diversity and was found under balancing selection [36].

Population and evolutionary genetics provide knowledge about how selective forces are networking in any species. This knowledge is useful for evaluating the impact of control and elimination measures on the population structure of vector parasites. Hence, the aim of the current longitudinal study was to investigate changes in P. vivax diversity and population structure over 19 years, based on analysis of pvama1 I-II , in samples taken in southern Mexico. The study was from 1993 to 2011, a period comprising control and pre-elimination phases.

Methods

Plasmodium vivax samples

A total of 288 samples (105 of whole blood and 183 of blood smeared on filter paper) had been previously obtained from symptomatic patients, for about two-decade period since 1993 [37,38,39]. The P. vivax diagnosis was carried out at the laboratory facility of the Regional Center for Research in Public Health (CRISP-INSP). Patients lived in Jurisdiction VII of the State of Chiapas in Mexico, which comprises an area of 4644.07 km2 of tropical and template regions, with the altitude ranging from sea level to foothills up to 4,000 m above sea level. The number of P. vivax samples analyzed per year were as follows: 1993, n = 12; 1994, n = 25; 1995, n = 8; 1996, n = 4; 1997, n = 18; 1998, n = 18; 1999, n = 21; 2000, n = 17; 2001, n = 11; 2002, n = 19; 2003, n = 22; 2004, n = 20; 2005, n = 21; 2008, n = 22; 2009, n = 20; 2010, n = 17; and 2011, n = 13. DNA extraction was carried out with the QIAamp ® DNA Blood Mini Kit (Qiagen, Redwood City, CA, USA), following the manufacturer’s instructions.

Pvama1 I-II amplification and sequencing

A DNA fragment comprising domains DI and DII (pvama1 I-II ) was amplified. The reaction mixture consisted of 5 μl of 10× reaction buffer (600 mM Tris-SO4 at pH 8.9 and 180 mM (NH4)2SO4, 2 μl of MgSO4 (2 mM), 1 μl of dNTPs (0.2 mM), 1.8 μl of each primer at 36 pM (PvPvama1F 5′-TCC AGC TGG AAG ATG TCC TG-3 ′ and Pvama1R 5′-CCG CCC TTT TCT CTA CAC AG-3′), 0.2 μl enzyme Platinum® TaqDNA Polymerase High Fidelity Invitrogen™ (1 U per reaction), and 2–4 μl of DNA. The final volume was adjusted to 50 μl. The PCR conditions were as follows: the first cycle at 95 °C for 5 min, followed by 35 cycles at 95 °C for 60 s, one cycle at 61 °C for 60 s, another at 72 °C for 75 s, and a final cycle at 72 °C for 10 min. The reactions were performed in a Thermal Cycler 2720 thermocycler (Applied Biosystems®, CA, USA). The samples that produced little or no visible PCR product were re-amplified with a nested PCR by using primers Pvama1F and Pvama1R2 5′-CGC AGG GAC ATT TGA TAC TCT CC-3′ and 1–2 μl of PCR product from the first PCR reaction.

The PCR-amplified products were separated by electrophoresis with agarose gels at 1%. For visualisation, 0.2 μg/ml ethidium bromide was added, and they were observed in an ultraviolet light chamber. The PCR products of approximately 1,100 bp were purified using the MinElute® PCR Purification Kit (Qiagen, Redwood City, CA, USA), according to the manufacturer’s instructions. The purified DNA samples were sequenced on an ABI 3730xl DNA Analyzer (ThermoFisher, Waltham, MA, USA) in the High-Performance Center in Genomics (Seattle, WA, USA).

The electropherograms of 24 samples showed double peaks, suggesting mixed genotype infections. Therefore, DNA fragments obtained from the first PCR were cloned by using the TOPO TA Cloning Kit with pCR® 2.1-TOPO (Invitrogen, Carlsbad, CA, USA). For each isolate, three to five different clones were sequenced.

The nucleotide sequences were revised by using the BioEdit version 5.0.9 Sequencing Alignment Editor Copyright © program and aligned with Clustal W [40], utilising the sequence of strain Sal I as a reference. The consensus sequences were deposited in the GenBank database (Accession numbers: KY094724–KY094901).

Genetic analyses

The population genetic analysis was performed for all samples (n = 213) and for subgroups (defined by time periods). DnaSP v5.1 was used to determine the number of polymorphisms (S), haplotypes (h), and synonymous (dS) and non-synonymous (dN) changes, as well as the minimum number of recombination events (Rm), the haplotype diversity (Hd), genetic diversity (indexes π and θ), and linkage disequilibrium (LD) given by index R2. Moreover, to ascertain whether the observed diversity departs from neutrality, the difference between the substitution rates of nonsynonymous (dN) and synonymous (dS) nucleotide changes (dN - dS) was evaluated within populations by using the MEGA v6.0 program [41]. This estimation was based on the maximum joint likelihood reconstruction of ancestral states under the models of Muse-Gaut [42] and substitution of codons [43]. In addition, data from the neutrality test McDonald-Kreitman was calculated [44] with DnaSP v5.1. As the outgroup, a Plasmodium cynomolgi sequence was used (GenBank X66099.1) [45].

Plasmodium vivax haplotype networks were constructed to investigate their genetic relationships, which change over time. First, all the sequences from the parasites in southern Mexico were included to prepare the “master” haplotype network with the TCS 2.21 program. This network was the template to prepare the layers of the temporal network, one per period. The networks were piled up to produce a three-dimensional effect, thus facilitating a comparison of the relationships of the haplotypes and their frequency. The relationship between haplotypes was resolved through mutational steps, which allowed for the assignment of existing haplotypes to a common ancestor of the population [46].

STRUCTURE (v2.3.4) analysis was implemented to explore how local haplotypes were grouped into subpopulations, and if possible to visualise changes in the frequency of the defined subpopulations over time. This technique employs a Bayesian approach to establish the probability of the number of populations in a sample of sequences. The analysis was carried out for 50,000 iterations, followed by 100,000 Markov Chain Monte Carlo (MCMC), and all runs were based on the model of admixture [47]. Twenty replicates run with a K-value that ranged from two to six predicted the optimum value of K, and the probability LNP [D / K] was calculated [48]. The fixation index (F ST) was determined between parasite subgroups of southern Mexico (defined by time periods) and for parasites from distinct geographic origins using the program DnaSP v5.1.

Homologous sequences of P. vivax from different geographic sites were found at NCBI, GenBank: Venezuela (VNZ), EU346015–EU346087 [33]; Iran (IR) JF682785–JF682790 (Unpublished), JX624732–JX624760 [32]; Sri Lanka (SLK), EF218679–EF218701 [28], India (IND) EU282774–EU282822 [29], EF025187–EF025197 [49]; Papua New Guinea (PNG), KC702458, KC702402–KC702503 [22]; Thailand (THL) FJ784891–FJ785121 [30] and South Korea (SK) KM230319–KM230384 [31].

Results

Plasmodium vivax ama1 I-II from southern Mexico, 1993–2011

Polymorphism and MGI

One hundred and seventy-eight pvama1 I-II sequences comprising a fragment of 702 bp were obtained from 176 parasite isolates; eighty-four percent (n = 89) of whole blood samples (period 2002–2011), and only 47.5% (n = 87) of blood impregnated onto filter paper. The pvama1 I-II gene fragment was obtained from all samples from 1995 (n = 8) but none of the samples from 1996 (n = 4). For other years the proportion of samples with a high-quality sequence varied between 17 and 58%. Failure to amplify parasite DNA or to obtain a good quality sequence could be due to low parasitemia, and the quality of DNA preserved onto filter paper.

The electropherograms suggested only one haplotype (assigned as a single genotype infection) for 158 isolates, while 18 samples (10.33%) had more than one peak in at least one nucleotide position (suggesting more than one genotype per sample). The samples were re-sequenced and electropherograms reviewed by three examiners. The cloning resolved only for two samples, two haplotypes per each sample.

Thirty-five previously reported sequences (KP759815–KP759849) from 2006 and 2007 [36] were included in the analysis, which made 213 sequences.

Genetic diversity

The16 polymorphic sites were dimorphic. There were ten changes in the DI domain and six in the DII domain (Additional file 1: Table S1). The whole gene fragment (DI-DII) resolved fifteen haplotypes (H1-H15). Only one nucleotide change existed between H1 and H7, and one between H8 and H11. Haplotype H7 was also present in a single genotype infection. Seven polymorphisms were exclusive to southern Mexico, as previously described [36], and only two other nucleotide changes were detected in a single parasite from 2010 (T583A and A911G). Eight of 15 haplotypes were exclusive for southern Mexico (H1, H4, H7, H8, H9, H10, H14 and H15) (Additional file 2: Figure S1). At least one of the other haplotypes was previously reported in Venezuela [33], Iran [32], India [29, 49], Sri Lanka [28], Thailand [30] or South Korea [31] (see also Additional file 2: Figure S1).

For pvama1 I-II in southern Mexico, the overall nucleotide (π) and genetic diversity (θ) were low, being 0.0067 ± 0.0002 and 0.0038 ± 0.0009, respectively. These values are similar to those reported for a sample of 35 parasites from 2006–2007 [36], and above those calculated for parasites from Venezuela [33] or South Korea [31], both sites of low transmission [2]. Other geographic sites with high transmission had greater genetic diversity (0.0073–0.0101), including Iran [32], India [29, 49], Sri Lanka [28], PNG [22] and Thailand [30] (Additional file 3: Table S2). Parasites from southern Mexico had the lowest haplotype diversity (0.734), with their Rm value being similar to that found in parasites from Venezuela and South Korea (Additional file 3: Table S2).

Natural selection

The dN - dS value was positive (1.527; P > 0.05) for southern Mexico, as it was at the global level (0.806), in neither case reaching significance. There were similar values for other geographic sites, except SLK (2.154, P <0.05). However, the McDonald-Kreitman test was positive and highly significant at the global and local level for all geographic sites herein analysed (Table 1).

Table 1 Natural selection tests on pvama1 I-II among P. vivax populations of different geographic origins

Linkage disequilibrium

The LD coefficient R2 for southern Mexican parasites (0.5103) was moderate (Additional file 3: Table S2) compared to the elevated values at sites with the low transmission. Contrarily, R2 was low at sites showing high transmission.

F ST values

Plasmodium vivax from distinct geographical origins had from moderate to high values. Southern Mexican parasites had F ST values from 0.1267 to 0.4601 when compared to other regions. Between southern Mexico or Venezuela and Iran, moderate F ST values of 0.1267 and 0.1605 were obtained. Between southern Mexico and Venezuela, the value was higher (F ST = 0.2251, P < 0.001). The greatest F ST values were between South Korea and other sites (0.3654 to 0.6054). Moreover, F ST values between Iran, India and Sri Lanka were close to zero (Additional file 4: Table S3).

Temporal genetic analysis of P. vivax ama1 I-II

The number of sequences obtained per year varied from two to nineteen. Figure 1 shows the proportion of each haplotype found per year. There were important changes in haplotype frequency. For example, H8 was at a proportion of 0.667 in 2007, followed by its reduction to 0.158 in 2008. Contrarily, H1 increased to predominate in 2008 and 2009. In 2002, H5 was present at a higher frequency than H1 and H8. These sequences were grouped into six consecutive 3-year time periods, each including a minimum of 20 sequences: 1993–1995 (n = 21), 1997–1999 (n = 23), 2000–2002 (n = 22), 2003–2005 (n = 52), 2006–2008 (n = 54) and 2009–2011 (n = 41).

Fig. 1
figure 1

Yearly distribution of pvama1 I-II haplotypes in southern Mexico during 19 years. Of 213 DNA sequences, 15 haplotypes were resolved (denominated H1-H15). The proportion of each haplotype in each year is represented by the size of the oval (only those representing ≥ 0.10 are indicated). Some haplotypes were highly frequent in a particular period and existed in samples collected from different years. There were variations in haplotype frequency. It is possible that certain haplotypes alternated with one another. For instance, H1 was present for almost the entire study, having an increased frequency during 2008–2009 and in 2011. This haplotype seemed to alternate with H8, which had a greater proportion in 2007 and in 2010. No samples from 1996 were tested

The highest proportion of samples presumably carrying mixed genotype infections (MGI) was determined for 1993–1995 and 1997–1999, resulting in a value of 23.1 and 47.8%, respectively. The percentage of mixed genotype infections dropped to 4.5 and 2.4% for 2000–2002 and 2009–2011, respectively. In only two cases, the cloning resolved two haplotypes per sample: sample number 8998 included haplotypes H1 and H7, and sample number 2500 had H8 and H11. For another 12 isolates, only one haplotype per sample was resolved, while four other isolates were not possible to clone.

The genetic parameters varied over time. In the samples from 1993–1995, the nucleotide diversity was the lowest (π = 0.00526). Then in the 1996–1999 period there was an increase in this π value (0.00758) (Fig. 2a), which coincided with the significant rise in malaria cases from 1998 to 2000 (Fig. 2d), likely as the result of the climactic event “El Niño”. Subsequently, a reduction in both malaria cases and genetic diversity took place during 2000–2002, followed by a greater decline in malaria cases (although genetic diversity is remaining the same) during 2003–2005. Despite the rise in malaria cases during 2006–2008, the nucleotide diversity (π) continued to diminish through the 2009–2011 period (0.00591 ± 0.00007). The theta value was consistently much lower than pi. In contrast, the LD index R2 underwent a gradual increase from 0.4886 to 0.9614 (Fig. 2b). The dN - dS value was positive and gradually decreased from 1.971 to 1.313, but was not significant (Fig. 2c). Figure 2d shows the number of confirmed malaria cases per year, compiled by the Sanitary Jurisdiction VII in the State of Chiapas, Mexico. The number of cases was high during 1998–2000, and then a sudden drop was registered. A second peak was observed in the 2006–2009 period, followed by a sustained decline after that. In this Jurisdiction, only 101 malaria cases were reported in 2011.

Fig. 2
figure 2

Temporal genetic changes for P. vivax ama1 I-II in southern Mexico from 1993 to 2011. a Nucleotide diversity (π) and genetic diversity (θ). b Linkage disequilibrium given by R2. c Natural selection by dN-dS values. d Annual incidence of malaria according to confirmed cases in the 1993–2011 period within Jurisdiction VII of Chiapas

Haplotype networks and temporal changes in pvama1I-II allelic frequencies

The master network showed that the 15 haplotypes from southern Mexico were separated by one to twelve mutational steps, and two sampled haplotypes by one to four mutational steps. The most frequent haplotypes, H1 and H8 at 42.9 and 24%, respectively, were highly divergent, being separated by at least nine mutational steps. Haplotype H1 was two mutational steps from H5 (13.6%) and two to eleven mutational steps from other low-frequency haplotypes (mostly detected in a single isolate). On the other hand, H8 was four mutational steps from H2 (4.2%) and H4 (8.5%), and two mutational steps from H15, which was found in a single isolate. The network also revealed two loops, suggesting recombination events between haplotypes H2, H3, H8 and H11, as well as between haplotypes H3, H5, H6, H7 and H14 (Fig. 3a).

Fig. 3
figure 3

Temporal haplotype network of P. vivax ama1 I-II in southern Mexico. a The master network for the entire period (1993–2011). b Haplotype networks in a structure of six layers; the master network was used as a template and each layer represented a three-year period. Each color corresponds to a haplotype (H1-H15) and the size of the ovals is proportional to the number of isolates which share the same haplotype (except for H1 and H8; 0.5 magnification). Solid lines connect the haplotypes; black dots represent the number of mutational steps between haplotypes (including those sampled or extinct). Only the haplotypes detected in each time layer are shown in color

Seven of 15 pvama1 I-II haplotypes (46.7%) were present in more than one period. Another eight haplotypes were unique, and six of these were observed from 1998 to 2001, the period with the maximum number of reported malaria cases (Fig. 3b). There were five divergent haplotypes (H2, H4, H5, H7 and H10) from the beginning to the end of the study. Though H8 was abundant from 2002 on, it was not encountered in the samples during 1993–1995. Haplotype H5, closely related to H1 and H7 and detected in the 1993–1995 period, was not found immediately after the climactic event “El Niño” (1997–2001) but was very frequently observed from 2002 to 2008. H2 and H4, which were separated from H8 by five mutational steps and from each other by nine mutational steps, had similar distribution overtime. Based on the haplotype network for global parasites, with 104 varieties, southern Mexico has three lineages that are closely related to haplotypes from Venezuela, India, Iran and South Korea (Additional file 2: Figure S1).

Structure analysis

The 15 haplotypes from southern Mexico were grouped into three subpopulations (P1, P2 and P3) (Fig. 4a). P2 comprised ~ 57%, followed by P3 at ~ 24% and P1 at ~ 10.8%. P2, which included haplotype H1, was the most abundant and widely distributed subpopulation over time. P3 was present from 1997 to 2011. The frequency of P1 fluctuated greatly, and was the only one to undergo a drastic reduction. Its ~ 50% frequency in the 90’s fell to ~ 30% during 2003–2005 and then dropped sharply to ~ 4% in the 2006–2008 and 2009–2011 periods. In the 1990s, interestingly, P1 and P2 were the two most frequent populations (Fig. 4b). About 8% of individuals showed mixed ancestry. One mixed type consisted of P2 and P3 ancestry and was found in several, which were more frequent from 1993 to 1999 (Fig. 4c). A moderate structure was determined in Mexico compared to the global level (Additional file 5: Figure S2).

Fig. 4
figure 4

Population structure of P. vivax ama1 I-II in southern Mexico. The subpopulations are represented by different colors (K = 3): red (P1), green (P2) and blue (P3). A vertical bar represents each subject. a Resolution of the three subpopulations resulting from Bayesian cluster analysis. b Individuals are chronologically ordered, and periods are indicated. c The frequency of individuals with mixed ancestry is represented; one frequent type of admixed genotype (P2/P3*) is shown in black and other admixed types in grey

The highest F ST values were detected between 1993–1995 and 1997–1999, and between 1997–1999 and 2000–2011 (F ST = 0.167 and 0.153, respectively; P < 0.01 in both cases). The sample size was similar during these periods (Table 2).

Table 2 F ST values based on ama1 I-II between P. vivax populations defined by 3-year periods in southern Mexico

Discussion

Due to its significant reduction in malaria cases, Mexico entered the pre-elimination phase in 2009 [50]. One of the most affected areas is Southern Mexico, a region sharing a border with Central America. The complex malaria epidemiology of this region is illustrated in the current study. For 19 years (1993–2011), the increases and decreases in P. vivax cases were accompanied by changes in the parameters of population genetics. The gradual reduction in P. vivax cases after 1999–2000 went with a reduction in the genetic diversity and the rate of mixed genotype infections, as well as an increase in linkage disequilibrium, alterations in the haplotype and subpopulation frequency, and variations in natural selection.

The 213 P. vivax sequences detected over 19 years represented a much higher number of pvama1 I-II haplotypes than that found in each of the three-year periods. Additionally, the most frequent haplotypes were present either during all or most of the study. In the majority of cases, the genetic diversity varied over time and coincided with the fluctuations in malaria cases. The greatest number of haplotypes and the highest genetic diversity (π = 0.00758), observed during 1997–1999, but the period of pre-elimination (2009–2011) witnessed the second lowest diversity (π = 0.00591) in a subgroup, and this was based on a comparable or higher number of samples than those from previous years, which indicates that this reduced value is very likely to be accurate. These results suggest that the effective population size of P. vivax was relatively small in southern Mexico. A previous study using microsatellites also showed that changes in the annual number of cases were accompanied by modifications in the nucleotide diversity [51].

The current study was conducted with only a few samples per year, as the focus was on the trends rather than yearly changes. The rise in malaria cases during 2006–2009 seemed not to affect the nucleotide diversity, this being similar to the previous period (2003–2005). Moreover, no genetic difference between the parasites from these two periods was suggested by the F ST value, meaning that parasite structure was little affected.

The differentiation values can increase due to the divergence of the remaining parasites (as suggested by Volkman) [52], in regions entering the elimination phase [7], or in periods of re-emergence [13]. It was herein observed that the pre-elimination parasites were highly differentiated from parasites of uppermost transmission (1997–1999).

The detection of multiple genotype infections (MGIs) in regions with various levels of P. vivax transmission is common, but its rate rises in a high transmission setting [53]. In this study, there was a low-moderate proportion of MGIs and admixture (based on pvama1), as reported previously when using microsatellites in southern Mexico [12]. During the control phase (especially between 1993 and 1999), either MGIs or admixture rates were higher than in the elimination phase, the latter being a time when the number of MGIs and samples with admixture became negligible. This coincides with the extremely high linkage disequilibrium present in this phase.

Though the MGI rate was the greatest during the control phase, the presence of several haplotypes over time suggest that clonal transmission of P. vivax has been occurring in southern Mexico. The differential mosquito susceptibility to P. vivax strains and the geographic isolation of certain strains [12, 34] might have reduced the chances of recombination between divergent haplotypes in the region. If recombination occurred, it was between identical or very closely related genotypes, as reported by others [54]. The clonal transmission of human Plasmodium species occurs in both high [7, 20] and low transmission settings, the latter observed in Korea when using pvmsp1 and pvcsp [55] and in Brazil with P. falciparum when employing msp1 and msp2 [56, 57].

After malaria was re-introduced in Korea, the increase in genetic diversity was the result of parasite migrations rather than the generation of new hybrid parasites [13]. Although southern Mexico shares a border with Central America, and human migration flows from South/Central to North America, the introduction of new divergent parasites was not herein demonstrated. Like Central America, the circumsporozoite vk210 and vk247 genotypes have been detected since the early 1990s in southern Mexico [58]. Consequently, the haplotype networks suggest that P. vivax ama1 I-II haplotypes are closely related and have recently expanded due to local adaptations [12]. This idea is supported by previous findings with msp1 icb5-6 [35] and more recently with the genomic SNP analysis revealing that parasites from southern Mexico have recent common ancestry [59]. Such adapted, highly frequent and persistent haplotypes that are exclusive to southern Mexico might have spread as a consequence of climactic events that impeded access to communities, either by roads or rivers, thus creating impediments to the delivery of control measures.

Genes coding for diverse antigens (e.g. pvama1) are present in many low-frequency haplotypes that fluctuate by genetic drift, especially in populations facing a drastic demographic reduction. It cannot be ruled out that haplotypes coming from relapse cases may influence the observed changes in allele frequencies. Hypnozoite genetic reservoirs contribute to maintaining the malaria genetic pool over time [60]. Up to 1998 in southern Mexico, P. vivax was treated for five days with primaquine [61], which is not highly effective in preventing relapse episodes [62]. Subsequently, intermittent single doses administered to confirmed cases [60] did not avoid symptomatic and asymptomatic relapse episodes [33, 36]. In addition, gametocytes were harvested from most P. vivax infected samples (primary and relapse infections) [51, 63]. Another factor involved in the local development of P. vivax strains is the influence of climatic events, such as “El Niño” and Hurricane Stan (2005), which might have contributed to the resilience of malaria transmission despite control measures. The reports on the number of malaria cases suggest that it took about three and four years to get malaria back to a downward trend after each of these events.

Although there was a reduced haplotype pool in the pre-elimination phase, the most divergent haplotypes representing the three P. vivax lineages or subpopulations persisted. The yearly frequencies of P. vivax haplotypes underwent expansion and alternation, and the latter case was more notorious among more frequent and divergent haplotypes (e.g. H1 vs H8, H4 or H5). This information may be relevant for vaccine development.

There is evidence that balancing selection has maintained the polymorphism of pvama1 I-II [33, 36]. The haplotype alternation herein observed could have been induced by specific antibody responses [64]. It is interesting that dN - dS values decreased along with the decline in the number of malaria cases. Perhaps the reduction in immune individuals diminished selective pressure. Robust evidence exists that antibodies against PvAMA1 can inhibit the invasion of parasites into reticulocytes of human hosts [65, 66], and that antibody response to PvAMA1 are naturally acquired by human populations in distinct endemic regions. The influence of the antibody response has been reported in South America [67,68,69], the Mideast [69] and the Indian subcontinent [70, 71]. However, this does not explain the steady decrease in dN - dS values when malaria cases were most numerous, during 1997–1999 and 2006–2008.

Conclusion

In summary, anti-malarial actions have rendered a gradual reduction in the number malaria cases and caused alterations in the genetic structure of the parasite in the hypoendemic region of southern Mexico, despite the fluctuations in these parameters caused by climactic events and parasite and vector characteristics. Hence, parasite elimination from the region might be possible in the near future.

Abbreviations

MGI:

Mixed genotype infections

pvama1 I-II :

P. vivax apical membrane antigen 1 gene domains I and II

References

  1. Guerra CA, Howes RE, Patil AP, Gething PW, Van Boeckel TP, Temperley WH, et al. The international limits and population at risk of Plasmodium vivax transmission in 2009. PLoS Negl Trop Dis. 2010;4(8):e774.

    Article  PubMed  PubMed Central  Google Scholar 

  2. WHO. World Malaria Report 2015. Geneva: World Health Organization; 2015. http://www.who.int/malaria/publications/world-malaria-report-2015/en/. Accessed 16 April 2017.

    Google Scholar 

  3. Mexican Ministry of Health: Epidemiological bulletin. http://www.epidemiologia.salud.gob.mx/dgae/boletin/intd_boletin.html (2014). Accessed 16 April 2017

  4. PAHO. Interactive Malaria Statistics. Washington: Pan American Health Organization; 2014. http://www.paho.org/hq/index.php?option=com_content&view=article&id=2632%3A2010-interactive-malaria-statistics&catid=1617%3Amalaria-statistics-maps&Itemid=2130&lang=en. Accessed 16 Apr 2017.

    Google Scholar 

  5. Janse CJ, van der Klooster PF, van der Kaay HJ, van der Ploeg M, Overdulve JP. DNA synthesis in Plasmodium berghei during asexual and sexual development. Mol Biochem Parasitol. 1986;20(2):173–82.

    Article  CAS  PubMed  Google Scholar 

  6. Neafsey DE, Galinsky K, Jiang RH, Young L, Sykes SM, Saif S, et al. The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum. Nat Genet. 2012;44(9):1046–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Koepfli C, Rodrigues PT, Antao T, Orjuela-Sánchez P, Van den Eede P, Gamboa D, et al. Plasmodium vivax diversity and population structure across four continents. PLoS Negl Trop Dis. 2015;9(6):e0003872.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Chevillon C, de Meeûs T, McCoy KD. Population Genetics and Molecular Epidemiology of Infectious Diseases. In: Morand S, Beaudeau F, Cabaret J, editors. New Frontiers of Molecular Epidemiology of Infectious Diseases. Dordrecht: Springer Netherlands; 2012. p. 45–76.

  9. Meirmans PG. The trouble with isolation by distance. Mol Ecol. 2012;21(12):2839–46.

    Article  PubMed  Google Scholar 

  10. Travis J. The interplay of population dynamics and the evolutionary process. Philos Trans R Soc Lond B Biol Sci. 1990;330(1257):253–9.

    Article  CAS  PubMed  Google Scholar 

  11. Carlton JM, Das A, Escalante AA. Genomics, population genetics and evolutionary history of Plasmodium vivax. Adv Parasitol. 2013;81:203–22.

    Article  PubMed  Google Scholar 

  12. Joy DA, Gonzalez-Ceron L, Carlton JM, Gueyes A, Fay M, McCutchan TF, et al. Local adaptation and vector-mediated population structure in Plasmodium vivax malaria. Mol Biol Evol. 2008;25(6):1245–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Iwagami M, Hwang S-Y, Kim S-H, Park S-J, Lee G-Y, Matsumoto-Takahashi ELA, et al. Microsatellite DNA analysis revealed a drastic genetic change of Plasmodium vivax population in the Republic of Korea during 2002 and 2003. PLoS Negl Trop Dis. 2013;7(10):e2522.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Batista CL, Barbosa S, Da Silva BM, Viana SA, Ferreira MU. Genetic diversity of Plasmodium vivax over time and space: a community-based study in rural Amazonia. Parasitology. 2014;142(2):374–84.

    Article  PubMed  Google Scholar 

  15. Delgado-Ratto C, Soto-Calle VE, Van den Eede P, Gamboa D, Rosas A, Abatih EN, et al. Population structure and spatio-temporal transmission dynamics of Plasmodium vivax after radical cure treatment in a rural village of the Peruvian Amazon. Malar J. 2014;13:8.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Culleton R, Coban C, Zeyrek FY, Cravo P, Kaneko A, Randrianarivelojosia M, et al. The origins of African Plasmodium vivax; insights from mitochondrial genome sequencing. PLoS One. 2011;6(12):e29137.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Miao M, Yang Z, Patch H, Huang Y, Escalante AA, Cui L. Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion. BMC Evol Biol. 2012;12:22.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Taylor JE, Pacheco MA, Bacon DJ, Beg MA, Machado RL, Fairhurst RM, et al. The evolutionary history of Plasmodium vivax as inferred from mitochondrial genomes: Parasite genetic diversity in the Americas. Mol Biol Evol. 2013;30(9):2050–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Pacheco MA, Poe AC, Collins WE, Lal AA, Tanabe K, Kariuki SK, et al. A comparative study of the genetic diversity of the 42 kDa fragment of the merozoite surface protein 1 in Plasmodium falciparum and P. vivax. Infect Genet Evol. 2007;7(2):180–7.

    Article  CAS  PubMed  Google Scholar 

  20. Chenet SM, Tapia LL, Escalante AA, Durand S, Lucas C, Bacon DJ. Genetic diversity and population structure of genes encoding vaccine candidate antigens of Plasmodium vivax. Malar J. 2012;11:68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ju H, Kang JM, Moon SU, Kim JY, Lee HW, Lin K, et al. Genetic polymorphism and natural selection of Duffy binding protein of Plasmodium vivax Myanmar isolates. Malar J. 2012;11:60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. 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(10):e2506.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Barry AE, Waltmann A, Koepfli C, Barnadas C, Mueller I. Uncovering the transmission dynamics of Plasmodium vivax using population genetics. Pathogens and Global Health. 2015;109(3):142–52.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Cowman AF, Berry D, Baum J. The cellular and molecular basis for malaria parasite invasion of the human red blood cell. J Cell Biol. 2012;198(6):961–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Riglar DT, Richard D, Wilson DW, Boyle MJ, Dekiwadia C, Turnbull L, et al. Super-resolution dissection of coordinated events during malaria parasite invasion of the human erythrocyte. Cell Host Microbe. 2011;9(1):9–20.

    Article  CAS  PubMed  Google Scholar 

  26. Hodder AN, Crewther PE, Matthew ML, Reid GE, Moritz RL, Simpson RJ, et al. The disulfide bond structure of Plasmodium apical membrane antigen-1. J Biol Chem. 1996;271(46):29446–52.

    Article  CAS  PubMed  Google Scholar 

  27. Pizarro JC, Normand BV-L, Chesne-Seck M-L, Collins CR, Withers-Martinez C, Hackett F, et al. Crystal structure of the malaria vaccine candidate apical membrane antigen 1. Science. 2005;308(5720):408–11.

  28. 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 (PvAMA1) locus in a Sri Lankan population. Mol Biol Evol. 2007;24(4):939–47.

    Article  CAS  PubMed  Google Scholar 

  29. Thakur A, Alam MT, Bora H, Kaur P, Sharma YD. Plasmodium vivax: Sequence polymorphism and effect of natural selection at apical membrane antigen 1 (PvAMA1) among Indian population. Gene. 2008;419(1–2):35–42.

    Article  CAS  PubMed  Google Scholar 

  30. 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(6):1295–300.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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(1):455.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Zakeri S, Sadeghi H, Mehrizi AA, Djadid ND. Population genetic structure and polymorphism analysis of gene encoding apical membrane antigen-1 (AMA1) of Iranian Plasmodium vivax wild isolates. Acta Trop. 2013;126(3):269–79.

    Article  CAS  PubMed  Google Scholar 

  33. 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(10):e3366.

    Article  PubMed  PubMed Central  Google Scholar 

  34. González-Cerón L, Alvarado-Delgado A, Martínez-Barnetche J, Rodríguez MH, Ovilla-Muñoz M, Pérez F, et al. Sequence variation of ookinete surface proteins Pvs25 and Pvs28 of Plasmodium vivax isolates from southern Mexico and their association to local anophelines infectivity. Infect Genet Evol. 2010;10(5):645–54.

    Article  PubMed  Google Scholar 

  35. Cerritos R, González-Cerón L, Nettel JA, Wegier A. Genetic structure of Plasmodium vivax using the merozoite surface protein 1 icb5-6 fragment reveals new hybrid haplotypes in southern Mexico. Malar J. 2014;13:35.

    Article  PubMed  PubMed Central  Google Scholar 

  36. González-Cerón L, Cerritos R, Corzo-Mancilla J, Santillán F. Diversity and evolutionary genetics of the three major Plasmodium vivax merozoite genes participating in reticulocyte invasion in southern Mexico. Parasit Vectors. 2015;8(1):651.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Gonzalez-Ceron L, Rodriguez MH, Wirtz RA, Sina BJ, Palomeque OL, Nettel JA, et al. Plasmodium vivax: a monoclonal antibody recognized a circumsporozoite protein precursor on the sporozoite surface. Exp Parasitol. 1998;90(3):203–11.

    Article  CAS  PubMed  Google Scholar 

  38. Rodriguez MH, Gonzalez-Ceron L, Hernandez JE, Nettel JA, Villarreal C, Kain KC, et al. Different prevalences of Plasmodium vivax phenotypes VK210 and VK247 associated with the distribution of Anopheles albimanus and Anopheles pseudopunctipennis in Mexico. Am J Trop Med Hyg. 2000;62(1):122–7.

    CAS  PubMed  Google Scholar 

  39. González-Cerón L, Martinez-Barnetche J, Montero-Solís C, Santillán F, Soto AM, Rodríguez MH, et al. Molecular epidemiology of Plasmodium vivax in Latin America: polymorphism and evolutionary relationships of the circumsporozoite gene. Malar J. 2013;12:243.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Muse SV, Gaut BS. A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. Mol Biol Evol. 1994;11(5):715–24.

    CAS  PubMed  Google Scholar 

  43. Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981;17(6):368–76.

    Article  CAS  PubMed  Google Scholar 

  44. McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991;351(6328):652–4.

  45. Dutta S, Malhotra P, Chauhan VS. Sequence analysis of apical membrane antigen 1 (AMA-l) of Plasmodium cynomolgi bastianelli. Mol Biochem Parasitol. 1995;73(1–2):267–70.

    Article  CAS  PubMed  Google Scholar 

  46. Templeton AR, Sing CF. A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. IV. Nested analyses with cladogram uncertainty and recombination. Genetics. 1993;134(2):659–69.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Earl DA, VonHoldt BM. Structure Harvester: a website and program for visualizing Structure output and implementing the Evanno method - Technical Note. Conserv Genet Resour. 2012;4(2):359–61.

    Article  Google Scholar 

  49. 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(3):252–6.

    Article  CAS  PubMed  Google Scholar 

  50. WHO. World Malaria Report 2016. Geneva: World Health Organization; 2016. http://www.who.int/malaria/publications/world-malaria-report-2016/report/en/. Accessed 16 April 2017.

    Google Scholar 

  51. Gonzalez-Ceron L, Mu J, Santillán F, Joy D, Sandoval MA, Camas G, et al. Molecular and epidemiological characterization of Plasmodium vivax recurrent infections in southern Mexico. Parasit Vectors. 2013;1(6):109.

    Article  Google Scholar 

  52. Volkman SK, Neafsey DE, Schaffner SF, Park DJ, Wirth DF. Harnessing genomics and genome biology to understand malaria biology. Nat Rev Genet. 2012;13(5):315–28.

    Article  CAS  PubMed  Google Scholar 

  53. Havryliuk T, Ferreira MU. A closer look at multiple-clone Plasmodium vivax infections: detection methods, prevalence and consequences. Mem Inst Oswaldo Cruz. 2009;104(1):67–73.

    Article  PubMed  Google Scholar 

  54. Escalante AA, Cornejo OE, Rojas A, Udhayakumar V, Lal AA. Assessing the effect of natural selection in malaria parasites. Trends Parasitol. 2004;20(8):388–95.

    Article  PubMed  Google Scholar 

  55. Choi YK, Choi KM, Park MH, Lee EG, Kim YJ, Lee BC, et al. Rapid dissemination of newly introduced Plasmodium vivax genotypes in South Korea. Am J Trop Med Hyg. 2010;82(3):426–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Ferreira MU, Ribeiro WL, Tonon AP, Kawamoto F, Rich SM. Sequence diversity and evolution of the malaria vaccine candidate merozoite surface protein-1 (MSP-1) of Plasmodium falciparum. Gene. 2003;304:65–75.

    Article  CAS  PubMed  Google Scholar 

  57. Tonon AP, Hoffmann EH, Silveira LA, Ribeiro AG, Gonçalves CR, Ribolla PE, et al. Plasmodium falciparum: sequence diversity and antibody recognition of the merozoite surface protein-2 (MSP-2) in Brazilian Amazonia. Exp Parasitol. 2004;108(3–4):114–25.

    Article  CAS  PubMed  Google Scholar 

  58. Kain KC, Brown AE, Kyle H, Wirtz RA, Jay S, Rodriguez MH, et al. Circumsporozoite genotyping of global isolates of Plasmodium vivax from dried blood specimens. J Clin Microbiol. 1992;30:1863–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Hupalo DN, Luo Z, Melnikov A, Sutton PL, Rogov P, Escalante A, Vallejo AF, et al. Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax. Nat Genet. 2016;48(8):953–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Imwong M, Snounou G, Pukrittayakamee S, Tanomsing N, Kim JR, Nandy A, et al. Relapses of Plasmodium vivax infection usually result from activation of heterologous hypnozoites. J Infect Dis. 2007;195(7):927–33.

    Article  CAS  PubMed  Google Scholar 

  61. Mexican Ministry of Health: Norma Oficial Mexicana NOM-032-SSA2-2002. http://www.salud.gob.mx/unidades/cdi/nom/032ssa202.html (2002). Accessed 16 April 2017

  62. Gogtay NJ, Desai S, Kamtekar KD, Kadam VS, Dalvi SS, Kshirsagar NA. Efficacies of 5- and 14-days primaquine regimens in the prevention of relapses in Plasmodium vivax infections. Ann Trop Med Parasitol. 1999;93(8):809–12.

    Article  CAS  PubMed  Google Scholar 

  63. Gonzalez-Ceron L, Rodriguez MH, Sandoval MA, Santillan F, Galindo-Virgen S, Betanzos AF, et al. Effectiveness of combined chloroquine and primaquine treatment in 14 days versus intermittent single dose regimen, in an open, non-randomized, clinical trial, to eliminate Plasmodium vivax in southern Mexico. Malar J. 2015;14:426.

  64. Bueno LL, Lobo FP, Morais CG, Mourão LC, de Ávila RA, Soares IS, et al. Identification of a highly antigenic linear B cell epitope within Plasmodium vivax apical membrane antigen 1 (AMA-1). PLoS One. 2011;6(6):e21289.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Vicentin EC, Françoso KS, Rocha MV, Iourtov D, dos Santos FL, Kubrusly FS, et al. Invasion-inhibitory antibodies elicited by immunization with Plasmodium vivax apical membrane antigen-1 expressed in Pichia pastoris yeast. Infect Immun. 2014;82(3):1296–307.

  66. Arévalo-Herrera M, Chitnis C, Herrera S. Current status of Plasmodium vivax vaccine. Hum Vaccin. 2010;6(1):124–32.

  67. Múfalo 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(12–13):1266–73.

    Article  PubMed  Google Scholar 

  68. Rodrigues MH, Rodrigues KM, Oliveira TR, Cômodo 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(2):185–92.

    Article  CAS  PubMed  Google Scholar 

  69. Salavatifar M, Zakeri S, Abouie Mehrizi A, Mirkhazemi S, Dinparast DN. Evaluation of naturally acquired antibody responses to two variant forms of Plasmodium vivax apical membrane antigen-1 in individuals living in areas of low and unstable malaria transmission of Iran. Arch Iran Med. 2015;18(12):834–43.

    PubMed  Google Scholar 

  70. Seth RK, Bhat AA, Rao DN, Biswas S. Acquired immune response to defined Plasmodium vivax antigens in individuals residing in northern India. Microbes Infect. 2010;12(3):199–206.

    Article  CAS  PubMed  Google Scholar 

  71. Wickramarachchi T, Premaratne PH, Perera KL, Bandara S, Kocken CH, Thomas AW, et al. Natural human antibody responses to Plasmodium vivax apical membrane antigen 1 under low transmission and unstable malaria conditions in Sri Lanka. Infect Immun. 2006;74(1):798–801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Alejandro Flores-Alanis is a doctoral student from Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM) and received fellowship 347268 from CONACYT. Thanks to Olga L. Palomeque at the Regional Center for Research in Public Health in Mexico for her technical assistance.

Funding

This work was funded by CONACyT-Mexico project CB-2009-01-131247.

Availability of data and material

The datasets supporting the conclusions of this article are included within the article and its additional files. The DNA sequences and the data related were deposited in the GenBank database with accession numbers KY094724–KY094901.

Authors’ contributions

AFA, run molecular techniques, analysis of data and its interpretation, drafted the manuscript. LGC and RC, participated in study design, data analysis and interpretation, and were major contributors in writing the manuscript. CX, drafted the manuscript. FSV and MAS, performed molecular techniques. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was approved by the Ethical Committee of the National Institute of Public Health in Mexico (INSP). P. vivax blood samples had been obtained from patients living in the southern Chiapas, Mexico.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Lilia González-Cerón or René Cerritos.

Additional files

Additional file 1: Table S1.

Polymorphism and haplotypes defined for P. vivax ama1 I-II in southern Mexican parasites. Sixteen variable sites were detected, and 15 haplotypes (H1-H15) were resolved. (DOCX 21 kb)

Additional file 2: Figure S1.

Global haplotype network of pvama1I-II. Each circle corresponds to a haplotype (104 haplotypes) and each color represent the parasite origin, the size of the circle is proportional to the number of isolates that share the same haplotype. Single haplotypes were not considered. Solid lines connect the haplotypes, while short and perpendicular lines on the solid lines represent the number of mutational steps between two haplotypes. Small black circles indicate haplotypes not sampled or extinct. (PDF 144 kb)

Additional file 3: Table S2.

Parameters of genetic diversity and recombination for pvama1 I-II in parasite populations of different geographic origin. (DOCX 24 kb)

Additional file 4: Table S3.

F ST values between P. vivax populations of different geographic origin. (DOCX 23 kb)

Additional file 5: Figure S2.

Population structure of P. vivax based on pvama1I-II. Colors (red, blue and green) represent the resulting populations (K = 3) from Bayesian clustering analysis. A vertical bar represents each individual. (PDF 197 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Flores-Alanis, A., González-Cerón, L., Santillán, F. et al. Temporal genetic changes in Plasmodium vivax apical membrane antigen 1 over 19 years of transmission in southern Mexico. Parasites Vectors 10, 217 (2017). https://doi.org/10.1186/s13071-017-2156-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-017-2156-y

Keywords