Skip to main content

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

Abstract

Background

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

Methods

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

Results

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

Conclusions

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

Background

Malaria has been a major public health problem worldwide. In 2021, an estimated 247 million cases of malaria occurred, and 619,000 deaths were reported globally. In the Americas, malaria is primarily distributed in the Southern hemisphere and Plasmodium vivax is the predominant parasite, representing 64% of all malaria cases [1]. Ecuador has effectively managed the spread of malaria in recent years. The diagnosed prevalence of the disease has decreased by approximately 99% in the last 15 years. Transmission of malaria is mostly periodic in the coastal region and perennial in the Amazon, mainly rural and affects mostly men in working age [2]. Periodic outbreaks occur in endemic regions and since 2015 the number of reported malaria cases increased each year, reaching 2,190 in 2021, out of which 88% corresponded to P. vivax [2]. However, it is likely that the real number of infected individuals is underestimated since there are continuous undocumented reports of vivax malaria, and several studies have shown that asymptomatic and submicroscopic infections are common in Ecuadorian endemic communities [3, 4].

Plasmodium vivax merozoites preferentially infect human reticulocytes in a dynamic multiple-step process progressing rapidly from initial attachment to apical orientation and junction formation, then invasion and membrane sealing [5]. Merozoite surface protein 1 (MSP1) is the most abundant protein on the parasite surface and is implicated in the invasive merozoite’s initial attachment [6]. The mature form of MSP1 has undergone proteolytic processing and exists on the merozoite surface as a complex of four peptide fragments of 83, 30, 38 and 42 kDa [7,8,9]. At the time of invasion, the 42 kDa fragment is further processed in two more fragments of 33 and 19 kDa. The latter, MSP119 (coded by pvmsp-1-19) remains on the merozoite surface and is carried into the invaded reticulocyte [10, 11]. The EGF-like domains and disulfide binding pattern of MSP119 are highly conserved and several studies demonstrated that it is highly immunogenic [12, 13]. Additionally, it has been shown that antibodies against MSP119 can significantly inhibit parasite growth in vitro and consequently anti-MSP1 immunity is associated with immune protection against malaria [14,15,16].

In the next steps of invasion, merozoites orient their apical end towards the reticulocyte surface membrane that leads to junction formation. This process involves proteins released from the organelles of the apical complex, including Duffy binding protein (DBP) and Reticulocyte binding protein 1a (RBP1a). These proteins are synthesized during the final stages of asexual development and then sequestered in the micronemes until transported to the surface of the parasite in the moments of invasion [17, 18]. Region II (DBPII) contains the receptor-binding motifs necessary for adherence to its cognate receptor, the Duffy antigen receptor for chemokines (DARC) [19,20,21]. Although less immunogenic than MSP1, antibodies directed against PvDBPII reduce binding efficiency to DARC and inhibit invasion of P. vivax merozoites into reticulocytes [22]. RBP1a has been suggested to play an important role in binding and re-orientation of merozoites when invading reticulocytes [23]. Although PvRBP1a has been well characterized in the last few years, the specific binding site is still under discussion and the binding motif appears to be between residues 361 and 599, corresponding to pvrbp1a-2 gene fragment [23, 26].

In junction formation, merozoites form an irreversible attachment as a tight junction with the reticulocyte, which is mediated by another microneme protein, Apical Membrane Antigen 1 (AMA1). Like DBP, AMA1 is transported to the merozoite surface after trafficking through the micronemes where it is processed prior to being transported to the parasite surface during invasion [27]. Once this antigen is on the surface, it can bind to RON2 (Rhoptry Neck Protein 2), which is previously synthesized in the rhoptries of the same merozoite, and then inserted into the surface of reticulocytes to fulfill the self-anchoring function [28, 29]. The extracellular region of AMA1 consists of three structural domains defined by a highly conserved pattern of disulfide bounds. Analysis of the three-dimensional structure of this ligand, revealed that domains I and II belong to the PAN module superfamily, which is involved in binding to receptors in reticulocytes [30]. In addition, several studies suggest that the immunogenic N-terminal extracellular domains I and II that bind RON2 are the targets of inhibitory antibodies, making it an important vaccine candidate against P. vivax [31].

In the final phase of invasion, several factors stimulate the actin-myosin motor to power the junction to move from the anterior to the posterior end of the merozoite, bringing the parasite into the reticulocyte while a parasitophorous vacuole is formed around it. Finally, the membrane seals to enclose the parasite inside a vacuole in the reticulocyte [32].

The knowledge about strain variation of these critical merozoite invasion ligands is important since these are targets of protective immunity and potential vaccine candidates. There are some insights about the genetic diversity of dbpII in P. vivax populations from Colombia [33] and Brazil [34], but the data about allelic variation of these ligands is generally limited in the Americas and the general levels of geographic diversity remain poorly understood. In Venezuela a comparative study of ama1 diversity among Plasmodium falciparum and P. vivax was carried out with Amazonian isolates [35]. Also, a recent analysis of ama1 sequences in Brazilian P. vivax isolates determined the implications of haplotype diversity in the immune response [36]. However, only one study, carried out in parasites from southern Mexico, analysed the main genes involved in the invasive process together (pvmsp1-42, pvdbpII and pvama1). The findings revealed some differences in the genetic diversity may occur by region and possible evolutionary forces may influence each gene, depending on its role in reticulocyte invasion [37].

To address this important knowledge gap of these vaccine candidates, studies are needed to provide essential information about the diversity of these genes, their distribution in endemic countries, as well as the influence of evolutionary forces in the generation and maintenance of the genetic variability. Though a relatively small country, the geography of Ecuador is ideally suited to analyse parasite allelic diversity of geographically distinct regions, since the malaria endemic areas are separated by the Andes Mountains, where there are no mosquito vectors for transmission, into a western Pacific region and an Amazonian region.

The aim of this study was to determine the allelic diversity of four P. vivax merozoite invasion proteins considered potential candidate vaccines (pvmsp1-19, pvdbpII, pvrbp1a-2 and pvama1), and how natural selection influences this diversity in Ecuador. These findings will help contribute to the feasibility of use of these genes as potential targets in future vaccine development for the region and the country.

Methods

The study protocol was approved by the Ethical Review Committee of Pontificia Universidad Católica del Ecuador Approval number (#CBE-016-2013). Written informed consent (IC) was provided by study participants and/or their legal guardians for sample collection.

Plasmodium vivax blood samples

A total of 87 samples from 2012 to 2015 were analysed from individuals diagnosed with P. vivax infection by microscopy and PCR [38]. Samples were grouped by their geographic origin: Coast (n = 27); Western Amazon (n = 11); and Eastern Amazon (n = 49) (Fig. 1). Total genomic DNA (gDNA) was extracted from 200 ul of infected blood samples using QIAamp DNA minikit (Qiagen, USA) and resuspended in 60 ul, following the manufacturer’s instructions.

Fig. 1
figure 1

Geographical distribution of the localities assigned for the study. A total of 87 P. vivax samples collected between 2012–2015 were used. Localities were assigned depending on the geographical position. Eastern Amazon (green), western Amazon (orange) and coast (purple)

PCR amplification and DNA sequencing

A pvmsp1-19 gene fragment of approximately 326 bp was amplified by PCR using primers: pvmsp1-19F 5’-ATGAGCTCCGAGCACACATGTATAG-3’ and pvmsp1-19R 5’-ATGCACAGGAGGAAAAGCAACATG-3’. The final volume of 37.5 µL contained: 1X GoTaq Colorless Mastermix (Promega, Madison, USA), 25 mM of MgCl2, 10 µM of each primer and 3 µL DNA. The PCR cycle was 94 ºC for 2 min, followed by 35 cycles of 94 ºC for 30 s, 62 ºC for 30 s and 72 ºC for 90 s, and a final extension at 72 ºC for 5 min.

A pvdbpII gene fragment of approximately 1200 bp was amplified using the primers: pvdbpIIF 5’-GATAAAACTGGGGAGGAAAAAGAT-3’ and pvdbpIIR 5’-CTTATCGGATTTGAATTGGTGGC-3’ [37]. The PCR was prepared as above, and the PCR cycle was at 94 ºC for 2 min, followed by 35 cycles of 94 ºC for 30 s, 62 ºC for 30 s and 72 ºC for 90 s, and a final extension at 72 ºC for 5 min.

A pvrbp1a-2 gene fragment of approximately 1064 bp was amplified using the primers: pvrbp1a-2F 5’-AATACAATGCGCAAGAATTTTATATC-3’ and pvrbp1a-2R 5’-GTGAAGTACACACTTGATTTCCC-3’. The PCR was prepared as above, and the PCR cycle was at 94 ºC for 2 min, followed by 35 cycles of 94 ºC for 30 s, 54.7 ºC for 30 s and 72 ºC for 90 s, and a final extension at 72 ºC for 5 min.

The pvama1 CDS of approximately 1600 bp was amplified by a nested PCR using the primers: nest 1 pvama1_F16 5’-GCGGTTACTTCCACCCC-3’ and pvama1_R613 5’-GCGTGGTGTGGGAGGCCC-3’, and nest 2 pvama1_F29 5’- GCAAACCAAATCGCTGCC-3’ and pvama1_R598 5’- GCCTCGGGGTCGAGCATCTCGTC-3’. The PCR for nest 1 was prepared as 1X of GoTaq Colorless Mastermix (Promega, Madison, USA), 25 mM of MgCl2, 10 µM of each primer and 2 µL of genomic DNA; for a final PCR volume of 25 µL. The PCR cycle was 94 ºC for 3 min, followed by 40 cycles of 94 ºC for 1 min, 46 ºC for 1 min and 72 ºC for 2.5 min, and a final extension at 72 ºC for 10 min [35].

The PCR for nest 2 of pvama1 was prepared as 1X of GoTaq Colorless Mastermix (Promega, Madison, USA), 25 mM of MgCl2, 10 µM of each primer, and 3 µL of the previous PCR product; for a total PCR reaction volume of 37.5 µL. The PCR cycle was 94 ºC for 3 min, followed by 40 cycles of 94 ºC for 1 min, 57 ºC for 1 min and 72 ºC for 2.5 min, and a final extension at 72 ºC for 10 min [35].

The amplified products were purified using Illustra™ ExoProStar™ Enzymatic PCR and Sequencing Clean-up kit (GE Healthcare Life Sciences, Buckinghamshire, UK). A solution was prepared with of 2.5 µL of exonuclease, 25 µL of alkaline phosphatase and 472.5 µL of H2O, for a final volume of 500 µL. Afterwards, 5 µL of this solution was poured in each sample. Then, they were incubated at 37 ºC for 25 min and, subsequently, at 80 ºC for 20 min.

The purified products of pvmsp1-19, pvdbpII and pvrbp1a were sequenced using their own forward and reverse primers. For pvama1 the following primers:

pvama1_F142 5’-AGAATTCCAGCTGGAAGATG-3’, pvama1_F301 5’-CGTAAAAATTTAGGAAACGCC-3’, pvama1_F441 5’-CCCCTGCAGCATATATAAAGAC-3’, pvama1_R138 5’-TTCCACTTCTGCATCTTCCCC-3’, pvama1_R293 5’-ACATTTTTGCTCAAATACACC-3’ and pvama1_R442 5’-GGGAAATCCCGGTCTACTTC-3’ were used to encompass the whole sequence [35]. The sequencing was performed in a 4-capillary electrophoresis genetic analyzer (Applied Biosystems model 3130xl) at Macrogen (Seoul, KOR).

The quality of pherograms was verified manually and incomplete and low-quality sequences were not further analysed. Geneious v10.1.3 (Biomatters Limited, NZ) was used to assemble, align, and edit the raw sequences to obtain consensus sequences for each gene fragment.

Data analysis

The consensus sequences were aligned to the corresponding genes of Sal1, as follows: for pvmsp1-19, PVX_099980 (PlasmoDB); for pvdbpII, PVX_110810 (PlasmoDB); for pvrbp1a-2, AY501884 (GenBank) and for pvama1, PVX_092275 (PlasmoDB). DNA polymorphisms were analysed for the four gene fragments. The number of segregating sites (S), singletons (Si), parsimony sites (Pa), the average number of nucleotide differences (k), number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π), and their corresponding standard deviations (SD) were calculated with DnaSP v5.1 software [39].

The genetic relationship between haplotypes, of pvdbpII, pvrbp1a-2 and pvama1 and how they were distributed in the locations, were constructed using Median Joining Network algorithm in NETWORK v5.0.0.1 software.

Tajima's D and Fu & Lì’s D* and F* were calculated to test the neutral theory of evolution. Positive or negative significant values of both tests reject the neutral theory of evolution, suggesting the influence of natural selection [38]. The McDonald-Kreitman test was used to analyse synonymous and nonsynonymous variation. Under neutrality, the ratio of replacement to synonymous fixed substitutions between groups should be the same as the ratio of replacement to synonymous polymorphisms within groups [39].

Wright’s fixation statistics (FST) indexes were calculated to determine the degree of differentiation between the P. vivax from Ecuador and from other geographic origins, based on pvdbpII, pvrbp1a-2 and pvama1 invasion genes. FST values range from 0 to 1, high values indicate considerable degree of differentiation between parasite populations. The respective analyses were completed using the Arlequin v3.5.2.2 software.

A total of 638 sequences for pvdbpII from other geographic sites were obtained from the NCBI GenBank. For Brazil (BRZ), n = 122: EU812839-EU812960 [34]. For Colombia (COL), n = 17: U50575-U50591 [33]. For India (IND), n = 100: FJ491142-FJ49124 (Prajapati and Joshi, unpublished). For Iran (IRN), n = 19: EU860428-EU860438 [40], KF318358-KF318359 [41], and KF791921-KF791926 [42]. For South Korea (SK), n = 15: JN989472-JN989484 [43], and AF215737-AF215738 [44]. For Mexico (MEX), n = 35: KP759780-KP759814 [37]. For Papua Nueva Guinea (PNG), n = 200: AF289480-AF291096 [45], AY970837-AY970925 [46], and AF469515-AF469602 [47]. For Sri Lanka (SLK), n = 100: GU143914-GU144013 [48]. For Thailand (THL), n = 30: EF219451-EF379135 [49].

A total of 437 sequences were obtained for pvama1. For India (IND), n = 8: EF025187-EF025197 [50]. For Papua Nueva Guinea (PNG), n = 102: KC702402-KC702503 [51]. For Sri Lanka (SLK), n = 23: EF218679-EF218701 [52]. For Thailand (THL), n = 231: FJ784891-FJ785121 [53]. For Venezuela (VNZ), n = 73: EU346015-EU346087 [35].

Polymorphic residues were mapped on three-dimensional structural models for PvDBPII (4NUV) and PvAMA1 (1W81) using Discovery Studio Visualizer v16.1.0. B-cell epitope predictions for PvAMA1 were based on P. vivax Sal1 (PVX_092275) using BepiPred-2.0: Sequential B-Cell Epitope Predictor [54]. An epitope threshold of 0.5 with a length of at least 5 amino acids long was applied, as suggested in a previous study [55].

Results

Genetic polymorphism of  P. vivax  merozoite invasion genes

This study analysed P. vivax samples from two geographically and ecologically separate malaria endemic regions of Ecuador on both sides of the Andes Mountains. A pvmsp1-19 gene fragment of 294 bp (codons 1646–1743) was analysed from a total of 72 samples. Only one nonsynonymous change was detected at codon 1660 Arg/Trp (TGG –AGG), belonging to sample Pv066 (Eastern Amazon) (Table 1). A gene fragment of 909 bp (codons 224–525) was analysed from a total of 50 samples for pvdbpII. Sequence analysis using the Sal1 strain as reference showed 44 polymorphisms, 9 were synonymous and 35 nonsynonymous. Ten of these 35 nonsynonymous mutations have been reported frequently worldwide (R308S, K371E, D384G, E385K, K386N, R390H, N417K, L424I, W437R and I503K) [56]. In addition, only one trimorphic site was detected (I516F/M), whereas the others were dimorphic. Interestingly, sample Pv009 exhibited an insertion (CTA) that coded for a Leu (L475) and was the only sequence with a change in residue F306L. Diversity analysis indicated a total number of 28 haplotypes (H) among the 50 samples analysed, which represented a nucleotide diversity (π) of 0.00736, and a haplotype diversity (Hd) of 0.940 (Table 1).

Table 1 Comparison of P. vivax merozoite gene diversity among Ecuadorian samples, 2012–2015

A pvrbp1a-2 gene fragment of 852 bp (codons 351–634) was analysed from a total of 70 samples. Fifty-six polymorphic sites were detected, 9 were synonymous and 47 nonsynonymous. Five trimorphic sites were observed (N352Y/H, E360K/G, D388G/E, G599L/E and Q625R/P), while the rest were dimorphic. The analysis showed 31 haplotypes (H), which represented, a nucleotide diversity (π) of 0.00747 and a haplotype diversity (Hd) of 0.924 (Table 1).

The entire CDS of the pvama1 gene, with a total length of 1593 bp (codons 1-531), was analysed from 73 samples. Nucleotide sequence comparison with the Sal1 sequence revealed 33 SNPs, 6 of which were synonymous and 27 nonsynonymous. Furthermore, 18 of the 33 found in Ecuadorian isolates were frequent polymorphic residues found globally (G42V, R66K, D107A, R112T, K120R, N130K, N132D, L140I, A141E, E145A, E189K, K190Q, P210S, E277K, K352N, Q380K, L384P y E385D) [51]. Diversity indexes revealed a total number of 16 haplotypes (H), which represented a nucleotide diversity (π) of 0.00439, and a haplotype diversity (Hd) of 0.836 (Table 1).

Diversity and genetic differentiation of  P. vivax  invasion genes

There were 28 different haplotypes of pvdbpII based on the nucleotide sequence analysis of CDS for this receptor binding domain. Eastern Amazon was the locality that had most of the haplotypes identified (n = 17), followed by Coast (n = 8) and Western Amazon (n = 5) (Fig. 2). Remarkably, pvdbpII haplotypes were generally distinct among these locations as only 2 variants (dbpII-h1 and dbpII-h3) were shared. For pvrbp1a-2, 31 haplotypes were detected, most of them distributed in Eastern Amazon (n = 20), followed by Coast (n = 9) and Western Amazon (n = 5) (Fig. 2). Only one pvrbp1a-2 haplotype was shared within the three locations (rbpa-h5) and the variant (rbp1a-h3) was reported in both eastern and Western Amazon. A total of 16 haplotypes were identified. For pvama1, 16 haplotypes were detected. Most of the pvama1 haplotypes were located in the Eastern Amazon (n = 11), followed by Western Amazon (n = 5) and Coast (n = 3). Two variants (ama1-h2 and ama1-h12) were localized within Eastern and Western Amazon, while the haplotype ama1-h3 was shared between Coast and Western Amazon (Fig. 2).

Fig. 2
figure 2

Haplotype frequencies of Plasmodium vivax genes encoding merozoite surface proteins in Ecuador, 2012–2015. a pvdbpII haplotypes frequencies in the study sites. b pvrbp1a-2 haplotypes frequencies in the study sites. c pvama1 haplotypes frequencies in the study sites

Based on these data a haplotype network was constructed for pvdbpII, pvrbp1a-2 and pvama1. In general, all haplotypes belonging to the Eastern Amazon were more related to those from Western Amazon. Not unexpectedly, haplotypes localized in Coast exhibited less relationship with the variants from the Amazon Region. While haplotypes are mixed in different locations, most haplotypes are restricted to a particular location and most locations have related haplotypes (Fig. 3).

Fig. 3
figure 3

Haplotype networks of Plasmodium vivax genes encoding merozoite surface proteins in Ecuador, 2012–2015. The haplotype network for each marker is shown. Each circle represents one haplotype, and the circle size indicates frequency. Colours represent the geographical origin of the samples, Eastern Amazon (green), Western Amazon (orange) and Coast (purple)

In particular, analysis of the pvdbpII gene fragment showed the highest haplotype diversity in Western Amazon (0.933), while the highest value of nucleotide diversity was found in Eastern Amazon (0.00588). FST values showed moderate differentiation between parasites from the Eastern and Western Amazon (0.14257). In contrast, there was high differentiation between P. vivax from the Amazon Region and those from the Coast Region (Additional file 1: Table S1). In pvrbp1a-2, the highest haplotype diversity was found in Eastern Amazon (0.875), whereas nucleotide diversity value was higher in Western Amazon (0.00626). Moreover, FST values considered a great differentiation of the populations among three locations (FST>0.25). However, the FST value was lower (0.30473) when comparing the populations of the Eastern and Western Amazon (Additional file 2: Table S2).

In the case of pvama1, the Western Amazon locality showed the highest values of haplotype and nucleotide diversity, 0.893 and 0.0693 respectively. Similar to pvdbpII and pvrbp1a-2, FST values demonstrated that parasites from the Amazon regions had lower genetic differentiation among themselves (0.27113). Likewise, genetic differentiation was more elevated among populations from the Coast and Western Amazon (0.35536), and between Coast and Eastern Amazon (0.39884) (Additional file 3: Table S3).

Genetic comparison of  P. vivax  merozoites genes from Ecuador with geographic sites around the world

Six hundred and eighty eight sequences from different parts of the world were compiled for the comparative analysis of the ligand domain DBP region II (pvdbpII). The analysis of this region revealed that every country, except Mexico, had high levels of haplotype diversity (Hd), and ranged from 0.919 to 0.993. The major nucleotide diversity (π) was seen in Thailand with 0.01093 (Table 2). The largest haplotypic distance of Ecuadorian parasites measured by Fst values was with Iran (0.00936), India (0.05616) and Brazil (0.05981). In contrast, the closest distance was with samples from Colombia (Table 2).

Table 2 Diversity and genetic differentiation among diverse geographic populations of P. vivax using available sequences of the pvdbpII gene

A total of 510 sequences for the pvama1 gene from six different countries were included in the comparative analysis. The haplotype diversity (Hd) observed in this gene ranged from 0.909 to 0.996. Ecuador was the country with the lowest haplotype diversity (0.836) (Table 3). The lowest haplotype distance of Ecuadorian P. vivax was found with Venezuela (Fst = 0.12935), indicating moderate differentiation between both nations. In contrast, the highest Fst value was observed when compared with Sri Lanka (0.48394) (Table 3).

Table 3 Estimation of the diversity and genetic differentiation among diverse geographic populations of P. vivax using available sequences of the pvama1 gene

Natural selection in  pvdbpII  and  pvama1

Tests for neutrality were applied to determine if the allelic variants from Ecuadorian endemic locations may have resulted from potential selective forces. Fu & Li’s D* and F* neutrality tests, for the pvdbpII region revealed significant values (p < 0.05) of − 2.99684 and − 2.77579, respectively. Tajima´s D and McDonald-Kreitman test were not statistically significant, which indicated no departure from neutrality for the pvdpbII region (Table 4). On the other hand, values of Fu & Li’s D* (− 3.48451) and F* (− 1.50687) were highly significant for the pvrbp1a-2 gene fragment (p < 0.02), while Tajima´s D value was not significant (Table 4).

Table 4 Natural selection of P. vivax genes encoding merozoite proteins

In the pvama1 CDS, results for the Fu & Li’s D*, F*, and Tajima´s D were not significant. However, the result of the McDonald-Kreitman test revealed a highly significant value (p = 0.000007), which demonstrated a departure from neutrality in pvama1 gene of Ecuadorian P. vivax (Table 4).

To determine the most variable region of pvama1, genetic diversity was evaluated independently in the three domains. The analysis revealed that domain I (DI) was the most diverse in the whole CDS and domain III (D3) showed the lowest values (Additional file 4: Table S4). The values of Fu & Li D*, and F* neutrality tests and Tajima’s D in DI were nor significant. On the other hand, McDonald-Kreitman test resulted in a highly significant value (p = 0.009933) (Table 5). In D II, the values of Fu & Li D*, F* and Tajima’s were not significant, but they were for the McDonald-Kreitman test (p = 0.000707). The results of the McDonald-Kreitman test suggested departure from neutrality in both domains (Table 5).

Table 5 Natural selection on pvama1 from Ecuadorian samples, 2012–2015

Analysis of polymorphisms associated with neutralizing antibody epitopes in PvDBPII and PvAMA1

Recent studies using inhibitory monoclonal antibodies (mAbs) have shown conserved epitopes that elicit a broadly neutralizing response and may confer protective immunity against different strains of P. vivax [57, 58]. Therefore, the presence of polymorphic residues in epitopes of three broadly neutralizing murine antibodies (bNAbs) 2D10, 2H2 y 2C6, and of two human mAbs (053054 and 092096) in PvDBPII were evaluated. One insertion at position 430 (L430) and two amino acid changes (A435G and D440G) were observed within 2D10/2H2 epitopes (residues 413–441). Moreover, five substitutions (T468A, I471F/M, V472L, Y475S and Q480H) were present in 2C6 epitope (residues 465–485) (Fig. 4).

Fig. 4
figure 4

Three-dimensional structure of Plasmodium vivax DBPII showing the polymorphisms and the neutralizing mAb epitopes. a Solvent-accessible surface representation of two opposite sides of the PvDBPII three-dimensional structure. Broadly neutralizing epitopes of murine monoclonal antibodies 2D10 and 2H2 are shown in light red whereas the 2C6 epitope in light green. Epitopes of human neutralizing monoclonal antibodies (053054 and 092096) are shown in yellow. b Solvent-accessible surface representation of two opposite sides of the PvDBPII three-dimensional structure. Polymorphic residues are shown in blue and labelled

The discontinuous conformational epitope for mAb 053054 is composed of amino acids 264–281, 356–372 and 249, whereas 092096 binds to the epitope comprised of residues 270–289, 355–375, 249 and 219 [59]. According to this, four polymorphic residues (T359R, M361L and N372K and I374M) were found for both mAbs epitopes (Fig. 4).

Twenty-two predicted linear B-cell epitopes were found across PvAMA1, with 9 containing at least one polymorphic site. Among the 27 nonsynonymous substitutions found at PvAMA1, 15 sites are located within B-cell epitopes (Additional file 5: Table S5). Mutation sites inside predicted B cell epitopes were mapped on a three-dimensional structure of PvAMA1 (Fig. 5).

Fig. 5
figure 5

Three-dimensional structure of Plasmodium vivax AMA1 with the polymorphisms and the predicted B-cell epitopes mapped a Solvent-accessible surface representation of the active (left) and silent face (right) of the PvAMA1 three-dimensional structure. Polymorphic residues are coloured according to location: DI (cyan), DII (magenta) and DIII (orange). Hydrophobic ligand binding cleft residues are shown in dark blue. b Solvent-accessible surface representation of two opposite sides of the PvAMA1 three-dimensional structure. Residues of the linear B-cell epitopes are shown in light red. c Solvent-accessible surface representation of two opposite sides of the PvAMA1 three-dimensional structure. Polymorphic residues are shown in light blue

Discussion

In Ecuador, malaria has been a major public health problem over the last decades. Between 2001 and 2015 documented clinical cases decreased by more than 99%, leading the country towards the elimination phase. Nevertheless, several periodic outbreaks of malaria have occurred in the endemic regions of the country, and from 2015 the number of cases has increased, especially affecting some of the poorest inhabitants in Ecuador’s endemic areas. It is not yet certain whether the increasing incidence of vivax malaria is due to destabilization of local transmission patterns or new immunologically distinct parasite strains from different endemic regions.

Different biological and socioeconomic factors contribute, in several ways, to the maintenance of P. vivax diversity, including relapses (caused by hypnozoites), early gametocyte development, asymptomatic infections and the late diagnosis and treatment [53, 56]. However, the main mechanisms to generate and maintain diversity in the merozoite invasion genes occur from natural selection caused by the host immune responses and the diversity created by natural meiotic recombination among different variants of parasites inside the mosquito’s midgut [34, 60, 61].

In this study, diversity of invasion ligands was analyzed in Ecuadorian P. vivax. A novel substitution was found in pvmsp1-19. There were previously only three reports of mutations in this region: D1706E [62], K1709E [63] and N1692K [64]. Overall, pvmsp1-19 was highly conserved in Ecuadorian P. vivax parasites, which is consistent with previous studies reported worldwide [37, 63,64,65,66]. MSP119 is highly preserved [37, 64, 67]. It is presumed the highly conserved nature of MSP119 is due to an important functional role, especially in the food vacuole formation [68].

The analysis of pvdbpII revealed 35 nonsynonymous substitutions, of which 10 (R308S, K371E, D384G, E385K, K386N, R390H, N417K, L424I, W437R y I503K) have been widely reported in different endemic regions [56]. In addition, polymorphic residues in the positions N417K, W437R and I503K, found in most of the Ecuadorian samples, have been considered as targets of binding inhibitory antibodies (Fig. 4). This suggests that these mutations, alone or in combination, could affect the binding efficiency of the antibodies [69, 70]. Interestingly, residue F306L, which was found in sample Pv009 (Coast), had only been reported previously in isolates from Asia [40, 49, 71]. The same sample had a CTA codon insertion, which coded for a leucine between amino acid 474 and 475. This insertion has been found in samples from India, Iran, Indonesia, Thailand, and Brazil [42].

A surprisingly high level of genetic diversity was discovered in Ecuadorian parasites, and pvdbpII was the highest in terms of haplotype diversity (0.940). This finding is consistent with several studies in different endemic countries in which similar levels of diversity have been found [33, 34, 48, 71, 72]. The results indicate that pvdbpII is consistently polymorphic at a national and worldwide scale, although there was no evidence of departure from neutrality in this region. However, significant values in Fu & Li’s D* and F* could be suggesting selective pressure. This would imply that a mutation involved in a more efficient way to evade the host immune response occurred within this region, thereby increasing parasite fitness. Such mutations provide an adaptative advantage in individuals to increase frequency in the population. These results are similar to those obtained in Mexico by González-Cerón and collaborators in 2015, although studies in other endemic regions revealed evidence of balancing selection in pvdbpII ligand domain that is considered to be the target of anti-DBP protective immunity [33, 48, 49, 71,72,73].

PvDBPII inhibitory mAbs epitopes in Ecuadorian parasites revealed a lack of high-frequency polymorphisms. This is consistent with a previous study, in which DBPII sequences that represent global diversity were analyzed for retention of the epitopes for 2D10, 2H2, and 2C6 antibodies [57, 59]. Together, 2D10, 2H2, and 2C6 recognize broadly conserved epitopes and play a role as bNAbs that will be effective against different P. vivax strains [57]. Only four amino acid changes (T359R, M361L and N372K and I374M) were identified in human mAbs epitopes (053054 and 092096) (Fig. 4). A comparative analysis demonstrates that mutations T359R, N372K and I374M individually and in combination do not affect antibody binding to DBPII.

Some epitopes accumulate a greater number of polymorphisms in order to create immunologically-distinct B-cell epitopes, allowing the parasite to evade the acquired immune responses and lead to successful invasion of the reticulocytes. On the other hand, there are regions evolving neutrally. The contrast of selective pressures in pvdbpII has been previously reported for region IV [34, 74]. This may imply a balance between the high functional restriction, maintaining invariable structures due to the important role that DBP protein plays in reticulocyte binding and the search for diversity as a response against the pressure of the host’s immune system [74].

The analysis of pvrbp1a-2 gene fragment revealed 9 synonymous substitutions and 47 nonsynonymous. Interestingly, sequences from four isolates (Pv009, Pv034, Pv057 and Pv066) were identical to the reference parasite Sal 1. pvrbp1a-2 had the highest nucleotide diversity (0.00747) (Table 1). Previous studies in pvrbp1a have shown that the region between amino acids 435 and 777 is the most polymorphic of the whole protein, suggesting it as a target of antibody response [75, 76].

Neutrality tests carried out in this gene fragment showed highly significant values for Fu & Li’s D* and F*, but not significant for Tajima’s D value (Table 4). A mutation could have increased rapidly its frequency within the population and led to a decrease in variability in this fragment. In addition, McDonald-Kreitman test could not be performed to corroborate this result due to P. knowlesi pvrbp1a orthologous gene corresponded to a pseudogene highly degenerated (ψpknbp1) possibly due to an ancient recombination event [77]. The results in this study suggest that pvdbpII, pvrbp1a-2 went through a selective sweep process as an adaptative advantage.

The analysis of the whole pvama1 gene showed a total of 33 polymorphic sites, 6 were synonymous and 27 nonsynonymous (Fig. 5). Of these, 18 have been widely reported in many endemic countries [51]. Besides, 15 polymorphisms have been identified in predicted linear B-cell epitopes and might play an important role in immune evasion process (Additional file 5: Table S5). This is consistent with a previous study in which residues encountered at positions E145A, P210S, R249H, K352N and N445D have been overlapping B-cell epitopes [78]. These findings suggested that variations in amino acids could affect the protein structure, causing changes in charge and polarity, and helping the parasite escape from the host immune response [79]. In the same way, mutations Q380K and L384P, found in this study, have been previously suggested to be implicated in the PvAMA1 protein structure modification [61, 78].

Furthermore, the major number of polymorphic residues were on the active face of the PvAMA1 ectodomain, principally within domain I and II (Fig. 5). This suggests that this side of the protein is more exposed to the immune response [51]. In contrast, only four polymorphic residues (G42V, T44S, R66K and S70F) were located on the other side of this protein (silent face), which seems to have no relevance in the invasion process (Fig. 5) [79]. In addition, several polymorphic residues were grouped forming shield-type structure, which maintain preserved the RON2 gap union site. This could be a parasite’s mechanism to divert the immune response, accumulating mutations in non-relevant regions and thus keep highly conserved relevant structures for the invasive process [79] (Fig. 5).

Although values of Fu & Li’s D* and F* and Tajima’s D were non-significant, a slightly departure from neutrality was shown. Nevertheless, McDonald-Kreitman test determined a highly significant value (p = 0.000007), which evidenced the presence of balancing selection in this gene (Tables 4 and 5). These results were consistent with several investigations that found evidence of balancing selection, principally, in domains I [35, 51, 80, 81] and II [31, 61, 82]. Thus, the selective pressure by the human immune system causes that parasites keep mutations and maintain diversity in both domains to evade the immune responses.

The haplotype arrangement of pvdbpII, pvrbp1a-2 and pvama1 revealed a major number of variants in Eastern Amazon (Fig. 2). Despite the low and unstable transmission in endemic areas of Ecuador, several periodic outbreaks have been recently reported in the east of the country. This event could lead to an increase in parasite circulation and, therefore a high level of recombination. Most samples come from a border area where there is common cross border movement, explaining a higher diversity of P. vivax in this locality. It has been shown that genetic diversity, particularly in pvdbpII, is not related to malaria endemic levels, since similar patterns have been observed in areas with low and unstable transmission, such as parts of the Brazilian Amazon and Papua New Guinea [56].

Only a few haplotypes were shared among localities and most of them showed a greater genetic relationship according to the locality where they belong. Moreover, the variants found in the Western Amazon were more related to those from the Eastern Amazon, and both had less genetic similarity with the haplotypes from the Coast (Fig. 3). Similarly, FST analyses exhibited a high difference between parasites from the Coast Region and those from the Amazon Region (0.355–0.556). (Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3). These results are in accordance with previous reports using microsatellites in the same samples [83] that revealed important difference among Amazonian and Coastal parasites due to the geographical distance and the presence of geological barriers such as the Andes Mountain range. Although the diversity found in invasion ligand genes (pvdbpII and pvama1) in Ecuador was high, when comparing it to other parts of the world, the country had one of the lowest values (Tables 2 and 3).

Analysis of global pvdbpII sequences revealed low differentiation of parasite populations among countries, finding the same haplotypes in different countries, which is consistent with a previous study by De Sousa et al. [56], in which a global analysis of pvdbpII showed low differentiation of parasite populations among countries. Moreover, haplotypes found were all present in the eight different countries.

These findings confirm that merozoite antigenic diversity, rather than being influenced by geographic barriers, is more likely dependent on intrinsic characteristics of human populations, principally immune response [84]. Moreover, due to the differences in the endemicity spectrum, both recombination and natural selection (product of host immune system) seem to be differentially influenced between distinct geographical areas [56], which suggests that both mechanisms play a determinant role in the haplotype structure.

Conclusions

The differences in diversity in the main P. vivax genes involved in reticulocyte invasion depend on the role that each play in this process. Possibly, the prolonged exposure of PvDBPII and PvRBP1a-2 to the host immune system contributed to maintain high diversity. PvAMA1, which seemed to be less exposed, has a moderate diversity. On the contrary, PvMSP1-19 despite of being present from the initial attachment, is not fully exposed and, therefore, is a conserved region. pvdbpII and pvama1 genetic diversity found in Ecuadorian P. vivax confirmed that diversity in these genes is not related to levels of endemicity. The major diversity of P. vivax in Ecuador was identified within the Amazonian localities possibly because of parasite migration from Peru in addition to locally circulating parasites and is in concordance with the increase in the number of infections in this area of the country. Natural selection behaves differently in each gene. Even more, acts differentially through them. pvdbpII and pvrbp1a-2 were possibly affected by a selective sweeping process, whereas in pvama1 the influence of balancing selection was evidenced, mainly in domains I and II .

Availability of data and materials

The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

P. vivax :

Plasmodium vivax

pvmsp1-19 :

Plasmodium vivax merozoite surface protein 1 carboxyl fragment of 19 kD

pvdbpII :

Plasmodium vivax duffy binding protein region II

pvrbp1a-2 :

Plasmodium vivax reticulocyte binding protein 1a region 2

pvama1 :

Plasmodium vivax apical membrane antigen 1

DARC:

Duffy antigen chemokine receptor

RON2:

Rhoptry neck protein 2

RBC:

Red blood cell

BRZ:

Brazil

COL:

Colombia

IND:

India

IRN:

Iran

SK:

South Korea

MEX:

Mexico

PNG:

Papua New Guinea

SLK:

Sri Lanka

THL:

Thailand

VNZ:

Venezuela

References

  1. WHO. World Malaria Report 2022. Geneva, World Health Organization, 2023. https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2022. Accessed 13 Jan 2023.

  2. Ministerio de Salud Pública del Ecuador (MSP). Gacetas Vectoriales SE-51-2021. 2021. https://www.salud.gob.ec/wp-content/uploads/2022/01/GACETA-GENERAL-VECTORIALES-SE-51.pdf. Accessed 16 Jan 2022.

  3. Sáenz FE, Arévalo-Cortés A, Valenzuela G, Vallejo AF, Castellanos A, Poveda-Loayza A, et al. Malaria epidemiology in low-endemicity areas of the northern coast of Ecuador: high prevalence of asymptomatic infections. Malar J. 2017;16:300.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Veloz-Haro E. Duffy blood group phenotypes/genotypes and their association with malaria prevalence in four communities of northwest Ecuador. (Disertación de pregrado). Ecuador: Pontificia Universidad Católica del Ecuador, Quito; 2015.

  5. Weiss GE, Gilson PR, Taechalertpaisarn T, Tham WH, de Jong NW, Harvey KL, et al. Revealing the sequence and resulting cellular morphology of receptor-ligand interactions during Plasmodium falciparum invasion of erythrocytes. PLoS Pathog. 2015;11:e1004670.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Gilson PR, Nebl T, Vukcevic D, Moritz RL, Sargeant T, Speed TP, et al. Identification and stoichiometry of glycosylphosphatidylinositol-anchored membrane proteins of the human malaria parasite Plasmodium falciparum. Mol Cell Proteom. 2006;5:1286–99.

    Article  CAS  Google Scholar 

  7. Stafford WH, Günder B, Harris A, Heidrich HG, Holder AA, Blackman MJ. A 22 kDa protein associated with the Plasmodium falciparum merozoite surface protein-1 complex. Mol Biochem Parasit. 1996;80:159–69.

    Article  CAS  Google Scholar 

  8. Trucco C, Fernandez-Reyes D, Howell S, Stafford WH, Scott-Finnigan TJ, Grainger M, et al. The merozoite surface protein 6 gene codes for a 36 kDa protein associated with the Plasmodium falciparum merozoite surface protein-1 complex. Mol Biochem Parasit. 2001;112:91–101.

    Article  CAS  Google Scholar 

  9. Pachebat JA, Kadekoppala M, Grainger M, Dluzewski AR, Gunaratne RS, Scott-Finnigan TJ, et al. Extensive proteolytic processing of the malaria parasite merozoite surface protein 7 during biosynthesis and parasite release from erythrocytes. Mol Biochem Parasit. 2007;151:59–69.

    Article  CAS  Google Scholar 

  10. Blackman MJ, Heidrich HG, Donachie S, McBride JS, Holder AA. A single fragment of a malaria merozoite surface protein remains on the parasite during red cell invasion and is the target of invasion-inhibiting antibodies. J Exp Med. 1990;172:379–82.

    Article  CAS  PubMed  Google Scholar 

  11. Blackman MJ, Holder AA. Secondary processing of the Plasmodium falciparum merozoite surface protein-1 (MSP1) by a calcium-dependent membrane-bound serine protease: shedding of MSP133 as a noncovalently associated complex with other fragments of the MSP1. Mol Biochem Parasit. 1992;50:307–15.

    Article  CAS  Google Scholar 

  12. Sachdeva S, Ahmad G, Malhotra P, Mukherjee P, Chauhan VS. Comparison of immunogenicities of recombinant Plasmodium vivax merozoite surface protein 1 19-and 42-kiloDalton fragments expressed in Escherichia coli. Infect Immun. 2004;72:5775–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Fernandez-Becerra C, Sanz S, Brucet M, Stanisic DI, Alves FP, Camargo, et al. Naturally-acquired humoral immune responses against the N-and C-termini of the Plasmodium vivax MSP1 protein in endemic regions of Brazil and Papua New Guinea using a multiplex assay. Malar J. 2010;9:29.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Egan AF, Burghaus P, Druilhe P, Holder AA, Riley EM. Human antibodies to the 19 kDa C-terminal fragment of Plasmodium falciparum merozoite surface protein 1 inhibit parasite growth in vitro. Parasite Immunol. 1999;21:133–9.

    Article  CAS  PubMed  Google Scholar 

  15. John CC, O’Donnell RA, Sumba PO, Moormann AM, de Koning-Ward TF, King, et al. Evidence that invasion-inhibitory antibodies specific for the 19-kDa fragment of merozoite surface protein-1 (MSP-119) can play a protective role against blood-stage Plasmodium falciparum infection in individuals in a malaria endemic area of Africa. J Immunol. 2004;173:666–72.

    Article  CAS  PubMed  Google Scholar 

  16. Okech BA, Corran PH, Todd J, Joynson-Hicks A, Uthaipibull C, Egwang, et al. Fine specificity of serum antibodies to Plasmodium falciparum merozoite surface protein, PfMSP-119, predicts protection from malaria infection and high-density parasitemia. Infect Immun. 2004;72:1557–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Chitnis CE, Sharma A. Targeting the Plasmodium vivax Duffy-binding protein. Trends Parasitol. 2008;24:29–34.

    Article  CAS  PubMed  Google Scholar 

  18. Galinski MR, Medina CC, Ingravallo P, Barnwell JW. A reticulocyte-binding protein complex of Plasmodium vivax merozoites. Cell. 1992;69:1213–26.

    Article  CAS  PubMed  Google Scholar 

  19. Adams JH, Hudson DE, Torii M, Ward GE, Wellems TE, Aikawa M, et al. The Duffy receptor family of Plasmodium knowlesi is located within the micronemes of invasive malaria merozoites. Cell. 1990;63:141–53.

    Article  CAS  PubMed  Google Scholar 

  20. Adams JH, Sim BK, Dolan SA, Fang X, Kaslow DC. Miller LH.A family of erythrocyte binding proteins of malaria parasites. Proc Natl Acad Sci USA. 1992;89:7085–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Singh SK, Hora R, Belrhali H, Chitnis CE, Sharma A. Structural basis for Duffy recognition by the malaria parasite duffy-binding-like domain. Nature. 2006;439:741–4.

    Article  CAS  PubMed  Google Scholar 

  22. Grimberg BT, Udomsangpetch R, Xainli J, McHenry A, PanichakulT, Sattabongkot J, et al. Plasmodium vivax invasion of human erythrocytes inhibited by antibodies directed against the Duffy binding protein. PLoS Med. 2007;4:e337.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Han JH, Lee SK, Wang B, Muh F, Nyunt MH, Na S, et al. Identification of a reticulocyte-specific binding domain of Plasmodium vivax reticulocyte-binding protein 1 that is homologous to the PfRh4 erythrocyte-binding domain. Sci Rep. 2016;6: 26993.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Gupta ED, Anand G, Singh H, Chaddha K, Bharti PK, Singh, et al. Naturally acquired human antibodies against reticulocyte-binding domains of Plasmodium vivax proteins, PvRBP2c and PvRBP1a, exhibit binding-inhibitory activity. J Infect Dis. 2017;215:1558–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ntumngia FB, Thomson-Luque R, Galusic S, Frato G, Frischmann S, Peabody DS, et al. Identification and immunological characterization of the ligand domain of Plasmodium vivax reticulocyte binding protein 1a. J Infect Dis. 2018;218:1110–8.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Gupta S, Singh S, Popovici J, Roesch C, Shakri AR, Guillotte-Blisnick M, et al. Targeting a reticulocyte binding protein and Duffy binding protein to inhibit reticulocyte invasion by Plasmodium vivax. Sci Rep. 2018;8:10511.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 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 Parasit. 1994;67:59–68.

    Article  CAS  Google Scholar 

  28. Besteiro S, Dubremetz JF, Lebrun M. The moving junction of apicomplexan parasites: a key structure for invasion. Cell Microbiol. 2011;13:797–805.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Pizarro JC, Vulliez-Le Normand B, Chesne-Seck ML, Collins CR, Withers-Martinez C, Hackett F, et al. Crystal structure of the malaria vaccine candidate apical membrane antigen 1. Science. 2005;308:408–11.

    Article  CAS  PubMed  Google Scholar 

  31. 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:1266–73.

    Article  PubMed  Google Scholar 

  32. 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:9–20.

    Article  CAS  PubMed  Google Scholar 

  33. Ampudia E, Patarroyo MA, Patarroyo ME, Murillo LA. Genetic polymorphism of the Duffy receptor binding domain of Plasmodium vivax in colombian wild isolates. Mol Biochem Parasit. 1996;78:269–72.

    Article  CAS  Google Scholar 

  34. Sousa TN, Tarazona-Santos EM, Wilson DJ, Madureira AP, Falcão PR, Fontes, et al. Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations. Malar J. 2010;9:334.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Bittencourt NC, Silva ABIEd Virgili NS, Schappo AP, Gervásio JHDB, Pimenta TS, et al. Plasmodium vivax AMA1: implications of distinct haplotypes for immune response. PLoS Negl Trop Dis. 2020;14:e0008471.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. 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:651.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Snounou G. Detection and identification of the four malaria parasite species infecting humans by PCR amplification. In: Clapp JP, editor. Species diagnostics protocols: PCR and other nucleic acid methods. Humana Press: New Jersey; 1996. p. 263–91.

    Google Scholar 

  39. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.

    Article  CAS  PubMed  Google Scholar 

  40. Babaeekho L, Zakeri S, Djadid ND. Genetic mapping of the Duffy binding protein (DBP) ligand domain of Plasmodium vivax from unstable malaria region in the Middle East. Am J Trop Med Hyg. 2009;80:112–8.

    Article  CAS  Google Scholar 

  41. Valizadeh V, Zakeri S, Mehrizi AA, Djadid ND. Non-allele specific antibody responses to genetically distinct variant forms of Plasmodium vivax Duffy binding protein (PvDBP-II) in Iranians exposed to seasonal malaria transmission. Acta Trop. 2014;136:89–100.

    Article  CAS  PubMed  Google Scholar 

  42. Valizadeh V, Zakeri S, Mehrizi AA, Djadid ND. Population genetics and natural selection in the gene encoding the Duffy binding protein II in iranian Plasmodium vivax wild isolates. Infect Genet Evol. 2014;21:424–35.

    Article  CAS  PubMed  Google Scholar 

  43. Ju HL, Kang JM, Moon SU, Bahk YY, Cho PY, Sohn WM, et al. Genetic diversity and natural selection of Duffy binding protein of Plasmodium vivax Korean isolates. Acta Trop. 2013;125:67–74.

    Article  CAS  PubMed  Google Scholar 

  44. Kho WG, Chung JY, Sim EJ, Kim DW, Chung WC. Analysis of polymorphic regions of Plasmodium vivax Duffy binding protein of korean isolates. Korean J Parasitol. 2001;39: 143.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Xainli J, Adams JH, King CL. The erythrocyte binding motif of Plasmodium vivax Duffy binding protein is highly polymorphic and functionally conserved in isolates from Papua New Guinea. Mol Biochem Parasit. 2000;111:253–60.

    Article  CAS  Google Scholar 

  46. Cole-Tobian JL, Biasor M, King CL. High complexity of Plasmodium vivax infections in Papua New guinean children. Am J Trop Med Hyg. 2005;73:626–33.

    Article  PubMed  Google Scholar 

  47. Cole-Tobian JL, Cortés A, Baisor M, Kastens W, Xainli J, Bockarie M, et al. Age-acquired immunity to a Plasmodium vivax invasion ligand, the Duffy binding protein. J Infect Dis. 2002;186:531–9.

    Article  CAS  PubMed  Google Scholar 

  48. Premaratne PH, Aravinda BR, Escalante AA, Udagama PV. Genetic diversity of Plasmodium vivax Duffy binding protein II (PvDBPII) under unstable transmission and low intensity malaria in Sri Lanka. Infect Genet Evol. 2011;11:1327–39.

    Article  CAS  PubMed  Google Scholar 

  49. Gosi P, Khusmith S, Khalambaheti T, Lanar DE, Schaecher KE, Fukuda MM, et al. Polymorphism patterns in Duffy-binding protein among thai Plasmodium vivax isolates. Malar J. 2008;7: 112.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 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.

    Article  CAS  PubMed  Google Scholar 

  51. 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 (pv AMA1). PLoS Negl Trop Dis. 2013;7: e2506.

    Article  PubMed  PubMed Central  Google Scholar 

  52. 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.

    Article  CAS  PubMed  Google Scholar 

  53. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Jespersen MC, Peters B, Nielsen M, Marcatili P. BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. 2017;45:W24-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Bittencourt NC, Silva ABIED, Virgili NS, Schappo AP, Gervásio JHD, Pimenta TS, et al. Plasmodium vivax AMA1: implications of distinct haplotypes for immune response. PLoS Negl Trop Dis. 2020;14:e0008471.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. De Sousa TN, Carvalho LH, de Brito CFA. Worldwide genetic variability of the Duffy binding protein: insights into Plasmodium vivax vaccine development. PLoS ONE. 2011;6:e22944.

    Article  Google Scholar 

  57. Chen E, Salinas ND, Huang Y, Ntumngia F, Plasencia MD, Gross ML, Tolia NH. Broadly neutralizing epitopes in the Plasmodium vivax vaccine candidate Duffy binding protein. Proc Natl Acad Sci USA. 2016;113:6277–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Carias LL, Dechavanne S, Nicolete VC, Sreng S, Suon S, Amaratunga C. Identification and characterization of functional human monoclonal antibodies to Plasmodium vivax duffy-binding protein. J Immunol Res. 2019;202:2648–60.

    CAS  Google Scholar 

  59. Urusova D, Carias L, Huang Y, Nicolete VC, Popovici J, Roesch C, Tolia NH. Structural basis for neutralization of Plasmodium vivax by naturally acquired human antibodies that target DBP. Nat Microbiol. 2019;4:1486–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Escalante AA, Grebert HM, Chaiyaroj SC, Magris M, Biswas S, Nahlen BL, et al. Polymorphism in the gene encoding the apical membrane antigen-1 (AMA-1) of Plasmodium falciparum. X. Asembo Bay Cohort Project. Mol Biochem Parasit. 2001;113:279–87.

    Article  CAS  Google Scholar 

  61. 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.

    Article  CAS  PubMed  Google Scholar 

  62. Zeyrek FY, Tachibana SI, Yuksel F, Doni N, Palacpac N, Arisue N, et al. Limited polymorphism of the Plasmodium vivax merozoite surface protein 1 gene in isolates from Turkey. Am J Trop Med Hyg. 2010;83:1230–7.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Pasay MC, Cheng Q, Rzepczyk C, Saul A. Dimorphism of the C terminus of the Plasmodium vivax merozoite surface protein 1. Mol Biochem Parasit. 1995;70:217–9.

    Article  CAS  Google Scholar 

  64. Kang JM, Ju HL, Kang YM, Lee DH, Moon SU, Sohn WM, et al. Genetic polymorphism and natural selection in the C-terminal 42 kDa region of merozoite surface protein-1 among Plasmodium vivax Korean isolates. Malar J. 2012;11: 206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Soares IS, Barnwell JW, Ferreira MU, Da Cunha MG, Laurino JP, Castilho BA, et al. A Plasmodium vivax vaccine candidate displays limited allele polymorphism, which does not restrict recognition by antibodies. Mol Med. 1999;5:459–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Putaporntip C, Jongwutiwes S, Sakihama N, Ferreira MU, Kho WG, Kaneko A, et al. Mosaic organization and heterogeneity in frequency of allelic recombination of the Plasmodium vivax merozoite surface protein-1 locus. Proc Natl Acad Sci USA. 2002;99:16348–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Huaman MC, Martin LB, Malkin E, Narum DL, Miller LH, Mahanty S, et al. Ex vivo cytokine and memory T cell responses to the 42-kDa fragment of Plasmodium falciparum merozoite surface protein-1 in vaccinated volunteers. J Immunol. 2008;180:1451–61.

    Article  CAS  PubMed  Google Scholar 

  68. Dluzewski AR, Ling IT, Hopkins JM, Grainger M, Margos G, Mitchell GH, et al. Formation of the food vacuole in Plasmodium falciparum: a potential role for the 19 kDa fragment of merozoite surface protein 1 (MSP119). PLoS ONE. 2008;3:e3085.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Van Buskirk KM, Tobian JLC, Baisor M, Sevova ES, Bockarie M, King CL, et al. Antigenic drift in the ligand domain of Plasmodium vivax Duffy binding protein confers resistance to inhibitory antibodies. J Infect Dis. 2004;190:1556–62.

    Article  Google Scholar 

  70. McHenry AM, Barnes SJ, Ntumngia FB, King CL, Adams JH. Determination of the molecular basis for a limited dimorphism, N417K, in the Plasmodium vivax Duffy-binding protein. PLoS ONE. 2011;6: e20192.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Ju HL, 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 

  72. Ju HL, Kang JM, Moon SU, Bahk YY, Cho PY, Sohn WM, et al. Genetic diversity and natural selection of Duffy binding protein of Plasmodium vivax Korean isolates. Acta Trop. 2013;125:67–74.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  74. Martinez P, Suarez CF, Cardenas PP, Patarroyo MA. Plasmodium vivax Duffy binding protein: a modular evolutionary proposal. Parasitology. 2004;128:353–66.

    Article  CAS  PubMed  Google Scholar 

  75. Caraballo JR, Delgado G, Rodriguez R, Patarroyo MA. The antigenicity of a Plasmodium vivax reticulocyte binding protein-1 (PvRBP1) recombinant fragment in humans and its immunogenicity and protection studies in Aotus monkeys. Vaccine. 2007;25:3713–21.

    Article  Google Scholar 

  76. Ferreira AR, Singh B, Cabrera-Mora M, De Souza ACM, Marques MTQ, Porto, et al. Evaluation of naturally acquired IgG antibodies to a chimeric and non-chimeric recombinant species of Plasmodium vivax reticulocyte binding protein-1: lack of association with HLA-DRB1*/DQB1* in malaria exposed individuals from the brazilian Amazon. PLoS ONE. 2014;9:e105828.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Meyer EV, Semenya AA, Okenu DM, Dluzewski AR, Bannister LH, Barnwell JW, et al. The reticulocyte binding-like proteins of P. knowlesi locate to the micronemes of merozoites and define two new members of this invasion ligand family. Mol Biochem Parasit. 2009;165:111–21.

    Article  CAS  Google Scholar 

  78. 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.

    Article  CAS  PubMed  Google Scholar 

  79. Bai T, Becker M, Gupta A, Strike P, Murphy VJ, Anders RF, et al. Structure of AMA1 from Plasmodium falciparum reveals a clustering of polymorphisms that surround a conserved hydrophobic pocket. Proc Natl Acad Sci USA. 2005;102:12736–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Arnott A, Wapling J, Mueller I, Ramsland PA, Siba PM, Reeder JC, et al. Distinct patterns of diversity, population structure and evolution in the AMA1 genes of sympatric Plasmodium falciparum and Plasmodium vivax populations of Papua New Guinea from an area of similarly high transmission. Malar J. 2014;13: 233.

    Article  PubMed  PubMed Central  Google Scholar 

  81. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Gentil F, Bargieri DY, Leite JA, Francoso KS, Patricio MB, Espindola N, 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.

    Article  CAS  PubMed  Google Scholar 

  83. Vélez-Álvarez E. Caracterización poblacional de Plasmodium vivax en Ecuador mediante microsatélites y análisis de resistencia de genes de resistencia frente a antimaláricos. (Disertación de pregrado). ESPE. Universidad de las Fuerzas Armadas, Quito, Ecuador. 2015.

  84. Barry AE, Schultz L, Buckee CO, Reeder JC. Contrasting population structures of the genes encoding ten leading vaccine-candidate antigens of the human malaria parasite, Plasmodium falciparum. PLoS ONE. 2009;4: e8497.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank the participants from all the communities and the members of SNEM and the Ministry of Health of Ecuador who participated in field work. In particular we thank Luis Enrique Castro (MOH) and Roberto Montoya (PAHO) for their involvement in the training of the personnel for sample collection.

Funding

This research was funded by the World Academy of Sciences (Grant Number 16–158 RG/BIO/LD), Pontificia Universidad Católica del Ecuador (Grants M131416 and N131416).

Author information

Authors and Affiliations

Authors

Contributions

FES, JHA and AN conceived and designed the study, AN, FN, YG and FES wrote the manuscript; AN and FN performed the experiments; FES, AN and JHA analyzed the data. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Fabián E. Sáenz.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Ethical Review Committee of Pontificia Universidad Católica del Ecuador Approval number (#CBE-016-2013). Written informed consent (IC) was provided by study participants and/or their legal guardians for sample collection.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

 Diversity and genetic differentiation in pvdbpII among Ecuadorian locations.

Additional file 2: Table S2.

 Diversity and genetic differentiation in pvrbp1a-2 among Ecuadorian locations.

Additional file 3: Table S3.

 Diversity and genetic differentiation in pvama1 among Ecuadorian locations.

Additional file 4: Table S4.

 Genetic diversity among pvama1 domains in Ecuadorian samples.

Additional file 5: Table S5.

 Polymorphisms associated with predicted linear B-cell epitopes in PvAMA1 among Ecuadorian samples.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Núñez, A., Ntumngia, F., Guerra, Y. et al. Genetic diversity and natural selection of Plasmodium vivax reticulocyte invasion genes in Ecuador. Malar J 22, 225 (2023). https://doi.org/10.1186/s12936-023-04640-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-023-04640-0

Keywords

  • Plasmodium vivax
  • Ecuador
  • Genetic diversity
  • Merozoite invasion
  • Reticulocytes
  • Natural selection