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.

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.

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.