Spatial variation in genetic diversity and natural selection on the thrombospondin-related adhesive protein locus of Plasmodium vivax (PvTRAP)

Genetic Diversity and Natural Selection of Plasmodium vivax Merozoite Surface Protein 8 in Global Populations

Plasmodium vivax Merozoite Surface Protein 8 (PvMSP8) is a promising candidate target for the development of multi-component vaccines. Therefore, determining the genetic variation pattern of Pvmsp8 is essential in providing a reference for the rational design of the P. vivax malaria vaccines. This study delves into the genetic characteristics of the Pvmsp8 gene, specifically focusing on samples from the China-Myanmar border (CMB) region, and contrasts these findings with broader global patterns. The study uncovers that Pvmsp8 exhibits a notable level of conservation across different populations, with limited polymorphisms and relatively low nucleotide diversity (0.00023-0.00120). This conservation contrasts starkly with the high polymorphisms found in other P. vivax antigens such as Pvmsp1. A total of 25 haplotypes and 14 amino acid mutation sites were identified in the global populations, and all mutation sites were confined to non-functional regions. The study also notes that most CMB Pvmsp8 haplotypes are shared among Burmese, Cambodian, Thai, and Vietnamese populations, indicating less geographical variance, but differ notably from those found in Pacific island regions or the Panama. The findings underscore the importance of considering regional genetic diversity in P. vivax when developing targeted malaria vaccines. Non departure from neutral evolution were found by Tajima's D test, however, statistically significant differences were observed between the kn/ks rates. The study's findings are crucial in understanding the evolution and population structure of the Pvmsp8 gene, particularly during regional malaria elimination efforts. The highly conserved nature of Pvmsp8, combined with the lack of mutations in its functional domain, presents it as a promising candidate for developing a broad and effective P. vivax vaccine. This research thus lays a foundation for the rational development of multivalent malaria vaccines targeting this genetically stable antigen.

Genetic Diversity and Population Genetic Analysis of Plasmodium falciparum Thrombospondin Related Anonymous Protein (TRAP) in Clinical Samples from Saudi Arabia

The thrombospondin related anonymous protein (TRAP) is considered one of the most important pre-erythrocytic vaccine targets. Earlier population genetic studies revealed the TRAP gene to be under strong balancing natural selection. This study is the first attempt to analyze genetic diversity, natural selection, phylogeography and population structure in 199 clinical samples from Saudi Arabia using the full-length PfTRAP gene. We found the rate of nonsynonymous substitutions to be significantly higher than that of synonymous substitutions in the clinical samples, indicating a strong positive or diversifying selection for the full-length gene and the Von Willebrand factor (VWF). The nucleotide diversity was found to be π~0.00789 for the full-length gene; however, higher nucleotide diversity was observed for the VWF compared to the thrombospondin repeat region (TSP). Deduction of the amino acid sequence alignment of the PNP repeat region in the Saudi samples revealed six genotypes characterized by tripeptide repeat motifs (PNP, ANP, ENP and SNP). Haplotype network, population structure and population differentiation analyses indicated four distinct sub-populations in spite of the low geographical distance between the sampling sites. Our results suggest the likeliness of independent parasite evolution, creating opportunities for further adaptation, including host transition, and making malaria control even more challenging.

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.

Insights into the molecular diversity of Plasmodium vivax merozoite surface protein-3γ (pvmsp3γ), a polymorphic member in the msp3 multi-gene family

Plasmodium vivax merozoite surface protein 3 (PvMSP3) is encoded by a multi-gene family. Of these, PvMSP3α, PvMSP3β and PvMSP3γ, are considered to be vaccine targets. Despite comprehensive analyses of PvMSP3α and PvMSP3β, little is known about structural and sequence diversity in PvMSP3γ. Analysis of 118 complete pvmsp3γ sequences from diverse endemic areas of Thailand and 9 reported sequences has shown 86 distinct haplotypes. Based on variation in insert domains, pvmsp3γ can be classified into 3 types, i.e. Belem, Salvador I and NR520. Imperfect nucleotide repeats were found in six regions of the gene; none encoded tandem amino acid repeats. Predicted coiled-coil heptad repeats were abundant in the protein and displayed variation in length and location. Interspersed phase shifts occurred in the heptad arrays that may have an impact on protein structure. Polymorphism in pvmsp3γ seems to be generated by intragenic recombination and driven by natural selection. Most P. vivax isolates in Thailand exhibit population structure, suggesting limited gene flow across endemic areas. Phylogenetic analysis has suggested that insert domains could have been subsequently acquired during the evolution of pvmsp3γ. Sequence and structural diversity of PvMSP3γ may complicate vaccine design due to alteration in predicted immunogenic epitopes among variants.

Population genetic structure analysis of thrombospondin-related adhesive protein (TRAP) as a vaccine candidate antigen in worldwide Plasmodium falciparum isolates

Antigenic diversity is a major concern in malaria vaccine development that requires to be considered in developing a malaria vaccine. Plasmodium falciparum thrombospondin-related adhesive protein (PfTRAP) is a leading malaria vaccine candidate antigen. In the current study, we investigated the level of genetic diversity and natural selection of pftrap sequences in P. falciparum isolates from Iran (n = 47). The gene diversity of Iranian pftrap sequences was also compared to available global pftrap sequences deposited in the GenBank or PlasmoDB databases (n = 220). Comparison of Iranian PfTRAP sequences with T9/96 reference sequence showed the presence of 35 amino acid changes in 32 positions and a limited variation in repeat sequences, leading to 13 distinct haplotypes. The overall nucleotide diversity (π) for the ectodomain of Iranian pftrap sequences was 0.00444 ± 0.00043, with the highest diversity in Domain IV. Alignment comparison of global PfTRAP sequences with T9/96 reference sequence indicated 96 amino acid replacements as well as extensive variable repeat sequences (9-23 repeats), which led to 192 haplotypes. Among the global isolates, the lowest nucleotide diversity was detected in French Guianan (0.00428 ± 0.00163) and Iranian (0.00444 ± 0.00043) pftrap sequences, and the most variation was observed in domains II and IV in all populations. The dN-dS value displayed the evidence of positive selection due to recombination and immune system pressure. The Fst analysis revealed a gene flow between African populations; however, genetic differentiation observed between Iranian and other populations probably was due to gene flow barriers. Both conserved and variable epitopes were predicted in B- and T-cell epitopes of PfTRAP antigen. The obtained results from this study could be helpful for developing a PfTRAP-based malaria vaccine.

Polymorphism and natural selection in the merozoite surface protein 3F2 (PVX_97710) locus of Plasmodium vivax among field isolates

Plasmodium vivax, the chronic relapsing human malaria parasite with the most widespread distribution, possesses proteins associated with the merozoite surface that could be targets for host immune responses and potential vaccine candidates. Of these, the merozoite surface protein 3 of P. vivax (PvMSP3) is an attractive vaccine target as well as a genetic marker for epidemiological surveillance. PvMSP3 comprises a group of protein members encoded by a multigene family. Although some protein members, i.e. PvMSP3α and PvMSP3β, have been targets for molecular and immunological investigations, the most abundantly expressed protein member during late asexual erythrocytic stages, PvMSP3F2 (PVX_97710), remains unexplored. To address domain organization and evolution of this locus, the complete coding sequences of 31 P. vivax isolates from diverse malaria endemic areas of Thailand were analyzed and compared with 10 previously reported sequences. Results revealed that all PvMSP3F2 sequences differed but could be divided into 5 repeat-containing domains flanked by 6 non-repeat domains. Repeat domains II and IV at the 5' portion and domain X at the 3' portion exhibited extensive sequence and length variation whereas repeat domains VI and VIII located at the central region were relatively conserved. Despite a repertoire of PvMSP3F2 variants, predicted coiled-coil tertiary structure and predicted B-cell epitopes seem to be maintained. Evidence of intragenic recombination has been detected among field isolates in Thailand that could enhance sequence diversity at this locus. Non-repeat domains I and IX located at the 5' end and at the 3' portion, respectively, seem to have evolved under purifying selection. Evidence of positive selection was found in non-repeat domains III, V and VII where a number of predicted HLA class I epitopes were identified. Amino acid substitutions in these predicted epitopes could alter predicted peptide binding affinity or abolish peptide epitope property, suggesting that polymorphism in these epitopes conferred host immune evasion. Further studies on PvMSP3F2 are warranted, particularly on interaction with host immune system and the potential role of this PvMSP3 protein member as a vaccine target.

Diversity and natural selection on the thrombospondin-related adhesive protein (TRAP) gene of Plasmodium knowlesi in Malaysia

Background

Plasmodium knowlesi a parasite of the macaques is currently the most common cause of human malaria in Malaysia. The thrombospondin-related adhesive protein (TRAP) gene is pre-erythrocytic stage antigen. It is a well-characterized vaccine candidate in Plasmodium vivax and Plasmodium falciparum, however, no study has been done in the orthologous gene of P. knowlesi. This study investigates nucleotide diversity, haplotypes, natural selection and population differentiation of full-length pktrap genes in clinical samples from Malaysia.

Methods

Forty full-length pktrap sequences from clinical isolates of Malaysia along with the reference H-strain were downloaded from published databases. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. McDonald–Kreitman test was conducted using P. vivax and Plasmodium coatneyi as ortholog sequence in DnaSP 5.10 software. Population genetic differentiation index (F_ST) of parasite populations was determined using Arlequin v3.5. Phylogenetic relationships between trap ortholog genes were determined using MEGA 5.0 software.

Results

Comparison of 40 full-length pktrap sequences along with the H-strain identified 74 SNPs (53 non-synonymous and 21 synonymous substitutions) resulting in 29 haplotypes. Analysis of the full-length gene showed that the nucleotide diversity was lower compared to its nearest ortholog pvtrap. Domain-wise analysis indicated that the proline/asparagine rich region had higher nucleotide diversity compared to the von Willebrand factor domain and the thrombospondin-type-1 domain. McDonald–Kreitman test identified that the ratio of the number of nonsynonymous to synonymous polymorphic sites within P. knowlesi was significantly higher than that of the number of nonsynonymous to synonymous fixed sites between P. knowlesi and P. vivax. The von Willebrand factor domain also indicated balancing selection using MK test, however, it did not give significant results when tested with P. coatneyi as an outgroup. Phylogenetic analysis of full-length genes identified three distinct sub-clusters of P. knowlesi, one originating from Peninsular Malaysia and two originating from Malaysian Borneo. High population differentiation values was observed within samples from Peninsular Malaysia and Malaysian Borneo.

Conclusions

This study is the first to report on the genetic diversity and natural selection of full-length pktrap. Low level of genetic diversity was found across the full-length gene of pktrap. Balancing selection of the von Willebrand factor domain indicated that TRAP could be a target in inducing immune response against P. knowlesi infections. However, higher number of samples would be necessary to further confirm the findings.

Electronic supplementary material

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

Naturally acquired IgG antibodies to thrombospondin-related anonymous protein of Plasmodium vivax (PvTRAP) in Thailand predominantly elicit immunological cross-reactivity

BACKGROUND

Thrombospondin-related anonymous protein (TRAP) is a prime candidate for a malaria vaccine. Antibodies to Plasmodium vivax TRAP (PvTRAP) occur upon natural infection while specific antigenic domains remain to be addressed.

METHODS

The PvTRAP sequences were determined from 73 P. vivax isolates from Tak and Ubon Ratchathani provinces collected in 2013. The recombinant proteins representing four variants each for domain II (A domain) and domain IV (thrombospondin repeat region) of PvTRAP circulating in these areas were used as antigens in enzyme-linked immunosorbent assay against 246 serum samples from P. vivax-infected patients in both provinces collected during 2013 and 2014.

RESULTS

The prevalence of total IgG antibodies to at least one variant antigen of domain II and domain IV was 63.8% and 71.5%, respectively. Differential IgG antibody responses to these variant antigens of each domain were observed. Total IgG antibody responses to the variant antigens of each domain upon pairwise comparisons were highly correlated, suggesting immunological cross-reactivity in the majority of serum samples. A smaller proportion of serum samples contained non-cross-reactive antibodies to variants of each domain; particularly domain II in which amino acid differences significantly influenced antibody recognition. Previous malaria exposure positively affected antibody responses to domain IV. Positive seroconversion and rising antibody titres occurred within a few weeks after resolution of infections.

CONCLUSIONS

Both domains II and IV are targets of naturally acquired IgG antibodies. Despite sequence variation in these domains, most antibody responses were cross-reactive. A cross-sectional evaluation of antibodies to PvTRAP during acute infection could underestimate the seroprevalence.

Polymorphism in merozoite surface protein-7E of Plasmodium vivax in Thailand: Natural selection related to protein secondary structure

Merozoite surface protein 7 (MSP-7) is a multigene family expressed during malaria blood-stage infection. MSP-7 forms complex with MSP-1 prior to merozoite egress from erythrocytes, and could affect merozoite invasion of erythrocytes. To characterize sequence variation in the orthologue in P. vivax (PvMSP-7), a gene member encoding PvMSP-7E was analyzed among 92 Thai isolates collected from 3 major endemic areas of Thailand (Northwest: Tak, Northeast: Ubon Ratchathani, and South: Yala and Narathiwat provinces). In total, 52 distinct haplotypes were found to circulate in these areas. Although population structure based on this locus was observed between each endemic area, no genetic differentiation occurred between populations collected from different periods in the same endemic area, suggesting spatial but not temporal genetic variation. Sequence microheterogeneity in both N- and C- terminal regions was predicted to display 4 and 6 α-helical domains, respectively. Signals of purifying selection were observed in α-helices II-X, suggesting structural or functional constraint in these domains. By contrast, α-helix-I spanning the putative signal peptide was under positive selection, in which amino acid substitutions could alter predicted CD4+ T helper cell epitopes. The central region of PvMSP-7E comprised the 5’-trimorphic and the 3’-dimorphic subregions. Positive selection was identified in the 3’ dimorphic subregion of the central domain. A consensus of intrinsically unstructured or disordered protein was predicted to encompass the entire central domain that contained a number of putative B cell epitopes and putative protein binding regions. Evidences of intragenic recombination were more common in the central region than the remainders of the gene. These results suggest that the extent of sequence variation, recombination events and selective pressures in the PvMSP-7E locus seem to be differentially affected by protein secondary structure.

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.

Plasmodium vivax Cell Traversal Protein for Ookinetes and Sporozoites (PvCelTOS) gene sequence and potential epitopes are highly conserved among isolates from different regions of Brazilian Amazon

The Plasmodium vivax Cell-traversal protein for ookinetes and sporozoites (PvCelTOS) plays an important role in the traversal of host cells. Although essential to PvCelTOS progress as a vaccine candidate, its genetic diversity remains uncharted. Therefore, we investigated the PvCelTOS genetic polymorphism in 119 field isolates from five different regions of Brazilian Amazon (Manaus, Novo Repartimento, Porto Velho, Plácido de Castro and Oiapoque). Moreover, we also evaluated the potential impact of non-synonymous mutations found in the predicted structure and epitopes of PvCelTOS. The field isolates showed high similarity (99.3% of bp) with the reference Sal-1 strain, presenting only four Single-Nucleotide Polymorphisms (SNP) at positions 24A, 28A, 109A and 352C. The frequency of synonymous C109A (82%) was higher than all others (p<0.0001). However, the non-synonymous G28A and G352C were observed in 9.2% and 11.7% isolates. The great majority of the isolates (79.8%) revealed complete amino acid sequence homology with Sal-1, 10.9% presented complete homology with Brazil I and two undescribed PvCelTOS sequences were observed in 9.2% field isolates. Concerning the prediction analysis, the N-terminal substitution (Gly10Ser) was predicted to be within a B-cell epitope (PvCelTOS Accession Nos. AB194053.1) and exposed at the protein surface, while the Val118Leu substitution was not a predicted epitope. Therefore, our data suggest that although G28A SNP might interfere in potential B-cell epitopes at PvCelTOS N-terminal region the gene sequence is highly conserved among the isolates from different geographic regions, which is an important feature to be taken into account when evaluating its potential as a vaccine candidate.

Cell-traversal protein for ookinetes and sporozoites (CelTOS) presents a pivotal role in the cell traversal of host cells in mosquito and vertebrate hosts. For this reason, it has been considered a potential novel alternative for a vaccine against malaria caused by P. falciparum. However, little is known about its orthologous P. vivax CelTOS. Although the genetic diversity of this protein could be a limiting factor for acquisition of immunity and present implications for an effective vaccine development, it has never been explored. Thus, considering that the epidemiology of malaria in Brazil presents variable transmission rates and the knowledge on the genetic polymorphism of PvCelTOS remains unknown, we aimed to identify the pvceltos gene in isolates from five different regions of the Brazilian Amazon and to study the potential impacts of the genetic diversity of PvCelTOS in protein structures and predicted epitopes. Our findings indicate that PvCelTOS is an extremely conserved protein, presenting only four SNPs in the entire sequences of field isolates from Brazilian Amazon. The two non-synonymous mutations found in our field isolates presented no significant effect on the protein structure and a very low impact on potential T and B-cell epitopes indicated by our epitope prediction. Collectively, our data suggest that the small need to avoid the immune recognition by the human host and its importance on the parasite’s survival and transmission reflects a very conservative profile of pvceltos gene in field samples from Brazil and other endemic areas worldwide.

A comparative study on worldwide genetic diversity and population structure analysis of Plasmodium vivax thrombospondin-related adhesive protein (PvTRAP) and its implications for the vivax vaccine design

Plasmodium vivax thrombospondin-related anonymous protein (PvTRAP) is a promising malaria vaccine candidate; however, it exhibits sequence heterogeneity. Therefore, to design a broadly protective vivax vaccine, it is essential to have adequate information on signatures of selection and geospatial genetic diversity of global PvTRAP. For this purpose, 50 Iranian pvtrap were sequenced and compared with related available global sequences in GenBank. The nucleotide sequence analysis of Iranian pvtrap in comparison with the Sal-1 sequence showed the occurrence of 15 SNPs, and all sites were dimorphic. In total, 12 amino acid substitutions were detected and 2 of which were novel, resulting in 10 haplotypes that 8 of them were not reported in any other geographic regions. In comparison with global population, haplotype and nucleotide diversities were lowest in South Korean populations while higher levels of diversities were observed in Thai and Brazilian P. vivax populations. All 12 amino acid replacements in ectodomain of Iranian PvTRAP were distributed in predicted either B- or T-cells epitope as well as intrinsically unstructured/disordered regions (IURs). The present results revealed that observing the relatively low-level diversity in PvTRAP protein might actually be selected by immune response. In summary, the present analysis in parallel to the limited available published data has shown that genetic diversity in the global pvtrap exhibits low-level diversity and geographic variation. These results are of practical significance for the strategic development and deployment of control measures in particular for development of PvTRAP-based malaria vaccine.