Tweet
Syndicated from PubMed RSS Feeds
Using non-uniform read distribution models to improve isoform expression inference in RNA-Seq.
Bioinformatics. 2010 Dec 17;
Authors: Wu Z, Wang X, Zhang X
MOTIVATION: RNA-Seq technology based on next-generation sequencing provides the unprecedented ability of studying transcriptomes at high resolution and accuracy, and the potential of measuring expression of multiple isoforms from the same gene at high precision. Solved by maximum likelihood estimation, isoform expression can be inferred in RNA-Seq using statistical models based on the assumption that sequenced reads are distributed uniformly along transcripts. Modification of the model is needed when considering situations where RNA-Seq data do not follow uniform distribution. RESULTS: We proposed two curves, the global bias curve (GBC) and the local bias curves (LBCs), to describe the non-uniformity of read distributions for all genes in a transcriptome and for each gene, respectively. Incorporating the bias curves into the uniform read distribution model, we introduced non-uniform read distribution models to infer isoform expression levels. On a series of systematic simulation studies, the proposed models outperform the original model in recovering major isoforms and the expression ratio of alternative isoforms. We also applied the new model to real RNA-Seq data sets and found that its inferences on expression ratios of alternative isoforms are more reasonable. The experiments indicate that incorporating non-uniform read distribution information can improve the accuracy in modeling and inferring isoform expression in RNA-Seq. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
PMID: 21169371 [PubMed - as supplied by publisher]
More...
Using non-uniform read distribution models to improve isoform expression inference in RNA-Seq.
Bioinformatics. 2010 Dec 17;
Authors: Wu Z, Wang X, Zhang X
MOTIVATION: RNA-Seq technology based on next-generation sequencing provides the unprecedented ability of studying transcriptomes at high resolution and accuracy, and the potential of measuring expression of multiple isoforms from the same gene at high precision. Solved by maximum likelihood estimation, isoform expression can be inferred in RNA-Seq using statistical models based on the assumption that sequenced reads are distributed uniformly along transcripts. Modification of the model is needed when considering situations where RNA-Seq data do not follow uniform distribution. RESULTS: We proposed two curves, the global bias curve (GBC) and the local bias curves (LBCs), to describe the non-uniformity of read distributions for all genes in a transcriptome and for each gene, respectively. Incorporating the bias curves into the uniform read distribution model, we introduced non-uniform read distribution models to infer isoform expression levels. On a series of systematic simulation studies, the proposed models outperform the original model in recovering major isoforms and the expression ratio of alternative isoforms. We also applied the new model to real RNA-Seq data sets and found that its inferences on expression ratios of alternative isoforms are more reasonable. The experiments indicate that incorporating non-uniform read distribution information can improve the accuracy in modeling and inferring isoform expression in RNA-Seq. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
PMID: 21169371 [PubMed - as supplied by publisher]
More...