Normalization of RNA-sequencing data from samples with varying mRNA levels.

TitleNormalization of RNA-sequencing data from samples with varying mRNA levels.
Publication TypeJournal Article
Year of Publication2014
AuthorsAanes, H, Winata, CL, Moen, LF, Østrup, O, Mathavan, S, Collas, P, Rognes, T, Aleström, P
JournalPLoS One
Volume9
Issue2
Paginatione89158
Date Published2014
ISSN1932-6203
KeywordsAnimals, Base Sequence, Gene Expression, Gene Expression Profiling, Polymerase Chain Reaction, RNA, Messenger, Sequence Analysis, RNA, Zebrafish
Abstract

Methods for normalization of RNA-sequencing gene expression data commonly assume equal total expression between compared samples. In contrast, scenarios of global gene expression shifts are many and increasing. Here we compare the performance of three normalization methods when polyA(+) RNA content fluctuates significantly during zebrafish early developmental stages. As a benchmark we have used reverse transcription-quantitative PCR. The results show that reads per kilobase per million (RPKM) and trimmed mean of M-values (TMM) normalization systematically leads to biased gene expression estimates. Biological scaling normalization (BSN), designed to handle differences in total expression, showed improved accuracy compared to the two other methods in estimating transcript level dynamics. The results have implications for past and future studies using RNA-sequencing on samples with different levels of total or polyA(+) RNA.

DOI10.1371/journal.pone.0089158
Alternate JournalPLoS ONE
Citation Key32
PubMed ID24586560
PubMed Central IDPMC3934880