title: The impacts of read length and transcriptome complexity for de novo assembly: a simulation study.
abstract: Transcriptome assembly using RNA-seq data - particularly in non-model organisms has been dramatically improved, but only recently have the pre-assembly procedures, such as sequencing depth and error correction, been studied. Increasing read length is viewed as a crucial condition to further improve transcriptome assembly, but it is unknown whether the read length really matters. In addition, though many assembly tools are available now, it is unclear whether the existing assemblers perform well enough for all data with different transcriptome complexities. In this paper, we studied these two open problems using two high-performing assemblers, Velvet/Oases and Trinity, on several simulated datasets of human, mouse and S.cerevisiae. The results suggest that (1) the read length of paired reads does not matter once it exceeds a certain threshold, and interestingly, the threshold is distinct in different organisms; (2) the quality of de novo assembly decreases sharply with the increase of transcriptome complexity, all existing de novo assemblers tend to corrupt whenever the genes contain a large number of alternative splicing events.
authors.affiliation: School of Mathematics, Shandong University, Jinan, Shandong, China.
authors.forename: Zheng
authors.initials: Z
authors.lastname: Chang
authors.affiliation: School of Mathematics, Shandong University, Jinan, Shandong, China.
authors.forename: Zhenjia
authors.initials: Z
authors.lastname: Wang
authors.affiliation: School of Mathematics, Shandong University, Jinan, Shandong, China.
authors.forename: Guojun
authors.initials: G
authors.lastname: Li
mesh_terms: D017398:Alternative Splicing / Q000235:genetics; D000818:Animals; D020869:Gene Expression Profiling / Q000379:methods*; D059014:High-Throughput Nucleotide Sequencing / Q000379:methods*; D006801:Humans; D051379:Mice; D008962:Models, Theoretical; D012441:Saccharomyces cerevisiae / Q000235:genetics; D017423:Sequence Analysis, RNA / Q000379:methods*
references.citation: Nat Genet. 2008 Dec;40(12):1413-5
references.pmid: 18978789
references.citation: Bioinformatics. 2011 Sep 15;27(18):2518-28
references.pmid: 21775302
references.citation: Nat Biotechnol. 2010 May;28(5):511-5
references.pmid: 20436464
references.citation: Nat Biotechnol. 2011 May 15;29(7):644-52
references.pmid: 21572440
references.citation: Nat Rev Genet. 2010 Jan;11(1):31-46
references.pmid: 19997069
references.citation: Methods. 2009 Jul;48(3):249-57
references.pmid: 19336255
references.citation: PeerJ. 2013 Jul 23;1:e113
references.pmid: 23904992
references.citation: Cell Mol Life Sci. 2010 Feb;67(4):569-79
references.pmid: 19859660
references.citation: Nature. 2008 Nov 27;456(7221):470-6
references.pmid: 18978772
references.citation: Genome Res. 2008 May;18(5):821-9
references.pmid: 18349386
references.citation: Bioinformatics. 2009 Nov 1;25(21):2872-7
references.pmid: 19528083
references.citation: Genome Res. 2008 Sep;18(9):1509-17
references.pmid: 18550803
references.citation: Nat Rev Genet. 2009 Jan;10(1):57-63
references.pmid: 19015660
references.citation: BMC Genomics. 2013 Mar 12;14:167
references.pmid: 23496952
references.citation: Nat Biotechnol. 2010 May;28(5):503-10
references.pmid: 20436462
references.citation: Genome Res. 2002 Apr;12(4):656-64
references.pmid: 11932250
references.citation: Nat Rev Genet. 2011 Sep 07;12(10):671-82
references.pmid: 21897427
references.citation: Genome Res. 2009 Feb;19(2):336-46
references.pmid: 19056694
references.citation: BMC Bioinformatics. 2011 Dec 14;12 Suppl 14:S2
references.pmid: 22373417
references.citation: Bioinformatics. 2009 Dec 1;25(23):3056-9
references.pmid: 19762346
references.citation: Bioinformatics. 2012 Apr 15;28(8):1086-92
references.pmid: 22368243
references.citation: Nat Rev Genet. 2011 Feb;12(2):87-98
references.pmid: 21191423