METHODS AND RESULTS: The library of 16S rDNA V3-V4 hypervariable regions of gut microbiota was amplified and sequenced using Illumina MiSeq. The sequencing data were analyzed using Quantitative Insights into Microbial Ecology (QIIME) pipeline and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt). The most abundant bacterial phyla in both wild and captive T. tambroides were Firmicutes, Proteobacteria, Fusobacteria and Bacteroidetes. Cetobacterium spp., Peptostreptococcaceae family, Bacteroides spp., Phosphate solubilizing bacteria PSB-M-3, and Vibrio spp. were five most abundant OTU in wild T. tambroides as compared to Cetobacterium spp., Citrobacter spp., Aeromonadaceae family, Peptostreptococcaceae family and Turicibacter spp. in captive T. tambroides.
CONCLUSION: In this study, the specimens of the wild T. tambroides contain more diverse gut microbiota than of the captive ones. The results suggested that Cetobacterium spp. is one of the core microbiota in guts of T. tambroides. Besides, high abundant Bacteroides spp., Citrobacter spp., Turicibacter spp., and Bacillus spp. may provide important functions in T. tambroides guts.
SIGNIFICANCE AND IMPACT OF THE STUDY: The results of this study provide significant information of T. tambroides gut microbiota for further understanding of their physiological functions including growth and disease resistance.
Method: Thirty-five full-length pk41 sequences from clinical isolates of Malaysia along with four laboratory lines (along with H-strain) were downloaded from public databases. For comparative analysis between species, orthologous P41 genes from P. falciparum, P. vivax, P. coatneyi and P. cynomolgi were also downloaded. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. Phylogenetic relationships between Pk41 genes were determined using MEGA 5.0 software.
Results: Analysis of 39 full-length pk41 sequences along with the H-strain identified 36 SNPs (20 non-synonymous and 16 synonymous substitutions) resulting in 31 haplotypes. Nucleotide diversity across the full-length gene was low and was similar to its ortholog in P. vivax; pv41. Domain-wise amino acid analysis of the two s48/45 domains indicated low level of polymorphisms for both the domains, and the glutamic acid rich region had extensive size variations. In the central domain, upstream to the glutamate rich region, a unique two to six (K-E)n repeat region was identified within the clinical isolates. Overall, the pk41 genes were indicative of negative/purifying selection due to functional constraints. Domain-wise analysis of the s48/45 domains also indicated purifying selection. However, analysis of Tajima's D across the genes identified non-synonymous SNPs in the s48/45 domain II with high positive values indicating possible epitope binding regions. All the 6-cysteine residues within the s48/45 domains were conserved within the clinical isolates indicating functional conservation of these regions. Phylogenetic analysis of full-length pk41 genes indicated geographical clustering and identified three subpopulations of P. knowlesi; one originating in the laboratory lines and two originating from Sarawak, Malaysian Borneo.
Conclusion: This is the first study to report on the polymorphism and natural selection of pk41 genes from clinical isolates of Malaysia. The results reveal that there is low level of polymorphism in both s48/45 domains, indicating that this antigen could be a potential vaccine target. However, genetic and molecular immunology studies involving higher number of samples from various parts of Malaysia would be necessary to validate this antigen's candidacy as a vaccine target for P. knowlesi.