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

In Silico Analysis Reveals High Levels of Genetic Diversity of Plasmodium knowlesi Cell Traversal Protein for Ookinetes and Sporozoites (PkCelTOS) in Clinical Samples

The cell-traversal protein for ookinetes and sporozoites (CelTOS), expressed on the surface of ookinetes and sporozoitesin Plasmodium species, is a promising malaria vaccine candidate. CelTOS is essential for parasite invasion into mosquito midgut and human hepatocytes, thereby contributing to malaria transmission and disease pathogenesis. This study explores the genetic diversity, polymorphisms, haplotypes, natural selection, phylogenetic analysis, and epitope prediction in the full-length Plasmodium knowlesi CelTOS gene in clinical samples from Sarawak, Malaysian Borneo, and long-term laboratory strains from Peninsular Malaysia and the Philippines. Our analysis revealed a high level of genetic variation in the PkCelTOS gene, with a nucleotide diversity of π ~ 0.021, which was skewed towards the 3' end of the gene. This level of diversity is double that observed in PfCelTOS and 20 times that observed in PvCelTOS from worldwide clinical samples. Tests of natural selection revealed evidence for positive selection within clinical samples. Phylogenetic analysis of the amino acid sequence of PkCelTOS revealed the presence of two distinct groups, although no geographical clustering was observed. Epitope prediction analysis identified two potential epitopes (96AQLKATA102 and 124TIKPPRIKED133) using the IEDB server and one epitope (125IKPPRIKED133) by Bcepred server on the C' terminal region of PkCelTOS protein. Both the servers predicted a common epitope region of nine amino acid length (IKPPRIKED) peptide, which can be studied in the future as a potential candidate for vaccine development. These findings shed light on the genetic diversity, polymorphism, haplotypes, and natural selection within PkCelTOS in clinical samples and provide insights about its future prospects as a potential candidate for P. knowlesi malaria vaccine development.

Inter-Population Genetic Diversity and Clustering of Merozoite Surface Protein-1 (pkmsp-1) of Plasmodium knowlesi Isolates from Malaysia and Thailand

The genetic diversity of pkmsp-1 of Malaysian Plasmodium knowlesi isolates was studied recently. However, the study only included three relatively older strains from Peninsular Malaysia and focused mainly on the conserved blocks of this gene. In this study, the full-length pkmsp-1 sequence of recent P. knowlesi isolates from Peninsular Malaysia was characterized, along with Malaysian Borneo and Thailand pkmsp-1 sequences that were retrieved from GenBank. Genomic DNA of P. knowlesi was extracted from human blood specimens and the pkmsp-1 gene was PCR-amplified, cloned, and sequenced. The sequences were analysed for genetic diversity, departure from neutrality, and geographical clustering. The pkmsp-1 gene was found to be under purifying/negative selection and grouped into three clusters via a neighbour-joining tree and neighbour net inferences. Of the four polymorphic blocks in pkmsp-1, block IV, was most polymorphic, with the highest insertion-deletion (indel) sites. Two allelic families were identified in block IV, thereby highlighting the importance of this block as a promising genotyping marker for the multiplicity of infection study of P. knowlesi malaria. A single locus marker may provide an alternate, simpler method to type P. knowlesi in a population.

Genetic Diversity and Population Genetic Structure Analysis of Plasmodium knowlesi Thrombospondin-Related Apical Merozoite Protein (TRAMP) in Clinical Samples

The simian malaria parasite Plasmodium knowlesi causes a high number of zoonotic infections in Malaysia. The thrombospondin-related apical merozoite protein (TRAMP) is an essential ligand for binding to the erythrocyte cell surface, whereby it facilitates the invasion. This study is the first attempt to determine the genetic diversity, phylogeography, natural selection and population structure from 97 full-length PkTRAMP gene sequences originating from Malaysia. We found low levels of nucleotide diversity (π~0.0065) for the full-length gene despite samples originating from geographically separated regions (i.e., Peninsular Malaysia and Malaysian Borneo). The rate of synonymous substitutions was significantly higher than that of non-synonymous substitutions, indicating a purifying selection for the full-length gene within the clinical samples. The population genetic analysis revealed that the parasite population is undergoing a significant population expansion. The analysis of the amino acid sequence alignment of 97 PkTRAMP sequences identified 15 haplotypes, of which a major shared haplotype was noted Hap 1 (n = 68, Sarawak; n = 34, Sabah; n = 12, Peninsular Malaysia; n = 22). The phylogenetic analysis using DNA sequences identified two clusters that separated due to geographical distance and three mixed clusters with samples from both Peninsular Malaysia and Malaysian Borneo. Population structure analyses indicated two distinct sub-populations (K = 2). Our findings point to the potential for independent parasite evolution, which could make zoonotic malaria control and elimination even more challenging.

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.

Plasmodium knowlesi: the game changer for malaria eradication

Plasmodium knowlesi is a zoonotic malaria parasite that has gained increasing medical interest over the past two decades. This zoonotic parasitic infection is prevalent in Southeast Asia and causes many cases with fulminant pathology. Despite several biogeographical restrictions that limit its distribution, knowlesi malaria cases have been reported in different parts of the world due to travelling and tourism activities. Here, breakthroughs and key information generated from recent (over the past five years, but not limited to) studies conducted on P. knowlesi were reviewed, and the knowledge gap in various research aspects that need to be filled was discussed. Besides, challenges and strategies required to control and eradicate human malaria with this emerging and potentially fatal zoonosis were described.

Is there evidence of sustained human-mosquito-human transmission of the zoonotic malaria Plasmodium knowlesi? A systematic literature review

Background

The zoonotic malaria parasite Plasmodium knowlesi has emerged across Southeast Asia and is now the main cause of malaria in humans in Malaysia. A critical priority for P. knowlesi surveillance and control is understanding whether transmission is entirely zoonotic or is also occurring through human-mosquito-human transmission.

Methods

A systematic literature review was performed to evaluate existing evidence which refutes or supports the occurrence of sustained human-mosquito-human transmission of P. knowlesi. Possible evidence categories and study types which would support or refute non-zoonotic transmission were identified and ranked. A literature search was conducted on Medline, EMBASE and Web of Science using a broad search strategy to identify any possible published literature. Results were synthesized using the Synthesis Without Meta-analysis (SWiM) framework, using vote counting to combine the evidence within specific categories.

Results

Of an initial 7,299 studies screened, 131 studies were included within this review: 87 studies of P. knowlesi prevalence in humans, 14 studies in non-human primates, 13 studies in mosquitoes, and 29 studies with direct evidence refuting or supporting non-zoonotic transmission. Overall, the evidence showed that human-mosquito-human transmission is biologically possible, but there is limited evidence of widespread occurrence in endemic areas. Specific areas of research were identified that require further attention, notably quantitative analyses of potential transmission dynamics, epidemiological and entomological surveys, and ecological studies into the sylvatic cycle of the disease.

Conclusion

There are key questions about P. knowlesi that remain within the areas of research that require more attention. These questions have significant implications for malaria elimination and eradication programs. This paper considers limited but varied research and provides a methodological framework for assessing the likelihood of different transmission patterns for emerging zoonotic diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04110-z.

Identification, Mapping, and Genetic Diversity of Novel Conserved Cross-Species Epitopes of RhopH2 in Plasmodium knowlesi With Plasmodium vivax

Malaria is a major public health concern, and any tangible intervention during the pre-elimination phase can result in a significant reduction in infection rates. Recent studies have reported that antigens producing cross-protective immunity can play an important role as vaccines and halt malaria transmission in different endemic regions. In this study, we studied the genetic diversity, natural selection, and discovered novel conserved epitopes of a high molecular weight rhoptry protein 2 (RhopH2) in clinical samples of Plasmodium knowlesi and Plasmodium vivax cross-protective domains, which has been proven to produce cross-protective immunity in both species. We found low levels of nucleotide diversity (P. knowlesi; π ~ 0.0093, SNPs = 49 and P. vivax π ~ 0.0014, SNPs = 23) in P. knowlesi (n = 40) and P. vivax (n = 65) samples in the PkRhopH2 cross-protective domain. Strong purifying selection was observed for both species (P. knowlesi; dS - dN = 2.41, p < 0.009, P. vivax; dS - dN = 1.58, p < 0.050). In silico epitope prediction in P. knowlesi identified 10 potential epitopes, of which 7 epitopes were 100% conserved within clinical samples. Of these epitopes, an epitope with 10 amino acids (QNSKHFKKEK) was found to be fully conserved within all P. knowlesi and P. vivax clinical samples and 80%–90% conservation within simian malaria ortholog species, i.e., P. coatneyi and P. cynomolgi. Phylogenetic analysis of the PkRhopH2 cross-protective domain showed geographical clustering, and three subpopulations of P. knowlesi were identified of which two subpopulations originated from Sarawak, Malaysian Borneo, and one comprised only the laboratory lines from Peninsular Malaysia. This study suggests that RhopH2 could be an excellent target for cross-protective vaccine development with potential for outwitting strain as well as species-specific immunity. However, more detailed studies on genetic diversity using more clinical samples from both species as well as the functional role of antibodies specific to the novel conserved epitope identified in this study can be explored for protection against infection.

Diversity pattern of Plasmodium knowlesi merozoite surface protein 4 (MSP4) in natural population of Malaysia

Human infections due to the monkey malaria parasite Plasmodium knowlesi are increasingly being reported from Malaysia. The parasite causes high parasitaemia, severe and fatal malaria in humans thus there is a need for urgent measures for its control. The MSP4 is a potential vaccine candidate, which is well studied in Plasmodium falciparum and Plasmodium vivax; however, no study has been conducted in the orthologous gene of P. knowlesi. In this study, we investigated the level of polymorphisms, haplotypes, natural selection and population structure of full-length pkmsp4 in 32 clinical samples from Malaysian Borneo along with 4 lab-adapted strains. We found low levels of polymorphism across the gene with exon I showing higher diversity than the exon II. The C- terminal epidermal growth factor (EGF) domains and GPI-anchored region within exon II were mostly conserved with only 2 non-synonymous substitutions. Although 21 amino acid haplotypes were found, the frequency of mutation at the majority of the polymorphic positions was low. We found evidence of negative selection at the exon II of the gene indicating existence of functional constraints. Phylogenetic haplotype network analysis identified shared haplotypes and indicated geographical clustering of samples originating from Peninsular Malaysia and Malaysian Borneo. High population differentiation values were observed within parasite populations originating from Malaysian Borneo (Kapit, Sarikei and Betong) and laboratory-adapted strains obtained from Peninsular Malaysia and Philippines indicating distinct population structure. This is the first study to genetically characterize the full-length msp4 gene from clinical isolates of P. knowlesi from Malaysia and thus would be very useful for future rational vaccine studies. Further studies with higher number of samples and functional characterization of the protein will be necessary.

Plasmodium knowlesi clinical isolates from Malaysia show extensive diversity and strong differential selection pressure at the merozoite surface protein 7D (MSP7D)

BACKGROUND

The high proportion of human cases due to the simian malaria parasite Plasmodium knowlesi in Malaysia is a cause of concern, as they can be severe and even fatal. Merozoite surface protein 7 (MSP7) is a multigene family which forms a non-covalent complex with MSP-1 prior to receptor-ligand recognition in Plasmodium falciparum and thus an important antigen for vaccine development. However, no study has been done in any of the ortholog family members in P. knowlesi from clinical samples. This study investigates the level of polymorphism, haplotypes, and natural selection acting at the pkmsp-7D gene in clinical samples from Malaysia.

METHODS

Thirty-six full-length pkmsp7D gene sequences (along with the reference H-strain: PKNH_1266000) obtained from clinical isolates of Malaysia, which were orthologous to pvmsp7H (PVX_082680) were downloaded from public databases. Population genetic, evolutionary and phylogenetic analyses were performed to determine the level of genetic diversity, polymorphism, recombination and natural selection.

RESULTS

Analysis of 36 full-length pkmsp7D sequences identified 147 SNPs (91 non-synonymous and 56 synonymous substitutions). Nucleotide diversity across the full-length gene was higher than its ortholog in Plasmodium vivax (msp7H). Region-wise analysis of the gene indicated that the nucleotide diversity at the central region was very high (π = 0.14) compared to the 5' and 3' regions. Most hyper-variable SNPs were detected at the central domain. Multiple test for natural selection indicated the central region was under strong positive natural selection however, the 5' and 3' regions were under negative/purifying selection. Evidence of intragenic recombination were detected at the central region of the gene. Phylogenetic analysis using full-length msp7D genes indicated there was no geographical clustering of parasite population.

CONCLUSIONS

High genetic diversity with hyper-variable SNPs and strong evidence of positive natural selection at the central region of MSP7D indicated exposure of the region to host immune pressure. Negative selection at the 5' and the 3' regions of MSP7D might be because of functional constraints at the unexposed regions during the merozoite invasion process of P. knowlesi. No evidence of geographical clustering among the clinical isolates from Malaysia indicated uniform selection pressure in all populations. These findings highlight the further evaluation of the regions and functional characterization of the protein as a potential blood stage vaccine candidate for P. knowlesi.

The Plasmodium knowlesi Pk41 surface protein diversity, natural selection, sub population and geographical clustering: a 6-cysteine protein family member

Introduction

The zoonotic malaria parasite Plasmodium knowlesi has currently become the most dominant form of infection in humans in Malaysia and is an emerging infectious disease in most Southeast Asian countries. The P41 is a merozoite surface protein belonging to the 6-cysteine family and is a well-characterized vaccine candidate in P. vivax and P. falciparum; however, no study has been done in the orthologous gene of P. knowlesi. This study investigates the level of polymorphism, haplotypes and natural selection of pk41 genes in clinical isolates from Malaysia.

Method

Thirty-five full-length pk41 sequences from clinical isolates of Malaysia along with four laboratory lines (along with H-strain) were downloaded from public databases. For comparative analysis between species, orthologous P41 genes from P. falciparum, P. vivax, P. coatneyi and P. cynomolgi were also downloaded. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. Phylogenetic relationships between Pk41 genes were determined using MEGA 5.0 software.

Results

Analysis of 39 full-length pk41 sequences along with the H-strain identified 36 SNPs (20 non-synonymous and 16 synonymous substitutions) resulting in 31 haplotypes. Nucleotide diversity across the full-length gene was low and was similar to its ortholog in P. vivax; pv41. Domain-wise amino acid analysis of the two s48/45 domains indicated low level of polymorphisms for both the domains, and the glutamic acid rich region had extensive size variations. In the central domain, upstream to the glutamate rich region, a unique two to six (K-E)_n repeat region was identified within the clinical isolates. Overall, the pk41 genes were indicative of negative/purifying selection due to functional constraints. Domain-wise analysis of the s48/45 domains also indicated purifying selection. However, analysis of Tajima’s D across the genes identified non-synonymous SNPs in the s48/45 domain II with high positive values indicating possible epitope binding regions. All the 6-cysteine residues within the s48/45 domains were conserved within the clinical isolates indicating functional conservation of these regions. Phylogenetic analysis of full-length pk41 genes indicated geographical clustering and identified three subpopulations of P. knowlesi; one originating in the laboratory lines and two originating from Sarawak, Malaysian Borneo.

Conclusion

This is the first study to report on the polymorphism and natural selection of pk41 genes from clinical isolates of Malaysia. The results reveal that there is low level of polymorphism in both s48/45 domains, indicating that this antigen could be a potential vaccine target. However, genetic and molecular immunology studies involving higher number of samples from various parts of Malaysia would be necessary to validate this antigen’s candidacy as a vaccine target for P. knowlesi.

Within-population genetic diversity and population structure of Plasmodium knowlesi merozoite surface protein 1 gene from geographically distinct regions of Malaysia and Thailand

BACKGROUND

The C-terminal 42 kDa domain of Plasmodium knowlesi merozoite surface protein 1 (PkMSP1) is a potential asexual blood-stage vaccine candidate, however, only a limited number of clinical isolates have been analysed from Malaysia and no inter-country comparative diversity study has been conducted. In the present study, nucleotide diversity, haplotypes and natural selection levels of pkmsp1 in clinical samples from geographically distinct regions of Malaysia and Thailand were investigated. The overall population structure of the parasite from the region was determined.

METHODS

Eleven full-length pkmsp1 sequences obtained from clinical isolates of Malaysia along with the H-strain were downloaded from the database for domain wise characterization of pkmsp1 gene. Additionally, 76 pkmsp-1₄₂ sequences from Thailand and Malaysia were downloaded from the database for intra and inter-population analysis. DnaSP 5.10 and MEGA 5.0 software were used to determine genetic diversity, polymorphism, haplotypes and natural selection. Genealogical relationships were determined using haplotype network tree in NETWORK software v5.0. Population genetic differentiation index (F_ST) of parasites were analysed using Arlequin v3.5.

RESULTS

Sequence analysis of 11 full-length pkmsp1 sequences along with the H-strain identified 477 (8.4%) polymorphic sites, of which 107 were singleton sites. The overall diversity observed in the full-length genes were high in comparison to its ortholog pvmsp1 and the 4 variable domains showed extensive size variations. The nucleotide diversity was low towards the pkmsp1-42 compared to the conserved domains. The 19 kDa domain was less diverse and completely conserved among isolates from Malaysian Borneo. The nucleotide diversity of isolates from Peninsular Malaysia and Thailand were higher than Malaysian Borneo. Network analysis of pkmsp1-42 haplotypes showed geographical clustering of the isolates from Malaysian Borneo and grouping of isolates from Peninsular Malaysia and Thailand. Population differentiation analysis indicated high F_ST values between parasite populations originating from Malaysian Borneo, Peninsular Malaysia and Thailand attributing to geographical distance. Moderate genetic differentiation was observed for parasite populations from Thailand and Peninsular Malaysia. Evidence of population expansion and purifying selection were observed in all conserved domains with strongest selection within the pkmsp1-42 domain.

CONCLUSIONS

This study is the first to report on inter country genetic diversity and population structure of P. knowlesi based on msp1. Strong evidence of negative selection was observed in the 42 kDa domain, indicating functional constrains. Geographical clustering of P. knowlesi and moderate to high genetic differentiation values between populations identified in this study highlights the importance of further evaluation using larger number of clinical samples from Southeast Asian countries.