METHODS: We conducted molecular detection, genetic characterization, and Bayesian time-scale evolution analyses of NiV using pooled Pteropid bat roost urine samples from an outbreak area in 2012 and archived RNA samples from NiV case patients identified during 2012-2018 in Bangladesh.
RESULTS: NiV-RNA was detected in 19% (38/456) of bat roost urine samples and among them; nine N gene sequences were recovered. We also retrieved sequences from 53% (21 out of 39) of archived RNA samples from patients. Phylogenetic analysis revealed that all Bangladeshi strains belonged to NiV-BD genotype and had an evolutionary rate of 4.64 × 10-4 substitutions/site/year. The analyses suggested that the strains of NiV-BD genotype diverged during 1995 and formed two sublineages.
CONCLUSION: This analysis provides further evidence that the NiV strains of the Malaysian and Bangladesh genotypes diverged recently and continue to evolve. More extensive surveillance of NiV in bats and human will be helpful to explore strain diversity and virulence potential to infect humans through direct or person-to-person virus transmission.
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.