Global Population Structure of the Genes Encoding the Malaria Vaccine Candidate, Plasmodium vivax Apical Membrane Antigen 1 (PvAMA1)

Potent AMA1-specific human monoclonal antibody against P. vivax Pre-erythrocytic and Blood Stages

New therapeutics are necessary for preventing Plasmodium vivax malaria due to easy transmissibility and dormancy in the liver that increases the clinical burden due to recurrent relapse. We isolated 12 Pv Apical Membrane Antigen 1 (PvAMA1) specific human monoclonal antibodies from Peripheral Blood Mononuclear Cells of a Pv exposed individual. PvAMA1 is essential for sporozoite and merozoite invasion, making it a unique therapeutic target. HumAb 826827 blocked the invasion of human erythrocytes using Pv clinical isolates and inhibited sporozoite invasion of human hepatocytes in vitro (IC50 of 0.3 to 3.7 ug/mL). It also significantly reduced liver infection of chimeric FRG humHep mice in vivo. The crystal structure of rPvAMA1 bound to 826827 shows that 826827 partially occupies the highly conserved hydrophobic groove in PvAMA1 that binds its known receptor, RON2. We have isolated a potent humAb that is isolate transcendent, blocks both pre erythrocytic and blood stage infection, and could be a new therapy for Pv.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Plasmodium GPI-anchored micronemal antigen is essential for parasite transmission through the mosquito host

Plasmodium parasites, the eukaryotic pathogens that cause malaria, feature three distinct invasive forms tailored to the host environment they must navigate and invade for life cycle progression. One conserved feature of these invasive forms is the micronemes, apically oriented secretory organelles involved in egress, motility, adhesion, and invasion. Here we investigate the role of GPI-anchored micronemal antigen (GAMA), which shows a micronemal localization in all zoite forms of the rodent-infecting species Plasmodium berghei. ∆GAMA parasites are severely defective for invasion of the mosquito midgut. Once formed, oocysts develop normally, however, sporozoites are unable to egress and exhibit defective motility. Epitope-tagging of GAMA revealed tight temporal expression late during sporogony and showed that GAMA is shed during sporozoite gliding motility in a similar manner to circumsporozoite protein. Complementation of P. berghei knockout parasites with full-length P. falciparum GAMA partially restored infectivity to mosquitoes, indicating conservation of function across Plasmodium species. A suite of parasites with GAMA expressed under the promoters of CTRP, CAP380, and TRAP, further confirmed the involvement of GAMA in midgut infection, motility, and vertebrate infection. These data show GAMA's involvement in sporozoite motility, egress, and invasion, implicating GAMA as a regulator of microneme function.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Plasmodium vivax: the potential obstacles it presents to malaria elimination and eradication

Initiatives to eradicate malaria have a good impact on P. falciparum malaria worldwide. P. vivax, however, still presents significant difficulties. This is due to its unique biological traits, which, in comparison to P. falciparum, pose serious challenges for malaria elimination approaches. P. vivax's numerous distinctive characteristics and its ability to live for weeks to years in liver cells in its hypnozoite form, which may elude the human immune system and blood-stage therapy and offer protection during mosquito-free seasons. Many malaria patients are not fully treated because of contraindications to primaquine use in pregnant and nursing women and are still vulnerable to P. vivax relapses, although there are medications that could radical cure P. vivax. Additionally, due to CYP2D6's highly variable genetic polymorphism, the pharmacokinetics of primaquine may be impacted. Due to their inability to metabolize PQ, some CYP2D6 polymorphism alleles can cause patients to not respond to treatment. Tafenoquine offers a radical treatment in a single dose that overcomes the potentially serious problem of poor adherence to daily primaquine. Despite this benefit, hemolysis of the early erythrocytes continues in individuals with G6PD deficiency until all susceptible cells have been eliminated. Field techniques such as microscopy or rapid diagnostic tests (RDTs) miss the large number of submicroscopic and/or asymptomatic infections brought on by reticulocyte tropism and the low parasitemia levels that accompany it. Moreover, P. vivax gametocytes grow more quickly and are much more prevalent in the bloodstream. P. vivax populations also have a great deal of genetic variation throughout their genome, which ensures evolutionary fitness and boosts adaptation potential. Furthermore, P. vivax fully develops in the mosquito faster than P. falciparum. These characteristics contribute to parasite reservoirs in the human population and facilitate faster transmission. Overall, no genuine chance of eradication is predicted in the next few years unless new tools for lowering malaria transmission are developed (i.e., malaria elimination and eradication). The challenging characteristics of P. vivax that impede the elimination and eradication of malaria are thus discussed in this article.

Population genetic analyses inferred a limited genetic diversity across the pvama-1 DI domain among Plasmodium vivax isolates from Khyber Pakhtunkhwa regions of Pakistan

Background

Plasmodium vivax apical membrane antigen-1 (pvama-1) is an important vaccine candidate against Malaria. The genetic composition assessment of pvama-1 from wide-range geography is vital to plan the antigen based vaccine designing against Malaria.

Methods

The blood samples were collected from 84 P. vivax positive malaria patients from different districts of Khyber Pakhtunkhwa (KP) province of Pakistan. The highly polymorphic and immunogenic domain-I (DI) region of pvama-1 was PCR amplified and DNA sequenced. The QC based sequences raw data filtration was done using DNASTAR package. The downstream population genetic analyses were performed using MEGA4, DnaSP, Arlequin v3.5 and Network.5 resources.

Results

The analyses unveiled total 57 haplotypes of pvama-1 (DI) in KP samples with majorly prevalent H-14 and H-5 haplotypes. Pairwise comparative population genetics analyses identified limited to moderate genetic distinctions among the samples collected from different districts of KP, Pakistan. In context of worldwide available data, the KP samples depicted major genetic differentiation against the Korean samples with Fst = 0.40915 (P-value = 0.0001), while least distinction was observed against Indian and Iranian samples. The statistically significant negative values of Fu and Li’s D* and F* tests indicate the evidence of population expansion and directional positive selection signature. The slow LD decay across the nucleotide distance in KP isolates indicates low nucleotide diversity. In context of reference pvama-1 sequence, the KP samples were identified to have 09 novel non-synonymous single nucleotide polymorphisms (nsSNPs), including several trimorphic and tetramorphic substitutions. Few of these nsSNPs are mapped within the B-cell predicted epitopic motifs of the pvama-1, and possibly modulate the immune response mechanism.

Conclusion

Low genetic differentiation was observed across the pvama-1 DI among the P. vivax isolates acquired from widespread regions of KP province of Pakistan. The information may implicate in future vaccine designing strategies based on antigenic features of pvama-1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07798-1.

Genetic diversity of Plasmodium vivax reticulocyte binding protein 2b in global parasite populations

Background

Plasmodium vivax reticulocyte binding protein 2b (PvRBP2b) plays a critical role in parasite invasion of reticulocytes by binding the transferrin receptor 1. PvRBP2b is a vaccine candidate based on the negative correlation between antibody titers against PvRBP2b recombinant proteins and parasitemia and risk of vivax malaria. The aim of this study was to analyze the genetic diversity of the PvRBP2b gene in the global P. vivax populations.

Methods

Near full-length PvRBP2b nucleotide sequences (190–8349 bp) were obtained from 88 P. vivax isolates collected from the China–Myanmar border (n = 44) and Thailand (n = 44). An additional 224 PvRBP2b sequences were retrieved from genome sequences from parasite populations worldwide. The genetic diversity, neutral selection, haplotype distribution and genetic differentiation of PvRBP2b were examined.

Results

The genetic diversity of PvRBP2b was distributed unevenly, with peak diversity found in the reticulocyte binding region in the N-terminus. Neutrality analysis suggested that this region is subjected to balancing selection or population bottlenecks. Several amino acid variants were found in all or nearly all P. vivax endemic regions. However, the critical residues responsible for reticulocyte binding were highly conserved. There was substantial population differentiation according to the geographical separation. The distribution of haplotypes in the reticulocyte binding region varied among regions; even the two major haplotypes Hap_6 and Hap_8 were found in only five populations.

Conclusions

Our data show considerable genetic variations of PvRBPb in global parasite populations. The geographic divergence may pose a challenge to PvRBP2b-based vaccine development.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05296-6.

Untangling population structure and genetic diversity of reticulocyte binding protein 2b (PvRBP2b) erythrocytic stage vaccine candidate in worldwide Plasmodium vivax isolates

BACKGROUNDS

Plasmodium vivax is the predominant Plasmodium species distributed extensively in the Americas and Asia-Pacific areas. Encoded protein by Plasmodium vivax Reticulocyte Binding Proteins (PvRBPs) family member are of critical prominence to parasite invasion and have been considered the significant targets in development of malaria vaccine for the blood stage. As high genetic polymorphism of parasites may impede the effectiveness of vaccine development, more research to unraveling genetic polymorphism of pvrbp2b from various geographical regions seems indispensable to map the exact pattern of field isolates.

METHODOLOGY/PRINCIPAL FINDINGS

The aim of this study was to determine the sequences of Iranian pvrbp2b (nt: 502-1896) gene and then, to ascertain polymorphism of pvrbp2b gene, recombination, the level of genetic distances, evaluation of natural selection, and the prediction of B-cell epitopes of Iranian and global P. vivax isolates. Pvrbp2b partial gene was amplified and sequenced from 60 Iranian P. vivax isolates. Iranian pvrbp2b sequences as well as 95 published sequences from five countries were used to evaluate the genetic diversity and neutral evolution signature in worldwide scale. A total of 38 SNPs were identified among 60 Iranian pvrbp2b sequences (32 non-synonymous and 6 synonymous mutations), and 32 amino acid substitutions were observed in 29 positions as compared to Sal-1 sequence. Worldwide sequence analysis showed that 44 amino acid changes had occurred in 37 positions of which seven polymorphic sites had trimorphic mutations while the rest was dimorphic. The overall nucleotide diversity for Iranian isolates was 0.00431 ± 0.00091 while the level of nucleotide diversity was ranged from 0.00337 ± 0.00076 (Peru) to 0.00452 ± 0.00092 (Thailand) in global scale.

CONCLUSIONS/SIGNIFICANCE

Of amino acid substitutions, 12 replacements were located in the B-cell epitopes in which nine polymorphic sites were positioned in N-terminal and three polymorphic sites in predicted B-cell epitopes of C-terminal, signifying both variable and conserved epitopes for vaccine designing. Using the achieved outcome of the current investigation interrogate questions to the selection of conserved regions of pvrbp2b and understanding polymorphism and immune system pressure to pave a way for developing a vaccine based on PvRBP2b candidate antigen.

Genetic Diversity and Selection of Plasmodium vivax Apical Membrane Antigen-1 in China-Myanmar Border of Yunnan Province, China, 2009-2016

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is an important vaccine candidate for vivax malaria. However, antigenic variation within PvAMA-1 is a major obstacle to the design of a global protective malaria vaccine. In this study, we analyzed the genetic polymorphism and selection of the PvAMA-1 gene from 152 P. vivax isolates from imported cases to China, collected in the China–Myanmar border (CMB) area in Yunnan Province (YP) during 2009–2011 (n = 71) and 2014–2016 (n = 81), in comparison with PvAMA-1 gene information from Myanmar (n = 73), collected from public data. The overall nucleotide diversity of the PvAMA-1 gene from the 152 YP isolates was 0.007 with 76 haplotypes identified (Hd = 0.958). Results from the population structure suggested three groups among the YP and Myanmar isolates with optimized clusters value of K = 7. In addition, YP (2014–2016) isolates generally lacked some K components that were commonly found in YP (2009–2011) and Myanmar. Meanwhile, PvAMA-1 domain I is found to be the dominant target of positive diversifying selection and most mutation loci were found in this domain. The mutation frequencies of D107N/A, R112K/T, K120R, E145A, E277K, and R438H in PvAMA-1 were more than 70% in the YP isolates. In conclusion, high genetic diversity and positive selection were found in the PvAMA-1 gene from YP isolates, which are significant findings for the design and development of PvAMA-1-based malaria vaccine.

Genetic Polymorphism and Natural Selection of Apical Membrane Antigen-1 in Plasmodium falciparum Isolates from Vietnam

Apical membrane antigen-1 of Plasmodium falciparum (PfAMA-1) is a leading malaria vaccine candidate antigen. However, the genetic diversity of pfama-1 and associated antigenic variation in global P. falciparum field isolates are major hurdles to the design of an efficacious vaccine formulated with this antigen. Here, we analyzed the genetic structure and the natural selection of pfama-1 in the P. falciparum population of Vietnam. A total of 37 distinct haplotypes were found in 131 P. falciparum Vietnamese isolates. Most amino acid changes detected in Vietnamese pfama-1 were localized in the ectodomain, domains I, II, and III. Overall patterns of major amino acid changes in Vietnamese pfama-1 were similar to those of global pfama-1, but the frequencies of the amino acid changes slightly differed by country. Novel amino acid changes were also identified in Vietnamese pfama-1. Vietnamese pfama-1 revealed relatively lower genetic diversity than currently analyzed pfama-1 in other geographical regions, and suggested a distinct genetic differentiation pattern. Evidence for natural selection was detected in Vietnamese pfama-1, but it showed purifying selection unlike the global pfama-1 analyzed so far. Recombination events were also found in Vietnamese pfama-1. Major amino acid changes that were commonly identified in global pfama-1 were mainly localized to predicted B-cell epitopes, RBC-binding sites, and IUR regions. These results provide important information for understanding the genetic nature of the Vietnamese pfama-1 population, and have significant implications for the design of a vaccine based on PfAMA-1.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

Genetic diversity in two leading Plasmodium vivax malaria vaccine candidates AMA1 and MSP119 at three sites in India

Plasmodium vivax, a major contributor to the malaria burden in India, has the broadest geographic distribution and shows higher genetic diversity than P. falciparum. Here, we investigated the genetic diversity of two leading P. vivax vaccine candidate antigens, at three geographically diverse malaria-endemic regions in India. Pvama1 and Pvmsp119 partial coding sequences were generated from one hundred P. vivax isolates in India (Chennai n = 28, Nadiad n = 50 and Rourkela n = 22) and ~1100 published sequences from Asia, South America, North America, and Oceania regions included. These data were used to assess the genetic diversity and potential for vaccine candidacy of both antigens on a global scale. A total of 44 single nucleotide polymorphism (SNPs) were identified among 100 Indian Pvama1 sequences, including 10 synonymous and 34 nonsynonymous mutations. Nucleotide diversity was higher in Rourkela and Nadiad as compared to Chennai. Nucleotide diversity measures showed a strong balancing selection in Indian and global population for domain I of Pvama1, which suggests that it is a dominant target of the protective immune response. In contrast, the Pvmsp119 region showed highly conserved sequences in India and across the Oceania, South America, North America and Asia, demonstrating low genetic diversity in the global population when compared to Pvama1. Results suggest the possibility of including Pvmsp119 in a multivalent vaccine formulation against P. vivax infections. However, the high genetic diversity seen in Pvama1 would be more challenging for vaccine development.

Antibodies Against the Plasmodium vivax Apical Membrane Antigen 1 From the Belem Strain Share Common Epitopes Among Other Worldwide Variants

Malaria is a human parasitic disease distributed in many tropical countries and caused by various Plasmodium species. Plasmodium vivax has the largest geographical distribution of the Plasmodium species and is predominant in the Americas, including Brazil. Only a small number of P. vivax vaccine formulations have successfully reached clinical trials relative to their P. falciparum counterparts. One of the candidate antigens for a blood-stage P. vivax vaccine is apical membrane antigen 1 (PvAMA-1). Due to the worldwide distribution of Plasmodium parasites, a high degree of variability has been detected in this antigen sequence, representing a considerable challenge to the development of a universal vaccine against malaria. In this study, we evaluated how PvAMA-1 polymorphisms influence vaccine-derived immune responses in P. vivax malaria. To this end, we expressed 9 recombinant protein representatives of different PvAMA-1 allelic variants in the yeast Pichia pastoris: Belem, Chesson I, Sal-1, Indonesia XIX, SK0814, TC103, PNG_05_ESP, PNG_62_MU, and PNG_68_MAS. After protein expression and purification, we evaluated the breadth of the immune responses derived from malaria-exposed individuals from the Amazon region. From 611 serum samples of malaria-exposed individuals, 53.68% of them reacted against the PvAMA-1 Belem through ELISA. Positive samples were further tested against recombinant proteins representing the other PvAMA-1 allelic variants. Whereas Sal-1, Chesson I and SK0814 variants were highly recognized by tested serum samples, Indonesia XIX, TC103, PNG_05_ESP, PNG_62_MU, and PNG_68_MAS were only slightly recognized. Moreover, polyclonal sera derived from C57BL/6 mice immunized with the PvAMA-1 Belem protein predominantly recognized Belem, Sal-1, Chesson I, SK0814, and Indonesia XIX through ELISA. Last, ELISA-based competition assays demonstrated that a previous interaction between anti-Belem polyclonal serum and Sal-1, Chesson I, SK0814, or Indonesia XIX proteins could further inhibit antibody binding to the Belem variant. Our human and mouse data suggest the presence of common epitopes or cross-reactivity between Belem, Sal-1, Chesson I, and SK0814 variants. Although the PvAMA-1 Belem variant induces strain-transcendent antibodies, PvAMA-1 variants from Thailand and Papua New Guinea may need to be included in a universal vaccine formulation to achieve protection against P. vivax malaria.

Plasmodium vivax AMA1: Implications of distinct haplotypes for immune response

In Brazil, Plasmodium vivax infection accounts for around 80% of malaria cases. This infection has a substantial impact on the productivity of the local population as the course of the disease is usually prolonged and the development of acquired immunity in endemic areas takes several years. The recent emergence of drug-resistant strains has intensified research on alternative control methods such as vaccines. There is currently no effective available vaccine against malaria; however, numerous candidates have been studied in the past several years. One of the leading candidates is apical membrane antigen 1 (AMA1). This protein is involved in the invasion of Apicomplexa parasites into host cells, participating in the formation of a moving junction. Understanding how the genetic diversity of an antigen influences the immune response is highly important for vaccine development. In this study, we analyzed the diversity of AMA1 from Brazilian P. vivax isolates and 19 haplotypes of P. vivax were found. Among those sequences, 33 nonsynonymous PvAMA1 amino acid sites were identified, whereas 20 of these sites were determined to be located in predicted B-cell epitopes. Nonsynonymous mutations were evaluated for their influence on the immune recognition of these antigens. Two distinct haplotypes, 5 and 16, were expressed and evaluated for reactivity in individuals from northern Brazil. Both PvAMA1 variants were reactive. Moreover, the IgG antibody response to these two PvAMA1 variants was analyzed in an exposed but noninfected population from a P. vivax endemic area. Interestingly, over 40% of this population had antibodies recognizing both variants. These results have implications for the design of a vaccine based on a polymorphic antigen.

Plasmodium vivax is the most abundant Plasmodium species in Brazil. While this species has been neglected for many years, the recent emergence of drug-resistant strains and the absence of a vaccine intensified the efforts for a better control method. Naturally acquired immune response analysis is a useful tool for understanding the antigenicity of Plasmodium proteins and evaluating the potential of a vaccine candidate. In this study, the genetic variability of one of the leading P. vivax vaccine candidates (PvAMA1) was analyzed. Two distinct variants were expressed and the antibody response was evaluated in infected and noninfected individuals in the Brazilian Amazon. This improved understanding of the magnitude and dynamics of the antibody response will contribute to the knowledge of a vaccine candidate and open new perspectives in vivax malaria vaccine development.

Antibody responses within two leading Plasmodium vivax vaccine candidate antigens in three geographically diverse malaria-endemic regions of India

Background

Identifying highly immunogenic blood stage antigens which can work as target for naturally acquired antibodies in different eco-epidemiological settings is an important step for designing malaria vaccine. Blood stage proteins of Plasmodium vivax, apical membrane antigen-1 (PvAMA-1) and 19 kDa fragment of merozoite surface protein (PvMSP-1₁₉) are such promising vaccine candidate antigens. This study determined the naturally-acquired antibody response to PvAMA-1 and PvMSP-1₁₉ antigens in individuals living in three geographically diverse malaria endemic regions of India.

Methods

A total of 234 blood samples were collected from individuals living in three different eco-epidemiological settings, Chennai, Nadiad, and Rourkela of India. Indirect ELISA was performed to measure human IgG antibodies against recombinant PvAMA-1 and PvMSP-1₁₉ antigens. The difference in seroprevalence and factors associated with antibody responses at each site was statistically analysed.

Results

The overall seroprevalence was 40.6% for PvAMA-1 and 62.4% for PvMSP-1₁₉. Seroprevalence to PvAMA-1 was higher in Chennai (47%) followed by Nadiad (46.7%) and Rourkela (27.6%). For PvMSP-1₁₉, seroprevalence was higher in Chennai (80.3%) as compared to Nadiad (53.3%) and Rourkela (57.9%). Seroprevalence for both the antigens were found to be higher in Chennai where P. vivax is the dominant malaria species. In addition, heterogeneous antibody response was observed for PvAMA-1 and PvMSP-1₁₉ antigens at each of the study sites. Two factors, age and malaria positivity were significantly associated with seropositivity for both the antigens PvAMA-1 and PvMSP-1₁₉.

Conclusion

These data suggest that natural acquired antibody response is higher for PvMSP-1₁₉ antigen as compared to PvAMA-1 antigen in individuals living in three geographically diverse malaria endemic regions in India. PvMSP-1₁₉ appears to be highly immunogenic in Indian population and has great potential as a malaria vaccine candidate. The differences in immune response against vaccine candidate antigens in different endemic settings should be taken into account for development of asexual stage based P. vivax malaria vaccine, which in turn can enhance malaria control efforts.

BioStructMap: a Python tool for integration of protein structure and sequence-based features

Summary

A sliding window analysis over a protein or genomic sequence is commonly performed, and we present a Python tool, BioStructMap, that extends this concept to three-dimensional (3D) space, allowing the application of a 3D sliding window analysis over a protein structure. BioStructMap is easily extensible, allowing the user to apply custom functions to spatially aggregated data. BioStructMap also allows mapping of underlying genomic sequences to protein structures, allowing the user to perform genetic-based analysis over spatially linked codons-this has applications when selection pressures arise at the level of protein structure.

Availability and implementation

The Python BioStructMap package is available at https://github.com/andrewguy/biostructmap and released under the MIT License. An online server implementing standard functionality is available at https://biostructmap.burnet.edu.au.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic diversity and natural selection of transmission-blocking vaccine candidate antigens Pvs25 and Pvs28 in Plasmodium vivax Myanmar isolates

Transmission-blocking vaccines (TBVs) target the sexual stages of malarial parasites to interrupt or reduce the transmission cycle have been one of approaches to control malaria. Pvs25 and Pvs28 are the leading candidate antigens of TBVs against vivax malaria. In this study, genetic diversity and natural selection of the two TBV candidate genes in Plasmodium vivax Myanmar isolates were analyzed. The 62 Myanmar P. vivax isolates showed 9 and 19 different haplotypes for Pvs25 and Pvs28, respectively. The nucleotide diversity of Pvs28 was slightly higher than Pvs25, but not significant. Most amino acid substitutions observed in Myanmar Pvs25 and Pvs28 were concentrated at the EGF-2 and EGF-3 like domains. Major amino acid changes found in Myanmar Pvs25 and Pvs28 were similar to those reported in the global population, but novel amino acid substitutions were also identified. Negative selection was predicted in Myanmar Pvs25, whereas Pvs28 was under positive selection. Comparative analysis of global Pvs25 and Pvs28 suggests a substantial geographical difference between the Asian and American/African Pvs25 and Pvs28. The geographical genetic differentiation and the evidence for natural selection in global Pvs25 and Pvs28 suggest that the functional consequences of the observed polymorphism need to be considered for the development of effective TBVs based on the antigens.

Effect of low complexity regions within the PvMSP3α block II on the tertiary structure of the protein and implications to immune escape mechanisms

Background

Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is a promising vaccine candidate which has shown strong association with immunogenicity and protectiveness. Its use is however complicated by evolutionary plasticity features which enhance immune evasion. Low complexity regions (LCRs) provide plasticity in surface proteins of Plasmodium species, but its implication in vaccine design remain unexplored. Here population genetic, comparative phylogenetic and structural biology analysis was performed on the gene encoding PvMSP3α.

Results

Three LCRs were found in PvMSP3α block II. Both the predicted tertiary structure of the protein and the phylogenetic trees based on this region were influenced by the presence of the LCRs. The LCRs were mainly B cell epitopes within or adjacent. In addition a repeat motif mimicking one of the B cell epitopes was found within the PvMSP3a block II low complexity region. This particular B cell epitope also featured rampant alanine substitutions which might impair antibody binding.

Conclusion

The findings indicate that PvMSP3α block II possesses LCRs which might confer a strong phenotypic plasticity. The phenomenon of phenotypic plasticity and implication of LCRs in malaria immunology in general and vaccine candidate genes in particular merits further exploration.

Electronic supplementary material

The online version of this article (10.1186/s12900-019-0104-0) contains supplementary material, which is available to authorized users.

Limited genetic variability of Cytauxzoon felis apical membrane antigen-1 (ama1) from domestic cats and bobcats

Background

Cytauxzoon felis is a tick-transmitted apicomplexan that causes cytauxzoonosis in domestic cats (Felis catus). Even with intensive care, the mortality rate of acute cytauxzoonosis approaches 40% in domestic cats, while bobcats (Lynx rufus), the natural intermediate host of C. felis, remain clinically asymptomatic. However, multiple reports of domestic cats surviving acute disease without any treatment exist. One hypothesis for survival of these cats is infection with unique C. felis genotypes of lower pathogenicity. Prior studies have identified genetically distinct C. felis isolates containing polymorphisms within internal transcribed spacer regions (ITS) of the rRNA operon. However, these polymorphisms do not correlate with the clinical outcome of cytauxzoonosis, and so additional genetic markers are needed to test this hypothesis. We selected C. felis apical membrane antigen-1 (ama1) as a potential genetic marker of differential pathogenicity. AMA1 is a vaccine candidate for relatives of C. felis within Plasmodium spp.; however its historically high level of genetic polymorphism has resulted in escape from vaccine-induced immunity. While such diversity has hindered vaccine development, the expected polymorphism within the ama1 gene may be useful to evaluate population genetics.

Results

A 677 bp sequence of the C. felis ama1 gene was PCR-amplified from 84 domestic cats and 9 bobcats and demonstrated 99.9% sequence identity across all samples. A single nucleotide polymorphism (SNP) was identified in domestic cats and bobcats with evidence for co-infection with both genotypes identified in two domestic cats. The prevalence of the two genotypes varied with geographical distribution in domestic cats. Nucleotide diversity (π) and haplotype diversity (H) were calculated for C. felis ama1 and ama1 of related apicomplexans to assess genetic diversity. Based on these values (π = 0.00067 and H = 0.457), the diversity of the C. felis ama1 gene region analyzed is considerably lower than what is documented in related apicomplexans.

Conclusions

In surprising contrast to related apicomplexans, our results support that the sequence of the C. felis ama1 gene is highly conserved. While lack of genetic diversity limits utility of C. felis AMA1 as a genetic marker for clinical outcome, it supports further investigation as a vaccine candidate for cytauxzoonosis.

Genetic sequence characterization and naturally acquired immune response to Plasmodium vivax Rhoptry Neck Protein 2 (PvRON2)

Background

The genetic diversity of malaria antigens often results in allele variant-specific immunity, imposing a great challenge to vaccine development. Rhoptry Neck Protein 2 (PvRON2) is a blood-stage antigen that plays a key role during the erythrocyte invasion of Plasmodium vivax. This study investigates the genetic diversity of PvRON2 and the naturally acquired immune response to P. vivax isolates.

Results

Here, the genetic diversity of PvRON2_1828–2080 and the naturally acquired humoral immune response against PvRON2_1828–2080 in infected and non-infected individuals from a vivax malaria endemic area in Brazil was reported. The diversity analysis of PvRON2_1828–2080 revealed that the protein is conserved in isolates in Brazil and worldwide. A total of 18 (19%) patients had IgG antibodies to PvRON2_1828–2080. Additionally, the analysis of the antibody response in individuals who were not acutely infected with malaria, but had been infected with malaria in the past indicated that 32 patients (33%) exhibited an IgG immune response against PvRON2.

Conclusions

PvRON2 was conserved among the studied isolates. The presence of naturally acquired antibodies to this protein in the absence of the disease suggests that PvRON2 induces a long-term antibody response. These results indicate that PvRON2 is a potential malaria vaccine candidate.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2543-7) contains supplementary material, which is available to authorized users.

Structural patterns of selection and diversity for Plasmodium vivax antigens DBP and AMA1

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Proteome-wide mapping of immune features onto Plasmodium protein three-dimensional structures

Humoral immune responses against the malaria parasite are an important component of a protective immune response. Antibodies are often directed towards conformational epitopes, and the native structure of the antigenic region is usually critical for antibody recognition. We examined the structural features of various Plasmodium antigens that may impact on epitope location, by performing a comprehensive analysis of known and modelled structures from P. falciparum. Examining the location of known polymorphisms over all available structures, we observed a strong propensity for polymorphic residues to be exposed on the surface and to occur in particular secondary structure segments such as hydrogen-bonded turns. We also utilised established prediction algorithms for B-cell epitopes and MHC class II binding peptides, examining predicted epitopes in relation to known polymorphic sites within structured regions. Finally, we used the available structures to examine polymorphic hotspots and Tajima's D values using a spatial averaging approach. We identified a region of PfAMA1 involving both domains II and III under a high degree of balancing selection relative to the rest of the protein. In summary, we developed general methods for examining how sequence-based features relate to one another in three-dimensional space and applied these methods to key P. falciparum antigens.

Population genetic structure and natural selection of Plasmodium falciparum apical membrane antigen-1 in Myanmar isolates

Background

Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is one of leading blood stage malaria vaccine candidates. However, genetic variation and antigenic diversity identified in global PfAMA-1 are major hurdles in the development of an effective vaccine based on this antigen. In this study, genetic structure and the effect of natural selection of PfAMA-1 among Myanmar P. falciparum isolates were analysed.

Methods

Blood samples were collected from 58 Myanmar patients with falciparum malaria. Full-length PfAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. PfAMA-1 sequence of each isolate was sequenced. Polymorphic characteristics and effect of natural selection were analysed with using DNASTAR, MEGA4, and DnaSP programs. Polymorphic nature and natural selection in 459 global PfAMA-1 were also analysed.

Results

Thirty-seven different haplotypes of PfAMA-1 were identified in 58 Myanmar P. falciparum isolates. Most amino acid changes identified in Myanmar PfAMA-1 were found in domains I and III. Overall patterns of amino acid changes in Myanmar PfAMA-1 were similar to those in global PfAMA-1. However, frequencies of amino acid changes differed by country. Novel amino acid changes in Myanmar PfAMA-1 were also identified. Evidences for natural selection and recombination event were observed in global PfAMA-1. Among 51 commonly identified amino acid changes in global PfAMA-1 sequences, 43 were found in predicted RBC-binding sites, B-cell epitopes, or IUR regions.

Conclusions

Myanmar PfAMA-1 showed similar patterns of nucleotide diversity and amino acid polymorphisms compared to those of global PfAMA-1. Balancing natural selection and intragenic recombination across PfAMA-1 are likely to play major roles in generating genetic diversity in global PfAMA-1. Most common amino acid changes in global PfAMA-1 were located in predicted B-cell epitopes where high levels of nucleotide diversity and balancing natural selection were found. These results highlight the strong selective pressure of host immunity on the PfAMA-1 gene. These results have significant implications in understanding the nature of Myanmar PfAMA-1 along with global PfAMA-1. They also provide useful information for the development of effective malaria vaccine based on this antigen.

Analysis of genetic diversity and population structure of gene encoding cell-traversal protein for ookinetes and sporozoites (CelTOS) vaccine candidate antigen in global Plasmodium falciparum populations

Plasmodium falciparum cell-traversal protein for ookinetes and sporozoites (PfCelTOS) has been reported as one of the most attractive malaria vaccine candidate antigens. To design a broadly effective malaria vaccine based on this antigen, it is crucial to have adequate information on genetic diversity in global PfCelTOS. Therefore, the extent of sequence diversity at the full-length of the pfceltos was assessed among both natural P. falciparum isolates collected from Iran (n = 93) and from available global pfceltos sequence data retrieved from PlasmoDB database (n = 159). Also, recombination, natural selection, the degree of genetic differentiation as well as the predicted immunodominant regions in PfCelTOS were analyzed. In total, 40 SNPs (including 1 synonymous and 39 non-synonymous) were detected in 34 positions, as compared to 3D7 sequence, which led to 66 distinct haplotypes with different frequencies. Among those haplotypes, 34 (51.5%, excluded from further analysis) were singleton haplotype and mostly detected among Senegalese parasite isolates. PfCelt-1 was found as predominant haplotype (32.6% total frequency) that was only detected in Iranian P. falciparum isolates. Nucleotide diversity was low in French Guiana (0.00236 ± 0.00203) and Iranian (0.00259 ± 0.00048) P. falciparum isolates in comparison with African populations. Evidence for positive selection by host immunity and intragenic recombination were detected that are two key factors responsible for gene evolution and genetic diversity of pfceltos gene. The results of Fst analysis and haplotype network revealed that PfCelTOS antigen displayed evident genetic structure between geographical parasite populations. In conclusion, the present analysis demonstrates that there is a limited antigenic diversity and geographic variation in global PfCelTOS, and this finding may be associated with the critical function of this antigen in cell traversal of the parasite in sporozoite and ookinete. Besides, most of the predicted B- and T-cell epitopes were located in the conserved region of the gene, but most of the amino acid replacements were located at the C-terminal region of PfCelTOS. The obtained results in this investigation could provide knowledge for better design of PfCelTOS-based malaria vaccine.

Substantial population structure of Plasmodium vivax in Thailand facilitates identification of the sources of residual transmission

Background

Plasmodium vivax transmission in Thailand has been substantially reduced over the past 10 years, yet it remains highly endemic along international borders. Understanding the genetic relationship of residual parasite populations can help track the origins of the parasites that are reintroduced into malaria-free regions within the country.

Methodology/Results

A total of 127 P. vivax isolates were genotyped from two western provinces (Tak and Kanchanaburi) and one eastern province (Ubon Ratchathani) of Thailand using 10 microsatellite markers. Genetic diversity was high, but recent clonal expansion was detected in all three provinces. Substantial population structure and genetic differentiation of parasites among provinces suggest limited gene flow among these sites. There was no haplotype sharing among the three sites, and a reduced panel of four microsatellite markers was sufficient to assign the parasites to their provincial origins.

Conclusion/Significance

Significant parasite genetic differentiation between provinces shows successful interruption of parasite spread within Thailand, but high diversity along international borders implies a substantial parasite population size in these regions. The provincial origin of P. vivax cases can be reliably determined by genotyping four microsatellite markers, which should be useful for monitoring parasite reintroduction after malaria elimination.

This study presents an updated view of the P. vivax populations along the Thai-Myanmar and the Thai-Cambodian borders. Genotyping of parasite samples collected after intensified malaria control demonstrated that despite the decline in overall transmission intensity, the genetic diversity of the P. vivax parasites remained high. Parasite populations from three border provinces showed clear genetic separation. This indicates successful interruption of parasite gene flow within Thailand, but suggests frequent parasite migration across international borders. From the analysis of 10 microsatellite markers, we further refined a set of four that are sufficient for distinguishing the provincial origins of these parasites, which should allow tracking of parasite introduction among these provinces.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

IgG antibodies to synthetic GPI are biomarkers of immune-status to both Plasmodium falciparum and Plasmodium vivax malaria in young children

Background

Further reduction in malaria prevalence and its eventual elimination would be greatly facilitated by the development of biomarkers of exposure and/or acquired immunity to malaria, as well as the deployment of effective vaccines against Plasmodium falciparum and Plasmodium vivax. A better understanding of the acquisition of immunity in naturally-exposed populations is essential for the identification of antigens useful as biomarkers, as well as to inform rational vaccine development.

Methods

ELISA was used to measure total IgG to a synthetic form of glycosylphosphatidylinositol from P. falciparum (PfGPI) in a cohort of 1–3 years old Papua New Guinea children with well-characterized individual differences in exposure to P. falciparum and P. vivax blood-stage infections. The relationship between IgG levels to PfGPI and measures of recent and past exposure to P. falciparum and P. vivax infections was investigated, as well as the association between antibody levels and prospective risk of clinical malaria over 16 months of follow-up.

Results

Total IgG levels to PfGPI were low in the young children tested. Antibody levels were higher in the presence of P. falciparum or P. vivax infections, but short-lived. High IgG levels were associated with higher risk of P. falciparum malaria (IRR 1.33–1.66, P = 0.008–0.027), suggesting that they are biomarkers of increased exposure to P. falciparum infections. Given the cross-reactive nature of antibodies to PfGPI, high IgG levels were also associated with reduced risk of P. vivax malaria (IRR 0.65–0.67, P = 0.039–0.044), indicating that these antibodies are also markers of acquired immunity to P. vivax.

Conclusions

This study highlights that in young children, IgG to PfGPI might be a useful marker of immune-status to both P. falciparum and P. vivax infections, and potentially useful to help malaria control programs to identify populations at-risk. Further functional studies are necessary to confirm the potential of PfGPI as a target for vaccine development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-2042-2) contains supplementary material, which is available to authorized users.

The Biology of Plasmodium vivax

Plasmodium vivax is the second most prevalent cause of malaria worldwide and the leading cause of malaria outside of Africa. Although infections are seldom fatal clinical disease can be debilitating and imposes significant health and economic impacts on affected populations. Estimates of transmission and prevalence intensity can be problematic because many episodes of vivax originate from hypnozoite stages in the liver that have remained dormant from previous infections by an unknown mechanism. Lack of treatment options to clear hypnozoites and the ability to infect mosquitoes before disease symptoms present represent major challenges for control and eradication of vivax malaria. Compounding these challenges is the unique biology of P. vivax and limited progress in development of experimental research tools, thereby hindering development of new drugs and vaccines. Renewed emphasis on vivax malaria research is beginning to make progress in overcoming some of these challenges.

Plasmodium vivax ligand-receptor interaction: PvAMA-1 domain I contains the minimal regions for specific interaction with CD71+ reticulocytes

The malarial parasite’s invasion is complex, active and coordinated, involving many low and high affinity interactions with receptors on target cell membrane. Proteomics analysis has described around 40 proteins in P. vivax which could be involved in reticulocyte invasion; few have been studied with the aim of elucidating how many of them establish specific interactions with their respective host cells. Given the importance of knowing which of the parasite’s protein regions are functionally important for invasion, minimum regions mediating specific interaction between Plasmodium vivax apical membrane antigen 1 (PvAMA-1) and its host cell were here elucidated. The region covering PvAMA-1 domains I and II (PvAMA-DI-II) specifically bound to the CD71⁺ red blood cell subpopulation. A 20 residue-long region (⁸¹EVENAKYRIPAGRCPVFGKG¹⁰⁰) located in domain I was capable of inhibiting PvAMA-DI-II recombinant protein binding to young reticulocytes (CD71⁺CD45⁻) and rosette formation. This conserved peptide specifically interacted with high affinity with reticulocytes (CD71⁺) through a neuraminidase- and chymotrypsin-treatment sensitive receptor. Such results showed that, despite AMA-1 having universal functions during late Plasmodium invasion stages, PvAMA-1 had reticulocyte-preferring binding regions, suggesting that P. vivax target cell selection is not just restricted to initial interactions but maintained throughout the erythrocyte invasion cycle, having important implications for designing a specific anti-P. vivax vaccine.

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

Background

Mexico advanced to the pre-elimination phase in 2009 due to a significant reduction in malaria cases, and since 2000, Plasmodium vivax is the only species transmitted. During the last two decades, malaria transmission has been mostly local and isolated to a few regions. It is important to gain further insights into the impact of control measures on the parasite population structure. Hence, the aim of the current study was to determine detailed changes in P. vivax genetic diversity and population structure based on analysing the gene that encodes the apical membrane antigen 1 (pvama1). This analysis covered from control to pre-elimination (1993–2011) in a hypo-endemic region in southern Mexico.

Results

The 213 pvama1_I-II sequences presently analysed were grouped into six periods of three years each. They showed low genetic diversity, with 15 haplotypes resolved. Among the DNA sequences, there was a gradual decrease in genetic diversity, the number of mixed genotype infections and the intensity of positive selection, in agreement with the parallel decline in malaria cases. At the same time, linkage disequilibrium (R²) increased. The three-dimensional haplotype network revealed that pvama1_I-II haplotypes were separated by 1–11 mutational steps, and between one another by 0–3 unsampled haplotypes. In the temporal network, seven haplotypes were detected in at least two of the six-time layers, and only four distinct haplotypes were evidenced in the pre-elimination phase. Structure analysis indicated that three subpopulations fluctuated over time. Only 8.5% of the samples had mixed ancestry. In the pre-elimination phase, subpopulation P1 was drastically reduced, and the admixture was absent.

Conclusions

The results suggest that P. vivax in southern Mexico evolved based on local adaptation into three “pseudoclonal” subpopulations that diversified at the regional level and persisted over time, although with varying frequency. Control measures and climate events influenced the number of malaria cases and the genetic structure. The sharp decrease in parasite diversity and other related genetic parameters during the pre-elimination phase suggests that malaria elimination is possible in the near future. These results are useful for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2156-y) contains supplementary material, which is available to authorized users.

Generation, characterization and immunogenicity of a novel chimeric recombinant protein based on Plasmodium vivax AMA-1 and MSP119

Plasmodium vivax is the most widely distributed malaria species and the most prevalent species of malaria in America and Asia. Vaccine development against P. vivax is considered a priority in the global program for the eradication of malaria. Earlier studies have characterized the Apical Membrane Antigen 1 (AMA-1) ectodomain and the C-terminal region (19kDa) of the Merozoite Surface Protein 1 (MSP-1) of P. vivax as immunodominant antigens. Based on this characterization, we designed a chimeric recombinant protein containing both merozoite immunodominant domains (PvAMA1₆₆-MSP1₁₉). The recombinant PvAMA1₆₆-MSP1₁₉ was successfully expressed in Pichia pastoris and used to immunize two different mouse strains (BALB/c and C57BL/6) in the presence of the Poly (I:C) as an adjuvant. Immunization with the chimeric protein induced high antibody titers against both proteins in both strains of mice as detected by ELISA. Antisera also recognized the native proteins expressed on the merozoites of mature P. vivax schizonts. Moreover, this antigen was able to induce IFN-gamma-secreting cells in C57BL/6 mice. These findings indicate that this novel yeast recombinant protein containing PvAMA1₆₆ and PvMSP1₁₉ is advantageous, because of improved antibody titers and cellular immune response. Therefore, this formulation should be further developed for pre-clinical trials in non-human primates as a potential candidate for a P. vivax vaccine.

Analysis of Pvama1 genes from China-Myanmar border reveals little regional genetic differentiation of Plasmodium vivax populations

Background

With the premise of diminishing parasite genetic diversity following the reduction of malaria incidence, the analysis of polymorphic antigenic markers may provide important information about the impact of malaria control on local parasite populations. Here we evaluated the genetic diversity of Plasmodium vivax apical membrane antigen 1 (Pvama1) gene in a parasite population from the China-Myanmar border and compared it with global P. vivax populations.

Methods

We performed evolutionary analysis to examine the genetic diversity, natural selection, and population differentiation of 73 Pvama1 sequences acquired from the China-Myanmar border as well as 615 publically available Pvama1 sequences from seven global P. vivax populations.

Results

A total of 308 Pvama1 haplotypes were identified among the global P. vivax isolates. The overall nucleotide diversity of Pvama1 gene among the 73 China-Myanmar border parasite isolates was 0.008 with 41 haplotypes being identified (Hd = 0.958). Domain I (DI) harbored the majority (26/33) of the polymorphic sites. The McDonald Kreitman test showed a significant positive selection across the ectodomain and the DI of Pvama1. The fixation index (F_ST) estimation between the China-Myanmar border, Thailand (0.01) and Myanmar (0.10) showed only slight geographical genetic differentiation. Notably, the Sal-I haplotype was not detected in any of the analyzed global isolates, whereas the Belem strain was restricted to the Thai population. The detected mutations are mapped outside the overlapped region of the predicted B-cell epitopes and intrinsically unstructured/disordered regions.

Conclusions

This study revealed high levels of genetic diversity of Pvama1 in the P. vivax parasite population from the China-Myanmar border with DI displaying stronger diversifying selection than other domains. There were low levels of population subdivision among parasite populations from the Greater Mekong Subregion.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1899-1) contains supplementary material, which is available to authorized users.

Diversity and evolutionary genetics of the three major Plasmodium vivax merozoite genes participating in reticulocyte invasion in southern Mexico

Background

Reported malaria cases in the Americas had been reduced to about one-half million by 2012. To advance towards elimination of this disease, it is necessary to gain insights into how the malaria parasite is evolving, including the emergence, spread and persistence of new haplotypes in affected regions. In here, the genetic diversity of the three major P. vivax merozoite genes was analyzed.

Methods

From P. vivax-infected blood samples obtained in southern Mexico (SMX) during 2006–2007, nucleotide sequences were achieved for: the 42 kDa carboxyl fragment of the merozoite surface protein-1 (msp1₄₂), domains I-II of the apical membrane antigen-1 (ama1_I-II), and domain II of the Duffy binding protein (dbp_II). Gene polymorphism was examined and haplotype networks were developed to depict parasite relationships in SMX. Then genetic diversity, recombination and natural selection were analyzed and the degree of differentiation was determined as F_ST values.

Results

The diversity of P. vivax merozoite genes in SMX was less than that of parasites from other geographic origins, with dbp_II < ama1_I-II < msp1₄₂. Ama1_I-II and msp1₄₂ exposed the more numerous haplotypes exclusive to SMX. While, all dbp_II haplotypes from SMX were separated from one to three mutational steps, the networks of ama1_I-II and msp1₄₂ were more complex; loops and numerous mutational steps were evidenced, likely due to recombination. Sings of local diversification were more evident for msp1₄₂. Sixteen combined haplotypes were determined; one of these haplotypes not detected in 2006 was highly frequent in 2007. The Rm value was higher for msp1₄₂than for ama1_I-II, being insignificant for dbp_II. The dN-dS value was highly significant for ama1_I-II and lesser so for dbp_II. The F_ST values were higher for dbp_II than msp1₄₂, and very low for ama1_I-II.

Conclusions

In SMX, P. vivax ama1_I-II, dbp_II and msp1₄₂ demonstrated limited diversity, and exhibited a differentiated parasite population. The results suggest that differential intensities of selective forces are operating on these gene fragments, and probably related to their timing, length of exposure and function during reticulocyte adhesion and invasion. Therefore, these finding are essential for mono and multivalent vaccine development and for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1266-7) contains supplementary material, which is available to authorized users.

Genetic polymorphism in domain I of the apical membrane antigen-1 among Plasmodium knowlesi clinical isolates from Peninsular Malaysia

The simian malaria parasite Plasmodium knowlesi is now recognized as a species that can cause human malaria. The first report of large scale human knowlesi malaria was in 2004 in Malaysia Borneo. Since then, hundreds of human knowlesi malaria cases have been reported in Southeast Asia. The present study investigates the genetic polymorphism of P. knowlesi DI domain of the apical membrane antigen-1 (AMA-1), a protein considered as a promising vaccine candidate for malaria. The DI domain of AMA-1 gene of P. knowlesi clinical isolates from Peninsular Malaysia was amplified by PCR, cloned into Escherichia coli, then sequenced and analysed. Ninety-seven DI domain sequences were obtained. Comparison at the nucleotide level against P. knowlesi strain H as reference sequence showed 21 synonymous and 25 nonsynonymous mutations. Nonetheless, nucleotide sequence analysis revealed low genetic diversity of the DI domain, and it was under purifying (negative) selection. At the amino acid level, 26 different haplotypes were identified and 2 were predominant haplotypes (H1, H2) with high frequencies. Phylogenetic analysis revealed that the 26 haplotypes could be clustered into 2 distinct groups (I and II). Members of the groups were basically derived from haplotypes H1 and H2, respectively.

Population genetic structure and natural selection of apical membrane antigen-1 in Plasmodium vivax Korean isolates

Background

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is a leading candidate antigen for blood stage malaria vaccine. However, antigenic variation is a major obstacle in the development of an effective vaccine based on this antigen. In this study, the genetic structure and the effect of natural selection of PvAMA-1 among Korean P. vivax isolates were analysed.

Methods

Blood samples were collected from 66 Korean patients with vivax malaria. The entire PvAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. The PvAMA-1 sequence of each isolate was sequenced and the polymorphic characteristics and effect of natural selection were analysed using the DNASTAR, MEGA4, and DnaSP programs.

Results

Thirty haplotypes of PvAMA-1, which were further classified into seven different clusters, were identified in the 66 Korean P. vivax isolates. Domain II was highly conserved among the sequences, but substantial nucleotide diversity was observed in domains I and III. The difference between the rates of non-synonymous and synonymous mutations suggested that the gene has evolved under natural selection. No strong evidence indicating balancing or positive selection on PvAMA-1 was identified. Recombination may also play a role in the resulting genetic diversity of PvAMA-1.

Conclusions

This study is the first comprehensive analysis of nucleotide diversity across the entire PvAMA-1 gene using a single population sample from Korea. Korean PvAMA-1 had limited genetic diversity compared to PvAMA-1 in global isolates. The overall pattern of genetic polymorphism of Korean PvAMA-1 differed from other global isolates and novel amino acid changes were also identified in Korean PvAMA-1. Evidences for natural selection and recombination event were observed, which is likely to play an important role in generating genetic diversity across the PvAMA-1. These results provide useful information for the understanding the population structure of P. vivax circulating in Korea and have important implications for the design of a vaccine incorporating PvAMA-1.

Designing malaria vaccines to circumvent antigen variability

Prospects for malaria eradication will be greatly enhanced by an effective vaccine, but parasite genetic diversity poses a major impediment to malaria vaccine efficacy. In recent pre-clinical and field trials, vaccines based on polymorphic Plasmodium falciparum antigens have shown efficacy only against homologous strains, raising the specter of allele-specific immunity such as that which plagues vaccines against influenza and HIV. The most advanced malaria vaccine, RTS,S, targets relatively conserved epitopes on the P. falciparum circumsporozoite protein. After more than 40 years of development and testing, RTS,S, has shown significant but modest efficacy against clinical malaria in phase 2 and 3 trials. Ongoing phase 2 studies of an irradiated sporozoite vaccine will ascertain whether the full protection against homologous experimental malaria challenge conferred by high doses of a whole organism vaccine can provide protection against diverse strains in the field. Here we review and evaluate approaches being taken to design broadly cross-protective malaria vaccines.

Malaria Molecular Epidemiology: Lessons from the International Centers of Excellence for Malaria Research Network

Molecular epidemiology leverages genetic information to study the risk factors that affect the frequency and distribution of malaria cases. This article describes molecular epidemiologic investigations currently being carried out by the International Centers of Excellence for Malaria Research (ICEMR) network in a variety of malaria-endemic settings. First, we discuss various novel approaches to understand malaria incidence and gametocytemia, focusing on Plasmodium falciparum and Plasmodium vivax. Second, we describe and compare different parasite genotyping methods commonly used in malaria epidemiology and population genetics. Finally, we discuss potential applications of molecular epidemiological tools and methods toward malaria control and elimination efforts.

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

Crystal structure of Plasmodium knowlesi apical membrane antigen 1 and its complex with an invasion-inhibitory monoclonal antibody

The malaria parasite Plasmodium knowlesi, previously associated only with infection of macaques, is now known to infect humans as well and has become a significant public health problem in Southeast Asia. This species should therefore be targeted in vaccine and therapeutic strategies against human malaria. Apical Membrane Antigen 1 (AMA1), which plays a role in Plasmodium merozoite invasion of the erythrocyte, is currently being pursued in human vaccine trials against P. falciparum. Recent vaccine trials in macaques using the P. knowlesi orthologue PkAMA1 have shown that it protects against infection by this parasite species and thus should be developed for human vaccination as well. Here, we present the crystal structure of Domains 1 and 2 of the PkAMA1 ectodomain, and of its complex with the invasion-inhibitory monoclonal antibody R31C2. The Domain 2 (D2) loop, which is displaced upon binding the Rhoptry Neck Protein 2 (RON2) receptor, makes significant contacts with the antibody. R31C2 inhibits binding of the Rhoptry Neck Protein 2 (RON2) receptor by steric blocking of the hydrophobic groove and by preventing the displacement of the D2 loop which is essential for exposing the complete binding site on AMA1. R31C2 recognizes a non-polymorphic epitope and should thus be cross-strain reactive. PkAMA1 is much less polymorphic than the P. falciparum and P. vivax orthologues. Unlike these two latter species, there are no polymorphic sites close to the RON2-binding site of PkAMA1, suggesting that P. knowlesi has not developed a mechanism of immune escape from the host’s humoral response to AMA1.

Low levels of polymorphisms and no evidence for diversifying selection on the Plasmodium knowlesi Apical Membrane Antigen 1 gene

Infection with Plasmodium knowlesi, a zoonotic primate malaria, is a growing human health problem in Southeast Asia. P. knowlesi is being used in malaria vaccine studies, and a number of proteins are being considered as candidate malaria vaccine antigens, including the Apical Membrane Antigen 1 (AMA1). In order to determine genetic diversity of the ama1 gene and to identify epitopes of AMA1 under strongest immune selection, the ama1 gene of 52 P. knowlesi isolates derived from human infections was sequenced. Sequence analysis of isolates from two geographically isolated regions in Sarawak showed that polymorphism in the protein is low compared to that of AMA1 of the major human malaria parasites, P. falciparum and P. vivax. Although the number of haplotypes was 27, the frequency of mutations at the majority of the polymorphic positions was low, and only six positions had a variance frequency higher than 10%. Only two positions had more than one alternative amino acid. Interestingly, three of the high-frequency polymorphic sites correspond to invariant sites in PfAMA1 or PvAMA1. Statistically significant differences in the quantity of three of the six high frequency mutations were observed between the two regions. These analyses suggest that the pkama1 gene is not under balancing selection, as observed for pfama1 and pvama1, and that the PkAMA1 protein is not a primary target for protective humoral immune responses in their reservoir macaque hosts, unlike PfAMA1 and PvAMA1 in humans. The low level of polymorphism justifies the development of a single allele PkAMA1-based vaccine.

Genetic diversity of Plasmodium Vivax in South of Iran: A cross-sectional study

Despite declining the number of malaria cases in Iran, increased prevalence of malaria is supposed to be due to migration from eastern neighboring countries of Iran, which are abundant in Plasmodium vivax (P. vivax). The circumsporozoite protein (CSP) of the P. vivax, is one of the candidate antigens for antimalaria vaccine. The diversity of P. vivax populations circulating in Iran has been investigated by using circumsporozoite protein (CSP) in this study. A hundred and eighteen blood samples were collected from patients diagnosed with P. vivax malaria from south of Iran during 2007-2008. All samples were analyzed by using nested PCR/ RFLP and 18 were sequenced. Genotyping of Pvcsp gene showed that VK210 type was predominant (95%) in south of Iran. Sequence analysis of Pvcsp gene revealed 6 distinct allelic variants in VK210 type. The present data indicate that there is some degree of genetic diversity among P. vivax populations in Hormozgan province of Iran. It seems that in neighbors of Iran, VK210 type is predominant, probably due to similar vector of malaria in these regions.

Strategies for designing and monitoring malaria vaccines targeting diverse antigens

After more than 50 years of intensive research and development, only one malaria vaccine candidate, “RTS,S,” has progressed to Phase 3 clinical trials. Despite only partial efficacy, this candidate is now forecast to become the first licensed malaria vaccine. Hence, more efficacious second-generation malaria vaccines that can significantly reduce transmission are urgently needed. This review will focus on a major obstacle hindering development of effective malaria vaccines: parasite antigenic diversity. Despite extensive genetic diversity in leading candidate antigens, vaccines have been and continue to be formulated using recombinant antigens representing only one or two strains. These vaccine strains represent only a small fraction of the diversity circulating in natural parasite populations, leading to escape of non-vaccine strains and challenging investigators’ abilities to measure strain-specific efficacy in vaccine trials. Novel strategies are needed to overcome antigenic diversity in order for vaccine development to succeed. Many studies have now cataloged the global diversity of leading Plasmodium falciparum and Plasmodium vivax vaccine antigens. In this review, we describe how population genetic approaches can be applied to this rich data source to predict the alleles that best represent antigenic diversity, polymorphisms that contribute to it, and to identify key polymorphisms associated with antigenic escape. We also suggest an approach to summarize the known global diversity of a given antigen to predict antigenic diversity, how to select variants that best represent the strains circulating in natural parasite populations and how to investigate the strain-specific efficacy of vaccine trials. Use of these strategies in the design and monitoring of vaccine trials will not only shed light on the contribution of genetic diversity to the antigenic diversity of malaria, but will also maximize the potential of future malaria vaccine candidates.

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

Background

As Plasmodium falciparum and Plasmodium vivax co-exist in most malaria-endemic regions outside sub-Saharan Africa, malaria control strategies in these areas must target both species in order to succeed. Population genetic analyses can predict the effectiveness of interventions including vaccines, by providing insight into patterns of diversity and evolution. The aim of this study was to investigate the population genetics of leading malaria vaccine candidate AMA1 in sympatric P. falciparum and P. vivax populations of Papua New Guinea (PNG), an area of similarly high prevalence (Pf = 22.3 to 38.8%, Pv = 15.3 to 31.8%).

Methods

A total of 72 Pfama1 and 102 Pvama1 sequences were collected from two distinct areas, Madang and Wosera, on the highly endemic PNG north coast.

Results

Despite a greater number of polymorphic sites in the AMA1 genes of P. falciparum (Madang = 52; Wosera = 56) compared to P. vivax (Madang = 36, Wosera = 34), the number of AMA1 haplotypes, haplotype diversity (Hd) and recombination (R) was far lower for P. falciparum (Madang = 12, Wosera = 20; Hd ≤0.92, R ≤45.8) than for P. vivax (Madang = 50, Wosera = 38; Hd = 0.99, R = ≤70.9). Balancing selection was detected only within domain I of AMA1 for P. vivax, and in both domains I and III for P. falciparum.

Conclusions

Higher diversity in the genes encoding P. vivax AMA1 than in P. falciparum AMA1 in this highly endemic area has important implications for development of AMA1-based vaccines in PNG and beyond. These results also suggest a smaller effective population size of P. falciparum compared to P. vivax, a finding that warrants further investigation. Differing patterns of selection on the AMA1 genes indicate that critical antigenic sites may differ between the species, highlighting the need for independent investigations of these two leading vaccine candidates.