Alignment-Free Sequence Comparison Based on Next Generation Sequencing Reads: Extended Abstract.
Kai Song,Jie Ren,Zhiyuan Zhai,Xuemei Liu,Minghua Deng,Fengzhu Sun
DOI: https://doi.org/10.1007/978-3-642-29627-7_29
2012-01-01
Abstract:Next generation sequencing (NGS) technologies have generated enormous amount of shotgun read data and assembly of the reads can be challenging, especially for organisms without template sequences. We study the power of genome comparison based on shotgun read data without assembly using three alignment-free sequence comparison statistics, \(D_2, D_2^*\), and \(D_2^S\), both theoretically and by simulations. Theoretical formulas for the power of detecting the relationship between two sequences related through a common motif model are derived. It is shown that both \(D_2^*\) and \(D_2^S\) outperform D 2 for detecting the relationship between two sequences based on NGS data. We then study the effects of length of the tuple, read length, coverage, and sequencing error on the power of \(D_2^*\) and \(D_2^S\). Finally, variations of these statistics, \(d_2, d_2^*\) and \(d_2^S\), respectively, are used to first cluster 5 mammalian species with known phylogenetic relationships and then cluster 13 tree species whose complete genome sequences are not available using NGS shotgun reads. The clustering results using \(d_2^S\) are consistent with biological knowledge for the 5 mammalian and 13 tree species, respectively. Thus, the statistic \(d_2^S\) provides a powerful alignment-free comparison tool to study the relationships among different organisms based on NGS read data without assembly.