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

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, D2, D2*, and D2S, 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 D2* and D2S outperform D2 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 D2* and D2S. Finally, variations of these statistics, d2, d2* and d2S, 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 d2S are consistent with biological knowledge for the 5 mammalian and 13 tree species, respectively. Thus, the statistic d2S provides a powerful alignment-free comparison tool to study the relationships among different organisms based on NGS read data without assembly.