Molecular characterization and geographical distribution of Zika virus worldwide from 1947 to 2022.

OBJECTIVES: We conducted molecular characterization, demonstrated the geographical distribution of Zika virus (ZIKV) circulating worldwide from 1947 to 2022 and explored the potential genetic recombination site in the Thailand ZIKV genomes. METHODS: We constructed phylogenetic trees based on ZIKV coding sequences (CDS) and determined the geographical distribution of the representative viruses by genetic relationship and timeline. We determined genetic recombination among ZIKV and between ZIKV and other flaviviruses using similarity plot and bootscan analyzes, together with the phylogeny encompassing the CDS and eight subgenomic regions. RESULTS: The phylogenetic trees comprising 717 CDS showed two distinct African and Asian lineages. ZIKV in the African lineage formed two sublineages, and ZIKV in the Asian lineage diversified into the Asian and American sublineages. The 1966 Malaysian isolate was designated the prototype of the Asian sublineage and formed a node of only one member, while the newer viruses formed a distinct node. We detected no genetic recombination in the Thailand ZIKV. CONCLUSION: Five Thailand isolates discovered in 2006 were the second oldest ZIKV after the Malaysian prototype. Our result suggested two independent routes of ZIKV spread from Southeast Asia to Micronesia in 2007 and French Polynesia in 2013 before further spreading to South American countries.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae.

BACKGROUND: Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed. RESULTS: The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis. CONCLUSIONS: In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

Single-Strand Conformation Polymorphism Fingerprint Method for Dictyostelids.

Dictyostelid social amoebae are a highly diverse group of eukaryotic soil microbes that are valuable resources for biological research. Genetic diversity study of these organisms solely relies on molecular phylogenetics of the SSU rDNA gene, which is not ideal for large-scale genetic diversity study. Here, we designed a set of PCR-single-strand conformation polymorphism (SSCP) primers and optimized the SSCP fingerprint method for the screening of dictyostelids. The optimized SSCP condition required gel purification of the SSCP amplicons followed by electrophoresis using a 9% polyacrylamide gel under 4°C. We also tested the optimized SSCP procedure with 73 Thai isolates of dictyostelid that had the SSU rDNA gene sequences published. The SSCP fingerprint patterns were related to the genus-level taxonomy of dictyostelids, but the fingerprint dendrogram did not reflect the deep phylogeny. This method is rapid, cost-effective, and suitable for large-scale sample screening as compared with the phylogenetic analysis of the SSU rDNA gene sequences.

Phylogenetic analysis of two new complete genomes of the hepatitis E virus (HEV) genotype 3 from Thailand.

Hepatitis E virus (HEV) is a causative agent of acute viral hepatitis globally. Evolutionary phylogeny classifies the HEV into eight genotypes that correlate with the viral transmission. Only four genotypes have been proven to be responsible for transmission in humans. However, there has been no report on the genomics and genotyping of HEV in Thailand during the past ten years. Here, we identified the genotype distributions of the Thai isolates of HEV and we sequenced two HEV genomes. We screened for 18 Thai isolates of HEV from Siriraj Hospital in Bangkok, from 2014-2016. The HEV genomes were sequenced from the serum and feces of a patient. The results showed that all Thai isolates of HEV were identified as genotype 3 (HEV-3). The ORF2 and genome phylogenies suggested two subgenotypes, called 3.1 and 3.2. The Thai isolates of HEV were frequently found in the subgenotype 3.1. The genome sequences of the two Thai isolates of HEV from the serum and fecal samples of the same patient showed 91% nucleotide similarity with the HEV genotype 3. Comparisons between the HEV genome and the ORF2 phylogenies illustrated that the ORF2 tree can be used to identify HEV genotypes, but it has less phylogenetic power for the HEV evolution. The two new genome sequences of HEV-3 from Thailand could contribute valuable information to the HEV genome study. (226 words).

Dysregulation of microRNA in cholangiocarcinoma identified through a meta-analysis of microRNA profiling.

BACKGROUND: In the past decades, the potential of microRNA (miRNA) in cancer diagnostics and prognostics has gained a lot of interests. In this study, a meta-analysis was conducted upon the pooled miRNA microarray data of cholangiocarcinoma (CCA). AIM: To identify differentially expressed (DE) miRNAs and perform functional analyses in order to gain insights to understanding miRNA-target interactions involved in tumorigenesis pathways of CCA. METHODS: Raw data from 8 CCA miRNA microarray datasets, consisting of 443 samples in total, were integrated and statistically analyzed to identify DE miRNAs via comparison of levels of miRNA expression between CCA and normal bile duct samples using t-tests (P < 0.001). The 10-fold cross validation was performed in order to increase the robustness of the t-test results. RESULTS: Our data showed 70 up-regulated and 48 down-regulated miRNAs in CCA. Gene Ontology and pathway enrichment analyses revealed that mRNA targets of DE miRNAs were significantly involved in several biological processes. The most prominent dysregulated pathways included phosphatidylinositol-3 kinases/Akt, mitogen-activated protein kinase and Ras signaling pathways. CONCLUSION: DE miRNAs found in our meta-analysis revealed dysregulation in major cancer pathways involved in the development of CCA. These results indicated the necessity of understanding the miRNA-target interactions and the significance of dysregulated miRNAs in terms of diagnostics and prognostics of cancers.

An optimized genomic VCF workflow for precise identification of Mycobacterium tuberculosis cluster from cross-platform whole genome sequencing data.

Whole-genome sequencing (WGS) data allow for an inference of Mycobacterium tuberculosis (Mtb) clusters by using a pairwise genetic distance of ≤12 single nucleotide polymorphisms (SNPs) as a threshold. However, a problem of discrepancies in numbers of SNPs and genetic distance measurement is a great concern when combining WGS data from different next generation sequencing (NGS) platforms. We performed SNP variant calling on WGS data of 9 multidrug-resistant (MDR-TB), 3 extensively drug-resistant tuberculosis (XDR-TB) and a standard M. tuberculosis strain H37Rv from an Illumina/NextSeq500 and an Ion Torrent PGM. Variant calls were obtained using four different common variant calling tools, including Genome Analysis Toolkit (GATK) HaplotypeCaller (GATK-VCF workflow), GATK HaplotypeCaller and GenotypeGVCFs (GATK-GVCF workflow), SAMtools, and VarScan 2. Cross-platform pairwise SNP differences, minimum spanning networks and average nucleotide identity (ANI) were analysed to measure performance of the variant calling tools. Minimum pairwise SNP differences ranged from 2 to 14 SNPs when using GVCF workflow while maximum pairwise SNP differences ranged from 7 to 158 SNPs when using VarScan 2. ANI comparison between SNPs data from NextSeq500 and PGM of MDR-TB and XDR-TB showed maximum ANI of 99.7% and 99.0%, respectively, with GVCF workflow while the other SNP calling results showed lower ANI in a range of 98.6% to 95.1%. In this study, we suggest that the GVCF workflow showed the best performing variant caller to avoid cross-platform pairwise SNP differences.

A novel Ancestral Beijing sublineage of Mycobacterium tuberculosis suggests the transition site to Modern Beijing sublineages.

Global Mycobacterium tuberculosis population comprises 7 major lineages. The Beijing strains, particularly the ones classified as Modern groups, have been found worldwide, frequently associated with drug resistance, younger ages, outbreaks and appear to be expanding. Here, we report analysis of whole genome sequences of 1170 M. tuberculosis isolates together with their patient profiles. Our samples belonged to Lineage 1-4 (L1-L4) with those of L1 and L2 being equally dominant. Phylogenetic analysis revealed several new or rare sublineages. Differential associations between sublineages of M. tuberculosis and patient profiles, including ages, ethnicity, HIV (human immunodeficiency virus) infection and drug resistance were demonstrated. The Ancestral Beijing strains and some sublineages of L4 were associated with ethnic minorities while L1 was more common in Thais. L2.2.1.Ancestral 4 surprisingly had a mutation that is typical of the Modern Beijing sublineages and was common in Akha and Lahu tribes who have migrated from Southern China in the last century. This may indicate that the evolutionary transition from the Ancestral to Modern Beijing sublineages might be gradual and occur in Southern China, where the presence of multiple ethnic groups might have allowed for the circulations of various co-evolving sublineages which ultimately lead to the emergence of the Modern Beijing strains.

Characterization of Flagellotropic, Chi-Like Salmonella Phages Isolated from Thai Poultry Farms.

Despite a wealth of knowledge on Salmonella phages worldwide, little is known about poultry-associated Salmonella phages from Thailand. Here, we isolated 108 phages from Thai poultry farms that infect Salmonellaenterica serovar Typhimurium. Phages STm101 and STm118 were identified as temperate Siphoviridae phages. Genome sequencing and analyses revealed these phages share approximately 96% nucleotide sequence similarity to phage SPN19, a member of the Chi-like virus genus. PCR amplification of the gene encoding capsid protein E of the Chi-like phage was positive for 50% of phage isolates, suggesting a predominance of this phage type among the sampled poultry farms. In addition to the flagella, two phages required the lipopolysaccharide to infect and lyse Salmonella. Furthermore, phylogenomic analysis demonstrated that phages STm101 and STm118 formed a monophyletic clade with phages isolated from Western countries, but not from closer isolated phages from Korea. However, further investigation and more phage isolates are required to investigate possible causes for this geographic distribution.

Infection of multiple Mycobacterium tuberculosis strains among tuberculosis/human immunodeficiency virus co-infected patients: A molecular study in Myanmar.

Background: Appearance of Mycobacterium tuberculosis (MTB) in the sputum of a tuberculosis (TB)/human immunodeficiency virus (HIV) co-infected patient under treatment may indicate either failure or new infection. This study aims to evaluate whether TB treatment failure among TB/HIV co-infected patients is a real failure. Methods: A prospective cohort study was conducted among 566 TB/HIV co-infected patients who started TB treatment in 12 townships in the upper Myanmar. Among the 566 participants, 16 (2.8%) resulted in treatment failure. We performed a molecular study using mycobacterial interspersed repetitive-unit-variable number of tandem repeat (MIRU-VNTR) genotyping for them. The MIRU-VNTR profiles were analyzed using the web server, MIRU-VNTRplus. All data were entered into EpiData version 3.1 and analyzed using R version 3.4.3. Results: Among 16 failure patients, seven had incomplete laboratory results. Of the nine remaining patients, nobody had exactly the same MIRU-VNTR pattern between the initial and final isolates. Four patients had persistent East-African Indian (EAI) lineages and one each had persistent Beijing lineage, changing from EAI to Beijing, from Beijing to EAI, NEW-1 to Beijing, and NEW-1 to X strains. Female patients have significantly larger genetic difference between MTB of the paired isolates than male patients (t-test, P = 0.04). Conclusion: Thus, in our study patients, infection of multiple MTB strains is a possible cause of TB treatment failure. Explanation for the association between gender and distance of genotypes from the initial to subsequent MTB infection needs further studies.

Transcriptomic analysis of male and female Schistosoma mekongi adult worms.

BACKGROUND: Schistosoma mekongi is one of five major causative agents of human schistosomiasis and is endemic to communities along the Mekong River in southern Lao People's Democratic Republic (Laos) and northern Cambodia. Sporadic cases of schistosomiasis have been reported in travelers and immigrants who have visited endemic areas. Schistosoma mekongi biology and molecular biology is poorly understood, and few S. mekongi gene and transcript sequences are available in public databases. RESULTS: Transcriptome sequencing (RNA-Seq) of male and female S. mekongi adult worms (a total of three biological replicates for each sex) were analyzed and the results demonstrated that approximately 304.9 and 363.3 million high-quality clean reads with quality Q30 (> 90%) were obtained from male and female adult worms, respectively. A total of 119,604 contigs were assembled with an average length of 1273 nt and an N50 of 2017 nt. From the contigs, 20,798 annotated protein sequences and 48,256 annotated transcript sequences were obtained using BLASTP and BLASTX searches against the UniProt Trematoda database. A total of 4658 and 3509 transcripts were predominantly expressed in male and female worms, respectively. Male-biased transcripts were mostly involved in structural organization while female-biased transcripts were typically involved in cell differentiation and egg production. Interestingly, pathway enrichment analysis suggested that genes involved in the phosphatidylinositol signaling pathway may play important roles in the cellular processes and reproductive systems of S. mekongi worms. CONCLUSIONS: We present comparative transcriptomic analyses of male and female S. mekongi adult worms, which provide a global view of the S. mekongi transcriptome as well as insights into differentially-expressed genes associated with each sex. This work provides valuable information and sequence resources for future studies of gene function and for ongoing whole genome sequencing efforts in S. mekongi.

Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates.

Tuberculosis presents a global health challenge. Mycobacterium tuberculosis is divided into several lineages, each with a different geographical distribution. M. tuberculosis lineage 1 (L1) is common in the high-burden areas in East Africa and Southeast Asia. Although the founder effect contributes significantly to the phylogeographic profile, co-evolution between the host and M. tuberculosis may also play a role. Here, we reported the genomic analysis of 480 L1 isolates from patients in northern Thailand. The studied bacterial population was genetically diverse, allowing the identification of a total of 18 sublineages distributed into three major clades. The majority of isolates belonged to L1.1 followed by L1.2.1 and L1.2.2. Comparison of the single nucleotide variant (SNV) phylogenetic tree and the clades defined by spoligotyping revealed some monophyletic clades representing EAI2_MNL, EAI2_NTM and EAI6_BGD1 spoligotypes. Our work demonstrates that ambiguity in spoligotype assignment could be partially resolved if the entire DR region is investigated. Using the information to map L1 diversity across Southeast Asia highlighted differences in the dominant strain-types in each individual country, despite extensive interactions between populations over time. This finding supported the hypothesis that there is co-evolution between the bacteria and the host, and have implications for tuberculosis disease control.

Molecular Phylogenetics: Concepts for a Newcomer.

Molecular phylogenetics is the study of evolutionary relationships among organisms using molecular sequence data. The aim of this review is to introduce the important terminology and general concepts of tree reconstruction to biologists who lack a strong background in the field of molecular evolution. Some modern phylogenetic programs are easy to use because of their user-friendly interfaces, but understanding the phylogenetic algorithms and substitution models, which are based on advanced statistics, is still important for the analysis and interpretation without a guide. Briefly, there are five general steps in carrying out a phylogenetic analysis: (1) sequence data preparation, (2) sequence alignment, (3) choosing a phylogenetic reconstruction method, (4) identification of the best tree, and (5) evaluating the tree. Concepts in this review enable biologists to grasp the basic ideas behind phylogenetic analysis and also help provide a sound basis for discussions with expert phylogeneticists.

Evolution of protein indels in plants, animals and fungi.

BACKGROUND: Insertions/deletions (indels) in protein sequences are useful as drug targets, protein structure predictors, species diagnostics and evolutionary markers. However there is limited understanding of indel evolutionary patterns. We sought to characterize indel patterns focusing first on the major groups of multicellular eukaryotes. RESULTS: Comparisons of complete proteomes from a taxonically broad set of primarily Metazoa, Fungi and Viridiplantae yielded 299 substantial (>250aa) universal, single-copy (in-paralog only) proteins, from which 901 simple (present/absent) and 3,806 complex (multistate) indels were extracted. Simple indels are mostly small (1-7aa) with a most frequent size class of 1aa. However, even these simple looking indels show a surprisingly high level of hidden homoplasy (multiple independent origins). Among the apparently homoplasy-free simple indels, we identify 69 potential clade-defining indels (CDIs) that may warrant closer examination. CDIs show a very uneven taxonomic distribution among Viridiplante (13 CDIs), Fungi (40 CDIs), and Metazoa (0 CDIs). An examination of singleton indels shows an excess of insertions over deletions in nearly all examined taxa. This excess averages 2.31 overall, with a maximum observed value of 7.5 fold. CONCLUSIONS: We find considerable potential for identifying taxon-marker indels using an automated pipeline. However, it appears that simple indels in universal proteins are too rare and homoplasy-rich to be used for pure indel-based phylogeny. The excess of insertions over deletions seen in nearly every genome and major group examined maybe useful in defining more realistic gap penalties for sequence alignment. This bias also suggests that insertions in highly conserved proteins experience less purifying selection than do deletions.

SeqFIRE: a web application for automated extraction of indel regions and conserved blocks from protein multiple sequence alignments.

Analyses of multiple sequence alignments generally focus on well-defined conserved sequence blocks, while the rest of the alignment is largely ignored or discarded. This is especially true in phylogenomics, where large multigene datasets are produced through automated pipelines. However, some of the most powerful phylogenetic markers have been found in the variable length regions of multiple alignments, particularly insertions/deletions (indels) in protein sequences. We have developed Sequence Feature and Indel Region Extractor (SeqFIRE) to enable the automated identification and extraction of indels from protein sequence alignments. The program can also extract conserved blocks and identify fast evolving sites using a combination of conservation and entropy. All major variables can be adjusted by the user, allowing them to identify the sets of variables most suited to a particular analysis or dataset. Thus, all major tasks in preparing an alignment for further analysis are combined in a single flexible and user-friendly program. The output includes a numbered list of indels, alignments in NEXUS format with indels annotated or removed and indel-only matrices. SeqFIRE is a user-friendly web application, freely available online at www.seqfire.org/.

Distribution and sequence analysis of virulence associated genes in Vibrio cholerae O1, O139 and non-O1/non-O139 isolates from Thailand.

Virulence-associated genes of Vibrio cholerae including O1, O139 and non-O1/non-O139 from an outbreak in Songkhla Province and sporadic cases occurred in Thailand during 1993 - 2002 were investigated. One hundred eighty-five V. cholerae strains were examined for the presence of virulence-associated genes including ctxA, tcpA, zot, toxR, toxS, toxT, and ace by polymerase chain reaction. DNA sequences of ctxA, tcpA, zot and toxR also were investigated in 8 selected isolates. Results showed that the virulence factors genes were distributed in the majority of V cholerae O1 biotype Inaba and Ogawa strains (85-97%). All 6 strains of the O139 harbored toxR, toxS and toxT whereas ctxA, tcpA, zot and ace were detected only 50-67%. Toxins genes found in non-O1/non-O139 strains ranged 8-30% except toxR (73.5%). Results of multiple sequence alignments among the isolates compared with V. cholerae O1 in database (embl M21249), showed that ctxA, tcpA and zot sequences in all 8 isolates were conserved, but base changes were found in toxR sequence. These molecular characteristics of V cholerae isolated from Thailand will provide detailed information for facilitating future studies on the development and design of appropriate vaccine providing protection against local strains.