A new method to compute K-mer frequencies and its application to annotate large repetitive plant genomes.
A new method to compute K-mer frequencies and its application to annotate large repetitive plant genomes.
复制标题
DOI:
10.1186/1471-2164-9-517
复制
发表时间:
2008-10-31
期刊:
影响因子:
4.4
通讯作者:
Ware D
中科院分区:
文献类型:
--
作者:
Kurtz S;Narechania A;Stein JC;Ware D
The challenges of accurate gene prediction and enumeration are further aggravated in large genomes that contain highly repetitive transposable elements (TEs). Yet TEs play a substantial role in genome evolution and are themselves an important subject of study. Repeat annotation, based on counting occurrences of k-mers, has been previously used to distinguish TEs from low-copy genic regions; but currently available software solutions are impractical due to high memory requirements or specialization for specific user-tasks. Here we introduce the Tallymer software, a flexible and memory-efficient collection of programs for k-mer counting and indexing of large sequence sets. Unlike previous methods, Tallymer is based on enhanced suffix arrays. This gives a much larger flexibility concerning the choice of the k-mer size. Tallymer can process large data sizes of several billion bases. We used it in a variety of applications to study the genomes of maize and other plant species. In particular, Tallymer was used to index a set of whole genome shotgun sequences from maize (B73) (total size 109 bp.). We analyzed k-mer frequencies for a wide range of k. At this low genome coverage (≈ 0.45×) highly repetitive 20-mers constituted 44% of the genome but represented only 1% of all possible k-mers. Similar low-complexity was seen in the repeat fractions of sorghum and rice. When applying our method to other maize data sets, High-C0t derived sequences showed the greatest enrichment for low-copy sequences. Among annotated TEs, the most highly repetitive were of the Ty3/gypsy class of retrotransposons, followed by the Ty1/copia class, and DNA transposons. Among expressed sequence tags (EST), a notable fraction contained high-copy k-mers, suggesting that transposons are still active in maize. Retrotransposons in Mo17 and McC cultivars were readily detected using the B73 20-mer frequency index, indicating their conservation despite extensive rearrangement across cultivars. Among one hundred annotated bacterial artificial chromosomes (BACs), k-mer frequency could be used to detect transposon-encoded genes with 92% sensitivity, compared to 96% using alignment-based repeat masking, while both methods showed 92% specificity. The Tallymer software was effective in a variety of applications to aid genome annotation in maize, despite limitations imposed by the relatively low coverage of sequence available. For more information on the software, see .
登录
查看更多内容
影响因子:
4.3
作者:
Li R;Ye J;Li S;Wang J;Han Y;Ye C;Wang J;Yang H;Yu J;Wong GK;Wang J
通讯作者:
Wang J
影响因子:
7.4
作者:
Haberer, G;Young, S;Messing, J
通讯作者:
Messing, J
影响因子:
64.8
作者:
Matsumoto, T;Wu, JZ;Sugiura, M
通讯作者:
Sugiura, M
影响因子:
64.8
作者:
Jaillon, Olivier;Aury, Jean-Marc;Wincker, Patrick
通讯作者:
Wincker, Patrick
DOI:
10.1073/pnas.0502923102
发表时间:
2005-06-21
影响因子:
11.1
作者:
Lai, JS;Li, YB;Dooner, HK
通讯作者:
Dooner, HK