Sparse linear modeling of next-generation mRNA sequencing (RNA-Seq) data for isoform discovery and abundance estimation

Proceedings of the National Academy of Sciences of the United States of America

Volumn 108, Issue 50, 2011, Pages 19867-19872

  Li, Jingyi Jessica ^a   Jiang, Ci Ren ^b   Brown, James B ^a   Huang, Haiyan ^a   Bickel, Peter J ^a  

^a UNIVERSITY OF CALIFORNIA   (United States)

^b Statistical and Applied Mathematical Sciences Institute USA   (United States)

Author keywords

Fragment length distribution; GC contents; mRNA isoform discovery; Penalized estimation; Single end vs. paired end sequencing

Indexed keywords

MESSENGER RNA;

ARTICLE; EXON; PRIORITY JOURNAL; PROBABILITY; RNA SEQUENCE; SAMPLING; STATISTICAL MODEL;

ANIMALS; COMPUTER SIMULATION; DATABASES, NUCLEIC ACID; DROSOPHILA MELANOGASTER; EXONS; LINEAR MODELS; PROTEIN ISOFORMS; RNA, MESSENGER; SEQUENCE ANALYSIS, RNA; SOFTWARE;

EID: 84055216924     PISSN: 00278424     EISSN: 10916490     Source Type: Journal    
DOI: 10.1073/pnas.1113972108     Document Type: Article

References (28)

1. Trapnell C, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28:511-515.
2. Dohm JC, Lottaz C, Borodina T, Himmelbauer H (2010) Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res 36:e105.
3. Hansen KD, Brenner SE, Dudoit S (2010) Biases in Illumina transcriptome sequencing caused by random hexamer priming. Nucleic Acids Res 38:e131.
4. Li J, Jiang H,WongWH(2010) Modeling non-uniformity in short-read rates in RNA-Seq data. Genome Biol 11:R50.
5. Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L (2011) Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol 12:R22.
6. modENCODE Consortium (2009) Unlocking the secrets of the genome. Nature 459:927-930.
7. modENCODE Consortium (2010) Identification of functional elements and regulatory circuits by Drosophila modENCODE. Science 330:1787-1797.
8. Gerstein MB, et al. (2010) Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project. Science 330:1775-1787.
9. Lee S, et al. (2011) Accurate quantification of transcriptome from RNA-Seq data by effective length normalization. Nucleic Acids Res 39:e9.
10. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5:621-628. (Pubitemid 351911867)
11. Feng J, et al. (2010) Inference of isoforms from short sequence reads. 14th Annual International Conference on Research in Computational Molecular Biology (RECOMB 2010), Lecture Notes on Computer Science, 6044 (Springer, Berlin/Heidelber), pp 138-157.
12. Jiang H, Wong WH (2009) Statistical inferences for isoform expression in RNA-Seq. Bioinformatics 25:1026-1032.
13. Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN (2010) RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics 26:493-500.
14. Richard H, et al. (2011) Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments. Nucleic Acids Res 38:e112.
15. Li W, Feng J, Jiang T (2011) IsoLasso: A LASSO regression approach to RNA-Seq based transcriptome assembly. 15th Annual International Conference on Research in Computational Molecular Biology (RECOMB 2011), Lecture Notes on Computer Science, 6577 (Springer, Berlin/Heidelber), pp 168-188.
16. Flicek P, et al. (2011) Ensembl 2011. Nucleic Acids Res 39(Suppl 1):D800-D806.
17. Fujita PA (2011) The UCSC Genome Browser database: Update 2011. Nucleic Acids Res 39(Suppl 1):D876-D882.
18. Frohman MA, Dush MK, Martin GR (1988) Rapid production of full-length cDNAs from rare transcripts: Amplification using a single gene-specific oligonucleotide primer. Proc Natl Acad Sci USA 85:8998-9002. (Pubitemid 19001832)
19. Shiraki T, et al. (2003) Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. Proc Natl Acad Sci USA 100:15776-15781. (Pubitemid 38021066)
20. Guttman M, et al. (2010) Ab initio reconstruction of cell typespecific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs. Nat Biotechnol 28:503-510.
21. Adams MD, et al. (2003) Complementary DNA sequencing: Expressed sequence tags and human genome project. Science 252:1651-1656.
22. Tibshirani R (1996) Regression shrinkage and selection via the lasso. J R Stat Soc Series B Stat Methodol 58:267-288.
23. Liu S, Lin L, Jiang P, Wang D, Xing Y (2011) A comparison of RNA-Seq and high-density exon array for detecting differential gene expression between closely related species. Nucleic Acids Res 39:578-588.
24. Xia Z, Wen J, Chang C, Zhou X (2011) NSMAP: A method for spliced isoforms identification and quantification from RNA-Seq. BMC Bioinformatics 12:162.
25. Meinshausen N, Bühlmann P (2010) Stability selection. J R Stat Soc Series B Stat Methodol 72:417-473.
26. Dahinden C, Parmigiani G, Emerick MC, Bhlmann P (2007) Penalized likelihood for sparse contingency tables with an application to full-length cDNA libraries. BMC Bioinformatics 8:1-11.
27. Goeman JJ (2010) Penalized: L1 (Lasso) and L2 (Ridge) penalized estimation in GLMs and in the Cox model. R package version 0.9-31., Available at http://cran.r-project.org/web/packages/penalized/.
28. Mullen KM, van Stokkum IHM (2010) nnls: The Lawson-Hanson algorithm for non-negative least squares (NNLS). R package version 1.3., Available at http://cran.r-project.org/web/packages/nnls/.