Abstract:BackgroundThe comparison of samples, or beta diversity, is one of the essential problems in ecological studies. Next generation sequencing (NGS) technologies make it possible to obtain large amounts of metagenomic and metatranscriptomic short read sequences across many microbial communities. De novo assembly of the short reads can be especially challenging because the number of genomes and their sequences are generally unknown and the coverage of each genome can be very low, where the traditional alignment-based sequence comparison methods cannot be used. Alignment-free approaches based on k-tuple frequencies, on the other hand, have yielded promising results for the comparison of metagenomic samples. However, it is not known if these approaches can be used for the comparison of metatranscriptome datasets and which dissimilarity measures perform the best.ResultsWe applied several beta diversity measures based on k-tuple frequencies to real metatranscriptomic datasets from pyrosequencing 454 and Illumina sequencing platforms to evaluate their effectiveness for the clustering of metatranscriptomic samples, including three d2-type dissimilarity measures, one dissimilarity measure in CVTree, one relative entropy based measure S2 and three classical 1p-norm distances. Results showed that the measure d2(S) can achieve superior performance on clustering metatranscriptomic samples into different groups under different sequencing depths for both 454 and Illumina datasets, recovering environmental gradients affecting microbial samples, classifying coexisting metagenomic and metatranscriptomic datasets, and being robust to sequencing errors. We also investigated the effects of tuple size and order of the background Markov model. A software pipeline to implement all the steps of analysis is built and is available at http://code.google.com/p/d2-tools/.ConclusionsThe k-tuple based sequence signature measures can effectively reveal major groups and gradient variation among metatranscriptomic samples from NGS reads. The d2(S) dissimilarity measure performs well in all application scenarios and its performance is robust with respect to tuple size and order of the Markov model.

The Application of CVTree in Structural Analysis of Microbial Communities by 454 Pyrosequencing

Comparison of Microbial Diversity Determined with the Same Variable Tag Sequence Extracted from Two Different PCR Amplicons

CVTree: a Phylogenetic Tree Reconstruction Tool Based on Whole Genomes

Polyphyly in 16S rRNA-based LVTree Versus Monophyly in Whole-genome-based CVTree.

Prokaryote Phylogeny Meets Taxonomy: an Exhaustive Comparison of Composition Vector Trees with Systematic Bacteriology

Tree-aggregated predictive modeling of microbiome data

Whole-genome based Archaea phylogeny and taxonomy: A composition vector approach

A phylogenetic scan test on Dirichlet-tree multinomial model for microbiome data

An Introduction to Next Generation Sequencing Bioinformatic Analysis in Gut Microbiome Studies

Comparison of Metatranscriptomic Samples Based on K-Tuple Frequencies

Composition-based Classification of Short Metagenomic Sequences Elucidates the Landscapes of Taxonomic and Functional Enrichment of Microorganisms

A Fungal Phylogeny Based on 82 Complete Genomes Using the Composition Vector Method.

Comparative and phylogenetic analysis of chloroplast genomes from ten species in Quercus section Cyclobalanopsis

A Model for Paired-Multinomial Data and Its Application to Analysis of Data on a Taxonomic Tree

Effect of k-tuple length on sample-comparison with high-throughput sequencing data

High-resolution Single-Molecule Long-Fragment Rrna Gene Amplicon Sequencing of Bacterial and Eukaryotic Microbial Communities

VEGA: Visual Comparison of Phylogenetic Trees for Evolutionary Genome Analysis (chinavis 2019)

Plastome‐based phylogeny improves community phylogenetics of subtropical forests in China

CRP-Tree: A phylogenetic association test for binary traits

Exploration of Carcinogenesis Based on Tree Models Using CGH Data

A multi-amplicon 16S rRNA sequencing and analysis method for improved taxonomic profiling of bacterial communities