IntAPT: integrated assembly of phenotype-specific transcripts from multiple RNA-seq profiles.
IntAPT: integrated assembly of phenotype-specific transcripts from multiple RNA-seq profiles.
复制标题
IntAPT:来自多个 RNA-seq 配置文件的表型特异性转录本的集成组装。
DOI:
10.1093/bioinformatics/btaa852
复制
发表时间:
2021
期刊:
影响因子:
--
通讯作者:
Xuan,Jianhua
中科院分区:
文献类型:
--
作者:
Shi,Xu;Neuwald,AndrewF;Wang,Xiao;Wang,Tian-Li;Hilakivi-Clarke,Leena;Clarke,Robert;Xuan,Jianhua
MotivationHigh-throughput RNA sequencing has revolutionized the scope and depth of transcriptome analysis. Accurate reconstruction of a phenotype-specific transcriptome is challenging due to the noise and variability of RNA-seq data. This requires computational identification of transcripts from multiple samples of the same phenotype, given the underlying consensus transcript structure.ResultsWe present a Bayesian method, integrated assembly of phenotype-specific transcripts (IntAPT), that identifies phenotype-specific isoforms from multiple RNA-seq profiles. IntAPT features a novel two-layer Bayesian model to capture the presence of isoforms at the group layer and to quantify the abundance of isoforms at the sample layer. A spike-and-slab prior is used to model the isoform expression and to enforce the sparsity of expressed isoforms. Dependencies between the existence of isoforms and their expression are modeled explicitly to facilitate parameter estimation. Model parameters are estimated iteratively using Gibbs sampling to infer the joint posterior distribution, from which the presence and abundance of isoforms can reliably be determined. Studies using both simulations and real datasets show that IntAPT consistently outperforms existing methods for the IntAPT. Experimental results demonstrate that, despite sequencing errors, IntAPT exhibits a robust performance among multiple samples, resulting in notably improved identification of expressed isoforms of low abundance.Availability and implementationThe IntAPT package is available at http://github.com/henryxushi/IntAPT.Supplementary informationSupplementary data are available atBioinformaticsonline.