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.

Local Alignment-Free Sequences Based on D2shepp Statistics

Alignment-free sequence comparison based on next-generation sequencing reads.

Alignment-Free Sequence Comparison Based on Next Generation Sequencing Reads: Extended Abstract.

Multiple Alignment-Free Sequence Comparison

The role of WHO in the control of nutritional anaemia.

An efficient parallel algorithm for multiple sequence similarities calculation using a low complexity method.

New Developments Of Alignment-Free Sequence Comparison: Measures, Statistics And Next-Generation Sequencing

Efficient Seeding for Error-Prone Sequences with SubseqHash2

Sequence Alignment As Hypothesis Testing

Similarity Evaluation of DNA Sequences Based on Frequent Patterns and Entropy

A New Statistic for Efficient Detection of Repetitive Sequences

A statistical physics perspective on alignment-independent protein sequence comparison

Similarity analysis of DNA sequences through local distribution of nucleotides in strategic neighborhood

[Alignment-Free Biomolecular Sequence Comparison Method].

An efficient Z-score algorithm for assessing sequence alignments

Comparison of Metatranscriptomic Samples Based on K-Tuple Frequencies

An Algorithm for Alignment-free Sequence Comparison using Logical Match

DNA sequence comparison by a novel probabilistic method

A Measure of Dna Sequence Similarity by Fourier Transform with Applications on Hierarchical Clustering

A simple method for finding related sequences by adding probabilities of alternative alignments

Parallel linear space algorithm for large-scale sequence alignment