The evolution and diversity of a low complexity vaccine candidate, merozoite surface protein 9 (MSP-9), in Plasmodium vivax and closely related species

Sequence Analysis of Pvama-1 among Plasmodium Vivax Isolates in Sistan-Baluchistan

Background

Apical Membrane antigen 1 (AMA-1) is an important membrane protein that presents in all Plasmodium species and participates in critical phases in the attraction of cells. In human, it is one of the most immunodominant antigens with a protective immune response simulation role Apical Membrane antigen 1 (AMA-1) is an important membrane protein which presents in all Plasmodium species and is located on the surface of merozoite and sporozoites that participates in critical phases in attraction of human red blood cells by merozoites and hepatocytes by sporozoites, so in human, it is one of the most immunodominant antigens with a protective immune response simulation role. Since extra information is necessary to lighten of AMA-1 scope, we equaled genetic variation in P.vivax AMA-1 from 40 Iranian isolates with those reported from the other malarious countries.

Methods

Blood samples were collected from 40 patients' positive of P.vivax, and genomic DNA was extracted from the blood. The nucleotide sequence for 446 amino acid (AA) residues (42-488 of PvAMA-1) of AMA-1 gene was amplified via PCR and then sequenced.

Result

A total of 24 different haplotypes were recognized between samples. No new haplotype was determined in this research that was reported previously in other regions of Iran and the world. We detected 37-point mutations at the nucleotide level in their sequences and showed 43 amino acid variations, at 37 positions in which 6 sites demonstrate trimorphic polymorphism, and the others were dimorphic.

Conclusion

Sequence analysis of the major haplotype showed 95% similarity with P.vivax Sal-1 AMA-1 gene and high level of allelic diversity at the domain I of PvAMA-1 among P. vivax isolates of Iran. Because PvAMA-1 is noticeable as vaccine candidate antigen, these documents provide valuable information for the development of malaria vaccine.

Immunoglobulin G responses to variant forms of Plasmodium vivax merozoite surface protein 9 upon natural infection in Thailand

Merozoite surface protein 9 (MSP9) constitutes a ligand complex involved in erythrocyte invasion by malarial merozoites and is a promising vaccine target. Plasmodium vivax MSP9 (PvMSP9) is immunogenic upon natural malaria exposure. To address whether sequence diversity in PvMSP9 among field isolates could affect natural antibody responses, the recombinant proteins representing two variants each for the N- and the C-terminal domains of PvMSP-9 were used as antigens to assess antibody reactivity among 246 P. vivax-infected patients’ sera from Tak and Ubon Ratchathani Provinces in Thailand. Results revealed that the seropositivity rates of IgG antibodies to the N-terminal antigens were higher than those to the C-terminal antigens (87.80% vs. 67.48%). Most seropositive sera were reactive to both variants, suggesting the presence of common epitopes. Variant-specific antibodies to the N- and the C-terminal antigens were detected in 15.85% and 16.70% of serum samples, respectively. These seropositivity rates were not significant difference between provinces. The seropositivity rates, levels and avidity of anti-PvMSP9 antibodies exhibited positive trends towards increasing malaria episodes. The IgG isotype responses to the N- and the C-terminal antigens were mainly IgG1 and IgG3. The profile of IgG responses may have implications for development of PvMSP9-based vaccine.

Structural diversity, natural selection and intragenic recombination in the Plasmodium vivax merozoite surface protein 9 locus in Thailand

The merozoite surface protein 9 (MSP9) of malarial parasite forms co-ligand complex with the 19 kDa fragment of merozoite surface protein 1 (MSP1) prior to erythrocyte invasion. Interruption of this process could hamper subsequent asexual erythrocytic development of malaria parasites; therefore, these proteins are considered potential vaccine candidates. In Plasmodium vivax, MSP9 (PvMSP9) contains both conserved and polymorphic repetitive domains that were immunogenic upon natural malaria exposure and conferred protection in vaccination studies in animal models. To investigate the extent of sequence diversity at this locus, 104 P. vivax isolates from 4 major malaria endemic areas of Thailand were analyzed. Results revealed that pvmsp9 contained 3 repeat domains (R1-R3) flanked by conserved domains. Repeat domains exhibit extensive sequence and length variation, in which 14, 39 and 16 haplotypes for domains R1-R3, respectively, circulated in this country. Sequence diversity in pvmsp9 among P. vivax isolates from each endemic area displayed population structure. The extent of sequence diversity in pvmsp9 isolates from the provinces of Tak, Chanthaburi, Ubon Ratchathani and Prachuap Khiri Khan in northwestern, eastern, northeastern and southwestern areas, respectively, was almost comparable and was remarkably higher than that from Yala/Narathiwat population in southern Thailand. Evidence for intragenic recombination in this locus was observed within each P. vivax population except among isolates from Yala and Narathiwat. Synonymous nucleotide diversity significantly exceeded nonsynonymous nucleotide diversity in domains R2 and R3, indicating purifying selection. However, micro-scale signatures of positive and negative selections occurred in both conserved and repeat domains, implying two opposing forces, probably from functional or structural constraint and host immune pressure, could have influenced diversity at this locus. The immunodominant T and B cell epitopes so far identified were invariant or highly conserved across isolates. Further analysis of global isolates is warranted for vaccine design based on this protein.

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.

Identifying Potential Plasmodium vivax Sporozoite Stage Vaccine Candidates: An Analysis of Genetic Diversity and Natural Selection

Parasite antigen genetic diversity represents a great obstacle when designing a vaccine against malaria caused by Plasmodium vivax. Selecting vaccine candidate antigens has been focused on those fulfilling a role in invasion and which are conserved, thus avoiding specific-allele immune responses. Most antigens described to date belong to the blood stage, thereby blocking parasite development within red blood cells, whilst studying antigens from other stages has been quite restricted. Antigens from different parasite stages are required for developing a completely effective vaccine; thus, pre-erythrocyte stage antigens able to block the first line of infection becoming established should also be taken into account. However, few antigens from this stage have been studied to date. Several P. falciparum sporozoite antigens are involved in invasion. Since 77% of genes are orthologous amongst Plasmodium parasites, P. vivax sporozoite antigen orthologs to those of P. falciparum might be present in its genome. Although these genes might have high genetic diversity, conserved functionally-relevant regions (ideal for vaccine development) could be predicted by comparing genetic diversity patterns and evolutionary rates. This study was thus aimed at searching for putative P. vivax sporozoite genes so as to analyse their genetic diversity for determining their potential as vaccine candidates. Several DNA sequence polymorphism estimators were computed at each locus. The evolutionary force (drift, selection and recombination) drawing the genetic diversity pattern observed was also determined by using tests based on polymorphism frequency spectrum as well as the type of intra- and inter-species substitutions. Likewise, recombination was assessed both indirectly and directly. The results showed that sporozoite genes were more conserved than merozoite genes evaluated to date. Putative domains implied in cell traversal, gliding motility and hepatocyte interaction had a negative selection signal, being conserved amongst different species in the genus. PvP52, PvP36, PvSPATR, PvPLP1, PvMCP1, PvTLP, PvCelTOS, and PvMB2 antigens or functionally restricted regions within them would thus seem promising vaccine candidates and could be used when designing a pre-erythrocyte and/or multi-stage vaccine against P. vivax to avoid allele-specific immune responses that could reduce vaccine efficacy.

Comparative analysis of low complexity regions in Plasmodia

Low complexity regions (LCRs) are a common feature shared by many genomes, but their evolutionary and functional significance remains mostly unknown. At the core of the uncertainty is a poor understanding of the mechanisms that regulate their retention in genomes, whether driven by natural selection or neutral evolution. Applying a comparative approach of LCRs to multiple strains and species is a powerful approach to identify patterns of conservation in these regions. Using this method, we investigate the evolutionary history of LCRs in the genus Plasmodium based on orthologous protein coding genes shared by 11 species and strains from primate and rodent-infecting pathogens. We find multiple lines of evidence in support of natural selection as a major evolutionary force shaping the composition and conservation of LCRs through time and signatures that their evolutionary paths are species specific. Our findings add a comparative analysis perspective to the debate on the evolution of LCRs and harness the power of sequence comparisons to identify potential functionally important LCR candidates.

Host Cell Tropism and Adaptation of Blood-Stage Malaria Parasites: Challenges for Malaria Elimination

Plasmodium falciparum and Plasmodium vivax account for most of the mortality and morbidity associated with malaria in humans. Research and control efforts have focused on infections caused by P. falciparum and P. vivax, but have neglected other malaria parasite species that infect humans. Additionally, many related malaria parasite species infect nonhuman primates (NHPs), and have the potential for transmission to humans. For malaria elimination, the varied and specific challenges of all of these Plasmodium species will need to be considered. Recent advances in molecular genetics and genomics have increased our knowledge of the prevalence and existing diversity of the human and NHP Plasmodium species. We are beginning to identify the extent of the reservoirs of each parasite species in humans and NHPs, revealing their origins as well as potential for adaptation in humans. Here, we focus on the red blood cell stage of human infection and the host cell tropism of each human Plasmodium species. Determinants of tropism are unique among malaria parasite species, presenting a complex challenge for malaria elimination.

Analysis of polymorphisms and selective pressures on ama1 gene in Plasmodium knowlesi isolates from Sabah, Malaysia

The apical membrane antigen-1 (AMA-1) of Plasmodium spp. is a merozoite surface antigen that is essential for the recognition and invasion of erythrocytes. Polymorphisms occurring in this surface antigen will cause major obstacles in developing effective malaria vaccines based on AMA-1. The objective of this study was to characterize ama1 gene in Plasmodium knowlesi isolates from Sabah. DNA was extracted from blood samples collected from Keningau, Kota Kinabalu and Kudat. The Pkama1 gene was amplified using nested PCR and subjected to bidirectional sequencing. Analysis of DNA sequence revealed that most of the nucleotide polymorphisms were synonymous and concentrated in domain I of PkAMA-1. Forteen haplotypes were identified based on amino acid variations and haplotype K5 was the most common haplotype. d_N/d_S ratios implied that purifying selection was prevalent in Pkama1 gene. Fu and Li's D and F values further provided evidence of negative selection acting on domain II of Pkama1. Lownucleotide diversitywas also detected for the Pkama1 sequences,which is similar to reports on Pkama1 from Peninsular Malaysia and Sarawak. The presence of purifying selection and low nucleotide diversity indicated that domain II of Pkama1 can be used as a target for vaccine development.

Evolution of the merozoite surface protein 7 (msp7) family in Plasmodium vivax and P. falciparum: A comparative approach

Malaria parasites (genus Plasmodium) are a diverse group found in many species of vertebrate hosts. These parasites invade red blood cells in a complex process comprising several proteins, many encoded by multigene families, one of which is merozoite surface protein 7 (msp7). In the case of Plasmodium vivax, the most geographically widespread human-infecting species, differences in the number of paralogs within multigene families have been previously explained, at least in part, as potential adaptations to the human host. To explore this in msp7, we studied its orthologs in closely related nonhuman primate parasites; investigating both paralog evolutionary history and genetic polymorphism. The emerging patterns were then compared with the human parasite Plasmodium falciparum. We found that the evolution of the msp7 family is consistent with a birth-and-death model, where duplications, pseudogenizations, and gene loss events are common. However, all paralogs in P. vivax and P. falciparum had orthologs in their closely related species in non-human primates indicating that the ancestors of those paralogs precede the events leading to their origins as human parasites. Thus, the number of paralogs cannot be explained as an adaptation to human hosts. Although there is no functional information for msp7 in P. vivax, we found evidence for purifying selection in the genetic polymorphism of some of its paralogs as well as their orthologs in closely related non-human primate parasites. We also found evidence indicating that a few of P. vivax's paralogs may have diverged from their orthologs in non-human primates by episodic positive selection. Hence, they may had been under selection when the lineage leading to P. vivax diverged from the Asian non-human primates and switched into Homininae. All these lines of evidence suggest that msp7 is functionally important in P. vivax.

Genetic diversity and natural selection in the rhoptry-associated protein 1 (RAP-1) of recent Plasmodium knowlesi clinical isolates from Malaysia

BACKGROUND

The Plasmodium rhoptry-associated protein 1 (RAP-1) plays a role in the formation of the parasitophorous vacuole following the parasite's invasion of red blood cells. Although there is some evidence that the protein is recognized by the host's immune system, study of Plasmodium falciparum RAP-1 (PfRAP-1) suggests that it is not under immune pressure. A previous study on five old (1953-1962) P. knowlesi strains suggested that RAP-1 has limited genetic polymorphism and might be under negative selection. In the present study, 30 recent P. knowlesi isolates were studied to obtain a better insight into the polymorphism and natural selection of PkRAP-1.

METHODS

Blood samples from 30 knowlesi malaria patients were used. These samples were collected between 2010 and 2014. The PkRAP-1 gene, which contains two exons, was amplified by PCR, cloned into Escherichia coli and sequenced. Genetic diversity and phylogenetic analyses were performed using MEGA6 and DnaSP ver. 5.10.00 programs.

RESULTS

Thirty PkRAP-1 sequences were obtained. The nucleotide diversity (π) of exons 1, 2 and the total coding region (0.00915, 0.01353 and 0.01298, respectively) were higher than those of the old strains. Further analysis revealed a lower rate of non-synonymous (dN) than synonymous (dS) mutations, suggesting negative (purifying) selection of PkRAP-1. Tajima's D test and Fu and Li's D test values were not significant. At the amino acid level, 22 haplotypes were established with haplotype H7 having the highest frequency (7/34, 20.5 %). In the phylogenetic analysis, two distinct haplotype groups were observed. The first group contained the majority of the haplotypes, whereas the second had fewer haplotypes.

CONCLUSIONS

The present study found higher genetic polymorphism in the PkRAP-1 gene than the polymorphism level reported in a previous study. This observation may stem from the difference in sample size between the present (n = 30) and the previous (n = 5) study. Synonymous and non-synonymous mutation analysis indicated purifying (negative) selection of the gene. The separation of PkRAP-1haplotypes into two groups provides further evidence to the postulation of two distinct P. knowlesi types or lineages.

In silico Identification and Validation of a Linear and Naturally Immunogenic B-Cell Epitope of the Plasmodium vivax Malaria Vaccine Candidate Merozoite Surface Protein-9

Synthetic peptide vaccines provide the advantages of safety, stability and low cost. The success of this approach is highly dependent on efficient epitope identification and synthetic strategies for efficacious delivery. In malaria, the Merozoite Surface Protein-9 of Plasmodium vivax (PvMSP9) has been considered a vaccine candidate based on the evidence that specific antibodies were able to inhibit merozoite invasion and recombinant proteins were highly immunogenic in mice and humans. However the identities of linear B-cell epitopes within PvMSP9 as targets of functional antibodies remain undefined. We used several publicly-available algorithms for in silico analyses and prediction of relevant B cell epitopes within PMSP9. We show that the tandem repeat sequence EAAPENAEPVHENA (PvMSP9_E795-A808) present at the C-terminal region is a promising target for antibodies, given its high combined score to be a linear epitope and located in a putative intrinsically unstructured region of the native protein. To confirm the predictive value of the computational approach, plasma samples from 545 naturally exposed individuals were screened for IgG reactivity against the recombinant PvMSP9-RIRII_729-972 and a synthetic peptide representing the predicted B cell epitope PvMSP9_E795-A808. 316 individuals (58%) were responders to the full repetitive region PvMSP9-RIRII, of which 177 (56%) also presented total IgG reactivity against the synthetic peptide, confirming it validity as a B cell epitope. The reactivity indexes of anti-PvMSP9-RIRII and anti-PvMSP9_E795-A808 antibodies were correlated. Interestingly, a potential role in the acquisition of protective immunity was associated with the linear epitope, since the IgG1 subclass against PvMSP9_E795-A808 was the prevalent subclass and this directly correlated with time elapsed since the last malaria episode; however this was not observed in the antibody responses against the full PvMSP9-RIRII. In conclusion, our findings identified and experimentally confirmed the potential of PvMSP9_E795-A808 as an immunogenic linear B cell epitope within the P. vivax malaria vaccine candidate PvMSP9 and support its inclusion in future subunit vaccines.

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.

Heterogeneous genetic diversity pattern in Plasmodium vivax genes encoding merozoite surface proteins (MSP) -7E, -7F and -7L

Background

The msp-7 gene has become differentially expanded in the Plasmodium genus; Plasmodium vivax has the highest copy number of this gene, several of which encode antigenic proteins in merozoites.

Methods

DNA sequences from thirty-six Colombian clinical isolates from P. vivax (pv) msp-7E, −7F and -7L genes were analysed for characterizing and studying the genetic diversity of these pvmsp-7 members which are expressed during the intra-erythrocyte stage; natural selection signals producing the variation pattern so observed were evaluated.

Results

The pvmsp-7E gene was highly polymorphic compared to pvmsp-7F and pvmsp-7L which were seen to have limited genetic diversity; pvmsp-7E polymorphism was seen to have been maintained by different types of positive selection. Even though these copies seemed to be species-specific duplications, a search in the Plasmodium cynomolgi genome (P. vivax sister taxon) showed that both species shared the whole msp-7 repertoire. This led to exploring the long-term effect of natural selection by comparing the orthologous sequences which led to finding signatures for lineage-specific positive selection.

Conclusions

The results confirmed that the P. vivax msp-7 family has a heterogeneous genetic diversity pattern; some members are highly conserved whilst others are highly diverse. The results suggested that the 3′-end of these genes encode MSP-7 proteins’ functional region whilst the central region of pvmsp-7E has evolved rapidly. The lineage-specific positive selection signals found suggested that mutations occurring in msp-7s genes during host switch may have succeeded in adapting the ancestral P. vivax parasite population to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-495) contains supplementary material, which is available to authorized users.

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.