METHODS: A total of 70 blood samples were collected from Macaca fascicularis dwelling in the forest of Hulu Selangor by the Department of Wildlife and National Parks Peninsular Malaysia, Kuala Lumpur, Malaysia. DNA was extracted using PureLink™ Genomic DNA Kits. Conventional and nested PCR were used to detect the genus and species of Plasmodium parasites respectively. In addition, phylogenetic analysis was carried out to confirm the species of Plasmodium parasites.

RESULTS: Thirty-five (50 %) of the 70 samples were positive for Plasmodium using genus-specific primers. These positive samples were then subjected to nested PCR targeting the 18S ribosomal RNA genes to detect all five simian malaria parasites: namely, P. knowlesi, Plasmodium inui, Plasmodium cynomolgi, Plasmodium fieldi, and Plasmodium coatneyi. All five species of simian malaria parasites were detected. Of these, P. inui was the predominant (65.7 %), followed by P. knowlesi (60 %), P. cynomolgi (51.4 %) P. coatneyi (45.7 %) and P. fieldi (2.9 %). A total of nine macaques had mono-infection with P. knowlesi (four), P. cynomolgi (two), P. coatneyi (two) and P. fieldi (one). Eleven of the macaques had dual infections while 12 had triple infections. Three macaques were infected with four species of Plasmodium. Molecular and phylogenetic analysis confirmed the five species of Plasmodium parasites.

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.

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.

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-142 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 (FST) 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 FST 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.

METHODS: A total of 36 full-length pkmsp1p sequences along with the reference H-strain and 40 C-terminal pkmsp1p sequences from clinical isolates of Malaysia were downloaded from published genomes. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 and MEGA 5.0 software. Genealogical relationships were determined using haplotype network tree in NETWORK software v5.0. Population genetic differentiation index (F ST ) and population structure of parasite was determined using Arlequin v3.5 and STRUCTURE v2.3.4 software.

RESULTS: Comparison of 36 full-length pkmsp1p sequences along with the H-strain identified 339 SNPs (175 non-synonymous and 164 synonymous substitutions). The nucleotide diversity across the full-length gene was low compared to its ortholog pvmsp1p. The nucleotide diversity was higher toward the N-terminal domains (pkmsp1p-83 and 30) compared to the C-terminal domains (pkmsp1p-38, 33 and 19). Phylogenetic analysis of full-length genes identified 2 distinct clusters of P. knowlesi from Malaysian Borneo. The 40 pkmsp1p-19 sequences showed low polymorphisms with 16 polymorphisms leading to 18 haplotypes. In total there were 10 synonymous and 6 non-synonymous substitutions and 12 cysteine residues were intact within the two EGF domains. Evidence of strong purifying selection was observed within the full-length sequences as well in all the domains. Shared haplotypes of 40 pkmsp1p-19 were identified within Malaysian Borneo haplotypes.