GC-Content Normalization for RNA-Seq Data.

[admin]

GC-Content Normalization for RNA-Seq Data.

BMC Bioinformatics. 2011 Dec 17;12(1):480

Authors: Risso D, Schwartz K, Sherlock G, Dudoit S

Abstract
ABSTRACT: BACKGROUND: Transcriptome sequencing (RNA-Seq) has become the assay of choice for high-throughput studies of gene expression. However, as is the case with microarrays, major technology-related artifacts and biases affect the resulting expression measures. Normalization is therefore essential to ensure accurate inference of expression levels and subsequent analyses thereof. RESULTS: We focus on biases related to GC-content and demonstrate the existence of strong sample-specific GC-content effects on RNA-Seq read counts, which can substantially bias differential expression analysis. We propose three simple within-lane gene-level GC-content normalization approaches and assess their performance on two different RNA-Seq datasets, involving different species and experimental designs. Our methods are compared to state-of-the-art normalization procedures in terms of bias and mean squared error for expression fold-change estimation and in terms of Type I error and p-value distributions for tests of differential expression. The exploratory data analysis and normalization methods proposed in this article are implemented in the open-source Bioconductor R package EDASeq. CONCLUSIONS: Our within-lane normalization procedures, followed by between-lane normalization, reduce GC-content bias and lead to more accurate estimates of expression fold-changes and tests of differential expression. Such results are crucial for the biological interpretation of RNA-Seq experiments, where downstream analyses can be sensitive to the supplied lists of genes.


PMID: 22177264 [PubMed - as supplied by publisher]



PubMed comprises more than 19 million citations for biomedical articles from MEDLINE and life science journals. This RSS feed searches for mentions of Bioconductor - the open source and open development software project for the analysis and comprehension of genomic data.
More from PubMed...