Abstract:MOTIVATIONNext-generation sequencing (NGS) technologies generate large amounts of short read data for many different organisms. The fact that NGS reads are generally short makes it challenging to assemble the reads and reconstruct the original genome sequence. For clustering genomes using such NGS data, word-count based alignment-free sequence comparison is a promising approach, but for this approach, the underlying expected word counts are essential.A plausible model for this underlying distribution of word counts is given through modeling the DNA sequence as a Markov chain (MC). For single long sequences, efficient statistics are available to estimate the order of MCs and the transition probability matrix for the sequences. As NGS data do not provide a single long sequence, inference methods on Markovian properties of sequences based on single long sequences cannot be directly used for NGS short read data.RESULTSHere we derive a normal approximation for such word counts. We also show that the traditional Chi-square statistic has an approximate gamma distribution ,: using the Lander-Waterman model for physical mapping. We propose several methods to estimate the order of the MC based on NGS reads and evaluate those using simulations. We illustrate the applications of our results by clustering genomic sequences of several vertebrate and tree species based on NGS reads using alignment-free sequence dissimilarity measures. We find that the estimated order of the MC has a considerable effect on the clustering results ,: and that the clustering results that use a N: MC of the estimated order give a plausible clustering of the species.AVAILABILITY AND IMPLEMENTATIONOur implementation of the statistics developed here is available as R package 'NGS.MC' at http://www-rcf.usc.edu/∼fsun/Programs/NGS-MC/NGS-MC.htmlCONTACTfsun@usc.eduSUPPLEMENTARY INFORMATIONSupplementary data are available at Bioinformatics online.

Reconsidering the significance of genomic word frequency

Scaling behaviors of CG clusters in coding and noncoding DNA sequences

Statistical properties and fractals of nucleotide clusters in DNA sequences

Long-range correlations in DNA sequences using 2D DNA walk based on pairs of sequential nucleotides

The mathematics of the genetic code reveal that frequency degeneracy leads to exponential scaling in the DNA codon distribution of Homo sapiens

Long-Tail Feature of DNA Words Over- and Under-Representation in Coding Sequences

Microbial genome as a fluctuating system: Distribution and correlation of coding sequence lengths

Range-Limited Heaps' Law for Functional DNA Words in the Human Genome

Comparative statistical analysis of bacteria genomes in "word" context

Comparing reverse complementary genomic words based on their distance distributions and frequencies

Universal power law behaviors in genomic sequences and evolutionary models

Statistics of DNA sequences: a low frequency analysis

Universal spectrum for DNA base CG frequency distribution in Takifugu rubripes (Puffer fish) genome

Statistical distributions and entropy considerations in gene codes

Inference of Markovian Properties of Molecular Sequences from NGS Data and Applications to Comparative Genomics

A Fourier Characteristic of Coding Sequences: Origins and a Non-Fourier Approximation.

Universal 1/f noise, cross-overs of scaling exponents, and chromosome specific patterns of GC content in DNA sequences of the human genome

The diversity of a distributed genome in bacterial populations

The Birth-Death-Mutation process: a new paradigm for fat tailed distributions

Rare-Event Sampling Analysis Uncovers the Fitness Landscape of the Genetic Code

Universality and Shannon entropy of codon usage