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.
RESULTS: Six EV-71 isolates from an outbreak in Malaysia, in 1997, were sequenced completely. These isolates were identified as EV-71 subgenotypes, B3, B4 and C2. A phylogenetic tree that correlated well with the present enterovirus classification scheme was established using these full genome sequences and all other available full genome sequences of EV-71 and human enterovirus A (HEV-A). Using the 5' UTR, P2 and P3 genomic regions, however, isolates of EV-71 subgenotypes B3 and C4 segregated away from other EV-71 subgenotypes into a cluster together with coxsackievirus A16 (CV-A16/G10) and EV-71 subgenotype C2 clustered with CV-A8. Results from the similarity plot analyses supported the clustering of these isolates with other HEV-A. In contrast, at the same genomic regions, a CV-A16 isolate, Tainan5079, clustered with EV-71. This suggests that amongst EV-71 and CV-A16, only the structural genes were conserved. The 3' end of the virus genome varied and consisted of sequences highly similar to various HEV-A viruses. Numerous recombination crossover breakpoints were identified within the non-structural genes of some of these newer EV-71 subgenotypes.
CONCLUSION: Phylogenetic evidence obtained from analyses of the full genome sequence supports the possible occurrence of inter-typic recombination involving EV-71 and various HEV-A, including CV-A16, the most common causal agent of HFMD. It is suggested that these recombination events played important roles in the emergence of the various EV-71 subgenotypes.