» Articles » PMID: 29731954

MSIQ: JOINT MODELING OF MULTIPLE RNA-SEQ SAMPLES FOR ACCURATE ISOFORM QUANTIFICATION

Overview
Journal Ann Appl Stat
Date 2018 May 8
PMID 29731954
Citations 3
Authors
Affiliations
Soon will be listed here.
Abstract

Next-generation RNA sequencing (RNA-seq) technology has been widely used to assess full-length RNA isoform abundance in a high-throughput manner. RNA-seq data offer insight into gene expression levels and transcriptome structures, enabling us to better understand the regulation of gene expression and fundamental biological processes. Accurate isoform quantification from RNA-seq data is challenging due to the information loss in sequencing experiments. A recent accumulation of multiple RNA-seq data sets from the same tissue or cell type provides new opportunities to improve the accuracy of isoform quantification. However, existing statistical or computational methods for multiple RNA-seq samples either pool the samples into one sample or assign equal weights to the samples when estimating isoform abundance. These methods ignore the possible heterogeneity in the quality of different samples and could result in biased and unrobust estimates. In this article, we develop a method, which we call "joint modeling of multiple RNA-seq samples for accurate isoform quantification" (MSIQ), for more accurate and robust isoform quantification by integrating multiple RNA-seq samples under a Bayesian framework. Our method aims to (1) identify a consistent group of samples with homogeneous quality and (2) improve isoform quantification accuracy by jointly modeling multiple RNA-seq samples by allowing for higher weights on the consistent group. We show that MSIQ provides a consistent estimator of isoform abundance, and we demonstrate the accuracy and effectiveness of MSIQ compared with alternative methods through simulation studies on genes. We justify MSIQ's advantages over existing approaches via application studies on real RNA-seq data from human embryonic stem cells, brain tissues, and the HepG2 immortalized cell line. We also perform a comprehensive analysis of how the isoform quantification accuracy would be affected by RNA-seq sample heterogeneity and different experimental protocols.

Citing Articles

scINSIGHT for interpreting single-cell gene expression from biologically heterogeneous data.

Qian K, Fu S, Li H, Li W Genome Biol. 2022; 23(1):82.

PMID: 35313930 PMC: 8935111. DOI: 10.1186/s13059-022-02649-3.


AIDE: annotation-assisted isoform discovery with high precision.

Li W, Li S, Tong X, Deng L, Shi H, Li J Genome Res. 2019; 29(12):2056-2072.

PMID: 31694868 PMC: 6886511. DOI: 10.1101/gr.251108.119.


Modeling and analysis of RNA-seq data: a review from a statistical perspective.

Li W, Li J Quant Biol. 2019; 6(3):195-209.

PMID: 31456901 PMC: 6711375. DOI: 10.1007/s40484-018-0144-7.

References
1.
Wu A, Neff N, Kalisky T, Dalerba P, Treutlein B, Rothenberg M . Quantitative assessment of single-cell RNA-sequencing methods. Nat Methods. 2013; 11(1):41-6. PMC: 4022966. DOI: 10.1038/nmeth.2694. View

2.
Steijger T, Abril J, Engstrom P, Kokocinski F, Hubbard T, Guigo R . Assessment of transcript reconstruction methods for RNA-seq. Nat Methods. 2013; 10(12):1177-84. PMC: 3851240. DOI: 10.1038/nmeth.2714. View

3.
Rosenbloom K, Armstrong J, Barber G, Casper J, Clawson H, Diekhans M . The UCSC Genome Browser database: 2015 update. Nucleic Acids Res. 2014; 43(Database issue):D670-81. PMC: 4383971. DOI: 10.1093/nar/gku1177. View

4.
Pruitt K, Brown G, Hiatt S, Thibaud-Nissen F, Astashyn A, Ermolaeva O . RefSeq: an update on mammalian reference sequences. Nucleic Acids Res. 2013; 42(Database issue):D756-63. PMC: 3965018. DOI: 10.1093/nar/gkt1114. View

5.
Katz Y, Wang E, Airoldi E, Burge C . Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010; 7(12):1009-15. PMC: 3037023. DOI: 10.1038/nmeth.1528. View