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.

Bayesian Approach to DNA Segmentation into Regions with Different Average Nucleotide Composition

Gene Prediction by the Noise-Assisted MEMD and Wavelet Transform for Identifying the Protein Coding Regions

New Stopping Criteria for Segmenting DNA Sequences

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

Genotype Sequence Segmentation: Handling Constraints and Noise

Bayesian Random Segmentation Models to Identify Shared Copy Number Aberrations for Array CGH Data

Analysis of DNA sequences through local distribution of nucleotides in strategic neighborhoods

Bayesian Integrative Region Segmentation in Spatially Resolved Transcriptomic Studies

Minimum Segmentation for Pan-genomic Founder Reconstruction in Linear Time

CompostBin: A DNA composition-based algorithm for binning environmental shotgun reads

Coding Region Prediction Based on a Universal DNA Sequence Representation Method

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

An advanced approach for DNA sequencing and similarities analysis on the basis of groupings of nucleotide bases

Partitioning gene expression data by data-driven Markov chain Monte Carlo

Extending the Recursive Jensen-Shannon Segmentation of Biological Sequences

Binning Metagenomic Contigs Using Contig Embedding and Decomposed Tetranucleotide Frequency

A Novel Approach to Clustering Genome Sequences Using Inter-nucleotide Covariance

DNA Hash Pooling and its Applications

Bayesian Functional Data Analysis over Dependent Regions and Its Application for Identification of Differentially Methylated Regions

Performance evaluation of DNA copy number segmentation methods

Partition Based Algorithms for Rearrangement Distances with Flexible Intergenic Regions