Study design.

<p>Thirty-three BXD lines plus the 2 parental strains and their reciprocal F1 progeny were phenotyped. Mice were submitted to either one of 2 experiments. In Experiment 1 (left), EEG/EMG signals and LMA were recorded under standard 12:12 h light–dark conditions (white and black bars under top-left panel) for 2 baseline days (B1, B2), a 6 h SD (red bar) from ZT0–6 (ZT0 = light onset), followed by 2 recovery days (R1, R2). The deep sleep-wake phenome consists of 341 sleep-wake state-, LMA-, and EEG-related phenotypes quantified in each mouse, among which time spent in NREM sleep (gray area spans mean maximum and minimum NREM sleep time among BXD lines, respectively, for consecutive 90 min intervals). Mice in Experiment 2 (right) were used to collect cortex, liver, and blood samples at ZT6. Half of the mice were challenged with an SD as in Experiment 1, the other half were left undisturbed and served as controls (labeled Ctr). Cortex and liver samples were used to quantify gene expression by RNA-seq, blood samples for a targeted analysis of 124 metabolites by LC/MS, or with FIA/MS. For <i>ph</i>QTLs, <i>m</i>QTLs, and <i>e</i>QTLs, a high-density genotype dataset (Genome; approximately 11,000 SNPs) was created, merging identified RNA-seq variants with a publicly available database (<a href="http://www.genenetwork.org/" target="_blank">www.genenetwork.org</a>). The entirety of the multilevel dataset was integrated in a systems genetics analysis to chart molecular pathways underlying the many facets of sleep and the EEG, using newly developed computational tools to interactively visualize the results and pathways, and to prioritize candidate genes. EEG/EMG, electroencephalography/electromyogram; <i>e</i>QTL, expression quantitative trait locus; FIA/MS, flow injection analysis/mass spectrometry; LC/MS, liquid chromatography/mass spectrometry; LMA, locomotor activity; <i>m</i>QTL, metabolic quantitative trait locus; NREM, non-REM; <i>ph</i>QTL, phenotypic quantitative trait locus; RNA-seq, RNA sequencing; SD, sleep deprivation; ZT, zeitgeber time.</p

EEG delta power in NREM sleep after SD is associated with <i>Kif16b</i> and <i>Wrn</i>.

<p>(A) NREM sleep EEG spectra in the first 3 h after SD (ZT6–9) for the 2 BXD lines that displayed the lowest and highest EEG activity in the fast delta frequency band (2.5–4.25 Hz, δ2; top, see panel E) and for the 2 BXD lines that displayed the smallest and largest increase (or gain) in EEG power in the slow delta band (1.0–2.25 Hz, δ1; bottom, see panel E). Spectra were “1/f-corrected” (and therefore not directly comparable to the values in panel E) for better visualization of activity in higher frequency bands (theta [5–9 Hz, θ], sigma [11–16 Hz, σ], beta [18–30 Hz, β], and slow [32–55 Hz, γ1] and fast gamma [55–80 Hz, γ2]). Subsequent analyses were performed without this correction. (B) QTL mapping and prioritization for δ2 power identified a significant association on chromosome 2 and <i>Kif16b</i> in cortex as top-ranked gene (top). For the δ1 increase after SD, we obtained a suggestive QTL on chromosome 8 and a significant prioritization score for the DNA-helicase <i>Wrn</i>. (C) Hiveplot visualization of network connections for the δ1 and δ2 power after SD (top-left panels) and the SD-induced increase in δ1 and δ2 power over baseline (bottom-left panels). Note the marked differences in the networks and QTLs regulating the expression of these 2 delta bands. Right hiveplots highlight <i>Kif16b</i> in the δ2 power–associated network (top), and <i>Wrn</i> in the network associated with the δ1 increase (bottom). Only <i>Kif16b</i> expression in the cortex was linked to the chromosome 2 <i>cis</i>-<i>e</i>QTL and was not associated with any metabolite. <i>Wrn</i> expression was significantly linked to the chromosome 8 <i>cis-e</i>QTL and to the long phosphatidylcholine, PC-ae-C38:5. (D) <i>Kif16b</i> is highly significantly down-regulated in cortex (left), while it remains unchanged in liver after SD (<i>p</i> = 0.15; not shown). Also, <i>Wrn</i> expression was strongly down-regulated by SD in cortex (right) and only marginally so, albeit significantly, in liver (<i>p</i> = 0.02; not shown). (E) Strain distribution patterns. BXD lines carrying a <i>B6-</i>allele on the chromosome 2–associated region showed higher δ2 power after SD (left) and a significantly higher <i>Kif16b</i> expression (<i>p</i> = 1.3e−15; second to left) than <i>D2-</i>allele carriers. <i>D2-</i>allele carriers of the chromosome 8–associated region showed a larger δ1 increase after SD (second to right) as well as a significantly larger decrease in <i>Wrn</i> expression after SD (right) than <i>B6-</i>allele carriers. For color-coding of genotypes, see <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.2005750#pbio.2005750.g004" target="_blank">Fig 4</a>. CPM, counts per million; Ctr, control; EEG, electroencephalography; <i>e</i>QTL, expression quantitative trait locus; FDR, false discovery rate; <i>Kif16b</i>, <i>Kinesin family member 16B</i>; NREM, non-REM; PC-ae, phosphatidylcholine acyl-alkyl; QTL, quantitative trait locus; SD, sleep deprivation; <i>Wrn</i>, <i>Werner syndrome RecQ like helicase</i>; ZT, zeitgeber time</p

Genetic diversity in the BXD panel greatly impacts behavioral, metabolic, and molecular traits.

<p>The phenome was divided into 3 phenotypic categories: (i) LMA, (ii) EEG features (labeled EEG), and (iii) sleep-wake state characteristics (labeled State), which were subdivided further (see <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.2005750#sec015" target="_blank">Materials and methods</a>). The 5 classes of metabolites and the gene expression represent intermediate molecular phenotypic categories. (A) Heritability for EEG/behavioral and metabolite phenotypes. Dots represent single phenotypes within each category and subcategory indicated along the x-axis. Red dots represent phenotypes recorded in baseline (labeled bsl; B1 and B2), blue in recovery (labeled rec; R1 and R2), purple during SD, and green dots refer to the recovery-to-baseline contrasts. Values represent narrow-sense heritability. (B) Overview of significant and highly suggestive (FDR < 0.1) QTLs obtained for all 341 EEG/behavioral phenotypes (<i>ph</i>QTLs: LMA in red, EEG in blue, and sleep-wake state in green) and 124 blood metabolite levels in baseline and recovery (<i>m</i>QTLs; purple). Note that overlap of neighboring QTLs renders color shading darker. (C) Venn diagram of genes under significant <i>cis-e</i>QTL effect in liver and cortex for the two experimental conditions (SD and controls [labeled Ctr]). EEG, electroencephalography; <i>e</i>QTL, expression quantitative trait locus; FDR, false discovery rate; LMA, locomotor activity; <i>m</i>QTL, metabolic quantitative trait locus; <i>ph</i>QTL, phenotypic quantitative trait locus; QTL, quantitative trait locus; SD, sleep deprivation.</p

NREM sleep gain in the first 6 h of the dark period after SD is associated with <i>Acot11</i>.

<p>(A) Time course of hourly values of time spent in NREM sleep in baseline, SD (red area), and recovery for the 2 BXD lines showing the largest (BXD70; green) and lowest (BXD83; blue) NREM sleep gain during ZT12–18 (left). NREM sleep gain during 4 consecutive 6 h intervals during recovery compared to corresponding baseline intervals shows that in the recovery dark period (gray area), BXD83 mice did not accumulate extra NREM sleep, while BXD70 mice gained 88 min (middle). Strain distribution of ZT12–18 NREM sleep gain (right). B6-allele carriers compensated less for NREM sleep lost during SD than D2-allele carriers. For color-coding, see <a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.2005750#pbio.2005750.g004" target="_blank">Fig 4</a>. (B) Hiveplots for NREM sleep gain in 4 six-hour recovery intervals after the end of SD at ZT6. Compared to the other 3 intervals, NREM sleep gain was strongly associated with a number of metabolites during the second 6 h interval, i.e., ZT12–18. (C) NREM sleep gain during ZT12–18 mapped to a significant QTL on chromosome 4, explaining 45% of the total phenotypic variance (top left). PC-ae-C38:2 mapped suggestively to the same region (top right). Prioritization of liver transcripts for both phenotypes yielded <i>Acot11</i> as top-ranked, significant gene (bottom). (D) Hiveplot for the ZT12–18 NREM sleep gain, highlighting <i>Acot11</i>. <i>Acot11</i> was positively correlated with several phosphatidylcholines and to <i>Ovgp1</i> expression in the cortex. (E) Allelic effect of the chromosome 4–associated region on <i>Acot11</i> expression and PC-ae-C38:2 levels in the BXDs. <i>Acot11</i> expression in liver after SD was under a strong <i>e</i>QTL effect (<i>p</i> = 1.6e−13) with <i>B6-</i>allele carriers showing a higher <i>Acot11</i> expression than <i>D2</i>-allele carriers. <i>B6-</i>allele carriers also showed higher PC-ae-C38:2 levels after SD. (F) Both <i>Acot11</i> and PC-ae-C38:2 levels changed after SD. <i>Acot11</i> in liver and PC-ae-C38:2 in blood were significantly down-regulated. In the cortex, <i>Acot11</i> was, however, significantly up-regulated, and the chromosome 4–associated region did not modulate cortical <i>Acot11</i> expression. (G) Mice carrying 1 or 2 KO alleles for <i>Acot11</i> displayed less extra NREM sleep during recovery. In contrast to the BXD panel, this difference was present in the second (ZT18–24, right) and not during the first (ZT12–18, left panel) 6 h of the recovery dark period. <i>Acot11</i>, <i>acyl-CoA thioesterase 11</i>; CPM, counts per million; Ctr, control; <i>e</i>QTL, expression quantitative trait locus; KO, knockout; NREM, non-REM; PC-ae, phosphatidylcholine acyl-alkyl; QTL, quantitative trait locus; SD, sleep deprivation; ZT, zeitgeber time</p

Profound effects of SD on transcriptome, metabolome, and phenome.

<p>EEG/behavioral phenotypes, metabolites, and transcripts are organized into 3 “columns” (from left to right). Top 3 panels show the SD response (recovery/baseline fold change). Bottom 3 panels depict examples of allelic effects on the SD responses, with color-coding indicating the presence of a C57BL/6J or DBA/2J haplotype under the mapped QTL peaks (B6: gray for BXD and black for parental; D2: light brown for BXD and dark brown for parental). White bars mark the F1s and hatched bars strain in which haplotype could not be unambiguously determined. (A) Phenotypic changes after SD. The top significantly changed phenotype was the increase in NREM sleep EEG delta power (1–4 Hz) after SD (far-left blue data point). The most up-regulated phenotype was time spent in REM sleep during the first 6 h of darkness (ZT12–18) after SD (highest green data point). (B) Metabolite changes after SD. Most amino acids (blue) were down-regulated and most sphingolipids (brown) up-regulated after SD. The acylcarnitines C18:1 and C18:2 (highest red dots) increased the most. Vertical red line: significant threshold (FDR-adjusted <i>p</i>-value = 0.05). (C) DE analysis (SD/Ctr) for cortex (left) and liver (right). Genes were sorted according to their ranked <i>p</i>-value along the x-axis. Significantly affected transcripts in red (FDR-adjusted <i>p</i>-value < 0.05), nonsignificant results in black. Blue dots indicate 78 genes considered core molecular components of the sleep homeostatic response in the cortex [<a href="http://www.plosbiology.org/article/info:doi/10.1371/journal.pbio.2005750#pbio.2005750.ref034" target="_blank">34</a>]. Note that no low fold change threshold was applied. (D-F) Examples of genetically driven EEG/behavioral, metabolic, and transcriptional responses to SD, respectively. See text for details. <i>Arc</i>, <i>activity-regulated cytoskeletal-associated protein</i>; Ctr, control; DE, differential gene expression; EEG, electroencephalography; <i>Egr2</i>, <i>early growth response 2</i>; <i>Fam107a</i>, <i>family with sequence similarity 107</i>, <i>A</i>; FDR, false discovery rate; LMA, locomotor activity; <i>Mlycd</i>, <i>malonyl-CoA decarboxylase</i>; NREM, non-REM; <i>Plin4</i>, <i>Perilipin 4; Pla2g4e</i>, <i>phospholipase A2</i>, <i>group IVE</i>; QTL, quantitative trait locus; SD, sleep deprivation; <i>Ttll8</i>, <i>tubulin tyrosine ligase-like family 8</i>; ZT, zeitgeber time.</p

Changes in the frequency of theta oscillation during REM sleep after SD are associated with <i>Cyp4a32</i>.

<p>(A) Spectral profiles of the REM sleep EEG for 2 strains displaying an opposite shift in the frequency of theta oscillations after SD relative to baseline. This shift was quantified by the decrease and increase in TPF for BXD61 and BXD101, respectively (see panel F). (B) Hiveplot for the SD-induced shift in TPF. (C) One significant QTL for the TPF shift was detected on chromosome 4 and 1 suggestive QTL on chromosome 8. Prioritization yielded <i>Cyp4a32</i> as the top-ranked significant gene, based on the significant <i>cis-e</i>QTL modifying its expression in liver and a predicted damaging variation (V314E). (D) Effects of SD and genotype on liver <i>Cyp4a32</i> expression. Carrying a <i>B6-</i>allele at the <i>Cyp4a32 cis-e</i>QTL–associated marker greatly decreased its expression. (E) Hiveplot for the SD-induced shift in TPF, highlighting <i>Cyp4a32</i>’s links to the amino acid Valine and the chromosome 4 <i>e</i>QTL marker. (F) Strain distribution patterns for TPF differences and liver <i>Cyp4a32</i> expression after SD. <i>B6-</i>allele carriers at the chromosome 4–associated region had lower <i>Cyp4a32</i> liver expression and a decrease in TPF after SD, while <i>D2-</i>carriers increase TPF and have higher <i>Cyp4a32</i> expression. CPM, counts per million; Ctr, control; <i>Cyp4a32</i>, <i>Cytochrome P450</i>, <i>family 4</i>, <i>subfamily a</i>, <i>polypeptide 32</i>; DE, differential expression; EEG, electroencephalography; <i>e</i>QTL, expression quantitative trait locus; lod, logarithm of odds ratio; QTL, quantitative trait locus; SD, sleep deprivation; TPF, theta-peak frequency</p