Removal of unwanted variation reveals novel patterns of gene expression linked to sleep homeostasis in murine cortex

Abstract

BackgroundWhy we sleep is still one of the most perplexing mysteries in biology. Strong evidence indicates that sleep is necessary for normal brain function and that sleep need is a tightly regulated process. Surprisingly, molecular mechanisms that determine sleep need are incompletely described. Moreover, very little is known about transcriptional changes that specifically accompany the accumulation and discharge of sleep need. Several studies have characterized differential gene expression changes following sleep deprivation. Much less is known, however, about changes in gene expression during the compensatory response to sleep deprivation (i.e. recovery sleep).ResultsIn this study we present a comprehensive analysis of the effects of sleep deprivation and subsequent recovery sleep on gene expression in the mouse cortex. We used a non-traditional analytical method for normalization of genome-wide gene expression data, Removal of Unwanted Variation (RUV). RUV improves detection of differential gene expression following sleep deprivation. We also show that RUV normalization is crucial to the discovery of differentially expressed genes associated with recovery sleep. Our analysis indicates that the majority of transcripts upregulated by sleep deprivation require 6 h of recovery sleep to return to baseline levels, while the majority of downregulated transcripts return to baseline levels within 1–3 h. We also find that transcripts that change rapidly during recovery (i.e. within 3 h) do so on average with a time constant that is similar to the time constant for the discharge of sleep need.ConclusionsWe demonstrate that proper data normalization is essential to identify changes in gene expression that are specifically linked to sleep deprivation and recovery sleep. Our results provide the first evidence that recovery sleep is comprised of two waves of transcriptional regulation that occur at different times and affect functionally distinct classes of genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3065-8) contains supplementary material, which is available to authorized users.

A complete statistical model for calibration of RNA-seq counts using external spike-ins and maximum likelihood theory

“…Our work follows up on and extends the work of [15, 16, 36, 43, 44]. Our inference method is linear and global for each library, like that of [19], [36] and [45].…”

Persistent activation of interlinked type 2 airway epithelial gene networks in sputum-derived cells from aeroallergen-sensitized symptomatic asthmatics

“…We observed asymmetry in the number of up and downregulated genes in our analyses. Previous studies suggest this may occur when RUV adjustment is employed to remove the effect of unwanted variation in the data, but this adjustment was necessary to obtain a uniform p-value distribution 63 . Prior knowledge was employed to reconstruct the gene networks, which relies on data derived from experimental settings that may be far removed from the current study, and this is oversimplified given that genes can function in a context specific manner.…”

Primed to Sleep: The Dynamics of Synaptic Plasticity Across Brain States

“…PRPs transcript levels generally are maximal during waking and then decline during sleep (Figure 2 and Table 1). Several IEGs show a peak and decline in their cortical expression following SD and recovery sleep that parallels the time constant for the discharge of sleep homeostasis (as measured by SWA) (Gerstner et al, 2016). In addition, some of the more dramatic sleep phenotypes in mutant animals result from changes in membrane excitability [e.g., potassium channels in drosophila (Cirelli et al, 2005; Bushey et al, 2007; Koh et al, 2008; Wu et al, 2014) and mice (Espinosa et al, 2004; Douglas et al, 2007)].…”