Identifying differentially expressed transcripts from RNA-seq data with biological variation
Abstract
Motivation: High-throughput sequencing enables expression analysis at the level of individual transcripts. The analysis of transcriptome expression levels and differential expression estimation requires a probabilistic approach to properly account for ambiguity caused by shared exons and finite read sampling as well as the intrinsic biological variance of transcript expression. Results: We present BitSeq (Bayesian Inference of Transcripts from Sequencing data), a Bayesian approach for estimation of transcript expression level from RNA-seq experiments. Inferred relative expression is represented by Markov chain Monte Carlo (MCMC) samples from the posterior probability distribution of a generative model of the read data. We propose a novel method for differential expression analysis across replicates which propagates uncertainty from the sample-level model while modelling biological variance using an expression-level-dependent prior. We demonstrate the advantages of our method using simulated data as well as an RNA-seq dataset with technical and biological replication for both studied conditions. Availability: The implementation of the transcriptome expression estimation and differential expression analysis, BitSeq, has been written in C++.
Cite
@article{arxiv.1109.0863,
title = {Identifying differentially expressed transcripts from RNA-seq data with biological variation},
author = {Peter Glaus and Antti Honkela and Magnus Rattray},
journal= {arXiv preprint arXiv:1109.0863},
year = {2012}
}
Comments
12 pages, 6 figures in main text; 11 pages, 5 figures in supplementary information (included in the same file)