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

Genetic differentiation of Plasmodium vivax duffy binding protein in Ethiopia and comparison with other geographical isolates

BACKGROUND

Plasmodium vivax Duffy binding protein (PvDBP) is a merozoite surface protein located in the micronemes of P. vivax. The invasion of human reticulocytes by P. vivax merozoites depends on the parasite DBP binding domain engaging Duffy Antigen Receptor for Chemokine (DARC) on these red blood cells (RBCs). PvDBPII shows high genetic diversity which is a major challenge to its use in the development of a vaccine against vivax malaria.

METHODS

A cross-sectional study was conducted from February 2021 to September 2022 in five study sites across Ethiopia. A total of 58 blood samples confirmed positive for P. vivax by polymerase chain reaction (PCR) were included in the study to determine PvDBPII genetic diversity. PvDBPII were amplified using primers designed from reference sequence of P. vivax Sal I strain. Assembling of sequences was done using Geneious Prime version 2023.2.1. Alignment and phylogenetic tree constructions using MEGA version 10.1.1. Nucleotide diversity and haplotype diversity were analysed using DnaSP version 6.12.03, and haplotype network was generated with PopART version 1.7.

RESULTS

The mean age of the participants was 25 years, 5 (8.6%) participants were Duffy negatives. From the 58 PvDBPII sequences, seven haplotypes based on nucleotide differences at 8 positions were identified. Nucleotide diversity and haplotype diversity were 0.00267 ± 0.00023 and 0.731 ± 0.036, respectively. Among the five study sites, the highest numbers of haplotypes were identified in Arbaminch with six different haplotypes while only two haplotypes were identified in Gambella. The phylogenetic tree based on PvDBPII revealed that parasites of different study sites shared similar genetic clusters with few exceptions. Globally, a total of 39 haplotypes were identified from 223 PvDBPII sequences representing different geographical isolates obtained from NCBI archive. The nucleotide and haplotype diversity were 0.00373 and 0.845 ± 0.015, respectively. The haplotype prevalence ranged from 0.45% to 27.3%. Two haplotypes were shared among isolates from all geographical areas of the globe.

CONCLUSIONS

PvDBPII of the Ethiopian P. vivax isolates showed low nucleotide but high haplotype diversity, this pattern of genetic variability suggests that the population may have undergone a recent expansion. Among the Ethiopian P. vivax isolates, almost half of the sequences were identical to the Sal-I reference sequence. However, there were unique haplotypes observed in the Ethiopian isolates, which does not share with isolates from other geographical areas. There were two haplotypes that were common among populations across the globe. Categorizing population haplotype frequency can help to determine common haplotypes for designing an effective blood-stage vaccine which will have a significant role for the control and elimination of P. vivax.

Genetic diversity and molecular evolution of Plasmodium vivax Duffy Binding Protein and Merozoite Surface Protein-1 in northwestern Thailand

The local diversity and population structure of malaria parasites vary across different regions of the world, reflecting variations in transmission intensity, host immunity, and vector species. This study aimed to use amplicon sequencing to investigate the genotypic patterns and population structure of P. vivax isolates from a highly endemic province of Thailand in recent years. Amplicon deep sequencing was performed on 70 samples for the 42-kDa region of pvmsp1 and domain II of pvdbp. Unique haplotypes were identified and a network constructed to illustrate genetic relatedness in northwestern Thailand. Based on this dataset of 70 samples collected between 2015 and 2021, 16 and 40 unique haplotypes were identified in pvdbpII and pvmsp142kDa, respectively. Nucleotide diversity was higher in pvmsp142kDa than in pvdbpII (π = 0.027 and 0.012), as was haplotype diversity (Hd = 0.962 and 0.849). pvmsp142kDa also showed a higher recombination rate and higher levels of genetic differentiation (Fst) in northwestern Thailand versus other regions (0.2761-0.4881). These data together suggested that the genetic diversity of P. vivax in northwestern Thailand at these two studied loci evolved under a balancing selection, most likely host immunity. The lower genetic diversity of pvdbpII may reflect its stronger functional constrain. In addition, despite the balancing selection, a decrease in genetic diversity was observed. Hd of pvdbpII decreased from 0.874 in 2015-2016 to 0.778 in 2018-2021; π of pvmsp142kDa decreased from 0.030 to 0.022 over the same period. Thus, the control activities must have had a strong impact on the parasite population size. The findings from this study provide an understanding of P. vivax population structure and the evolutionary force on vaccine candidates. They also established a new baseline for tracking future changes in P. vivax diversity in the most malarious area of Thailand.

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.

Spatiotemporal Changes in Plasmodium vivax msp142 Haplotypes in Southern Mexico: From the Control to the Pre-Elimination Phase

For 20 years, Plasmodium vivax has been the only prevalent malaria species in Mexico, and cases have declined significantly and continuously. Spatiotemporal genetic studies can be helpful for understanding parasite dynamics and developing strategies to weaken malaria transmission, thus facilitating the elimination of the parasite. The aim of the current contribution was to analyze P. vivax-infected blood samples from patients in southern Mexico during the control (1993–2007) and pre-elimination phases (2008–2011). Nucleotide and haplotype changes in the pvmsp142 fragment were evaluated over time. The majority of multiple genotype infections occurred in the 1990s, when the 198 single nucleotide sequences exhibited 57 segregating sites, 64 mutations, and 17 haplotypes. Nucleotide and genetic diversity parameters showed subtle fluctuations from across time, in contrast to the reduced haplotype diversity and the increase in the R2 index and Tajima’s D value from 2008 to 2011. The haplotype network consisted of four haplogroups, the geographical distribution of which varied slightly over time. Haplogroup-specific B-cell epitopes were predicted. Since only high-frequency and divergent haplotypes persisted, there was a contraction of the parasite population. Given that 84% of haplotypes were exclusive to Mesoamerica, P. vivax flow is likely circumscribed to this region, representing important information for parasite surveillance.

Plasmodium vivax Cysteine-Rich Protective Antigen Polymorphism at Exon-1 Shows Recombination and Signatures of Balancing Selection

Plasmodium vivax Cysteine-Rich Protective Antigen (CyRPA) is a merozoite protein participating in the parasite invasion of human reticulocytes. During natural P. vivax infection, antibody responses against PvCyRPA have been detected. In children, low anti-CyRPA antibody titers correlated with clinical protection, which suggests this protein as a potential vaccine candidate. This work analyzed the genetic and amino acid diversity of pvcyrpa in Mexican and global parasites. Consensus coding sequences of pvcyrpa were obtained from seven isolates. Other sequences were extracted from a repository. Maximum likelihood phylogenetic trees, genetic diversity parameters, linkage disequilibrium (LD), and neutrality tests were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. In 22 sequences from Southern Mexico, two synonymous and 21 nonsynonymous mutations defined nine private haplotypes. These parasites had the highest LD-R² index and the lowest nucleotide diversity compared to isolates from South America or Asia. The nucleotide diversity and Tajima's D values varied across the coding gene. The exon-1 sequence had greater diversity and Rm values than those of exon-2. Exon-1 had significant positive values for Tajima's D, β-α values, and for the Z (HA: dN > dS) and MK tests. These patterns were similar for parasites of different origin. The polymorphic amino acid residues at PvCyRPA resembled the conformational B-cell peptides reported in PfCyRPA. Diversity at pvcyrpa exon-1 is caused by mutation and recombination. This seems to be maintained by balancing selection, likely due to selective immune pressure, all of which merit further study.

Exploration of Plasmodium vivax merozoite surface proteins 1 and 7 genetic diversity in Brazilian Amazon and Rio de Janeiro Atlantic Forest

Plasmodium vivax merozoite surface proteins (PvMSP) 1 and 7 are considered vaccine targets. Genetic diversity knowledge is crucial to assess their potential as immunogens and to provide insights about population structure in different epidemiological contexts. Here, we investigate the variability of pvmsp-1₄₂, pvmsp-7E, and pvmsp-7F genes in 227 samples from the Brazilian Amazon (BA) and Rio de Janeiro Atlantic Forest (AF). pvmsp-1₄₂ has 63 polymorphisms - 57 nonsynonymous - generating a nucleotide diversity of π = 0.009 in AF, and π = 0.018 in BA. In pvmsp-7E, 134 polymorphisms - 103 nonsynonymous - generate the nucleotide diversity of π = 0.027 in AF, and π = 0.042 in BA. The pvmsp-7F has only two SNPs - A610G and A1054T -, with nucleotide diversity of π = 0.0004 in AF, and π = 0.0007 in BA. The haplotype diversity of pvmsp-1₄₂, pvmsp-7E, and pvmsp-7F genes is 0.997, 1.00, and 0.649, respectively. None of the pvmsp-1₄₂ or pvmsp-7E sequences are identical to the Salvador 1 strain's sequence. Conversely, most of pvmsp-7F sequences (94/48%) are identical to Sal-1. We evaluated eight B-cell epitopes in pvmsp-7E, four of them showed higher nucleotide diversity compared to pvmsp-7E's epitopes. Positive selection was detected in pvmsp-1₄₂, pvmsp-7E central region, and pvmsp-7F with Tajima's D. In pvmsp-7E, the significant nucleotide and haplotype diversities with low genetic differentiation, could be indicative of balancing selection. The genetic differentiation of pvmsp-1₄₂ (0.315) and pvmsp-7F (0.354) genes between AF and BA regions is significant, which is not the case for pvmsp-7E (0.193). We conclude that pvmsp-1₄₂ and pvmsp-7E have great genetic diversity even in AF region, an enclosure area with deficient transmission levels of P. vivax zoonotic malaria. In both Brazilian regions, pvmsp-1₁₉, pvmsp-7E, and pvmsp-7F are conserved, most likely due to their roles in parasite survival, and could be considered potential targets for a "blood-stage vaccine".

Global distribution of single amino acid polymorphisms in Plasmodium vivax Duffy-binding-like domain and implications for vaccine development efforts

Plasmodium vivax (Pv) malaria continues to be geographically widespread with approximately 15 million worldwide cases annually. Along with other proteins, Duffy-binding proteins (DBPs) are used by plasmodium for RBC invasion and the parasite-encoded receptor binding regions lie in their Duffy-binding-like (DBL) domains-thus making it a prime vaccine candidate. This study explores the sequence diversity in PvDBL globally, with an emphasis on India as it remains a major contributor to the global Pv malaria burden. Based on 1358 PvDBL protein sequences available in NCBI, we identified 140 polymorphic sites within 315 residues of PvDBL. Alarmingly, country-wise mapping of SAAPs from field isolates revealed varied and distinct polymorphic profiles for different nations. We report here 31 polymorphic residue positions in the global SAAP profile, most of which map to the PvDBL subdomain 2 (α1-α6). A distinct clustering of SAAPs distal to the DARC-binding sites is indicative of immune evasive strategies by the parasite. Analyses of PvDBL-neutralizing antibody complexes revealed that between 24% and 54% of interface residues are polymorphic. This work provides a framework to recce and expand the polymorphic space coverage in PvDBLs as this has direct implications for vaccine development studies. It also emphasizes the significance of surveying global SAAP distributions before or alongside the identification of vaccine candidates.

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.

Genetic characterization of Plasmodium vivax in the Kyrgyz Republic

At the end of 2016, Kyrgyz Republic was certified by the World Health Organization as a malaria-free country, while only a decade ago this disease posed a serious health threat. The progress achieved by Kyrgyz Republic provides a unique example of tertian (Plasmodium vivax) malaria elimination. This success was based on an integrated approach, including measures for the treatment of infected people and disease prevention, vector control and the development of an effective national epidemiological surveillance system. Lower P. vivax msp-1, msp-3α, csp and dbp_II genes polymorphism was revealed in Kyrgyz Republic in compare with that in Tajikistan. Molecular characterization of the causative agent found that P. vivax populations in Kyrgyz Republic was comprised by several lineages, highly divergent in the south-western and genetically homogeneous in the northern regions of Kyrgyz Republic, d. Such profile in the northern regions was compatible with several recent introductions rather than a long-term endemic circulation of the parasite. A low level of genetic variability suggested that the parasitic systems of tertian malaria, were not adapted, which, along with other factors, largely determined the possibility of malaria elimination in northern Kyrgyz Republic. Other determinants included environmental, social, and epidemiological factors that limited the spread of malaria. South-western Kyrgyz Republic, a region with a high level of interstate migration, requires considerable attention to prevent the spread of malaria.

Diversity pattern of Duffy binding protein sequence among Duffy-negatives and Duffy-positives in Sudan

Background

Vivax malaria is a leading public health concern worldwide. Due to the high prevalence of Duffy-negative blood group population, Plasmodium vivax in Africa historically is less attributable and remains a neglected disease. The interaction between Duffy binding protein and its cognate receptor, Duffy antigen receptor for chemokine plays a key role in the invasion of red blood cells and serves as a novel vaccine candidate against P. vivax. However, the polymorphic nature of P. vivax Duffy binding protein (DBP), particularly N-terminal cysteine-rich region (PvDBPII), represents a major obstacle for the successful design of a DBP-based vaccine to enable global protection. In this study, the level of pvdbpII sequence variations, Duffy blood group genotypes, number of haplotypes circulating, and the natural selection at pvdbpII in Sudan isolates were analysed and the implication in terms of DBP-based vaccine design was discussed.

Methods

Forty-two P. vivax-infected blood samples were collected from patients from different areas of Sudan during 2014–2016. For Duffy blood group genotyping, the fragment that indicates GATA-1 transcription factor binding site of the FY gene (− 33T > C) was amplified by PCR and sequenced by direct sequencing. The region II flanking pvdbpII was PCR amplified and sequenced by direct sequencing. The genetic diversity and natural selection of pvdbpII were done using DnaSP ver 5.0 and MEGA ver 5.0 programs. Based on predominant, non-synonymous, single nucleotide polymorphisms (SNPs), prevalence of Sudanese haplotypes was assessed in global isolates.

Results

Twenty SNPs (14 non-synonymous and 6 synonymous) were identified in pvdbpII among the 42 Sudan P. vivax isolates. Sequence analysis revealed that 11 different PvDBP haplotypes exist in Sudan P. vivax isolates and the region has evolved under positive selection. Among the identified PvDBP haplotypes five PvDBP haplotypes were shared among Duffy-negative as well as Duffy-positive individuals. The high selective pressure was mainly found on the known B cell epitopes (H3) of pvdbpII. Comparison of Sudanese haplotypes, based on 10 predominant non-synonymous SNPs with 10 malaria-endemic countries, demonstrated that Sudanese haplotypes were prevalent in most endemic countries.

Conclusion

This is the first pvdbp genetic diversity study from an African country. Sudanese isolates display high haplotype diversity and the gene is under selective pressure. Haplotype analysis indicated that Sudanese haplotypes are a representative sample of the global population. However, studies with a large number of samples are needed. These findings would be valuable for the development of PvDBP-based malaria vaccine.

Electronic supplementary material

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

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.