Sleep is a ubiquitous state in animals (Anafi et al., 2019) that is controlled by an ensemble of nuclei and their brain-wide interactions (Pace-Schott and Hobson, 2002; Sakai, 2020; Saper et al., 2010). In mammals (Cirelli, 2009), sleep is broadly comprised of two stages: non-rapid eye movement (NREM or slow-wave sleep) and rapid eye movement (REM or paradoxical sleep) (Weber and Dan, 2016). Each of these states is associated with distinct patterns of electrical activity in the brain. During NREM sleep, there are broad-band increases in the local field potential (LFP) power in the cortex that are modulated at <1 Hz (Amzica and Steriade, 1998; Datta and Maclean, 2007; Saper and Fuller, 2017). During REM sleep, gamma band power (nominally 30–100 Hz) is elevated in the cortex (Cantero et al., 2004; Le Van Quyen et al., 2010). In the hippocampus, REM sleep is also associated with a marked increase in power in the theta band (nominally 4–10 Hz) in rodents (Montgomery et al., 2008; Sullivan et al., 2014).

While the dynamics of neural activity in the cortex and other brain structures during sleep are well characterized, the cerebrovascular manifestations of sleep are less clear. Pioneering studies using positron-emission tomography (PET) or ¹³³Xenon in humans have suggested that cerebral blood flow (CBF) and metabolism is reduced during NREM sleep as compared with those during awake levels, and increased above said levels during REM sleep (Braun et al., 1997; Townsend et al., 1973), though the temporal and spatial resolutions of these techniques are poor. The degree to which CBF changes during the different sleep states appears dependent upon brain region (Madsen et al., 1991; Maquet and Phillips, 1998), complicating the interpretation of functional connectivity studies looking at correlations between brain regions where subjects may be transitioning among arousal states (Gu et al., 2019). Several fMRI studies have shown significant alterations in hemodynamic signals and functional connectivity during NREM sleep (Boly et al., 2012; Dang-Vu et al., 2008; Fukunaga et al., 2006; Horovitz et al., 2008; Larson-Prior et al., 2009; Mitra et al., 2015), suggesting that the blood oxygen level dependent (BOLD) signal changes during sleep. Because BOLD signals are generated by a complicated interplay of cerebral metabolism and changes in blood flow and volume (Kim and Ogawa, 2012), the vascular basis of these changes and their relation to neural activity are not well understood. Although functional ultrasound measures indicate cerebral blood volume rises during sleep (likely due to arterial dilation Mateo et al., 2017) and is correlated with hippocampal theta and gamma band power (Bergel et al., 2018), the relationship of cortical vascular changes to local cortical neural activity is unknown.

Understanding the vascular basis of hemodynamic signals during sleep is relevant to many aspects of brain health and function. First, BOLD signal changes during sleep are associated with movement of cerebrospinal fluid (CSF) (Fultz et al., 2019), and the movement of CSF is thought to play an important role in maintaining brain health (Simon and Iliff, 2016; Tarasoff-Conway et al., 2015; Xie et al., 2013). Elucidating the vascular changes associated with these fluid movements would help resolve the actual drivers of fluid movement. Secondly, there is accumulating evidence that arousal state transitions drive large hemodynamic changes, both in animals performing tasks (Cardoso et al., 2019), and in humans and animals undergoing ‘resting-state’ studies (Chang et al., 2016; Liu, 2016; Liu et al., 2018; Tagliazucchi and Laufs, 2014). If the hemodynamic signals during the periods of altered arousal are large or correlated enough, the activity during the sleep states could dominate the functional connectivity signal. Complicating mechanistic studies in mice is the fact that head-fixed mice do not close their eyes during NREM and REM sleep (Yüzgeç et al., 2018), meaning that without careful monitoring or a task it is possible that many studies examining ‘resting-state’ correlations in head-fixed mice were compounded by sleep.

Here we measured behavior, neural activity, blood volume and arteriole dilations from head-fixed mice during the awake state and NREM and REM sleep. We found that arteriole dilations and blood volume changes during NREM and REM sleep could be two to five times larger than those occurring in the awake animal. The correlations between neural activity and hemodynamic signals was greatly increased during NREM sleep, and the functional connectivity between interhemispheric regions of somatosensory cortex also increased.

We used intrinsic optical signal (IOS) (Huo et al., 2014; Winder et al., 2017) (14 mice, nine males) and 2-photon microscopy (Drew et al., 2011; Echagarruga et al., 2020; Shih et al., 2012a) (six mice, two males) in concert with electrophysiology to measure neural activity (Buzsáki et al., 2012; Harris et al., 2016) from the whisker representation of somatosensory cortex and the CA1 region of the hippocampus in un-anesthetized, head-fixed C57BL6/J mice (Figure 1A) during the light cycle. After mice were habituated to head-fixation, data was acquired from each mouse for 5–7 days. We obtained 357.2 total hours of data from these mice (mean: 23 ± 5.1 hr per mouse from IOS mice Supplementary file 1), 5.8 ± 2.0 hr per mouse from 2-photon imaged mice (Supplementary file 1). All experiments were performed during the animal’s light cycle. We tracked whisker position (Kleinfeld and Deschênes, 2011; O'Connor et al., 2010), body movement, and nuchal muscle EMG (Datta and Maclean, 2007; Veasey et al., 2000), as spontaneous ‘fidgeting’ behaviors drive a substantial portion of neural activity and hemodynamic signals in the awake mouse (Drew et al., 2019; Musall et al., 2019; Stringer et al., 2019; Winder et al., 2017), and these measures can be used to determine the arousal state of the animal. We recorded neural activity differentially across stereotrodes from the whisker representation of the somatosensory cortex and the hippocampus to reject non-local electrical signals (Buzsáki et al., 2012; Kajikawa and Schroeder, 2011; Nicholson and Freeman, 1975). For IOS data, we used a bootstrap aggregation random forest to determine the arousal state from these behavioral measures and hippocampal and cortical LFPs (see Materials and methods), categorizing every non-overlapping five second interval into one of three categories: rfc-Awake, rfc-NREM, or rfc-REM, where ‘rfc’ denotes the arousal state assigned via an automated random forest classifier. rfc-Awake periods were further characterized into awake rest (defined as periods lasting longer than 10 s that lack whisker and body movement), awake whisking (defined as bouts of whisking lasting between 2 and 5 s in duration), and awake stimulation (directed air puffs to contralateral whiskers). Each of these three awake behaviors were manually identified in rfc-Awake data to exclude transitional arousal states. Periods of contiguous NREM sleep and contiguous REM sleep were defined as events lasting at least 30 s and 60 s, respectively (Supplementary file 1). The 2-photon microscopy data was manually scored in a fashion similar to rfc-IOS data (Manual(m)-Awake, m-NREM, m-REM) and then further subdivided into awake rest, awake whisking, contiguous NREM, and contiguous REM (no whisker was stimulation presented during 2-photon experiments) (Supplementary file 1). As awake behaviors can be much shorter in duration than NREM and REM events, we used different durations when categorizing each arousal state. We chose the specific durations for each arousal state based on the typical minimum duration of each event, allowing us to catch the neural and vascular dynamics of each behavior. All reported values are mean ± standard deviation, unless otherwise indicated.

Figure 1 with 9 supplements see all

Download asset Open asset

Sleep drives larger fluctuations than awake behaviors

We first examined how arousal state affected hemodynamic signals using intrinsic optical signal imaging (Huo et al., 2014; Sirotin and Das, 2009; Vazquez et al., 2014), which detects changes in total hemoglobin [HbT] from changes in reflectance. Periods of awake rest without whisking or stimulation were set as the zero baseline (see Methods). Increases in blood volume (vasodilation) cause decreases in reflectance, which are converted into hemoglobin changes (∆[HbT]) using the Beer-Lambert law (Ma et al., 2016). These hemodynamic measurements were done bilaterally through a polished and reinforced thinned-skull window (Figure 1A; Drew et al., 2010a; Shih et al., 2012b) encompassing the whisker-related section of the somatosensory cortex (Petersen, 2007). The LFP was recorded from the vibrissa cortex within the window (Figure 1B,C) to provide a direct measurement of neural activity from whisker-related portions of somatosensory cortex from both hemispheres (Winder et al., 2017). To assist in arousal state classification, we also recorded the CA1 hippocampal LFP (Figure 1B,C), nuchal muscle electromyography (EMG), motion of the whiskers, body motion, and heart rate. Hemodynamic measurements were taken from a region of interest within a 1 mm diameter circle centered on the pixels that showed the highest cross-correlation between reflectance and gamma band power (1–2 s lag) (Winder et al., 2017) during the first 15–60 min of data (see Figure 1—figure supplement 1A–C). This period contains not just rest, but also whisking and some whisker stimulation. This differs from Winder et al., 2017, where only resting periods, lack whisking and stimulation, were used to define the ROI. These pixels corresponded to a region putatively within the vibrissa cortex and were consistent across imaging days. This is consistent with the hemodynamic point-spread function having full-width at half-max of several hundred microns (Vazquez et al., 2014) and the hyperemic response being conducted over several hundred microns (Rungta et al., 2018). The spatial scale of the measured hemodynamic signal is very close to that measured by the electrodes.

At the beginning of each day’s imaging session, we stimulated the vibrissae on either side with a brief puff of air, which drove canonical neural and vascular responses (Figure 1D). Consistent with previous work (Winder et al., 2017), gamma band power in the somatosensory cortex (30–100 Hz) increased by 78.1 ± 66.7%, followed by a 16.8 ± 5.1 µM increase in [HbT] (corresponding to a −2.4 ± 0.7% decrease in reflectance, ∆R/R, see Figure 1—figure supplement 2). Volitional, awake whisking led to an increase in gamma band power and [HbT] that increased with the duration of the whisking event. Brief (0.5–2 s), moderate (2–5 s), and extended (>5 s) whisking events lead to an increase of 5 ± 4.8%, 9.7 ± 7%, and 16.5 ± 10.8% in gamma band power in the somatosensory cortex. These changes led to an increase in [HbT] of 2.5 ± 1.4 µM, 7.1 ± 3.3 µM, 12.1 ± 5.6 µM respectively (see Figure 1—figure supplement 3).

In comparison to the sensory-evoked responses, we found much larger changes in [HbT] associated with sleep. NREM sleep is characterized by low EMG activity, lack of whisker and body movement, and pronounced power in the low-frequency bands of the cortical LFP (Figure 1E–J, Videos 1–3; Scammell et al., 2017; Vyazovskiy and Harris, 2013). REM sleep is characterized by neck muscle atonia, sporadic whisker movement, and increased theta band power in the hippocampus (Figure 1E–J; Hobson and Pace-Schott, 2002; Walker and Stickgold, 2004). During contiguous NREM sleep, we observed large oscillations in total hemoglobin, up to 87.3 ± 9.9 µM [HbT] in peak-to-peak amplitude, compared to 32.3 ± 4.4 µM during awake rest (generalized linear mixed-effects (GLME), p<9 × 10⁻³²). During contiguous REM bouts, there was a prolonged (>30 s) increase of [HbT] of up to 107.7 ± 13.3 µM, greatly surpassing the maximum of 17 ± 3.2 µM during awake rest (GLME, p<3.5 × 10⁻⁵⁵, see Figure 1—figure supplement 4A,B). Examples of sleep-related changes in blood volume are shown in Figure 1E–J and Figure 1—figure supplement 5, Figure 1—figure supplement 6, Figure 1—figure supplement 7, and Figure 1—figure supplement 8. Note that as previously observed (Yüzgeç et al., 2018), the eyes of the mouse are open during both REM and NREM sleep states (Figure 1G). These results show that sleep in head-fixed mice is associated with large increase in blood volume, particularly during REM sleep.

Video 1

Download asset

Video 2

Download asset

Video 3

Download asset

Mice regularly enter into NREM and REM sleep during head-fixation

We then asked how often head-fixed mice sleep, and how the probability of sleep changes over time from the start of head-fixation. A hypnogram for a single mouse over 6 days with 5 s resolution is shown in Figure 2A. The white lines denote brief breaks in recording while data is saved. From this example, it is clear that the mouse has many periods of NREM and REM sleep interspersed with awake. Periods of REM sleep canonically follow NREM sleep (Saper et al., 2010). REM periods are typically followed by awakening, though REM sleep can be followed by NREM periods. The mean percentage of each classified state is shown in Figure 2B, and the breakdown for individual animals in Figure 2C. We noted no clear difference in (rfc-)Awake:NREM:REM sleep ratios between males and females, so they were pooled for other analyses. Plotting the probability of finding the mouse in each of the states as a function of time throughout the imaging session (Figure 2D) shows that as time goes by the mouse is more likely to be asleep. As awake mice typically whisk every ~10 s, we quantified the probability that the mouse had fallen asleep after a given period lacking whisking and body movement (‘Rest’) across all animals. We found that during only ~50% of ‘resting’ events lasting 10–15 s were the mice awake for the whole event, and with longer ‘resting’ events showing even lower probability of wakefulness (Figure 2E). This result is reminiscent of studies in humans showing the probability of being awake falls rapidly with time during a resting-state fMRI scan (Tagliazucchi and Laufs, 2014), though human sleep/wake behavior is much less fragmented than in the mouse, and the transition times are correspondingly longer. EMG activity (Figure 2F) was much lower during sleep than the awake state (see Figure 2—figure supplement 1A: rfc-NREM: p<3.2 × 10⁻¹², rfc-REM: p<1.7 × 10⁻²¹, GLME), though the amount of whisking (quantified as the variance in the whisker angle) was much smaller during rfc-REM and rfc-NREM sleep than during the awake state (Figure 2G, see Figure 2—figure supplement 1B: rfc-NREM: p<1.3 × 10⁻¹², rfc-REM: p<0.003, GLME). The heart rate was quantified by finding the peak in the intrinsic signal power spectrum between 5–15 Hz (Huo et al., 2015a; Huo et al., 2015b; Winder et al., 2017). Heart rate was lowest during rfc-NREM sleep (Figure 2H, see Figure 2—figure supplement 1C: rfc-NREM: p<1.8 × 10⁻¹³, rfc-REM: p<0.004, GLME), though heart rate during contiguous REM was comparable to heart rate in the awake, resting mouse, and the heart rate was elevated during awake whisking (Figure 2I, Awake Whisk: p<1.1 × 10⁻¹⁰, contiguous NREM: p<3.8 × 10⁻¹⁴, contiguous REM: p<0.0009, GLME). These observations demonstrate the prevalence of sleep in ‘resting-state’ data in head-fixed mice and show how unimodal measures, such as whisker movement or heart rate alone, are insufficient to detect these sleep states.

Figure 2 with 2 supplements see all

Download asset Open asset

Sleep drives arteriole dilatations much larger than those seen in the awake brain

The intrinsic optical signal contains contributions from arteries, veins, and capillaries (Huo et al., 2015a; Huo et al., 2015b; Zhang et al., 2019), so to better understand how arterioles changes contribute to the signal, we used two-photon microscopy (Shih et al., 2012a) to image pial and penetrating arterioles (25 and 4 arterioles respectively, from six mice) (Figure 3A). We recorded hippocampal and cortical LFP contralateral to the imaging window (Figure 3B,C). Two-photon imaging in awake mice has shown that sensory stimulation and locomotion drive dilations of arteries of approximately 20% above the baseline diameter (Drew et al., 2011; Echagarruga et al., 2020; Gao and Drew, 2016; Huo et al., 2015a; Huo et al., 2015b) and whisking drives dilations of 5–10% (Drew et al., 2020). We observed a similar dilation with spontaneous whisking as a function of whisking duration. Brief (0.5–2 s), moderate (2–5 s), and extended (>5 s) whisking events led to an increase in vessel diameter of 0.8 ± 0.9%, 8.3 ± 4.9%, and 10.9 ± 4.3% respectively (Figure 3D, see Figure 3—figure supplement 1). Examples of sleep-related changes in arteriole diameter is shown in Figure 3E–I and Figure 3—figure supplement 2, Figure 3—figure supplement 3, Figure 3—figure supplement 4, and Figure 3—figure supplement 5. In comparison to the awake state, arteriole diameter during NREM sleep follows a low-frequency dilation/constriction with peak dilations that can exceed those seen during moderate whisking (Figure 3D). During REM sleep, the arterioles slowly dilate over tens of seconds, and can reach peak dilations in excess of 50% of the baseline diameter during periods of awake rest. The dilation amplitudes during NREM and REM sleep dwarf those seen in the awake animals. Penetrating arterioles had similar sleep-wake dynamics as the pial arterioles (see Figure 3—figure supplement 2 vs Figure 3—figure supplement 3) and were combined into a single group for all analysis. These results show that both NREM and REM sleep cause pronounced cortical arteriole dilations. Contiguous NREM was associated with peak-to-peak diameter changes of 38 ± 15.8%, compared to 16.6 ± 4% seen during awake rest (GLME, p<3.6 × 10⁻¹⁰). Contiguous REM drove a prolonged ramp-up to peak dilation of 49.9 ± 9.1% compared to a peak dilation of 7.4 ± 2.6% during the awake resting baseline (GLME, p<2 × 10⁻²⁴, see Figure 1—figure supplement 4C,D). These large dilations were followed by a pronounced and rapid constriction upon waking.

Figure 3 with 5 supplements see all

Download asset Open asset

Hemodynamic changes associated with transitions between arousal states

To quantify the dynamics of sleep-related changes in blood volume and arteriole diameter, we looked at the dynamics of these signals, as well as LFP and EMG signals, aligned to the arousal state transition (Figure 4A). We used transitions between two arousal states, where time in each arousal state was at least 30 s in duration. Mammals will typically progress through the Awake-NREM-REM-Awake pattern of the sleep cycle (Saper et al., 2010), even if the awake periods between the end of a REM event and the initiation of the subsequent NREM period are very brief. The transition from the rfc-Awake state into rfc-NREM (Figure 4B) shows an increase in total hemoglobin from the baseline of very low [HbT] in the awake state, up to 30 µM over the course of 30 s. During this time, the EMG power decreases with similar temporal dynamics. The LFP power in the vibrissa cortex shows increased power in the delta band [1–4 Hz] of around 300%. The transition from rfc-NREM into the rfc-Awake state (Figure 4C) was largely a temporally reversed version of the awake to NREM transition, where [HbT], LFP, and EMG signals all quickly return to baseline values as the animal wakes up. During the transition from rfc-NREM to rfc-REM (Figure 4D) there was a slow increase in total hemoglobin from 45 µM up to 75 µM. During the NREM-REM transition, the muscles become atonic and the EMG power decreased even more. The theta band [4–10 Hz] power in the hippocampal LFP increased by approximately 300%. The transition from rfc-REM into the rfc-Awake state (Figure 4E) was the largest in terms of the magnitude of the total hemoglobin change, as the large blood volume increase seen during REM rapidly reverses (within seconds) as the animal wakes up. Transitions from awake to REM, as well as REM to NREM are possible, but are much less common, and did not occur often enough to quantify reliably. The amplitude and temporal dynamics of ∆[HbT] transitions between each arousal state were consistent across imaging days (see Figure 4—figure supplement 1).

Figure 4 with 1 supplement see all

Download asset Open asset

When looking at single arterioles during the transition from NREM to REM, arterioles will on average go from an approximate 20% dilation, up to approximately 40% dilation after a minute or so in REM (Figure 4F). Upon waking from REM, arterioles constrict back to the baseline diameter within a few seconds (Figure 4G). When the arterial diameter and ∆[HbT] are plotted together, they show very similar temporal dynamics during the NREM to REM transition and during the REM to awake transition (Figure 4F,G). The close match in dynamics suggests that arteriole dilations are a significant driver of the blood volume changes during sleep, as is seen in the awake brain (Huo et al., 2015a; Rungta et al., 2018).

Cortical hemodynamic signals increase during NREM and REM sleep

We quantitatively compared hemodynamic signals during different arousal states. For this quantification, we used awake resting events ≥ 10 s in duration, awake whisking events 2–5 s in duration, brief stimulations of the whiskers, contiguous NREM sleep events ≥ 30 s in duration, and contiguous REM sleep events ≥ 60 s in duration. The average ∆[HbT] during each day’s awake resting condition was set as zero (n = 14 mice, 28 hemispheres). Awake whisking events between 2 and 5 s in duration caused a slight increase in [HbT] of 3.1 ± 2.5 µM (GLME, p<0.11). The increase caused by whisking and fidgeting, which are the primary drivers of resting-state neural and hemodynamic signals in the awake mouse (Drew et al., 2019; Drew et al., 2020; Musall et al., 2019; Stringer et al., 2019; Winder et al., 2017), were dwarfed by those that occurred during contiguous NREM sleep (32.2 ± 11.2 µM, GLME, p<1 × 10⁻³⁴) and during contiguous REM sleep (77.1 ± 11.9 µM, GLME, p<1 × 10⁻⁷⁶). These sleep-driven changes were much larger than sensory-evoked changes generated by contralateral whisker stimulation (12.9 ± 5.6 µM, GLME, p<4.9 × 10⁻¹⁰). The probability distribution of ∆[HbT] during each arousal state is shown in Figure 5D. There was a clear and pronounced separation in the hemodynamic signals between the sleep and awake states, as well as between contiguous NREM and REM sleep. However, the increases in [HbT] during sleep were smaller than those caused by the vasodilator isoflurane (Flynn et al., 1992; Gao et al., 2015; Figure 5—figure supplement 1), showing the vessels are not maximally dilated during sleep.

Figure 5 with 1 supplement see all

Download asset Open asset

The average dilation of arterioles during each arousal state followed the same trend as the ∆[HbT] data (Figure 5B). Awake whisking events spurred an average diameter increase of 4.2 ± 3.2% (GLME, p<0.0006, n = 6 mice, 29 arterioles), with contiguous NREM sleep dilating 9.4 ± 9.3% (GLME, p<8.7 × 10⁻¹⁹, n = 6 mice, 21 arterioles) and contiguous REM sleep reaching dilations of 32.7 ± 8.6% (GLME, p<5 × 10⁻³⁴, n = 5 mice, 10 arterioles). The probability distribution of arteriole ∆D/D during each arousal state with data taken from all animal’s events is shown in Figure 5E Recordings of changes in blood flow measured using laser Doppler flowmetry (∆Q/Q) followed the same trend (Figure 5C) (n = 8 mice, one mouse was excluded due to a weak signal), with the largest flow increases during sleep. Awake volitional whisking events caused negligible increases of 0.4 ± 2.6% (GLME,  p<0.91). Flow changes were more prominent during contiguous NREM sleep, 14.4 ± 9.1% (GLME, p<0.0004), and even more-so during contiguous REM sleep, 28.3 ± 12.8% (GLME, p<1.5 × 10⁻⁸). The probability distribution of ∆Q/Q during each arousal state with data taken from all animal’s events is shown in Figure 5F. These measurements of blood volume, arteriole diameter, and flow all show a consistent and pronounced trend of vasodilation in the somatosensory cortex during sleep which far surpassed those seen during volitional behaviors and sensory-evoked stimulation in the awake animal.

Neurovascular coupling is strongest during NREM sleep

It has previously been shown that neurovascular coupling is similar across awake arousal states (Winder et al., 2017) and that both spontaneous and sensory-evoked hemodynamics are most strongly correlated with gamma band power and multi-unit average (MUA – a measure of local spiking activity) in the absence of overt stimulation (Mateo et al., 2017; Schölvinck et al., 2010; Shmuel and Leopold, 2008; Winder et al., 2017). To see if these relationships held true for neurovascular coupling in other arousal states, we looked at the relationship between the power in different frequency bands of the LFP or multi-unit activity and ∆[HbT] during periods of contiguous NREM and REM sleep. Correlations between the MUA and ∆[HbT] as well as the LFP and ∆[HbT] during awake rest showed a peak correlation of 0.24 ± 0.06 for the MUA and 0.16 ± 0.04 for the gamma band, consistent with previous work (Winder et al., 2017; Mateo et al., 2017). The hemodynamic response lagged the MUA by 1.16 ± 0.12 s and the gamma band power by 1.15 ± 0.13 s, respectively (Figure 6A, see Figure 6—figure supplement 1). The peak cross-correlation during contiguous NREM sleep was nearly double that of awake rest: MUA: 0.44 ± 0.06 (GLME, p<1.2 × 10⁻²³ compared with awake rest); gamma band: 0.3 ± 0.05 (GLME, p<8 × 10⁻²³ compared with awake rest). During NREM, there were similar dynamics in the lag of the hemodynamic response to the MUA as during awake rest (1.28 ± 0.11 s, GLME, p<0.08) and for the gamma band power (1.26 ± 0.12 s, GLME, p<0.15) (Figure 6B). The peak cross-correlations during contiguous REM sleep were significantly higher than those during awake rest for the MUA (0.32 ± 0.07, GLME, p<4.8 × 10⁻⁷), but not for gamma band power (0.15 ± 0.05; GLME, p<0.22). During REM there was a slight, but significant increase in the lag of the hemodynamic response to the MUA (1.49 ± 0.49 s, GLME, p<1.2 × 10⁻⁵) and for the gamma band power (1.4 ± 0.47 s; GLME, p<0.0016) (Figure 6C). However, the distribution of correlations across frequencies differed across states. In addition to a highly correlated gamma band, the correlations during contiguous NREM extended down into the beta band [13–30 Hz]. These results show that the correlation between neural activity and cerebral blood volume change markedly across states, though the hemodynamic lag differences are of order a few hundred milliseconds. We found weak to negative correlations for frequencies below 30 Hz in the awake state and during contiguous REM. In contrast, there were strong positive correlations between lower frequencies of LFP power and the [HbT] in contiguous NREM sleep. This suggests that the positive correlations between lower frequencies of the LFP and vasodilation seen in some studies may be due to the inclusion of NREM sleep in the ‘resting-state’.

Figure 6 with 1 supplement see all

Download asset Open asset

Neural and hemodynamic correlations across hemispheres increase during sleep

We then explored how the correlations and coherency of left and right hemisphere [HbT] and gamma band power were affected by arousal state. When quantifying signals in the frequency domain, the lowest resolvable frequency will be the inverse of the event duration. Because there were very few periods of awake rest greater than 10 s in duration (see Figure 2E), the lowest frequency for resting periods we could characterize was 0.1 Hz (Mitra and Pesaran, 1999). We thus examined all 15 min periods from our data without whisker stimulation that had at least 12 min (>80%) of rfc-Awake arousal state classifications as alert, as an extension of the awake resting and awake whisking arousal state into the lower frequencies. We also examined all 15 min periods without whisker stimulation that had at least 12 min (>80%) of rfc-NREM or rfc-REM arousal state classifications as asleep to extrapolate the contiguous NREM/REM arousal states into the lower frequencies. Note that the asleep condition contains some awake state and transitions, which may explain the lower power in the gamma band envelope above 0.1 Hz than the pure REM and NREM conditions. Lastly, we took every 15 min period without whisker stimulation, regardless of random forest classified arousal state, as all data, which has an arousal state distribution similar to Figure 2B and is more representative of the type of data obtained if behavioral monitoring and arousal state classification were not used.

The power spectrum of ∆[HbT], normalized to the peak power in the resting condition, is shown in Figure 7D (see Table 1, see Figure 7—figure supplement 1E,F). For [HbT], there was substantially more power in the lower frequencies for nearly all arousal states. The lower power of the ∆[HbT] power spectra at higher frequencies were consistent with the lowpass nature of the hemodynamic response (Drew, 2019; Silva et al., 2007; Winder et al., 2017; de Zwart et al., 2005). We examined the power in two frequency bands, 0.1 Hz (the ‘vasomotion’ frequency Mateo et al., 2017), and the ultra-low 0.01 Hz band. The power at 0.1 Hz was highest during contiguous NREM (awake rest: 0.9 ± 0.1 (A.U.); contiguous NREM: 6.7 ± 3.4, GLME vs. awake rest, p<7.6 × 10⁻¹⁹), and was also higher during contiguous REM than during awake rest (contiguous REM: 4.1 ± 2.6, GLME vs. awake rest, p<1.3 × 10⁻⁷). Extending into the ultra-low frequencies, power in the ∆[HbT] signal at 0.01 Hz was significantly higher in the asleep state than the alert state (176.2 ± 98.9 vs. 36.5 ± 28.6, GLME, p<5.3 × 10⁻¹⁰). Similar low-frequency dominated power spectra were seen in the diameters of arterioles, with the power at 0.1 Hz during contiguous NREM (4.5 ± 3.9, GLME, p<2.3 × 10⁻⁸) and contiguous REM (5.1 ± 5.7, GLME, p<2.7 × 10⁻⁷) also exceeding those seen during awake rest (0.9 ± 0.1). There was not enough single arteriole data that met the inclusion criteria for the asleep condition for 0.01 Hz (Figure 7G, Table 1, see Figure 7—figure supplement 1I,J).

Figure 7 with 3 supplements see all

Download asset Open asset

Table 1

Spectral Power at 0.1 Hz	Awake Rest	Cont. NREM	Cont. REM	Alert	Asleep	All Data
Gamma band Power (a.u.)	0.9 ± 0.1	13.3 ± 37.4 (p<0.01)	12.8 ± 23.4 (p<0.02)	1.3 ± 0.8 (p<0.93)	6.2 ± 8 (p<0.29)	3.6 ± 2.5 (p<0.58)
∆[HbT] Power (a.u.)	0.9 ± 0.1	6.7 ± 3.4 (p<7.6 × 10−19)	4.1 ± 2.6 (p<1.3 × 10−7)	2.2 ± 1.3 (p<0.03)	4.7 ± 2.3 (p<1.3 × 10−9)	3.8 ± 1.6 (p<1.3 × 10−6)
∆D/D Power (a.u.)	0.9 ± 0.1	4.5 ± 3.9 (p<2.7 × 10−10)	3.8 ± 3.7 (p<2.1 × 10−5)	3.4 ± 2.1 (p<9.3 × 10−7)		3.4 ± 1.9 (p<5.8 × 10−7)

1. Mean ± 1 standard deviation. p-values as a comparison to ‘Rest’. Shaded region indicates insufficient data for that arousal state.

Table 2

Spectral Power at 0.01 Hz	Alert	Asleep	All Data d
Gamma band Power (a.u.)	5 ± 3.2	34 ± 51 (p<0.001)	21.7 ± 20.6 (p<0.06)
∆[HbT] Power (a.u.)	36.5 ± 28.6	176.2 ± 98.9 (p<5.3 × 10−10)	128.9 ± 63.8 (p<7.9 × 10−6)
∆D/D Power (a.u.)	34.9 ± 25		45.6 ± 33.6 (p<0.06)

1. Mean ± 1 standard deviation. p-values as a comparison to ‘Alert’. Shaded region indicates insufficient data for that arousal state.

Previous functional imaging studies have observed marked increases in hemodynamic activity and functional connectivity between brain regions during periods of sleep (Boly et al., 2012; Dang-Vu et al., 2008; Fukunaga et al., 2006; Horovitz et al., 2008; Larson-Prior et al., 2009; Mitra et al., 2015), but how these hemodynamic changes relate to neural activity was unknown.

We then looked at gamma band power, which of all neural signals is most closely related to hemodynamic changes in the awake brain (Echagarruga et al., 2020; Mateo et al., 2017; Schölvinck et al., 2010; Winder et al., 2017). The fluctuations in the gamma band power had substantial power at all frequencies (Figure 7A, Table 1, see Figure 7—figure supplement 1A,B), consistent with previous recordings in primates (Leopold et al., 2003). Fluctuations in the gamma band power envelope at frequencies near 0.1 Hz have been proposed to drive vasomotion and underly functional connectivity (Mateo et al., 2017), and lower frequency fluctuations in the envelope (in the ~0.01 Hz range) could underly changes in connectivity at these frequencies (Chang and Glover, 2010). The power in the gamma band envelope at 0.1 Hz was substantially higher during sleep (awake rest: 0.9 ± 0.1; contiguous NREM sleep: 13.3 ± 37.4, GLME, p<0.01; contiguous REM Sleep: 12.8 ± 23.4, GLME, p<0.02). All power normalized by the peak in the power spectrum during rest. There was a similar increase in the gamma band power fluctuations at 0.01 Hz during sleep relative to the awake condition (alert: 5 ± 3.2; asleep: 34 ± 51, GLME, p<0.001). During sleep, the modulations of gamma band power at 0.1 and 0.01 Hz were much larger than during the awake state.

To quantify the relationship of neural activity or [HbT] between the left and right hemispheres, we looked at the magnitude of the coherence squared (MoC2) (Figure 7B, Tables 3 and 4, see Figure 7—figure supplement 1C,D). We used MoC2 because it tells us the amount of variance explained in one signal by the other for any given frequency, making MoC2 equivalent to the R² between two signals at any given frequency (Drew et al., 2020). The MoC2 between left and right hemisphere gamma band power at 0.1 Hz was relatively low. In the 0.01 Hz frequency band, left-right somatosensory cortex gamma band power MoC2 was higher. Interestingly, for all conditions, the MoC2 for the gamma band power was substantially lower than the MoC2 of [HbT] in the same frequency band (Tables 3 and 4).

Table 3

Coherence2 at 0.1 Hz	Awake Rest	Cont. NREM	Cont. REM	Alert	Asleep	All Data
Gamma band (Coherence2)	0.08 ± 0.06	0.26 ± 0.14 (p<4.8 × 10−8)	0.16 ± 0.1 (p<0.006)	0.13 ± 0.08 (p<0.04)	0.2 ± 0.11 (p<0.0002)	0.19 ± 0.09 (p<0.0004)
∆[HbT] (Coherence2)	0.61 ± 0.11	0.86 ± 0.05 (p<1.2 × 10−16)	0.77 ± 0.08 (p<5.7 × 10−10)	0.69 ± 0.14 (p<0.002)	0.82 ± 0.07 (p<2.9 × 10−13)	0.78 ± 0.06 (p<2.6 × 10−10)

1. Mean ± 1 standard deviation. p-values as a comparison to ‘Rest’.

Table 4

Coherence2 at 0.01 Hz	Alert	Asleep	All Data
Gamma band (Coherence2)	0.38 ± 0.27	0.69 ± 0.18 (p<7.2 × 10−6)	0.61 ± 0.21 (p<0.0002)
∆[HbT] Power (Coherence2)	0.82 ± 0.15	0.97 ± 0.02 (p<1.5 × 10−5)	0.95 ± 0.03 (p<7.5 × 10−5)

1. Mean ± 1 standard deviation. p-values as a comparison to ‘Alert’.

The MoC2 between left and right cortical somatosensory cortex ∆[HbT] was uniformly high (Figure 7E, Tables 3 and 4, see Figure 7—figure supplement 1G,H), though it was higher during contiguous NREM sleep (0.86 ± 0.05, GLME, p<1.2 × 10⁻¹⁶) and contiguous REM sleep (0.77 ± 0.08, GLME, p<5.7 × 10⁻¹⁰) as compared to awake rest (0.61 ± 0.11). low-frequency ∆[HbT] MoC2 during the asleep state reached near-unity at 0.01 Hz (0.97 ± 0.02, GLME, p<1.5 × 10⁻⁵), continuously surpassing those seen during the alert state (0.82 ± 0.15). To compare these results to functional connectivity measures, we evaluated the Pearson’s correlation coefficient between bilateral ∆[HbT] during each arousal state (Figure 7F, Table 5). The correlation between the left and right hemisphere ∆[HbT] was 0.74 ± 0.06 during the awake resting condition. They were elevated during awake whisking 2–5 s in duration (0.82 ± 0.09, GLME, p<4.8 × 10⁻⁸), during periods of contiguous NREM (0.9 ± 0.03, GLME, p<9.8 × 10⁻²⁰) and contiguous REM (0.91 ± 0.03, GLME, p<1.4 × 10⁻²¹). These correlations, however, were not fully explained by the correlations in underlying neural activity, which were much lower. This may be explained by a lower signal-to-noise ratio for neural signals. The correlation coefficients for the envelope of bilateral gamma band power [30–100 Hz] (Figure 7C, Table 5) were relatively small during rest (0.16 ± 0.07), and were not significantly elevated during awake whisking (0.15 ± 0.05, GLME, p<0.47) or during contiguous REM sleep (0.21 ± 0.13, GLME, p<0.11). However, correlations in gamma band power doubled during contiguous NREM sleep (0.3 ± 0.12, GLME, p<6.6 × 10⁻⁷). Analysis of the cross-hemisphere delta band [1–4 Hz], theta band [4–10 Hz], alpha band [10–13 Hz], and beta band [13–30 Hz] power showed similar trends as the gamma band, with elevated power, MoC2, and Pearson’s correlations during contiguous NREM/REM sleep compared to awake rest (see Figure 7—figure supplement 2, Figure 7—figure supplement 3, Supplementary file 1). The exception was the delta band, where the sleeping MoC2 and Pearson’s correlations were not statistically different than those during rest. These results show large changes in gamma band power during sleep, but relatively low cross-hemisphere correlations in neural activity. In contrast, the correlation coefficients for the left and right ∆[HbT] were uniformly high (Table 3), and significantly elevated by awake whisking. Contiguous NREM and REM sleep substantially elevated the ∆[HbT] correlations.

Table 5

Pearson’s Correlation Coef.	Awake Rest	Whisking	Cont. NREM	Cont. REM	Alert	Asleep	All Data
Gamma band (R)	0.16 ± 0.07	0.15 ± 0.05 (p<0.47)	0.3 ± 0.12 (p<6.1 × 10−7)	0.21 ± 0.13 (p<0.11)	0.26 ± 0.1 (p<0.0004)	0.35 ± 0.14 (p<4.7 × 10−10)	0.33 ± 0.11 (p<8.5 × 10−9)
∆[HbT] (R)	0.74 ± 0.06	0.82 ± 0.09 (p<3 × 10−8)	0.9 ± 0.03 (p<3.6 × 10−20)	0.91 ± 0.03 (p<4.9 × 10−22)	0.87 ± 0.07 (p<7.9 × 10−15)	0.96 ± 0.02 (p<8.1 × 10−28)	0.93 ± 0.02 (p<1.1 × 10−24)

1. Mean ± 1 standard deviation. p-values as a comparison to ‘Rest’.

All and all, these result show that sleep enhances the already strong cross-hemisphere coherency and correlation in ∆[HbT]. The increases in ∆[HbT] during sleep are larger than those during whisking/fidgeting behavior, which are the primary drivers of spontaneous hemodynamic signals in the awake mouse (Drew et al., 2019; Drew et al., 2020; Winder et al., 2017). Sleep increases the cross-hemisphere correlations and coherency of the LFP envelope as well, but the neural correlations are substantially lower than the hemodynamic correlations.

Impact of arousal state on total blood volume and neural-vascular coherence

To understand how well we can determine the arousal state (Awake/NREM/REM) from the changes in hemodynamic measurements (∆[HbT] or ∆D/D), we plotted the probability that the animals were in a given arousal state for a given amount of ∆[HbT] (Figure 8A) or ∆D/D (Figure 8B). Using only [HbT], the strength of the coupling between spontaneous neural activity and changes in blood volume is relevant to resting-state fMRI. Many studies have looked into this coupling (Ma et al., 2016; Mateo et al., 2017; Schölvinck et al., 2010; Shmuel and Leopold, 2008; Winder et al., 2017) but the measured strength of this coupling varies from study to study (Drew et al., 2020). One possible contribution to the varying strength of the coupling observed between neural and hemodynamic signals could be changes in the arousal state. When solely looking at changes in arteriole diameter (∆D/D), the variability in arteriole diameter during NREM sleep make it difficult to deduce arousal state. Fluctuations in arteriole diameter during NREM may appear similar to volitional whisking events, which occur every 10–15 s (Drew et al., 2020; Winder et al., 2017).

Figure 8 with 1 supplement see all

Download asset Open asset

The coherence between [HbT] and the gamma band power in various arousal states is shown in Figure 8C (other LFP bands are shown in Figure 8—figure supplement 1). Coherence at lower frequencies was highest in contiguous NREM sleep and for all data, and lowest during awake rest and during contiguous REM sleep. It may seem puzzling that the coherence is higher when taking all the data together than for any individual subset, but this can be explained by a restriction of range effect when grouping by arousal state. In Figure 8D, we show the ∆[HbT] versus the power in the gamma band for all of our ‘rfc’ data including stimulation. When all the data is considered together, it is clear that there was a robust relationship between neural activity and blood volume. However, as the level of neural activity and [HbT] were different across the arousal states, within any given state this relationship may not be so strong, again due to a restriction of range effect. All in all, we would like to emphasize the observation that correlations between ongoing neural activity and hemodynamic signals appear to be lowest during the awake state and highest during sleep (Figure 8e) and is most prominent in the lower frequencies.

Using optical imaging and electrophysiology in head-fixed mice, we found that sleep was associated with large increases in blood volume and arterial dilation in the somatosensory cortex. These vasodilations and increases in cerebral blood volume were very large, up to five times the size of those evoked by sensory stimulation or fidgeting behaviors like whisking. We found that the correlation and coherence between neural activity and blood volume were substantially stronger during NREM sleep than during any other arousal state, as were bilateral correlations and coherency. Because, at least under our imaging conditions, mice frequently and rapidly enter into sleep, the hemodynamic signals due to sleep will dominate over any awake signals due to the very large changes in neural activity and blood volume in the sleeping cortex. The net effect is that, without careful monitoring of arousal state, sleep-related hemodynamics and neural activity will dominate any ‘resting-state’ study in the un-anesthetized mouse.

We note that care should be taken in interpreting these results. The sleep patterns and depth may be altered by head-fixation. We only imaged cortical blood volume changes, and other areas of the brain may have different patterns of vasodilation (Braun et al., 1997; Townsend et al., 1973). We used electrophysiology to assay neural activity as it is used in humans (Buzsáki et al., 2012), and calcium indicators buffer intracellular calcium, can cause epilepsy, and do not report the majority of spikes even under optimal imaging conditions (McMahon and Jackson, 2018; Steinmetz et al., 2017; Theis et al., 2016). The activity of interneurons that express neuronal nitric oxide synthase (nNOS) can drive vasodilation without causing detectable changes in the LFP (Echagarruga et al., 2020; Krawchuk et al., 2020; Lee et al., 2020), so the overall neural activity captured in the LFP may not reflect the activity of neurons that drive hemodynamic responses. We note there may be a difference in sleep-state identification with local LFP recordings when compared with conventional EEG, which sum electrical signals over a larger area. Previous studies in rats have shown discrete regions of cortex going into periods of ‘localized sleep’ (Vyazovskiy et al., 2011), and in some cases, slow-wave oscillations can be seen during REM sleep in the cortex (Funk et al., 2016), either of which may reduce the accuracy of our arousal state classifications. We noted no arousal state discrepancies between left and right somatosensory cortices during manual sleep scoring, however, this does not exclude localized sleep in other cortical regions. The spatial resolution of IOS and our electrodes are of order 100 µm, and any vascular or neural changes on smaller spatial scales would be ‘blurred’, potentially resulting in lower estimates of neurovascular coupling strength. Additionally, as hemoglobin is the major absorber of both visible and infrared light in the brain, and changes in the level of hemoglobin can cause decreases in florescence signals (Haiss et al., 2009; Shen et al., 2012), sleep-related vasodilations attenuate signals from calcium indicators in neurons when visualized with either 1- or 2-photon methods. As there are multiple pathways by which neurons (and astrocytes) communicate with the vasculature (Attwell et al., 2010; Drew, 2019; Iadecola, 2017), it is likely that multiple mechanisms underly the vasodilation we observed during sleep (Özbay et al., 2018).

There are several implications for our work. First, it shows that it is critical to monitor arousal state, particularly in head-fixed mice, as they frequently fall asleep. Arousal state cannot be detected by simple monitoring of the face/eyes, as mice can sleep with their eyes open (Yüzgeç et al., 2018), and thus requires additional methods such as a combination of EMG/LFP/whisker tracking. If an animal falls asleep, this will result in larger neurovascular and bilateral correlations than in the awake state. There has been some disagreement as to how well hemodynamic signals track neural activity in the absence of any sensory stimulation (the ‘resting-state’) (Drew et al., 2020). With a few exceptions (e.g. Winder et al., 2017), these studies did not carefully monitor arousal state, and it could be that studies finding higher neurovascular correlations have episodes of sleep in them. Similar issues have been noted for human resting-state fMRI studies (Tagliazucchi and Laufs, 2014) and with non-human primates (Cardoso et al., 2019; Chang et al., 2016), so this is likely a ubiquitous problem with any experiment where the subject is unmotivated to stay awake. On a more fundamental level, these results show that the vasodilation seen in the awake brain is smaller than the hemodynamic changes seen during sleep. This vasodilation during sleep may serve to circulate cerebrospinal fluid (Aldea et al., 2019; Fultz et al., 2019; Kedarasetti et al., 2020; van Veluw et al., 2020) or some other physiological role.

Source data and code for generation of all figures can be found here: Code repository location:https://github.com/DrewLab/Turner_Gheres_Proctor_Drew_eLife2020 (copy archived at https://archive.softwareheritage.org/swh:1:rev:a8318fc5cb29b88504fe72f5d3d80867bd9791f2/) Data repository location: https://doi.org/10.5061/dryad.6hdr7sqz5.

1. 1. Drew P
   2. Turner K
   3. Gheres K
   4. Proctor E

   (2020) Dryad Digital Repository

   Neurovascular coupling and bilateral connectivity during NREM and REM sleep.

   https://doi.org/10.5061/dryad.6hdr7sqz5

1. 1. Adams MD
   2. Winder AT
   3. Blinder P
   4. Drew PJ

   (2018) The pial vasculature of the mouse develops according to a sensory-independent program

   Scientific Reports 8:1–12.

   https://doi.org/10.1038/s41598-018-27910-3

   • Google Scholar
2. 1. Aldea R
   2. Weller RO
   3. Wilcock DM
   4. Carare RO
   5. Richardson G

   (2019) Cerebrovascular smooth muscle cells as the drivers of intramural periarterial drainage of the brain

   Frontiers in Aging Neuroscience 11:1–17.

   https://doi.org/10.3389/fnagi.2019.00001

   • PubMed
   • Google Scholar
3. 1. Amzica F
   2. Steriade M

   (1998) Electrophysiological correlates of sleep Delta waves

   Electroencephalography and Clinical Neurophysiology 107:69–83.

   https://doi.org/10.1016/S0013-4694(98)00051-0

   • PubMed
   • Google Scholar
4. 1. Anafi RC
   2. Kayser MS
   3. Raizen DM

   (2019) Exploring phylogeny to find the function of sleep

   Nature Reviews Neuroscience 20:109–116.

   https://doi.org/10.1038/s41583-018-0098-9

   • PubMed
   • Google Scholar
5. 1. Attwell D
   2. Buchan AM
   3. Charpak S
   4. Lauritzen M
   5. Macvicar BA
   6. Newman EA

   (2010) Glial and neuronal control of brain blood flow

   Nature 468:232–243.

   https://doi.org/10.1038/nature09613

   • PubMed
   • Google Scholar
6. 1. Bergel A
   2. Deffieux T
   3. Demené C
   4. Tanter M
   5. Cohen I

   (2018) Local hippocampal fast gamma rhythms precede brain-wide hyperemic patterns during spontaneous rodent REM sleep

   Nature Communications 9:5364.

   https://doi.org/10.1038/s41467-018-07752-3

   • PubMed
   • Google Scholar
7. 1. Boly M
   2. Perlbarg V
   3. Marrelec G
   4. Schabus M
   5. Laureys S
   6. Doyon J
   7. Pélégrini-Issac M
   8. Maquet P
   9. Benali H

   (2012) Hierarchical clustering of brain activity during human nonrapid eye movement sleep

   PNAS 109:5856–5861.

   https://doi.org/10.1073/pnas.1111133109

   • PubMed
   • Google Scholar
8. 1. Braun AR
   2. Balkin TJ
   3. Wesenten NJ
   4. Carson RE
   5. Varga M
   6. Baldwin P
   7. Selbie S
   8. Belenky G
   9. Herscovitch P

   (1997) Regional cerebral blood flow throughout the sleep-wake cycle an H2(15)O PET study

   Brain 120:1173–1197.

   https://doi.org/10.1093/brain/120.7.1173

   • PubMed
   • Google Scholar
9. 1. Buzsáki G
   2. Anastassiou CA
   3. Koch C

   (2012) The origin of extracellular fields and currents--EEG, ECoG, LFP and spikes

   Nature Reviews Neuroscience 13:407–420.

   https://doi.org/10.1038/nrn3241

   • PubMed
   • Google Scholar
10. 1. Cantero JL
    2. Atienza M
    3. Madsen JR
    4. Stickgold R

    (2004) Gamma EEG dynamics in neocortex and Hippocampus during human wakefulness and sleep

    NeuroImage 22:1271–1280.

    https://doi.org/10.1016/j.neuroimage.2004.03.014

    • PubMed
    • Google Scholar
11. 1. Cardoso MMB
    2. Lima B
    3. Sirotin YB
    4. Das A

    (2019) Task-related hemodynamic responses are modulated by reward and task engagement

    PLOS Biology 17:e3000080.

    https://doi.org/10.1371/journal.pbio.3000080

    • PubMed
    • Google Scholar
12. 1. Chang C
    2. Leopold DA
    3. Schölvinck ML
    4. Mandelkow H
    5. Picchioni D
    6. Liu X
    7. Ye FQ
    8. Turchi JN
    9. Duyn JH

    (2016) Tracking brain arousal fluctuations with fMRI

    PNAS 113:4518–4523.

    https://doi.org/10.1073/pnas.1520613113

    • PubMed
    • Google Scholar
13. 1. Chang C
    2. Glover GH

    (2010) Time-frequency dynamics of resting-state brain connectivity measured with fMRI

    NeuroImage 50:81–98.

    https://doi.org/10.1016/j.neuroimage.2009.12.011

    • PubMed
    • Google Scholar
14. 1. Cirelli C

    (2009) The genetic and molecular regulation of sleep: from fruit flies to humans

    Nature Reviews Neuroscience 10:549–560.

    https://doi.org/10.1038/nrn2683

    • PubMed
    • Google Scholar
15. 1. Dang-Vu TT
    2. Schabus M
    3. Desseilles M
    4. Albouy G
    5. Boly M
    6. Darsaud A
    7. Gais S
    8. Rauchs G
    9. Sterpenich V
    10. Vandewalle G
    11. Carrier J
    12. Moonen G
    13. Balteau E
    14. Degueldre C
    15. Luxen A
    16. Phillips C
    17. Maquet P

    (2008) Spontaneous neural activity during human slow wave sleep

    PNAS 105:15160–15165.

    https://doi.org/10.1073/pnas.0801819105

    • PubMed
    • Google Scholar
16. 1. Datta S
    2. Maclean RR

    (2007) Neurobiological mechanisms for the regulation of mammalian sleep-wake behavior: reinterpretation of historical evidence and inclusion of contemporary cellular and molecular evidence

    Neuroscience & Biobehavioral Reviews 31:775–824.

    https://doi.org/10.1016/j.neubiorev.2007.02.004

    • PubMed
    • Google Scholar
17. 1. de Zwart JA
    2. Silva AC
    3. van Gelderen P
    4. Kellman P
    5. Fukunaga M
    6. Chu R
    7. Koretsky AP
    8. Frank JA
    9. Duyn JH

    (2005) Temporal dynamics of the BOLD fMRI impulse response

    NeuroImage 24:667–677.

    https://doi.org/10.1016/j.neuroimage.2004.09.013

    • Google Scholar
18. 1. Drew PJ
    2. Shih AY
    3. Driscoll JD
    4. Knutsen PM
    5. Blinder P
    6. Davalos D
    7. Akassoglou K
    8. Tsai PS
    9. Kleinfeld D

    (2010a) Chronic optical access through a polished and reinforced thinned skull

    Nature Methods 7:981–984.

    https://doi.org/10.1038/nmeth.1530

    • PubMed
    • Google Scholar
19. 1. Drew PJ
    2. Blinder P
    3. Cauwenberghs G
    4. Shih AY
    5. Kleinfeld D

    (2010b) Rapid determination of particle velocity from space-time images using the radon transform

    Journal of Computational Neuroscience 29:5–11.

    https://doi.org/10.1007/s10827-009-0159-1

    • PubMed
    • Google Scholar
20. 1. Drew PJ
    2. Shih AY
    3. Kleinfeld D

    (2011) Fluctuating and sensory-induced vasodynamics in rodent cortex extend arteriole capacity

    PNAS 108:8473–8478.

    https://doi.org/10.1073/pnas.1100428108

    • PubMed
    • Google Scholar
21. 1. Drew PJ

    (2019) Vascular and neural basis of the BOLD signal

    Current Opinion in Neurobiology 58:61–69.

    https://doi.org/10.1016/j.conb.2019.06.004

    • PubMed
    • Google Scholar
22. 1. Drew PJ
    2. Winder AT
    3. Zhang Q

    (2019) Twitches, blinks, and fidgets: important generators of ongoing neural activity

    The Neuroscientist 25:298–313.

    https://doi.org/10.1177/1073858418805427

    • PubMed
    • Google Scholar
23. 1. Drew PJ
    2. Mateo C
    3. Turner KL
    4. Yu X
    5. Kleinfeld D

    (2020) Ultra-slow oscillations in fMRI and Resting-State connectivity: neuronal and vascular contributions and technical confounds

    Neuron 107:782–804.

    https://doi.org/10.1016/j.neuron.2020.07.020

    • PubMed
    • Google Scholar
24. Software

    1. Drew PJ

    (2020) LabVIEW-DAQ VIs and Hardware, version VI

    GitHub.

    https://github.com/DrewLab/LabVIEW-DAQ
25. 1. Drew PJ
    2. Feldman DE

    (2009) Intrinsic signal imaging of deprivation-induced contraction of whisker representations in rat somatosensory cortex

    Cerebral Cortex 19:331–348.

    https://doi.org/10.1093/cercor/bhn085

    • PubMed
    • Google Scholar
26. 1. Echagarruga CT
    2. Gheres KW
    3. Norwood JN
    4. Drew PJ

    (2020) nNOS-expressing interneurons control basal and behaviorally evoked arterial dilation in somatosensory cortex of mice

    eLife 9:e60533.

    https://doi.org/10.7554/eLife.60533

    • PubMed
    • Google Scholar
27. 1. Flynn NM
    2. Buljubasic N
    3. Bosnjak ZJ
    4. Kampine JP

    (1992) Isoflurane produces endothelium-independent relaxation in canine middle cerebral arteries

    Anesthesiology 76:461–467.

    https://doi.org/10.1097/00000542-199203000-00021

    • PubMed
    • Google Scholar
28. Book

    1. Franklin KBJ
    2. Paxinos G

    (2007)

    The Mouse Brain in Stereotaxic Coordinates

    Academic Press.

    • Google Scholar
29. 1. Fukunaga M
    2. Horovitz SG
    3. van Gelderen P
    4. de Zwart JA
    5. Jansma JM
    6. Ikonomidou VN
    7. Chu R
    8. Deckers RH
    9. Leopold DA
    10. Duyn JH

    (2006) Large-amplitude, spatially correlated fluctuations in BOLD fMRI signals during extended rest and early sleep stages

    Magnetic Resonance Imaging 24:979–992.

    https://doi.org/10.1016/j.mri.2006.04.018

    • PubMed
    • Google Scholar
30. 1. Fultz NE
    2. Bonmassar G
    3. Setsompop K
    4. Stickgold RA
    5. Rosen BR
    6. Polimeni JR
    7. Lewis LD

    (2019) Coupled electrophysiological, Hemodynamic, and cerebrospinal fluid oscillations in human sleep

    Science 366:628–631.

    https://doi.org/10.1126/science.aax5440

    • PubMed
    • Google Scholar
31. 1. Funk CM
    2. Honjoh S
    3. Rodriguez AV
    4. Cirelli C
    5. Tononi G

    (2016) Local slow waves in superficial layers of primary cortical Areas during REM sleep

    Current Biology 26:396–403.

    https://doi.org/10.1016/j.cub.2015.11.062

    • PubMed
    • Google Scholar
32. 1. Gao Y-R
    2. Greene SE
    3. Drew PJ

    (2015) Mechanical restriction of intracortical vessel dilation by brain tissue sculpts the hemodynamic response

    NeuroImage 115:162–176.

    https://doi.org/10.1016/j.neuroimage.2015.04.054

    • Google Scholar
33. 1. Gao YR
    2. Drew PJ

    (2014) Determination of vessel cross-sectional area by thresholding in radon space

    Journal of Cerebral Blood Flow & Metabolism 34:1180–1187.

    https://doi.org/10.1038/jcbfm.2014.67

    • PubMed
    • Google Scholar
34. 1. Gao YR
    2. Drew PJ

    (2016) Effects of voluntary locomotion and calcitonin Gene-Related peptide on the dynamics of single dural vessels in awake mice

    The Journal of Neuroscience 36:2503–2516.

    https://doi.org/10.1523/JNEUROSCI.3665-15.2016

    • PubMed
    • Google Scholar
35. 1. Gu Y
    2. Han F
    3. Liu X

    (2019) Arousal contributions to Resting-State fMRI connectivity and dynamics

    Frontiers in Neuroscience 13:1–8.

    https://doi.org/10.3389/fnins.2019.01190

    • Google Scholar
36. 1. Haiss F
    2. Jolivet R
    3. Wyss MT
    4. Reichold J
    5. Braham NB
    6. Scheffold F
    7. Krafft MP
    8. Weber B

    (2009) Improved in vivo two-photon imaging after blood replacement by perfluorocarbon

    The Journal of Physiology 587:3153–3158.

    https://doi.org/10.1113/jphysiol.2009.169474

    • PubMed
    • Google Scholar
37. 1. Harris KD
    2. Quiroga RQ
    3. Freeman J
    4. Smith SL

    (2016) Improving data quality in neuronal population recordings

    Nature Neuroscience 19:1165–1174.

    https://doi.org/10.1038/nn.4365

    • PubMed
    • Google Scholar
38. 1. Hobson JA
    2. Pace-Schott EF

    (2002) The cognitive neuroscience of sleep: neuronal systems, consciousness and learning

    Nature Reviews Neuroscience 3:679–693.

    https://doi.org/10.1038/nrn915

    • PubMed
    • Google Scholar
39. 1. Horovitz SG
    2. Fukunaga M
    3. de Zwart JA
    4. van Gelderen P
    5. Fulton SC
    6. Balkin TJ
    7. Duyn JH

    (2008) Low frequency BOLD fluctuations during resting wakefulness and light sleep: a simultaneous EEG-fMRI study

    Human Brain Mapping 29:671–682.

    https://doi.org/10.1002/hbm.20428

    • PubMed
    • Google Scholar
40. 1. Huo BX
    2. Smith JB
    3. Drew PJ

    (2014) Neurovascular coupling and decoupling in the cortex during voluntary locomotion

    Journal of Neuroscience 34:10975–10981.

    https://doi.org/10.1523/JNEUROSCI.1369-14.2014

    • PubMed
    • Google Scholar
41. 1. Huo BX
    2. Gao YR
    3. Drew PJ

    (2015a) Quantitative separation of arterial and venous cerebral blood volume increases during voluntary locomotion

    NeuroImage 105:369–379.

    https://doi.org/10.1016/j.neuroimage.2014.10.030

    • PubMed
    • Google Scholar
42. 1. Huo BX
    2. Greene SE
    3. Drew PJ

    (2015b) Venous cerebral blood volume increase during voluntary locomotion reflects cardiovascular changes

    NeuroImage 118:301–312.

    https://doi.org/10.1016/j.neuroimage.2015.06.011

    • PubMed
    • Google Scholar
43. 1. Iadecola C

    (2017) The neurovascular unit coming of age: a journey through neurovascular coupling in health and disease

    Neuron 96:17–42.

    https://doi.org/10.1016/j.neuron.2017.07.030

    • PubMed
    • Google Scholar
44. 1. Kajikawa Y
    2. Schroeder CE

    (2011) How local is the local field potential?

    Neuron 72:847–858.

    https://doi.org/10.1016/j.neuron.2011.09.029

    • PubMed
    • Google Scholar
45. 1. Kedarasetti RT
    2. Turner KL
    3. Echagarruga C
    4. Gluckman BJ
    5. Drew PJ
    6. Costanzo F

    (2020) Functional hyperemia drives fluid exchange in the paravascular space

    Fluids and Barriers of the CNS 17:52.

    https://doi.org/10.1186/s12987-020-00214-3

    • PubMed
    • Google Scholar
46. 1. Kim SG
    2. Ogawa S

    (2012) Biophysical and physiological origins of blood oxygenation level-dependent fMRI signals

    Journal of Cerebral Blood Flow & Metabolism 32:1188–1206.

    https://doi.org/10.1038/jcbfm.2012.23

    • PubMed
    • Google Scholar
47. 1. Kleinfeld D
    2. Deschênes M

    (2011) Neuronal basis for object location in the vibrissa scanning sensorimotor system

    Neuron 72:455–468.

    https://doi.org/10.1016/j.neuron.2011.10.009

    • PubMed
    • Google Scholar
48. 1. Krawchuk MB
    2. Ruff CF
    3. Yang X
    4. Ross SE
    5. Vazquez AL

    (2020) Optogenetic assessment of VIP, PV, SOM and NOS inhibitory neuron activity and cerebral blood flow regulation in mouse somato-sensory cortex

    Journal of Cerebral Blood Flow & Metabolism 40:1427–1440.

    https://doi.org/10.1177/0271678X19870105

    • PubMed
    • Google Scholar
49. 1. Larson-Prior LJ
    2. Zempel JM
    3. Nolan TS
    4. Prior FW
    5. Snyder AZ
    6. Raichle ME

    (2009) Cortical network functional connectivity in the descent to sleep

    PNAS 106:4489–4494.

    https://doi.org/10.1073/pnas.0900924106

    • PubMed
    • Google Scholar
50. 1. Le Van Quyen M
    2. Staba R
    3. Bragin A
    4. Dickson C
    5. Valderrama M
    6. Fried I
    7. Engel J

    (2010) Large-scale microelectrode recordings of high-frequency gamma oscillations in human cortex during sleep

    Journal of Neuroscience 30:7770–7782.

    https://doi.org/10.1523/JNEUROSCI.5049-09.2010

    • PubMed
    • Google Scholar
51. 1. Lee L
    2. Boorman L
    3. Glendenning E
    4. Christmas C
    5. Sharp P
    6. Redgrave P
    7. Shabir O
    8. Bracci E
    9. Berwick J
    10. Howarth C

    (2020) Key aspects of neurovascular control mediated by specific populations of inhibitory cortical interneurons

    Cerebral Cortex 30:2452–2464.

    https://doi.org/10.1093/cercor/bhz251

    • PubMed
    • Google Scholar
52. 1. Leopold DA
    2. Murayama Y
    3. Logothetis NK

    (2003) Very slow activity fluctuations in monkey visual cortex: implications for functional brain imaging

    Cerebral Cortex 13:422–433.

    https://doi.org/10.1093/cercor/13.4.422

    • PubMed
    • Google Scholar
53. 1. Liu TT

    (2016) Noise contributions to the fMRI signal: an overview

    NeuroImage 143:141–151.

    https://doi.org/10.1016/j.neuroimage.2016.09.008

    • PubMed
    • Google Scholar
54. 1. Liu X
    2. de Zwart JA
    3. Schölvinck ML
    4. Chang C
    5. Ye FQ
    6. Leopold DA
    7. Duyn JH

    (2018) Subcortical evidence for a contribution of arousal to fMRI studies of brain activity

    Nature Communications 9:1–10.

    https://doi.org/10.1038/s41467-017-02815-3

    • Google Scholar
55. 1. Ma Y
    2. Shaik MA
    3. Kozberg MG
    4. Kim SH
    5. Portes JP
    6. Timerman D
    7. Hillman EMC

    (2016) Resting-state hemodynamics are spatiotemporally coupled to synchronized and symmetric neural activity in excitatory neurons

    PNAS 113:E8463–E8471.

    https://doi.org/10.1073/pnas.1525369113

    • Google Scholar
56. 1. Madsen PL
    2. Holm S
    3. Vorstrup S
    4. Friberg L
    5. Lassen NA
    6. Wildschiødtz G

    (1991) Human regional cerebral blood flow during rapid-eye-movement sleep

    Journal of Cerebral Blood Flow & Metabolism 11:502–507.

    https://doi.org/10.1038/jcbfm.1991.94

    • PubMed
    • Google Scholar
57. 1. Maquet P
    2. Phillips C

    (1998) Functional brain imaging of human sleep

    Journal of Sleep Research 7:42–47.

    https://doi.org/10.1046/j.1365-2869.7.s1.7.x

    • PubMed
    • Google Scholar
58. 1. Mateo C
    2. Knutsen PM
    3. Tsai PS
    4. Shih AY
    5. Kleinfeld D

    (2017) Entrainment of arteriole vasomotor fluctuations by neural activity is a basis of Blood-Oxygenation-Level-Dependent "Resting-State" Connectivity

    Neuron 96:936–948.

    https://doi.org/10.1016/j.neuron.2017.10.012

    • PubMed
    • Google Scholar
59. 1. McMahon SM
    2. Jackson MB

    (2018) An inconvenient truth: calcium sensors are calcium buffers

    Trends in Neurosciences 41:880–884.

    https://doi.org/10.1016/j.tins.2018.09.005

    • PubMed
    • Google Scholar
60. 1. Mitra A
    2. Snyder AZ
    3. Tagliazucchi E
    4. Laufs H
    5. Raichle ME

    (2015) Propagated infra-slow intrinsic brain activity reorganizes across wake and slow wave sleep

    eLife 4:e10781.

    https://doi.org/10.7554/eLife.10781

    • PubMed
    • Google Scholar
61. 1. Mitra PP
    2. Pesaran B

    (1999) Analysis of dynamic brain imaging data

    Biophysical Journal 76:691–708.

    https://doi.org/10.1016/S0006-3495(99)77236-X

    • PubMed
    • Google Scholar
62. 1. Montgomery SM
    2. Sirota A
    3. Buzsáki G

    (2008) Theta and gamma coordination of hippocampal networks during waking and rapid eye movement sleep

    Journal of Neuroscience 28:6731–6741.

    https://doi.org/10.1523/JNEUROSCI.1227-08.2008

    • PubMed
    • Google Scholar
63. 1. Musall S
    2. Kaufman MT
    3. Juavinett AL
    4. Gluf S
    5. Churchland AK

    (2019) Single-trial neural dynamics are dominated by richly varied movements

    Nature Neuroscience 22:1677–1686.

    https://doi.org/10.1038/s41593-019-0502-4

    • PubMed
    • Google Scholar
64. 1. Nicholson C
    2. Freeman JA

    (1975) Theory of current source-density analysis and determination of conductivity tensor for anuran cerebellum

    Journal of Neurophysiology 38:356–368.

    https://doi.org/10.1152/jn.1975.38.2.356

    • PubMed
    • Google Scholar
65. 1. O'Connor DH
    2. Clack NG
    3. Huber D
    4. Komiyama T
    5. Myers EW
    6. Svoboda K

    (2010) Vibrissa-based object localization in head-fixed mice

    Journal of Neuroscience 30:1947–1967.

    https://doi.org/10.1523/JNEUROSCI.3762-09.2010

    • PubMed
    • Google Scholar
66. 1. Özbay PS
    2. Chang C
    3. Picchioni D
    4. Mandelkow H
    5. Moehlman TM
    6. Chappel-Farley MG
    7. van Gelderen P
    8. de Zwart JA
    9. Duyn JH

    (2018) Contribution of systemic vascular effects to fMRI activity in white matter

    NeuroImage 176:541–549.

    https://doi.org/10.1016/j.neuroimage.2018.04.045

    • PubMed
    • Google Scholar
67. 1. Pace-Schott EF
    2. Hobson JA

    (2002) The neurobiology of sleep: genetics, cellular physiology and subcortical networks

    Nature Reviews Neuroscience 3:591–605.

    https://doi.org/10.1038/nrn895

    • PubMed
    • Google Scholar
68. 1. Peirson SN
    2. Brown LA
    3. Pothecary CA
    4. Benson LA
    5. Fisk AS

    (2018) Light and the laboratory mouse

    Journal of Neuroscience Methods 300:26–36.

    https://doi.org/10.1016/j.jneumeth.2017.04.007

    • PubMed
    • Google Scholar
69. 1. Petersen CC

    (2007) The functional organization of the barrel cortex

    Neuron 56:339–355.

    https://doi.org/10.1016/j.neuron.2007.09.017

    • PubMed
    • Google Scholar
70. 1. Pilorz V
    2. Tam SK
    3. Hughes S
    4. Pothecary CA
    5. Jagannath A
    6. Hankins MW
    7. Bannerman DM
    8. Lightman SL
    9. Vyazovskiy VV
    10. Nolan PM
    11. Foster RG
    12. Peirson SN

    (2016) Melanopsin regulates both Sleep-Promoting and Arousal-Promoting responses to light

    PLOS Biology 14:e1002482.

    https://doi.org/10.1371/journal.pbio.1002482

    • PubMed
    • Google Scholar
71. 1. Rungta RL
    2. Chaigneau E
    3. Osmanski BF
    4. Charpak S

    (2018) Vascular compartmentalization of functional hyperemia from the synapse to the pia

    Neuron 99:362–375.

    https://doi.org/10.1016/j.neuron.2018.06.012

    • PubMed
    • Google Scholar
72. 1. Sakai K

    (2020) What single‐unit recording studies tell Us about the basic mechanisms of sleep and wakefulness

    European Journal of Neuroscience 52:3507–3530.

    https://doi.org/10.1111/ejn.14485

    • Google Scholar
73. 1. Saper CB
    2. Fuller PM
    3. Pedersen NP
    4. Lu J
    5. Scammell TE

    (2010) Sleep state switching

    Neuron 68:1023–1042.

    https://doi.org/10.1016/j.neuron.2010.11.032

    • PubMed
    • Google Scholar
74. 1. Saper CB
    2. Fuller PM

    (2017) Wake-sleep circuitry: an overview

    Current Opinion in Neurobiology 44:186–192.

    https://doi.org/10.1016/j.conb.2017.03.021

    • PubMed
    • Google Scholar
75. 1. Scammell TE
    2. Arrigoni E
    3. Lipton JO

    (2017) Neural circuitry of wakefulness and sleep

    Neuron 93:747–765.

    https://doi.org/10.1016/j.neuron.2017.01.014

    • PubMed
    • Google Scholar
76. 1. Schölvinck ML
    2. Maier A
    3. Ye FQ
    4. Duyn JH
    5. Leopold DA

    (2010) Neural basis of global resting-state fMRI activity

    PNAS 107:10238–10243.

    https://doi.org/10.1073/pnas.0913110107

    • PubMed
    • Google Scholar
77. 1. Shen Z
    2. Lu Z
    3. Chhatbar PY
    4. O'Herron P
    5. Kara P

    (2012) An artery-specific fluorescent dye for studying neurovascular coupling

    Nature Methods 9:273–276.

    https://doi.org/10.1038/nmeth.1857

    • PubMed
    • Google Scholar
78. 1. Shih AY
    2. Driscoll JD
    3. Drew PJ
    4. Nishimura N
    5. Schaffer CB
    6. Kleinfeld D

    (2012a) Two-photon microscopy as a tool to study blood flow and neurovascular coupling in the rodent brain

    Journal of Cerebral Blood Flow & Metabolism 32:1277–1309.

    https://doi.org/10.1038/jcbfm.2011.196

    • PubMed
    • Google Scholar
79. 1. Shih AY
    2. Mateo C
    3. Drew PJ
    4. Tsai PS
    5. Kleinfeld D

    (2012b) A polished and reinforced Thinned-skull window for Long-term imaging of the mouse brain

    Journal of Visualized Experiments 61:1–6.

    https://doi.org/10.3791/3742

    • Google Scholar
80. 1. Shirey MJ
    2. Smith JB
    3. Kudlik DE
    4. Huo BX
    5. Greene SE
    6. Drew PJ

    (2015) Brief anesthesia, but not voluntary locomotion, significantly alters cortical temperature

    Journal of Neurophysiology 114:309–322.

    https://doi.org/10.1152/jn.00046.2015

    • PubMed
    • Google Scholar
81. 1. Shmuel A
    2. Leopold DA

    (2008) Neuronal correlates of spontaneous fluctuations in fMRI signals in monkey visual cortex: implications for functional connectivity at rest

    Human Brain Mapping 29:751–761.

    https://doi.org/10.1002/hbm.20580

    • PubMed
    • Google Scholar
82. 1. Silva AC
    2. Koretsky AP
    3. Duyn JH

    (2007) Functional MRI impulse response for BOLD and CBV contrast in rat somatosensory cortex

    Magnetic Resonance in Medicine 57:1110–1118.

    https://doi.org/10.1002/mrm.21246

    • PubMed
    • Google Scholar
83. 1. Simon MJ
    2. Iliff JJ

    (2016) Regulation of cerebrospinal fluid (CSF) flow in neurodegenerative, neurovascular and neuroinflammatory disease

    Biochimica Et Biophysica Acta (BBA) - Molecular Basis of Disease 1862:442–451.

    https://doi.org/10.1016/j.bbadis.2015.10.014

    • Google Scholar
84. 1. Sirotin YB
    2. Das A

    (2009) Anticipatory haemodynamic signals in sensory cortex not predicted by local neuronal activity

    Nature 457:475–479.

    https://doi.org/10.1038/nature07664

    • PubMed
    • Google Scholar
85. 1. Steinmetz NA
    2. Buetfering C
    3. Lecoq J
    4. Lee CR
    5. Peters AJ
    6. Jacobs EAK
    7. Coen P
    8. Ollerenshaw DR
    9. Valley MT
    10. de Vries SEJ
    11. Garrett M
    12. Zhuang J
    13. Groblewski PA
    14. Manavi S
    15. Miles J
    16. White C
    17. Lee E
    18. Griffin F
    19. Larkin JD
    20. Roll K
    21. Cross S
    22. Nguyen TV
    23. Larsen R
    24. Pendergraft J
    25. Daigle T
    26. Tasic B
    27. Thompson CL
    28. Waters J
    29. Olsen S
    30. Margolis DJ
    31. Zeng H
    32. Hausser M
    33. Carandini M
    34. Harris KD

    (2017) Aberrant cortical activity in multiple GCaMP6-Expressing transgenic mouse lines

    Eneuro 4:ENEURO.0207-17.2017.

    https://doi.org/10.1523/ENEURO.0207-17.2017

    • PubMed
    • Google Scholar
86. 1. Stringer C
    2. Pachitariu M
    3. Steinmetz N
    4. Reddy CB
    5. Carandini M
    6. Harris KD

    (2019) Spontaneous behaviors drive multidimensional, brainwide activity

    Science 364:eaav7893.

    https://doi.org/10.1126/science.aav7893

    • PubMed
    • Google Scholar
87. 1. Sullivan D
    2. Mizuseki K
    3. Sorgi A
    4. Buzsáki G

    (2014) Comparison of sleep spindles and theta oscillations in the Hippocampus

    Journal of Neuroscience 34:662–674.

    https://doi.org/10.1523/JNEUROSCI.0552-13.2014

    • PubMed
    • Google Scholar
88. 1. Tagliazucchi E
    2. Laufs H

    (2014) Decoding wakefulness levels from typical fMRI resting-state data reveals reliable drifts between wakefulness and sleep

    Neuron 82:695–708.

    https://doi.org/10.1016/j.neuron.2014.03.020

    • PubMed
    • Google Scholar
89. 1. Tarasoff-Conway JM
    2. Carare RO
    3. Osorio RS
    4. Glodzik L
    5. Butler T
    6. Fieremans E
    7. Axel L
    8. Rusinek H
    9. Nicholson C
    10. Zlokovic BV
    11. Frangione B
    12. Blennow K
    13. Ménard J
    14. Zetterberg H
    15. Wisniewski T
    16. de Leon MJ

    (2015) Clearance systems in the brain-implications for alzheimer disease

    Nature Reviews Neurology 11:457–470.

    https://doi.org/10.1038/nrneurol.2015.119

    • PubMed
    • Google Scholar
90. 1. Theis L
    2. Berens P
    3. Froudarakis E
    4. Reimer J
    5. Román Rosón M
    6. Baden T
    7. Euler T
    8. Tolias AS
    9. Bethge M

    (2016) Benchmarking spike rate inference in population calcium imaging

    Neuron 90:471–482.

    https://doi.org/10.1016/j.neuron.2016.04.014

    • PubMed
    • Google Scholar
91. 1. Townsend RE
    2. Prinz PN
    3. Obrist WD

    (1973) Human cerebral blood flow during sleep and waking

    Journal of Applied Physiology 35:620–625.

    https://doi.org/10.1152/jappl.1973.35.5.620

    • PubMed
    • Google Scholar
92. Software

    1. Turner KL
    2. Gheres KW
    3. Proctor EA
    4. Drew PJ

    (2020) Neurovascular coupling and bilateral connectivity during NREM and REM sleep, version 2.6

    GitHub.

    https://github.com/KL-Turner/Turner_Gheres_Proctor_Drew_eLife2020
93. 1. van Veluw SJ
    2. Hou SS
    3. Calvo-Rodriguez M
    4. Arbel-Ornath M
    5. Snyder AC
    6. Frosch MP
    7. Greenberg SM
    8. Bacskai BJ

    (2020) Vasomotion as a driving force for paravascular clearance in the awake mouse brain

    Neuron 105:549–561.

    https://doi.org/10.1016/j.neuron.2019.10.033

    • PubMed
    • Google Scholar
94. 1. Vazquez AL
    2. Fukuda M
    3. Crowley JC
    4. Kim SG

    (2014) Neural and hemodynamic responses elicited by forelimb- and photo-stimulation in channelrhodopsin-2 mice: insights into the hemodynamic point spread function

    Cerebral Cortex 24:2908–2919.

    https://doi.org/10.1093/cercor/bht147

    • PubMed
    • Google Scholar
95. 1. Veasey SC
    2. Valladares O
    3. Fenik P
    4. Kapfhamer D
    5. Sanford L
    6. Benington J
    7. Bucan M

    (2000) An automated system for recording and analysis of sleep in mice

    Sleep 23:1–16.

    https://doi.org/10.1093/sleep/23.8.1c

    • Google Scholar
96. 1. Vyazovskiy VV
    2. Olcese U
    3. Hanlon EC
    4. Nir Y
    5. Cirelli C
    6. Tononi G

    (2011) Local sleep in awake rats

    Nature 472:443–447.

    https://doi.org/10.1038/nature10009

    • Google Scholar
97. 1. Vyazovskiy VV
    2. Harris KD

    (2013) Sleep and the single neuron: the role of global slow oscillations in individual cell rest

    Nature Reviews Neuroscience 14:443–451.

    https://doi.org/10.1038/nrn3494

    • PubMed
    • Google Scholar
98. 1. Walker MP
    2. Stickgold R

    (2004) Sleep-dependent learning and memory consolidation

    Neuron 44:121–133.

    https://doi.org/10.1016/j.neuron.2004.08.031

    • PubMed
    • Google Scholar
99. 1. Weber F
    2. Dan Y

    (2016) Circuit-based interrogation of sleep control

    Nature 538:51–59.

    https://doi.org/10.1038/nature19773

    • PubMed
    • Google Scholar
100. 1. Winder AT
     2. Echagarruga C
     3. Zhang Q
     4. Drew PJ

     (2017) Weak correlations between hemodynamic signals and ongoing neural activity during the resting state

     Nature Neuroscience 20:1761–1769.

     https://doi.org/10.1038/s41593-017-0007-y

     • PubMed
     • Google Scholar
101. 1. Xie L
     2. Kang H
     3. Xu Q
     4. Chen MJ
     5. Liao Y
     6. Thiyagarajan M
     7. O'Donnell J
     8. Christensen DJ
     9. Nicholson C
     10. Iliff JJ
     11. Takano T
     12. Deane R
     13. Nedergaard M

     (2013) Sleep drives metabolite clearance from the adult brain

     Science 342:373–377.

     https://doi.org/10.1126/science.1241224

     • PubMed
     • Google Scholar
102. 1. Yüzgeç Ö
     2. Prsa M
     3. Zimmermann R
     4. Huber D

     (2018) Pupil size coupling to cortical states protects the stability of deep sleep via parasympathetic modulation

     Current Biology : CB 28:392–400.

     https://doi.org/10.1016/j.cub.2017.12.049

     • PubMed
     • Google Scholar
103. 1. Zhang Q
     2. Roche M
     3. Gheres KW
     4. Chaigneau E
     5. Kedarasetti RT
     6. Haselden WD
     7. Charpak S
     8. Drew PJ

     (2019) Cerebral oxygenation during locomotion is modulated by respiration

     Nature Communications 10:5515.

     https://doi.org/10.1038/s41467-019-13523-5

     • PubMed
     • Google Scholar

Author details

1. Kevin L Turner

   1. Department of Biomedical Engineering, The Pennsylvania State University, University Park, United States
   2. Center for Neural Engineering, The Pennsylvania State University, University Park, United States

   Contribution

   Conceptualization, Resources, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-3044-7079

2. Kyle W Gheres

   1. Center for Neural Engineering, The Pennsylvania State University, University Park, United States
   2. Graduate Program in Molecular, Cellular, and Integrative Biosciences, The Pennsylvania State University, University Park, United States

   Contribution

   Resources, Software, Methodology, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-7568-9023

3. Elizabeth A Proctor

   1. Department of Biomedical Engineering, The Pennsylvania State University, University Park, United States
   2. Center for Neural Engineering, The Pennsylvania State University, University Park, United States
   3. Department of Engineering Science and Mechanics, The Pennsylvania State University, University Park, United States
   4. Department of Neurosurgery, Penn State College of Medicine, Hershey, United States
   5. Department of Pharmacology, Penn State College of Medicine, Hershey, United States

   Contribution

   Formal analysis, Supervision, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7627-2198

4. Patrick J Drew

   1. Department of Biomedical Engineering, The Pennsylvania State University, University Park, United States
   2. Center for Neural Engineering, The Pennsylvania State University, University Park, United States
   3. Department of Engineering Science and Mechanics, The Pennsylvania State University, University Park, United States
   4. Department of Neurosurgery, Penn State College of Medicine, Hershey, United States

   Contribution

   Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Visualization, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   PJD17@PSU.EDU

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7483-7378

• 5,805

  views

• 793

  downloads

• 77

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.