title: On the importance of small changes in RNA expression.
abstract: The analysis of the differential expression of genes has been the key goal of many molecular biology methods for decades and will remain with us for decades to come. It constitutes a fundamental resource at our disposal for determining the relationship between products of transcription, biology and disease. The completed genome sequencing of many common species allowed microarrays and RNA sequencing (RNAseq) to become major tools in Systems Biology. However, we estimate that at least half of all experiments ignore transcripts that change less than some subjectively chosen threshold, typically around 2-3 fold. Here we show that a majority of the informative RNAs and differentially expressed transcripts can exhibit fold changes less than 2. We use highly quantitative single-molecule sequencing of total cellular RNA derived from a time course of inflammatory response, a process critical to a large number of diseases. Furthermore, we show that enrichment of biologically-relevant functions occurs even at very low fold changes in RNA levels. In addition, we show that most of the common statistical methods can reliably detect transcripts with low fold change when as few as 3 biological replicates are sequenced using single-molecule based RNAseq. In conclusion, given the prevalence of expression profiling in current research, the loss of data in half of all expression studies results in a significant, yet needless drain on the discovery process.
authors.affiliation: Immunovirology - Biogenisis Group, University of Antioquia, A.A. 1226, Medellin, Colombia; St. Laurent Institute, One Kendall Square, Cambridge, MA, USA. Electronic address: georgest98@yahoo.com.
authors.forename: Georges
authors.initials: G
authors.lastname: St Laurent
authors.forename: Dmitry
authors.initials: D
authors.lastname: Shtokalo
authors.forename: Michael R
authors.initials: MR
authors.lastname: Tackett
authors.forename: Zhaoqing
authors.initials: Z
authors.lastname: Yang
authors.forename: Yuri
authors.initials: Y
authors.lastname: Vyatkin
authors.forename: Patrice M
authors.initials: PM
authors.lastname: Milos
authors.forename: Bernd
authors.initials: B
authors.lastname: Seilheimer
authors.forename: Timothy A
authors.initials: TA
authors.lastname: McCaffrey
authors.forename: Philipp
authors.initials: P
authors.lastname: Kapranov
mesh_terms: D001483:Base Sequence; D020869:Gene Expression Profiling / Q000379:methods*; D006801:Humans; D007249:Inflammation / Q000235:genetics / Q000378:metabolism; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D012313:RNA / Q000235:genetics*; D017423:Sequence Analysis, RNA / Q000379:methods; D049490:Systems Biology; D059467:Transcriptome / Q000235:genetics
chemical_list: D012313:RNA
keywords: DE; Differential gene expression; FC; FDR; FISH; GAEA; GO; Gene Annotation Enrichment Analysis; Inflammation; LPS; MAQC; MicroArray quality control; RNAseq; RPK10M; SMS; Single molecule sequencing; Transcriptome; differential expression; false discovery rate; fluorescent in situ hybridization; fold change; gene ontology; lipopolysaccharide; reads per thousand bp of exon per 10million total informative reads per sample; single-molecule sequencing