METHODS: A total of 46 csp genes were subjected to polymerase chain reaction amplification. The genes were obtained from P. knowlesi isolates collected from different divisions of Sabah, Malaysian Borneo, and Peninsular Malaysia. The targeted gene fragments were cloned into a commercial vector and sequenced, and a phylogenetic tree was constructed while incorporating 168 csp sequences retrieved from the GenBank database. The genetic diversity and natural evolution of the csp sequences were analysed using MEGA6 and DnaSP ver. 5.10.01. A genealogical network of the csp haplotypes was generated using NETWORK ver. 4.6.1.3.
RESULTS: The phylogenetic analysis revealed indistinguishable clusters of P. knowlesi isolates across different geographic regions, including Malaysian Borneo and Peninsular Malaysia. Nucleotide analysis showed that the csp non-repeat regions of zoonotic P. knowlesi isolates obtained in this study underwent purifying selection with population expansion, which was supported by extensive haplotype sharing observed between humans and macaques. Novel variations were observed in the C-terminal non-repeat region of csp.
CONCLUSIONS: The csp non-repeat regions are relatively conserved and there is no distinct cluster of P. knowlesi isolates from Malaysian Borneo and Peninsular Malaysia. Distinctive variation data obtained in the C-terminal non-repeat region of csp could be beneficial for the design and development of vaccines to treat P. knowlesi.
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 (FST) 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.
METHODS: Phylogenetic analysis of locally acquired P. knowlesi infections, based on circumsporozoite, small subunit ribosomal ribonucleic acid (SSU rRNA), merozoite surface protein 1 and COX1 gene targets, was performed. The results were compared with the published sequences of regional isolates from Malaysia and Thailand.
RESULTS: Phylogenetic analysis of the circumsporozoite, SSU rRNA and merozoite surface protein 1 gene sequences for regional P. knowlesi isolates showed no obvious differentiation that could be attributed to their geographical origin. However, COX1 gene analysis showed that it was possible to differentiate between Singapore-acquired P. knowlesi infections and P. knowlesi infections from Peninsular Malaysia and Sarawak, Borneo, Malaysia.
CONCLUSION: The ability to differentiate between locally acquired P. knowlesi infections and imported P. knowlesi infections has important utility for the monitoring of P. knowlesi malaria control programmes in Singapore.
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.
METHODS: rPkMSP-119 protein was heterologously expressed using Expresso® Solubility and Expression Screening System and competent E. cloni® 10G cells according to protocol. Three rounds of biopanning were performed on purified rPkMSP-119 to identify binding peptides towards rPkMSP-119 using Ph.D.™-12 random phage display library. Binding sites of the identified peptides to PkMSP-119 were in silico predicted using the CABS-dock web server.
RESULTS: Four phage peptide variants that bound to PkMSP-119 were identified after three rounds of biopanning, namely Pkd1, Pkd2, Pkd3 and Pkd4. The sequences of both Pkd1 and Pkd2 consist of a large number of histidine residues. Pkd1 showed positive binding signal with 6.1× vs. BSA control. Docking results showed that Pkd1 and Pkd2 were ideal binding peptides for PkMSP-119 .
CONCLUSION: We identified two novel binding peptides of PkMSP-119 , Pkd1 (HFPFHHHKLRAH) and Pkd2 (HPMHMLHKRQHG), through phage display. They provide a valuable starting point for the development of novel therapeutics.