劉雪梅 臧翔 黃天來 楊哲，2 李文 葉宇中 胡珊

(1.華南理工大學(xué) 物理與光電學(xué)院， 廣東 廣州 510640； 2.中國工商銀行 廣州東城支行， 廣東 廣州 510100；3.中山大學(xué) 中山醫(yī)學(xué)院計(jì)算機(jī)中心， 廣東 廣州 510275)

1 序列非比對法原理

(1)

式中：CSi(W)為在序列Si中W出現(xiàn)的次數(shù)，W為序列中所有可能的k-tuple組成的集合；wi為k-tuple元素；A為字母表，是一組符號(hào)或字符，字母表中的元素組成序列.

(2)

式中，pW為W在序列中出現(xiàn)的概率.

(3)

2 病毒和宿主細(xì)胞的DNA序列非比對比較

圖1 k=5、Morder=1，2，3時(shí)和的AUC柱狀圖

2.2 最佳閾值的確定

圖2 k=5、Morder=1時(shí)的ROC曲線

表和統(tǒng)計(jì)打分值

2.3 最佳閾值應(yīng)用實(shí)例

3 結(jié)語

：

[1] WOMMACK K E，COLWELL R R.Virioplankton:viruses in aquatic ecosystems [J].Microbiol Mol Biol Rev，2000，64(1):69- 114．

[2] WEINBAUER M G.Ecology of prokaryotic viruses [J].FEMS Microbiol Rev，2004，28(2):127- 181．

[3] SUTTLE C A.Marine viruses-major players in the global ecosystem [J].Nat Rev Microbiol，2007，5(10):801- 812．

[4] CRAM J A，LI C X，NEEDHAM D M，et al.Cross-depth analysis of marine bacterial networks suggests downward propagation of temporal changes [J].ISME Journal，2015，9(12):2573- 2586．

[5] QI J，WANG B，HAO B I.Whole proteome prokaryote phylogeny without sequence alignment:aK-string composition approach [J].Journal of Molecular Evolution，2004，58(1):1- 11．

[6] DOMAZETLOSO M，HAUBOLD B.Alignment- free detection of local similarity among viral and bacterial genomes [J].Bioinformatics，2011，27(11):1466- 1472．

[7] REINERT G，CHEW D，SUN F，et al.Alignment-free sequence comparison(I):statistics and power [J].Journal of Computational Biology：A Journal of Computational Molecular Cell Biology，2009，16(12):1615- 1634．

[8] LIPPERT R A，HUANG H，WATERMAN M S.Distributional regimes for the number ofk-word matches between two random sequences [J].Proceedings of the National Academy of Sciences of the United States of America，2002，99(22):13980- 13989．

[9] WANG Y，LEI X，WANG S，et al.Effect ofk-tuple length on sample-comparison with high-through put sequencing data [J].Biochemical & Biophysical Research Communications，2015，469(4):1021- 1027．

[10] SONG K，REN J，LIU X M，et al.Alignment-free sequence comparison based on next-generation sequencing reads [J].Journal of Computational Biology，2012，20(2):64- 79．

[11] ZHAI Z，KU S Y，LUAN Y，et al.The power of detecting enriched patterns:an HMM approach [J].Journal of Computational Biology：A Journal of Computational Molecular Cell Biology，2010，17(4):581- 592．

[12] WANG Y，LIU L，CHEN L，et al.Comparison of meta-transcriptomic samples based onk-tuple frequencies [J].PLoS ONE，2014，9(1):e84348/1- 19．

[13] 劉雪梅，文德華，於黃忠，等.基于D2shepp統(tǒng)計(jì)法的非序列局部比對 [J].華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012，40(8):106- 110．

LIU Xue-mei，WEN De-hua，YU Huang-zhong，et al.Local alignment-free sequences based on D2shepp statistics [J].Journal of South China University of Technology(Natural Science Edition)，2012，40(8):106- 110．

[14] LIU X M，WAN L，LI J，et al.New powerful statistics for alignment-free sequence comparison under a pattern transfer model [J].Journal of Theoretical Biology，2011，284(1):106- 116．

[15] WAN L，REINERT G，SUN F，et al.Alignment-free sequence comparison(II):theoretical power of comparison statistics [J].Journal of Computational Biology，2010，17(11):1467- 1490．

[16] BAI J，KAI S，JIE R，et al.Comparison of metagenomics samples using sequence signatures [J].BMC Genomics，2012，13(1):138- 140．