- Open Access
Structural patterns of selection and diversity for Plasmodium vivax antigens DBP and AMA1
© The Author(s) 2018
- Received: 12 January 2018
- Accepted: 18 April 2018
- Published: 2 May 2018
Plasmodium vivax is a significant contributor to the global malaria burden, and a vaccine targeting vivax malaria is urgently needed. An understanding of the targets of functional immune responses during the course of natural infection will aid in the development of a vaccine. Antibodies play a key role in this process, with responses against particular epitopes leading to immune selection pressure on these epitopes. A number of techniques exist to estimate levels of immune selection pressure on particular epitopes, with a sliding window analysis often used to determine particular regions likely to be under immune pressure. However, such analysis neglects protein three-dimensional structural information. With this in mind, a newly developed tool, BioStructMap, was applied to two key antigens from Plasmodium vivax: PvAMA1 and PvDBP Region II. This tool incorporates structural information into tests of selection pressure.
Sequences from a number of populations were analysed, examining spatially-derived nucleotide diversity and Tajima’s D over protein structures for PvAMA1 and PvDBP. Structural patterns of nucleotide diversity were similar across all populations examined, with Domain I of PvAMA1 having the highest nucleotide diversity and displaying significant signatures of immune selection pressure (Tajima’s D > 0). Nucleotide diversity for PvDBP was highest bordering the dimerization and DARC-binding interface, although there was less evidence of immune selection pressure on PvDBP compared with PvAMA1. This study supports previous work that has identified Domain I as the main target of immune-mediated selection pressure for PvAMA1, and also supports studies that have identified functional epitopes within PvDBP Region II.
The BioStructMap tool was applied to leading vaccine candidates from P. vivax, to examine structural patterns of selection and diversity across a number of geographic populations. There were striking similarities in structural patterns of diversity across multiple populations. Furthermore, whilst regions of high diversity tended to surround conserved binding interfaces, a number of protein regions with very low diversity were also identified, and these may be useful targets for further vaccine development, given previous evidence of functional antibody responses against these regions.
- Plasmodium vivax
- Protein structure
- Immune selection
- Population genetics
Plasmodium vivax infected an estimated 8.55 million people in 2016 and is a significant contributor to global malaria morbidity, with the majority of P. vivax cases occurring within South-East Asia . There remains a significant need for a vaccine against P. vivax, and an understanding of the targets of natural immune responses following P. vivax infection is likely to aid such an effort. A key challenge in vaccine development is the identification of specific antigens and epitopes that are targets of protective antibody responses. It is possible to use population genetic data to identify regions of proteins that are under immune-mediated selection pressure, which gives rise to balancing selection within that protein region. Tajima’s D is one test statistic that is often used to identify departures from a neutral model of selection, and has been applied to malaria genes both on a per-gene basis [2, 3], or as a sliding window analysis along a gene [4–7]. A number of studies have previously examined P. vivax proteins such as apical membrane antigen 1 (PvAMA1) and Duffy-binding protein (PvDBP) for evidence of immune selection pressure using these approaches [4, 8–11]. However, because a sliding window analysis is typically performed over the linear gene sequence, it does not take into account the impact of the three-dimensional (3D) structural constraints of the protein in the calculation of selection pressures. A new method that allows incorporation of protein structural information into tests for selection pressure has recently been described , and has been applied here to two leading P. vivax vaccine candidates: PvAMA1 and PvDBP.
AMA1 is a type I transmembrane protein present in all Plasmodium species . It is localized to the parasite micronemes, and is released onto the surface of the merozoite prior to invasion of red blood cells . AMA1 binds to RON2 during the formation of the tight junction between parasite and host–cell membranes [15, 16] and is a target of protective immune responses [17–20]. The ectodomain of Plasmodium AMA1 proteins is divided into three domains, termed Domains I (DI), II (DII) and III (DIII) (Additional file 1) . DI is considered to be the most polymorphic, and is also the site of RON2 binding . RON2 binds a conserved hydrophobic cleft that is surrounded by a number of highly polymorphic regions, the most notable being the C1L loop, a surface exposed loop with high variability that is suggested to define strain-specificity in anti-AMA1 responses in Plasmodium falciparum infection . While DI is generally considered to be the most important for functional antibody responses, there is evidence that DII and DIII may also be targets of functional antibody responses in P. falciparum [8, 24]. A number of studies have investigated selection pressures on PvAMA1. Evidence for balancing selection within Domain I has been observed in a Venezuelan population , in a number of Papua New Guinean populations (Madang and Wosera, Madang and East Sepik) [4, 8], an Iranian population  and in a Peruvian population . In contrast, two other studies examining PvAMA1 sequences in isolates from Korea  and Myanmar  did not find any evidence of balancing selection, but instead observed evidence of recent population bottleneck and expansion in those populations.
PvDBP is an important micronemal protein that binds to the Duffy antigen/receptor for chemokines (DARC) on human reticulocytes during invasion [28, 29]. Whilst there is evidence that PvDBP is not absolutely essential for invasion of reticulocytes [30–32], Duffy-negative individuals are largely resistant to P. vivax infection, and hence PvDBP makes an attractive vaccine target . PvDBP is part of the erythrocyte-binding like (EBL) family of proteins, which include EBA175, EBA181 and EBA140 in P. falciparum . PvDBP is the sole EBL family protein in P. vivax . EBL family proteins are composed of a number of distinct domains, with Region II (RII) being a cysteine-rich Duffy-binding like (DBL) domain that is involved in binding to erythrocytes. EBL family proteins each recognize a different receptor via their respective DBL domains [29, 36–38]; PvDBP binds to DARC via its DBL domain (RII) . During this process two PvDBP molecules form a dimer around two DARC molecules [39, 40]. RII of PvDBP has been divided into a number of subdomains (subdomains 1–3)  (Additional file 2), and it is subdomain 2 that contains both the dimer interface and DARC binding residues . Immune responses against PvDBP have been associated with protection from clinical malaria in naturally exposed cohorts [42, 43] whilst antibodies against PvDBP RII epitopes have been found to inhibit both attachment of PvDBP RII to erythrocytes  and in vitro invasion of erythrocytes . With regards to immune selection pressure on PvDBP, a study of 100 Sri Lankan isolates found no evidence of significant selection pressure on this region using Tajima’s D, dN/dS or Fu and Li’s D and F statistics . Another study examining genetic diversity of PvDBP RII across multiple populations showed a significantly positive value of dN/dS in this region, suggesting that this region may be under immune selection pressure .
In this study, selection pressures on PvAMA1 and PvDBP Region II were examined in the context of protein structure, using a newly developed tool called BioStructMap . BioStructMap enables the application of a 3D sliding window over a protein structure. This allows incorporation of protein structural information into tests such as Tajima’s D or nucleotide diversity that are traditionally performed as a linear 2D sliding window over a protein or nucleotide sequence. A previous study identified a discontinuous region of PfAMA1 bordering DII and DIII that had a strong signature of balancing selection when considering spatially derived Tajima’s D . Given that other studies have identified DI of PvAMA1, rather than DII or DIII, as being under balancing selection, it was considered that incorporation of protein structural information might yield additional insights into other regions under immune selection pressure. Genomic sequences from a number of populations were analysed, and spatially-derived nucleotide diversity and Tajima’s D were examined using protein structural information for PvAMA1 and PvDBP Region II. Structural patterns of nucleotide diversity were similar across all populations examined, with Domain I of PvAMA1 having the highest nucleotide diversity and displaying significant signatures of balancing selection (Tajima’s D > 0). Nucleotide diversity for PvDBP was highest bordering the dimerization and DARC-binding interface, although there was less evidence of immune selection pressure on this antigen.
Reference sequences for Sal-1 PvDBP (PVX_110810) and PvAMA1 (PVX_092275) were obtained from PlasmoDB, v34 (http://www.plasmodb.org) . Genomic sequences from field isolates were extracted from GenBank, restricted to P. vivax, and with the condition that sequences had to cover > 95% of the structured domains examined here (i.e. partial fragments from these regions were excluded). Sequences from non-human hosts were excluded, as were sequences without a known geographic location described either in the associated literature or clearly annotated in the sequence record. Single isolate populations were also excluded. For PvDBP, bases 766–1659 were used, while for PvAMA1 bases 121–1422 were used.
Sequence polymorphism analysis
Sequences from all isolates were aligned using MUSCLE v3.8.31 , and alignments manually adjusted to minimize gaps. Maximum likelihood phylogenetic trees were constructed for AMA-1 and DBP sequences using IQ-TREE v184.108.40.206 . The ultrafast bootstrap estimation (UFBoot)  was used with 5000 bootstrap replicates and a best fit model was chosen according to the Bayesian Inference Criterion . Phylogenetic trees were visualized with iTOL .
Incorporation of structural information into sequence analysis
A Python package, BioStructMap, which allows for the application of a 3D sliding window over a protein structure has been previously described . BioStructMap was used to compute spatially derived Tajima’s D and nucleotide diversity (π) values for both PvAMA1 and PvDBP structures, with a radius of 15 Å for each window. BioStructMap was also used to calculate Normalized Shannon Entropy on a per-residue basis.
Protein structural models
Known PvDBP structures (i.e. 4NUU, 3RRC, 4YFS, 4NUV, 5F3J) all contain a number of unresolved residues, and to ensure complete coverage of DBP Region II in this analysis, a template-based model of P. vivax DBP was generated for use with BioStructMap. This model was created using ModPipe , an automated software pipeline that utilizes MODELLER for the generation of comparative protein structure models . The PDB structure 4NUU was used to generate this comparative model. The generated model is accessible via ModBase (https://modbase.compbio.ucsf.edu/; ModPipe model ID f7602e019fac5be4a79c4cca6751b392) .
The PvAMA1 model used has been previously described , and uses a chimeric template to generate a structural model of PvAMA1.
Comparing patterns of selection and diversity between populations
To compare structural patterns of nucleotide diversity and Tajima’s D between all populations considered in this study, Spearman’s rank correlation coefficient was computed for both PvAMA1 and PvDBP residue data between each pair of populations. Residues with missing data in one or both populations (i.e. Tajima’s D was undefined) were excluded from analysis for that pair of populations.
Data analysis, statistics and other software used
The majority of data analysis was performed using the Anaconda distribution of Python 3.5. Plotting was performed with the Python Matplotlib package, version 1.5.1 . Statistical analysis was performed using SciPy . Protein structures were visualized using PyMol .
Population structure of PvAMA1 and PvDBP sequences
Traditional measures of selection pressure and diversity for PvAMA1 and PvDBP
Population genetics parameters for PvAMA1 sequences from various geographic locations
Number of polymorphic sites
π (× 10−3)
Papua New Guinea (Madang)
Papua New Guinea (East Sepik)
Thailand (Tak; 1996)
Thailand (Tak; 2007)
Population genetics parameters for PvDBP sequences from various geographic locations
Number of polymorphic sites
π (× 10−3)
Papua New Guinea
For PvAMA1, nucleotide diversity (π) and mean number of pairwise differences (k) were highest in the three Thai populations examined, with a maximum π of 10.09 × 10−3 in the 2007 Tak Province, Thailand samples. Nucleotide diversity was lowest in the South Korean population (π = 5.90 × 10−3), which may be explained by a recent bottleneck in this population, limiting overall diversity; this is supported by the negative Tajima’s D value observed for PvAMA1 in this South Korean population (D = − 1.20). Tajima’s D is often used to identify regions under balancing selection (Tajima’s D > 0), which can be the result of immune selection pressure, and it is noted that for PvAMA1, only the Venezuelan population had a significantly positive Tajima’s D value (D = 2.12, p < 0.05), as per the confidence limits outlined in Tajima , although these limits may be overly conservative . Most other populations also had a positive Tajima’s D value, such as Papua New Guinea (Madang: D = 1.59; East Sepik: D = 1.09) and Thailand (Chantaburi: D = 1.22), with the sole exception of South Korea (discussed above).
When examining PvDBP, nucleotide diversity was highest in South Korea (π = 10.53 × 10−3), Bangkok, Thailand (π = 10.37 × 10−3) and Myanmar (π = 10.17 × 10−3). Interestingly, the South Korean population had a negative Tajima’s D value that was statistically significant, suggestive of a recent population bottleneck and expansion. This population also had the highest number of polymorphic sites (58 non-synonymous). This high variability is probably due to 2 or 3 divergent isolates (Fig. 1b), and haplotype diversity is relatively low for this population (Hd = 0.85). Tajima’s D values for PvDBP were mostly close to zero, with the highest value observed for a Colombian population (D = 1.72).
These observations agree with previous studies that have generally observed greater signatures of immune selection pressure on PvAMA1 as compared to PvDBP .
Diversity and selection on the PvAMA1 structure
Diversity and selection on the PvDBP structure
Amino acid mutations within binding interfaces
While binding interfaces were generally observed to have low nucleotide diversity, there were some polymorphic residues located within these binding interfaces for both PvAMA1 and PvDBP (Additional file 8). There were a total of 9 polymorphic residues within the PvAMA1:RON2 interface, although most of these were present at very low frequencies within the global population. The three most frequent polymorphisms within the PvAMA1:RON2 interface fall within a short stretch from residues 130–133 (N132D, N130K, D133N). The most frequent of these polymorphisms (N132D) is a relatively conservative change from asparagine to aspartic acid, and is the only one of the three residues whose side chains are directly involved in hydrogen bonds with RON2; the side chains from both N130 and D133 are only involved in intramolecular hydrogen bonds. It has previously been shown that a PvRON2 peptide binds well to PfAMA1 . With this in mind, it is noted that the region corresponding to N130–D133 in PvAMA1 is not conserved in PfAMA1, and that this region is also not involved in the binding between PvRON2 and PfAMA1 (PDB structure 5NQF) . These observations, coupled with the observation that a number of polymorphisms fall within this region, suggests that these residues are not a major determinant for binding to PvRON2 and are hence amenable to polymorphic variation.
There were 7 polymorphic residues observed within the PvDBP dimerization and DARC binding interface, although again most of these were present at low frequencies. The highest frequency polymorphism within this region is R263S, with a minor allele frequency of 20%. In the PDB structure 4NUU, the backbone of R263 is involved in hydrogen bonds to DARC, while the side chain is involved in a single intra-molecular hydrogen bond. This residue is also part of a loop which is disordered when not bound to DARC (loop 254–267) , suggesting a level of structural plasticity. The other polymorphic residue of note is T359R, which was identified as a major contact within the DARC binding interface by Batchelor et al. , although this polymorphism was observed at a relatively low minor allele frequency of 4.1% across all populations. However, this residue is not conserved between PvDBP and PkDBPα; PkDBPα has an arginine at this position and is also able to bind DARC, supporting the viability of the T359R mutation with regards to maintaining DARC binding ability.
Comparison between 3D and linear sliding window approaches
While a number of previous studies have utilized a linear sliding window approach when calculating Tajima’s D or nucleotide diversity [4–7, 65], this study is one of the first to utilize a spatially derived Tajima’s D. As such, the Tajima’s D values obtained using a 3D sliding window were compared with those obtained using a conventional linear sliding window, using a window size of 102 base pairs and a step size of 3 base pairs for the linear sliding window (Additional files 9, 10). In general there was broad correspondence between the two approaches for both PvDBP RII and PvAMA1, although there were some additional regions within PvAMA1 with Tajima’s D values above the threshold for significance when using a 3D sliding window approach that were not identified with a linear sliding window. However, these regions were generally structurally contiguous with other stretches of sequence that were above the threshold for significance when using a linear sliding window.
Comparison of spatially derived nucleotide diversity and Tajima’s D between populations
In this study, patterns of nucleotide diversity and selection were examined over the protein structures for the P. vivax antigens PvAMA1 and PvDBP. A number of major observations stand out from this work. Firstly, patterns of diversity on both PvAMA1 and PvDBP were remarkably similar across multiple geographic populations, despite phylogenetic trees for both PvAMA1 and PvDBP sequences suggesting a level of clustering according to geographic location. The only exception to this for PvAMA1 was the South Korean population which displayed evidence of a recent bottleneck and expansion. This similarity in patterns of diversity is important when trying to extend conclusions made from studies from single geographic locations to a worldwide population, and these observations suggest a universality with regards to major epitopes on these antigens. It is also interesting to note that highly polymorphic residues for both PvAMA1 and PvDBP tended to fall within regions surrounding, but not a part of, ligand binding interfaces. For AMA1, RON2 binds in a hydrophobic cleft in DI, and polymorphic residues fall on one side of this hydrophobic cleft, but generally not within residues that make contact with the RON2 peptide. Similarly with PvDBP, contact with DARC occurs primarily via subdomain 2, and the most polymorphic regions were near the DARC binding and dimerization interface. In contrast, the residues directly involved in the PvDBP binding interface were highly conserved. These results are suggestive of two things. Firstly, the residues that make up the key binding interfaces in these two antigens have limited capacity for polymorphic variation due to functional constraints, as has been previously suggested by other studies [39, 63]. This makes them attractive vaccine targets, as potential epitopes within these binding sites would have very limited antigenic diversity, and are also less likely to undergo extensive mutations to evade immune responses. Secondly, the high degree of polymorphism around these interfaces suggests that antibody responses that target these polymorphic sites are capable of inhibiting parasite invasion. This inhibition is likely the result of steric hindrance preventing receptor binding and/or dimerization. Future efforts could involve epitope focusing techniques [66, 67] to direct antibody responses to these key conserved interfaces.
For PvAMA1, we observed balancing selection primarily on DI in all populations examined. Additionally, while the 3D sliding window approach highlighted additional residues under balancing selection as compared to a linear sliding window approach, nearly all of these regions fell within DI. This agrees with a number of other studies in which PvAMA DI is the only domain found to be under significant balancing selection [4, 8–10]. This is in contrast to selection pressures observed on PfAMA1, in which both DI and DIII have been observed to be under balancing selection [5–9, 68]. Previous work of ours has applied spatially-derived Tajima’s D calculations over a PfAMA1 structure and identified strong balancing selection in a region bordering DII and DIII , lending further evidence to DIII being under immune selection pressure in P. falciparum but not P. vivax AMA1. The biological reasons for such a difference are unclear, as AMA1 has a conserved role between Plasmodium species, although it is possible that DIII of PvAMA1 is less immunogenic than the corresponding PfAMA1 domain due to structural or sequence differences between the two antigens.
Although individuals in malaria endemic areas develop antibodies to PvAMA1 [69–71], there are no comprehensive studies on how these antibodies interact with the different domains of PvAMA1. However, AMA1 is functionally conserved across Plasmodium species, and there is evidence that PvAMA1 is functionally equivalent in a P. falciparum transgenic line in which PfAMA1 is replaced by PvAMA1 . As such, comparisons can be drawn with antibody studies on PfAMA1. Dutta et al.  generated a panel of monoclonal antibodies (mAbs) to PfAMA1 and observed that their strain specificity and functional activity was determined by the diversity of the epitope sequence. The limited diversity in PvAMA1 DIII observed in this current study aligns with the observations that mAbs to PfAMA1 DIII were the most strain transcending. Similarly, we observed that the polymorphic face of PvAMA1 DI had the highest diversity, in line with the observation by Dutta et al. that mAbs that bound the polymorphic face of PfAMA1 DI were strain specific, compared to the others that bound the conserved face . Importantly, mAbs that bound to the conserved face of PfAMA1 still showed strong growth inhibitory activity suggesting that epitopes on the conserved face can be targets of neutralizing antibodies, despite being under less immune pressure than immunodominant polymorphic regions.
For PvDBP, one of the regions of high diversity across all populations was a previously identified epitope within subdomain 2 termed the DEK epitope . Others have observed that the DEK epitope is highly polymorphic and immunodominant . However, due to the polymorphic nature of this epitope, cross-strain specificity is a concern when creating a PvDBP RII-based vaccine. Recent work has characterized the location of several conserved epitopes within PvDBP RII that are the target of inhibitory mAbs 2D10, 2H2 and 2C6 [73, 74], and all of these conserved epitopes fall within subdomain III, which had the lowest overall nucleotide diversity in this current study. This highlights some of the limitations of using population-level genomic data for identification of functionally important targets of antibody responses—the possibility that conserved regions may contain potential epitopes that are the targets of inhibitory antibodies cannot be excluded.
Several attempts have been made to divert immune responses away from these highly polymorphic regions of PvDBP RII and towards conserved epitopes. One attempt to focus away from polymorphic regions involved mutating residues in the DEK epitope to reduce its immunogenicity, and these DEKnull mutants induce antibodies that bind PvDBP and can inhibit the interaction with Duffy Binding Ligand . More recently, further epitope focusing techniques have been employed with PvDBP RII, with one strategy involving mutation of all polymorphic residues to alanine, threonine or serine residues . This ‘DEKnull-2’ recombinant PvDBP RII construct was shown to elicit broadly neutralizing antibodies following mouse immunization, with some naturally exposed individuals also shown to recognize the conserved epitopes on this construct . Other recent efforts towards the development of a PvDBP RII vaccine include a Phase 1a trial of a prime-boost viral-vectored vaccine that demonstrated both safety and immunogenicity, with cross-strain inhibition demonstrated for the single heterologous strain tested .
Given the success using epitope focusing techniques for PvDBP RII, we suggest that such an approach could be applied to AMA1 to focus antibody responses towards conserved epitopes on the silent face of DI or within DIII. This would involve mutating major polymorphic residues to reduce the immunodominance of epitopes within the polymorphic face of DI; the most polymorphic residues identified in this study (Additional file 11) could serve as starting point for this work. Alternatively, it has been shown for PfAMA1 that immunization with multiple heterologous strains of PfAMA1 is capable of inducing strain transcending antibody responses , and this approach could also be applied to PvAMA1. Other epitope-focusing approaches also exist, including the use of small protein scaffolds to mimic native epitopes [66, 67], although these might be challenging given the discontinuous nature of many potential epitopes within AMA1. The approaches used in this work could also be applied to other antigens such as PvRBP2b, which has recently been identified as a ligand for reticulocyte invasion via binding to transferrin receptor 1 (TfR1) .
In this work, signatures of diversity and selection were identified on PvAMA1 and PvDBP, and it was shown that the regions of high diversity and balancing selection on the protein structure are remarkably similar across a number of populations. This suggests that dominant epitopes are the same across multiple human populations, which has positive implications for the development of a universal P. vivax vaccine. Furthermore, polymorphisms were observed to cluster around binding interfaces on both PvDBP and PvAMA1, suggesting a level of immune pressure on residues surrounding these key binding interfaces. Large regions with very low diversity were also identified for both antigens, and it is suggested that these areas may also be useful targets to focus on for further vaccine development, given previous evidence of functional antibody responses against these conserved regions.
AJG designed the study, performed the research, analysed the results and wrote the manuscript; PAR contributed the modelled PvAMA1 structure; VI, JSR and PAR provided critical interpretation of the data. All authors revised, commented and read the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
Funding was provided by the National Health and Medical Research Council (NHMRC) of Australia including an Early Career Research Fellowship (APP1037722) and Project Grant (APP1125788) to JSR; and an Australian Postgraduate Award to support AJG through Monash University. Burnet Institute received funding from the NHMRC Independent Research Institutes Infrastructure Support Scheme, and the Victorian State Government Operational Infrastructure Support Scheme.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis 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.
- WHO. World Malaria Report 2017. Geneva: World Health Organization; 2017. http://www.who.int/malaria/publications/world-malaria-report-2017/en/ Accessed 21 Dec 2017.
- Mobegi VA, Duffy CW, Amambua-Ngwa A, Loua KM, Laman E, Nwakanma DC, et al. Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity. Mol Biol Evol. 2014;31:1490–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Amambua-Ngwa A, Tetteh KKA, Manske M, Gomez-Escobar N, Stewart LB, Elizabeth Deerhake M, et al. Population genomic scan for candidate signatures of balancing selection to guide antigen characterization in malaria parasites. PLoS Genet. 2012;8:e1002992.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Osier FHA, Weedall GD, Verra F, Murungi L, Tetteh KKA, Bull P, et al. Allelic diversity and naturally acquired allele-specific antibody responses to Plasmodium falciparum apical membrane antigen 1 in Kenya. Infect Immun. 2010;78:4625–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Polley SD, Chokejindachai W, Conway DJ. Allele frequency-based analyses robustly map sequence sites under balancing selection in a malaria vaccine candidate antigen. Genetics. 2003;165:555–61.PubMedPubMed CentralGoogle Scholar
- Polley SD, Conway DJ. Strong diversifying selection on domains of the Plasmodium falciparum apical membrane antigen 1 gene. Genetics. 2001;158:1505–12.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Guy AJ, Irani V, Beeson JG, Webb B, Sali A, Richards JS, et al. Proteome-wide mapping of immune features onto Plasmodium protein three-dimensional structures. Sci Rep. 2018;8:4355.View ArticlePubMedPubMed CentralGoogle Scholar
- Peterson MG, Marshall VM, Smythe JA, Crewther PE, Lew A, Silva A, et al. Integral membrane protein located in the apical complex of Plasmodium falciparum. Mol Cell Biol. 1989;9:3151–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Healer J, Crawford S, Ralph S, McFadden G, Cowman AF. Independent translocation of two micronemal proteins in developing Plasmodium falciparum merozoites. Infect Immun. 2002;70:5751–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lamarque M, Besteiro S, Papoin J, Roques M, Normand BV-L, 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Richard D, MacRaild CA, Riglar DT, Chan J-A, Foley M, Baum J, et al. Interaction between Plasmodium falciparum apical membrane antigen 1 and the rhoptry neck protein complex defines a key step in the erythrocyte invasion process of malaria parasites. J Biol Chem. 2010;285:14815–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Srinivasan P, Beatty WL, Diouf A, Herrera R, Ambroggio X, Moch JK, et al. Binding of Plasmodium merozoite proteins RON2 and AMA1 triggers commitment to invasion. Proc Natl Acad Sci USA. 2011;108:13275–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Mugyenyi CK, Elliott SR, McCallum FJ, Anders RF, Marsh K, Beeson JG. Antibodies to polymorphic invasion-inhibitory and non-inhibitory epitopes of Plasmodium falciparum apical membrane antigen 1 in human malaria. PLoS ONE. 2013;8:e68304.View ArticlePubMedPubMed CentralGoogle Scholar
- Hodder AN, Crewther PE, Anders RF. Specificity of the protective antibody response to apical membrane antigen 1. Infect Immun. 2001;69:3286–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Fowkes FJI, Richards JS, Simpson JA, Beeson JG. The relationship between anti-merozoite antibodies and incidence of Plasmodium falciparum malaria: a systematic review and meta-analysis. PLoS Med. 2010;7:e1000218.View ArticlePubMedPubMed CentralGoogle Scholar
- Pizarro JC, Vulliez-Le Normand B, 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:408–11.View ArticlePubMedGoogle Scholar
- Vulliez-Le Normand B, Saul FA, Hoos S, Faber BW, Bentley GA. Cross-reactivity between apical membrane antigen 1 and rhoptry neck protein 2 in P. vivax and P. falciparum: a structural and binding study. PLoS ONE. 2017;12:e0183198.View ArticlePubMedPubMed CentralGoogle Scholar
- Dutta S, Dlugosz LS, Clayton JW, Pool CD, Haynes JD, Gasser RA 3rd, et al. Alanine mutagenesis of the primary antigenic escape residue cluster, c1, of apical membrane antigen 1. Infect Immun. 2010;78:661–71.View ArticlePubMedGoogle Scholar
- Dutta S, Dlugosz LS, Drew DR, Ge X, Ge X, Ababacar D, et al. Overcoming antigenic diversity by enhancing the immunogenicity of conserved epitopes on the malaria vaccine candidate apical membrane antigen-1. PLoS Pathog. 2013;9:e1003840.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Kang J-M, Lee J, Cho P-Y, Moon S-U, Ju H-L, 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Moon S-U, Na B-K, Kang J-M, Kim J-Y, Cho S-H, Park Y-K, et al. Genetic polymorphism and effect of natural selection at domain I of apical membrane antigen-1 (AMA-1) in Plasmodium vivax isolates from Myanmar. Acta Trop. 2010;114:71–5.View ArticlePubMedGoogle Scholar
- Wertheimer SP, Barnwell JW. Plasmodium vivax interaction with the human Duffy blood group glycoprotein: identification of a parasite receptor-like protein. Exp Parasitol. 1989;69:340–50.View ArticlePubMedGoogle Scholar
- Horuk R, Chitnis CE, Darbonne WC, Colby TJ, Rybicki A, Hadley TJ, et al. A receptor for the malarial parasite Plasmodium vivax: the erythrocyte chemokine receptor. Science. 1993;261:1182–4.View ArticlePubMedGoogle Scholar
- Cavasini CE, de Mattos LC, Couto AAD, Bonini-Domingos CR, Valencia SH, de Neiras WC, et al. Plasmodium vivax infection among Duffy antigen-negative individuals from the Brazilian Amazon region: an exception? Trans R Soc Trop Med Hyg. 2007;101:1042–4.View ArticlePubMedGoogle Scholar
- Ménard D, Barnadas C, Bouchier C, Henry-Halldin C, Gray LR, Ratsimbasoa A, et al. Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc Natl Acad Sci USA. 2010;107:5967–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Ntumngia FB, Thomson-Luque R, de Torres LM, Gunalan K, Carvalho LH, Adams JH. A novel erythrocyte binding protein of Plasmodium vivax suggests an alternate invasion pathway into Duffy-positive reticulocytes. MBio. 2016;7:e011261.View ArticleGoogle Scholar
- Ntumngia FB, Adams JH. Design and immunogenicity of a novel synthetic antigen based on the ligand domain of the Plasmodium vivax Duffy binding protein. Clin Vaccine Immunol. 2012;19:30–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Adams JH, Blair PL, Kaneko O, Peterson DS. An expanding EBL family of Plasmodium falciparum. Trends Parasitol. 2001;17:297–9.View ArticlePubMedGoogle Scholar
- Carlton JM, Adams JH, Silva JC, Bidwell SL, Lorenzi H, Caler E, et al. Comparative genomics of the neglected human malaria parasite Plasmodium vivax. Nature. 2008;455:757–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Sim BK, Chitnis CE, Wasniowska K, Hadley TJ, Miller LH. Receptor and ligand domains for invasion of erythrocytes by Plasmodium falciparum. Science. 1994;264:1941–4.View ArticlePubMedGoogle Scholar
- Maier AG, Duraisingh MT, Reeder JC, Patel SS, Kazura JW, Zimmerman PA, et al. Plasmodium falciparum erythrocyte invasion through glycophorin C and selection for Gerbich negativity in human populations. Nat Med. 2003;9:87–92.View ArticlePubMedGoogle Scholar
- Mayer DCG, Cofie J, Jiang L, Hartl DL, Tracy E, Kabat J, et al. Glycophorin B is the erythrocyte receptor of Plasmodium falciparum erythrocyte-binding ligand, EBL-1. Proc Natl Acad Sci USA. 2009;106:5348–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Batchelor JD, Malpede BM, Omattage NS, DeKoster GT, Henzler-Wildman KA, Tolia NH. Red blood cell invasion by Plasmodium vivax: structural basis for DBP engagement of DARC. PLoS Pathog. 2014;10:e1003869.View ArticlePubMedPubMed CentralGoogle Scholar
- Batchelor JD, Zahm JA, Tolia NH. Dimerization of Plasmodium vivax DBP is induced upon receptor binding and drives recognition of DARC. Nat Struct Mol Biol. 2011;18:908–14.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Cole-Tobian JL, Michon P, Biasor M, Richards JS, Beeson JG, Mueller I, et al. Strain-specific Duffy binding protein antibodies correlate with protection against infection with homologous compared to heterologous Plasmodium vivax strains in Papua New Guinean children. Infect Immun. 2009;77:4009–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Xainli J, Baisor M, Kastens W, Bockarie M, Adams JH, King CL. Age-dependent cellular immune responses to Plasmodium vivax Duffy binding protein in humans. J Immunol. 2002;169:3200–7.View ArticlePubMedGoogle Scholar
- Chootong P, Ntumngia FB, VanBuskirk KM, Xainli J, Cole-Tobian JL, Campbell CO, et al. Mapping epitopes of the Plasmodium vivax Duffy binding protein with naturally acquired inhibitory antibodies. Infect Immun. 2010;78:1089–95.View ArticlePubMedGoogle Scholar
- Grimberg BT, Udomsangpetch R, Xainli J, McHenry A, Panichakul T, Sattabongkot J, et al. Plasmodium vivax invasion of human erythrocytes inhibited by antibodies directed against the Duffy binding protein. PLoS Med. 2007;4:e337.View ArticlePubMedPubMed CentralGoogle Scholar
- Aurrecoechea C, Brestelli J, Brunk BP, Dommer J, Fischer S, Gajria B, et al. PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2009;37:D539–43.View ArticlePubMedGoogle Scholar
- Sukumaran J, Holder MT. DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010;26:1569–71.View ArticlePubMedGoogle Scholar
- Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1948;27:379–423.View ArticleGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.View ArticlePubMedGoogle Scholar
- Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Le SV. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35:518–22.View ArticlePubMedGoogle Scholar
- Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Eswar N, John B, Mirkovic N, Fiser A, Ilyin VA, Pieper U, et al. Tools for comparative protein structure modeling and analysis. Nucleic Acids Res. 2003;31:3375–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Sali A, Blundell TL. Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993;234:779–815.View ArticlePubMedGoogle Scholar
- Pieper U, Webb BM, Barkan DT, Schneidman-Duhovny D, Schlessinger A, Braberg H, et al. ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 2011;39:D465–74.View ArticlePubMedGoogle Scholar
- Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9:90–5.View ArticleGoogle Scholar
- van der Walt S, Colbert SC, Varoquaux G. The NumPy array: a structure for efficient numerical computation. Comput Sci Eng. 2011;13:22–30.View ArticleGoogle Scholar
- Schrödinger LLC. The PyMOL molecular graphics system, version 1.8. 2015.Google Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMedPubMed CentralGoogle Scholar
- Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133:693–709.PubMedPubMed CentralGoogle Scholar
- Faber BW, Kadir KA, Rodriguez-Garcia R, Remarque EJ, Saul FA, Normand BV-L, et al. Low levels of polymorphisms and no evidence for diversifying selection on the Plasmodium knowlesi apical membrane antigen 1 Gene. PLoS ONE. 2015;10:e0124400.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Takala SL, Coulibaly D, Thera MA, Batchelor AH, Cummings MP, Escalante AA, et al. Extreme polymorphism in a vaccine antigen and risk of clinical malaria: implications for vaccine development. Sci Transl Med. 2009;1:2ra5.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 Parasitol. 2001;113:279–87.View ArticlePubMedGoogle Scholar
- Oscherwitz J. The promise and challenge of epitope-focused vaccines. Hum Vaccin Immunother. 2016;12:2113–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Correia BE, Bates JT, Loomis RJ, Baneyx G, Carrico C, Jardine JG, et al. Proof of principle for epitope-focused vaccine design. Nature. 2014;507:201.View ArticlePubMedPubMed CentralGoogle Scholar
- Mehrizi AA, Sepehri M, Karimi F, Djadid ND, Zakeri S. Population genetics, sequence diversity and selection in the gene encoding the Plasmodium falciparum apical membrane antigen 1 in clinical isolates from the south-east of Iran. Infect Genet Evol. 2013;17:51–61.View ArticlePubMedGoogle Scholar
- Fowkes FJI, McGready R, Cross NJ, Hommel M, Simpson JA, Elliott SR, et al. New insights into acquisition, boosting, and longevity of immunity to malaria in pregnant women. J Infect Dis. 2012;206:1612–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Cutts JC, Powell R, Agius PA, Beeson JG, Simpson JA, Fowkes FJI. Immunological markers of Plasmodium vivax exposure and immunity: a systematic review and meta-analysis. BMC Med. 2014;12:150.View ArticlePubMedPubMed CentralGoogle Scholar
- McLean ARD, Boel M, McGready R, Ataide R, Drew D, Tsuboi T, et al. Antibody responses to Plasmodium falciparum and Plasmodium vivax and prospective risk of Plasmodium spp. infection postpartum. Am J Trop Med Hyg. 2017;96:1197–204.View ArticlePubMedPubMed CentralGoogle Scholar
- Drew DR, Sanders PR, Weiss G, Gilson PR, Crabb BS, Beeson JG. Functional conservation of the AMA1 host-cell invasion ligand between P. falciparum and P. vivax: a novel platform to accelerate vaccine and drug development. J Infect Dis. 2018;217:498–507.View ArticlePubMedGoogle Scholar
- Chen E, Salinas ND, Huang Y, Ntumngia F, Plasencia MD, Gross ML, et al. Broadly neutralizing epitopes in the Plasmodium vivax vaccine candidate Duffy Binding Protein. Proc Natl Acad Sci USA. 2016;113:6277–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Ntumngia FB, Schloegel J, Barnes SJ, McHenry AM, Singh S, King CL, et al. Conserved and variant epitopes of Plasmodium vivax Duffy binding protein as targets of inhibitory monoclonal antibodies. Infect Immun. 2012;80:1203–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ntumngia FB, Pires CV, Barnes SJ, George MT, Thomson-Luque R, Kano FS, et al. An engineered vaccine of the Plasmodium vivax Duffy binding protein enhances induction of broadly neutralizing antibodies. Sci Rep. 2017;7:13779.View ArticlePubMedPubMed CentralGoogle Scholar
- Payne RO, Silk SE, Elias SC, Milne KH, Rawlinson TA, Llewellyn D, et al. Human vaccination against Plasmodium vivax Duffy-binding protein induces strain-transcending antibodies. JCI Insight. 2017;2:93683 (Epub ahead of print).View ArticlePubMedGoogle Scholar
- Gruszczyk J, Kanjee U, Chan L-J, Menant S, Malleret B, Lim NTY, et al. Transferrin receptor 1 is a reticulocyte-specific receptor for Plasmodium vivax. Science. 2018;359:48–55.View ArticlePubMedGoogle Scholar