Estimation of Alternative Splicing Isoform Frequencies from RNA-Seq Data

Overview

Journal Algorithms Mol Biol

Publisher Biomed Central

Date 2011 Apr 21

PMID 21504602

Citations 82

Authors

Marius Nicolae

Serghei Mangul

Ion I Mandoiu

Alex Zelikovsky

Affiliations

Soon will be listed here.

Abstract

Background: Massively parallel whole transcriptome sequencing, commonly referred as RNA-Seq, is quickly becoming the technology of choice for gene expression profiling. However, due to the short read length delivered by current sequencing technologies, estimation of expression levels for alternative splicing gene isoforms remains challenging.

Results: In this paper we present a novel expectation-maximization algorithm for inference of isoform- and gene-specific expression levels from RNA-Seq data. Our algorithm, referred to as IsoEM, is based on disambiguating information provided by the distribution of insert sizes generated during sequencing library preparation, and takes advantage of base quality scores, strand and read pairing information when available. The open source Java implementation of IsoEM is freely available at http://dna.engr.uconn.edu/software/IsoEM/.

Conclusions: Empirical experiments on both synthetic and real RNA-Seq datasets show that IsoEM has scalable running time and outperforms existing methods of isoform and gene expression level estimation. Simulation experiments confirm previous findings that, for a fixed sequencing cost, using reads longer than 25-36 bases does not necessarily lead to better accuracy for estimating expression levels of annotated isoforms and genes.

Citing Articles

Oarfish: Enhanced probabilistic modeling leads to improved accuracy in long read transcriptome quantification.

Jousheghani Z, Patro R bioRxiv. 2024; .

PMID: 38464200 PMC: 10925290. DOI: 10.1101/2024.02.28.582591.

RNA-seq data science: From raw data to effective interpretation.

Deshpande D, Chhugani K, Chang Y, Karlsberg A, Loeffler C, Zhang J Front Genet. 2023; 14:997383.

PMID: 36999049 PMC: 10043755. DOI: 10.3389/fgene.2023.997383.

Comparison of Metagenomics and Metatranscriptomics Tools: A Guide to Making the Right Choice.

Terron-Camero L, Gordillo-Gonzalez F, Salas-Espejo E, Andres-Leon E Genes (Basel). 2022; 13(12).

PMID: 36553546 PMC: 9777648. DOI: 10.3390/genes13122280.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Baaijens J, Zulli A, Ott I, Nika I, van der Lugt M, Petrone M Genome Biol. 2022; 23(1):236.

PMID: 36348471 PMC: 9643916. DOI: 10.1186/s13059-022-02805-9.

T Cell Epitope Prediction and Its Application to Immunotherapy.

Schaap-Johansen A, Vujovic M, Borch A, Hadrup S, Marcatili P Front Immunol. 2021; 12:712488.

PMID: 34603286 PMC: 8479193. DOI: 10.3389/fimmu.2021.712488.

References

Temple G, Gerhard D, Rasooly R, Feingold E, Good P, Robinson C . The completion of the Mammalian Gene Collection (MGC). Genome Res. 2009; 19(12):2324-33. PMC: 2792178. DOI: 10.1101/gr.095976.109. View

Bloom J, Khan Z, Kruglyak L, Singh M, Caudy A . Measuring differential gene expression by short read sequencing: quantitative comparison to 2-channel gene expression microarrays. BMC Genomics. 2009; 10:221. PMC: 2686739. DOI: 10.1186/1471-2164-10-221. View

She Y, Hubbell E, Wang H . Resolving deconvolution ambiguity in gene alternative splicing. BMC Bioinformatics. 2009; 10:237. PMC: 2739860. DOI: 10.1186/1471-2105-10-237. View

Carninci P, Kasukawa T, Katayama S, Gough J, Frith M, Maeda N . The transcriptional landscape of the mammalian genome. Science. 2005; 309(5740):1559-63. DOI: 10.1126/science.1112014. View

Howard B, Heber S . Towards reliable isoform quantification using RNA-SEQ data. BMC Bioinformatics. 2010; 11 Suppl 3:S6. PMC: 2863065. DOI: 10.1186/1471-2105-11-S3-S6. View

Jackson B, Schnable P, Aluru S . Parallel short sequence assembly of transcriptomes. BMC Bioinformatics. 2009; 10 Suppl 1:S14. PMC: 2648799. DOI: 10.1186/1471-2105-10-S1-S14. View

Roberts A, Trapnell C, Donaghey J, Rinn J, Pachter L . Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011; 12(3):R22. PMC: 3129672. DOI: 10.1186/gb-2011-12-3-r22. View

Richard H, Schulz M, Sultan M, Nurnberger A, Schrinner S, Balzereit D . Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments. Nucleic Acids Res. 2010; 38(10):e112. PMC: 2879520. DOI: 10.1093/nar/gkq041. View

Hiller D, Jiang H, Xu W, Wong W . Identifiability of isoform deconvolution from junction arrays and RNA-Seq. Bioinformatics. 2009; 25(23):3056-9. PMC: 3167695. DOI: 10.1093/bioinformatics/btp544. View

10.

Sultan M, Schulz M, Richard H, Magen A, Klingenhoff A, Scherf M . A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008; 321(5891):956-60. DOI: 10.1126/science.1160342. View

11.

Langmead B, Trapnell C, Pop M, Salzberg S . Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009; 10(3):R25. PMC: 2690996. DOI: 10.1186/gb-2009-10-3-r25. View

12.

Oshlack A, Wakefield M . Transcript length bias in RNA-seq data confounds systems biology. Biol Direct. 2009; 4:14. PMC: 2678084. DOI: 10.1186/1745-6150-4-14. View

13.

Trapnell C, Williams B, Pertea G, Mortazavi A, Kwan G, van Baren M . Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010; 28(5):511-5. PMC: 3146043. DOI: 10.1038/nbt.1621. View

14.

Jiang H, Wong W . Statistical inferences for isoform expression in RNA-Seq. Bioinformatics. 2009; 25(8):1026-32. PMC: 2666817. DOI: 10.1093/bioinformatics/btp113. View

15.

Li B, Ruotti V, Stewart R, Thomson J, Dewey C . RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2009; 26(4):493-500. PMC: 2820677. DOI: 10.1093/bioinformatics/btp692. View

16.

Pasaniuc B, Zaitlen N, Halperin E . Accurate estimation of expression levels of homologous genes in RNA-seq experiments. J Comput Biol. 2011; 18(3):459-68. DOI: 10.1089/cmb.2010.0259. View

17.

Wang E, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C . Alternative isoform regulation in human tissue transcriptomes. Nature. 2008; 456(7221):470-6. PMC: 2593745. DOI: 10.1038/nature07509. View

18.

Anton M, Gorostiaga D, Guruceaga E, Segura V, Carmona-Saez P, Pascual-Montano A . SPACE: an algorithm to predict and quantify alternatively spliced isoforms using microarrays. Genome Biol. 2008; 9(2):R46. PMC: 2374713. DOI: 10.1186/gb-2008-9-2-r46. View

19.

Birol I, Jackman S, Nielsen C, Qian J, Varhol R, Stazyk G . De novo transcriptome assembly with ABySS. Bioinformatics. 2009; 25(21):2872-7. DOI: 10.1093/bioinformatics/btp367. View

20.

Griffith M, Griffith O, Mwenifumbo J, Goya R, Morrissy A, Morin R . Alternative expression analysis by RNA sequencing. Nat Methods. 2010; 7(10):843-7. DOI: 10.1038/nmeth.1503. View