Functional differences in cerebral activation between slow wave-coupled and uncoupled sleep spindles

Spindles are often temporally coupled to slow waves (SW). These SW-spindle complexes have been implicated in memory consolidation that involves transfer of information from the hippocampus to the neocortex. However, spindles and SW, which are characteristic of NREM sleep, can occur as part of this complex, or in isolation. It is not clear whether dissociable parts of the brain are recruited when coupled to SW vs. when spindles or SW occur in isolation. Here, we tested differences in cerebral activation time-locked to uncoupled spindles, uncoupled SW and coupled SW-spindle complexes using simultaneous EEG-fMRI. Consistent with the “active system model,” we hypothesized that brain activations time-locked to coupled SW-spindles would preferentially occur in brain areas known to be critical for sleep-dependent memory consolidation. Our results show that coupled spindles and uncoupled spindles recruit distinct parts of the brain. Specifically, we found that hippocampal activation during sleep is not uniquely related to spindles. Rather, this process is primarily driven by SWs and SW-spindle coupling. In addition, we show that SW-spindle coupling is critical in the activation of the putamen. Importantly, SW-spindle coupling specifically recruited frontal areas in comparison to uncoupled spindles, which may be critical for the hippocampal-neocortical dialogue that preferentially occurs during sleep.

Introduction

There is a large body of evidence indicating that the optimal transformation of newly acquired, labile memory traces into lasting and integrated memories is facilitated by sleep (Diekelmann and Born, 2010; Lewis and Durrant, 2011; Dudai et al., 2015; Witkowski et al., 2020). The process of memory consolidation involves transfer of information from the hippocampus to the neocortex (Wang and Morris, 2009; Eichenbaum, 2017; Takehara-Nishiuchi, 2020), and during sleep, this process is facilitated by hippocampal ripples. Ripples are involved in the repeated reactivation and replaying of newly-formed hippocampal-dependent memory traces (Buzsáki, 2015; Sadowski et al., 2016). Ripples, however, do not occur in isolation. Rather, they are part of slow wave–spindle–ripple complexes, whereby ripples are nested in the excitatory troughs of spindles, and spindles are nested preferentially, but not exclusively (Yordanova et al., 2017), in the excitatory troughs (“up-states”) of slow waves (Helfrich et al., 2019). The “active system model” of memory consolidation suggests that these three oscillations form a temporal hierarchy, which serves as an endogenous timing mechanism that supports the hippocampal-neocortical dialogue. This dialogue is required for the consolidation of new memories (Latchoumane et al., 2017; Klinzing et al., 2019). While it is not possible to record hippocampal ripples in humans without employing invasive techniques, slow wave-spindle (SW-SP) coupling can be measured from scalp recordings and is considered to be an index of memory reactivation during sleep (Schreiner et al., 2021).

Advances in simultaneous EEG-fMRI recording have enabled the measurement of regional brain activity and functional connectivity, time-locked precisely to spontaneous events during sleep, such as slow waves (Kaufmann et al., 2006; Dang-Vu et al., 2008; Betta et al., 2021) and spindles (Laufs et al., 2007; Schabus et al., 2007; Tyvaert et al., 2008; Andrade et al., 2011; Caporro et al., 2012; Fang et al., 2019, 2020), hence providing insights into systems-level changes in brain activity that reflect memory consolidation during sleep. These studies have shown that spindles do not only recruit areas necessary for spindle generation (e.g., thalamus), but interestingly, that they spontaneously recruit and reactivate brain areas known to support memory functions (e.g., putamen, hippocampus). At the cellular level, spindles are associated with large influxes of Ca²⁺, which facilitates neocortical long-term potentiation via synaptic plasticity (Sejnowski and Destexhe, 2000; Timofeev et al., 2002; Nevian and Sakmann, 2006). This effect is enhanced when spindles occur during SW up states (Niethard et al., 2018). Thus, the spaced and repeated action of SW-SP complexes appear to be an ideal mechanism to facilitate hippocampal-neocortical systems-level consolidation.

However, it is not clear whether these spontaneous activations occur during all spindle events, or specifically during coupled SW-SP complexes. The present study thus aimed to address this unanswered question. We also investigated whether coupled spindles and uncoupled (i.e., isolated) spindles are dissociable from one another in terms of their functional neuroanatomical substrates. Here, we expected to observe dissociable patterns of recruited brain activity in areas time-locked to SW-SP coupled events, as compared to uncoupled spindle events, or compared to uncoupled slow wave events. We hypothesize that brain activations time-locked to coupled SW-spindles would preferentially occur in brain areas known to be critical for the hippocampal-neocortical dialogue and the striatum, required for the consolidation of new memories.

Materials and methods

Participants

All participants were 20–35 years of age. Participants were initially screened for irregular sleep schedules (bedtime outside the hours of ∼22:00–24:00 h, wake time outside the hours of 07:00–09:00 h), left-handedness, shift work, and the use of medications known to affect sleep. Participants were also excluded from the study if they considered themselves a smoker, consumed >1–2 caffeinated beverages/day, consumed >7 alcoholic beverages/week, or had a history of chronic pain, seizures or head injury. Participants were required to refrain from recreational drug use and limit caffeine (to no more than 1 beverage in the A.M.) and alcohol intake at least 3 days prior to, and throughout participation in the study. Actigraphy and sleep logs were used to confirm the participants’ sleep and activity cycles throughout the study. To rule out participants with signs of depression or anxiety and ensure normal sleep-wake patterns, participants also completed the Beck Depression (Beck et al., 1988b) and Anxiety Inventories (Beck et al., 1988a) as well as the Sleep Disorders Questionnaire (Douglass et al., 1994).

A total of 35 participants met study inclusion criteria. A minimum of 30 coupled and 30 uncoupled spindle events were necessary for analysis purposes. An additional fourteen participants did not meet this criterion. In total, N = 21 participants (mean age 24, 13 females) were included in the final data analyses.

Ethics statement

All study procedures and methods adhered to the Declaration of Helsinki and were approved by the Western University Health Science research ethics board. All participants were given a letter with details of the study, provided informed consent, and were financially compensated for their participation.

Procedures

All participants who met the study inclusion criteria underwent an orientation session where they were given the study instructions, a sleep diary and an activity monitor to verify their sleep–wake cycle (Figure 1). A minimum of 1 week following the orientation session, participants completed the EEG-fMRI sleep recording night. Participants arrived at the sleep laboratory on the recording night at approximately 8 P.M. The scanning procedures began at 9 P.M., at which point, the EEG equipment was installed and configured. Following this, localizer scans, a T1 structural scan, and an 8-min eyes-closed awake resting-state scan were completed. EEG was acquired during the initial scanning procedures to confirm that participants remained awake. These procedures took up to one hour to complete. The sleep session (“lights out”) began at about 10 P.M., within the range of the participants’ habitual bedtime. After the sleep session, participants slept the rest of the night in the nearby sleep laboratory.

FIGURE 1

Figure 1. Study design. Participants underwent an initial screening to rule out any signs of sleep disorders, unusual sleep habits, or other health-related criteria and MRI compatibility. Eligible participants then visited the sleep laboratory for the orientation session at least 1 week before the EEG-fMRI sleep recording night, in which participants were given detailed instructions about the study procedure, the sleep diary, and an activity monitor. Finally, participants completed the EEG-fMRI sleep recording session beginning at 9 P.M, with lights out for the sleep session at 10 P.M. The sleep session ended at 12 A.M. Reproduced with permission from cerebral cortex.

Polysomnographic acquisition and analysis

Polysomnographic recording parameters

Polysomnographic (PSG) recordings were obtained using a 64-channel magnetic resonance (MR)-compatible EEG cap, which included one electrocardiogram (ECG) lead (Braincap MR, Easycap) and 64-channels of EEG recorded via two MR-compatible 32-channel amplifiers (Brainamp MR Plus, Brain Products GmbH). EEG recordings were referenced to FCz and digitized at 5,000 samples per second with a 500-nV resolution. The single ECG lead included in the EEG cap offers limited visualization of the r-peak of the QRS complex, needed for accurate ballistocardiograph (BCG) correction. Thus, three additional bipolar ECG recordings were taken using a MR-compatible 15-channel amplifier (Brainamp ExG MR, Brain Products GmbH). To reduce BCG artifact by up to 40% (Mullinger et al., 2011), participants were positioned in the MRI scanner so that they were shifted away from the isocenter of the magnetic field by 40 mm, making BCG correction more straightforward.

Data were analog filtered using a 500 Hz band-limiter low-pass filter and a 0.0159 Hz high-pass filter with a 10-sec time constant. Data were recorded with Brain Products Recorder Software, Version 1.x and transferred to the recording computer via fiber-optic cable and hardware synchronized to the scanner clock using the Brain Products “Sync Box” (Brain Products GmbH). As recommended in the literature (Mulert and Lemieux, 2010), the MRI sequence parameters were selected to ensure that the gradient artifact would be time stable, and that the lowest harmonic of the gradient artifact (18.52 Hz) would occur at the highest possible frequency and above the spindle band (11–16 Hz). Thus, the MR scan repetition time was set to 2,160 msec, matching a common multiple of the EEG sample time (0.2 msec), the product of the scanner clock precision (0.1 μsec), and the number of slices (40) used.

EEG preprocessing

EEG scanner artifacts were removed in several steps. First, an adaptative average template subtraction method (Allen et al., 2000) implemented in Brain Products Analyzer Software, Version 2.x was used and data was downsampled to 250 Hz. Next, r-peaks in the ECG were detected semi-automatically. Each r-peak was visually verified and, when necessary, manually adjusted to correct both false-positive and false-negative r-peak detections. This step is crucial for optimal BCG correction. Next, adaptive template subtraction (Allen et al., 1998) was used to remove BCG artifacts time-locked to the r-peak of the QRS complex. Once the MRI-related artifacts had been corrected, the data was visually inspected and the amplitude of the residual artifacts time-locked to the r-peaks were examined. If the peak of the maximum amplitude of the residual artifact exceeded 3 μV during the QRS complex (e.g., 0–600 msec), an independent component-analysis-based approach (Srivastava et al., 2005; Mantini et al., 2007) was applied to remove any remaining BCG residual artifacts. Lastly, the EEG was re-referenced to the averaged mastoids and a low-pass filter of 60 Hz was applied. Following the artifact correction, sleep stages were scored in accordance with standard criteria (Iber et al., 2007) using the “VisEd Marks” toolbox¹ for EEGLAB (Delorme and Makeig, 2004).

Slow wave detection

Slow waves were automatically detected from Fz, Cz, and Pz during movement artifact-free NREM sleep (N2 and SWS) via a period amplitude analysis detection algorithm² similar to that previously described (Bersagliere and Achermann, 2010), adapted for EEGlab (Delorme and Makeig, 2004), and written for MATLAB R2019b (The MathWorks Inc.). First, the EEG signal was band-pass filtered (32nd order Chebyshev Type 2 low-pass filter, 80 dB stopband attenuation, 2.15 Hz frequency cut-off; 64th-order Chebyshev type 2 high-pass filter, 80 dB stopband attenuation, 0.46 Hz frequency cut-off). The cut-off frequencies were selected to achieve minimal attenuation in the band of interest while keeping a good attenuation of the neighboring frequencies. The filters were applied in the forward and reverse directions to achieve zero-phase distortion. Next, half-waves were determined as negative or positive deflections between two consecutive zero crossings in the band-pass filtered signal for frequencies between 0.5 and 2 Hz. Only adjacent half-waves with a peak-to-peak amplitude higher than 75 μV and longer than 0.250 s were considered for the analysis. The latency of the negative peak of each slow wave was extracted.

Spindle detection

Sleep spindles were automatically detected from Fz, Cz, and Pz during movement artifact-free NREM sleep (N2 and SWS) using an established (Fogel et al., 2014; Albouy et al., 2015) and validated (Ray et al., 2015) method employing EEGlab-compatible (Delorme and Makeig, 2004) software³ written for MATLAB R2019b (The MathWorks Inc.). Detailed processing steps, procedures, and method validation are reported elsewhere (Ray et al., 2015). Briefly, the spindle data were extracted from movement artifact-free, NREM epochs. The detection method (Ray et al., 2015) used a complex demodulation transformation of the EEG signal with a bandwidth of 5 Hz centered about a carrier frequency of 13.5 Hz (i.e., 11–16 Hz) (Iber et al., 2007). The method employs and adaptive amplitude threshold at the 99th percentile on the transformed signal. Following automatic detection, detected spindles were visually verified by a single expert scorer. The variables of interest extracted from this method include spindle amplitude, duration, and density (number of spindles/minute of NREM sleep) for each participant. This method is consistent with several previous EEG-fMRI studies investigating brain activations time-locked to spindles (Laufs et al., 2007; Schabus et al., 2007; Tyvaert et al., 2008; Andrade et al., 2011; Bergmann et al., 2012; Caporro et al., 2012; Fogel et al., 2017a; Fang et al., 2019, 2020; Smith et al., 2020). From the spindle data, the onset of and peak of each spindle event were extracted.

Slow wave-spindle coupling

Using the slow wave negative peak latencies and the spindle peak latencies from Fz, Cz, and Pz, we performed coupling detection procedures using the approach originally developed and validated by Mölle et al. (2002, 2011), Clemens et al. (2011), Staresina et al. (2015), Baena et al. (2020), Denis et al. (2021) employing EEGlab-compatible (Delorme and Makeig, 2004) software written for MATLAB R2019b (The MathWorks Inc.). This procedure involved building a time window of 4 s around the negative peaks of the detected slow waves. The spindles were then marked as coupled SW-SP complexes when the spindle peak latency occurred within the 4-s time window of the slow wave peak detected in the same channel as the spindle event. Alternately, spindles were marked as uncoupled when the spindle peak latency occurred outside the 4-s time window. Due to the low number of coupled spindle events on each channel, detections from Fz, Cz, and Pz were merged and used in the final analyses. To address duplicate detection of the same slow wave event by the different channels, we removed events which latency difference was lower than the minimum duration threshold for the half wave (0.125 s). Lag was measured as the distance between the slow wave negative peak and the spindle onset. The lag variance was then calculated for each individual as a measure of coupling strength. Paired sampled t-test were used for the comparison of the SW-SP complexes, measured in time bins of 200 ms along the 4 s window.

Additionally, the phase of the bandpass-filtered slow-wave signal in radians at the spindle peak latency was computed. The mean direction of the phase angles for all coupled spindle events were determined using the CircStat toolbox (Berens, 2009). Hilbert transform was applied to extract the preferred phase of SW-SP coupling for each participant averaging all individual events preferred phases. Then, we perform uniformity tests (Rayleigh test) and uniformity using positive slow wave peaks as the predefined mean direction (V-test). From the SW-SP coupling data, and in addition to a list of SW-SP and uncoupled spindle events for each recording, the spindle peak oscillatory frequency and lag for the SW-SP events were extracted.

MRI acquisition and analysis

MRI recording parameters

A 3.0 Tesla Magnetom Prisma MRI system (Siemens) with a 64-channel head coil were used to obtain brain images. Structural T1-weighted MRI images were acquired using a 3-D MPRAGE sequence [repetition time (TR) = 2,300 msec, echo time = 2.98 msec, inversion time = 900 msec, flip angle = 9°, 176 slices, field of view = 256 mm× 256 mm, matrix size = 256 × 256 × 176, voxel size = 1 mm× 1 mm× 1 mm]. During the sleep session, multislice T2*-weighted fMRI images were acquired using an EPI sequence using axial slice orientation (TR = 2,160 msec, echo time = 30 msec, flip angle = 90°, 40 transverse slices, 3-mm slice thickness, 10% interslice gap, field of view = 220 mm× 220 mm, matrix size = 64 × 64 × 40, voxel size = 3.44 mm× 3.44 mm× 3 mm). Up to 2 h of sleep EEG-fMRI data was acquired among all participants during fMRI acquisition.

Image preprocessing

Functional images were preprocessed and analyzed using SPM12 (Welcome Department of Imaging Neuroscience)⁴ implemented in MATLAB R2015a (The Mathworks Inc.). fMRI images were corrected for slice acquisition time differences and realigned to correct head motion using rigid body transformation. A mean realigned image was then created from the resulting images. The structural T1 image was coregistered to this mean functional image using a rigid body transformation optimized to maximize the normalized mutual information between the two images. Coregistration parameters were then applied to the realigned BOLD time series. The co-registered structural images were segmented into gray matter, white matter, and cerebrospinal fluid. An average participant-based template was created using DARTEL in SPM12. All functional and anatomical images were spatially normalized using the resulting template, which was generated from the structural scans. Finally, spatial smoothing was applied to all functional images (Gaussian kernel, 8-mm FWHM).

Statistical analysis

First-level (within-subject) GLM

The onset and duration for each SW-SP and uncoupled spindle and slow wave events were identified from the EEG data and were considered events of interest. Friston-24 movement parameters (Friston et al., 1996), the mean white matter intensity, and the mean cerebral spinal fluid intensity for each participant were entered into the model as nuisance variables. A high-pass filter with a cut-off at 128 s was used to remove low frequency drifts from the time series. Brain activations time-locked to each event type (e.g., spindle, SW, and coupled SW-SP events) were estimated using the onset and duration of each event in a fixed-effects GLM using an event-related fMRI design. The BOLD time series data were modeled using the canonical hemodynamic response function (HRF). In addition, to account for the variability in the latency of the peak response and variability in the duration of the peak response, the temporal and dispersion basis function derivatives were also included in the model [this approach is commonly used for similar sleep EEG-fMRI data sets. See (Schabus et al., 2007; Dang-Vu et al., 2008; Urbain et al., 2013) for details]. Consequently, this approach yields an “informed basis set,” generating three contrast t-maps, one for the canonical HRF, one for the temporal derivative, and one for the dispersion derivative, for each participant. The resulting contrast images (HRF, temporal, and dispersion derivates) are then entered in a second-level (i.e., group-level), random effects analysis (Penny and Holmes, 2007). Because latency and dispersion derivative effects cannot be computed separately from their canonical HRF at the random effect level (Henson et al., 2002), we computed F tests to jointly test relevant contrasts of parameter estimates. This allowed us to identify cerebral regions in which BOLD response was modulated by the events of interest (e.g., spindles, slow waves, coupled events), using the HRF, and its onset latency and temporal dispersion derivatives (Friston et al., 1998).

Second-level (group) GLM

The first-level contrasts maps of coupled SW-SP and uncoupled spindles, and uncoupled SW events comprised of the informed basis set (contrasts for the HRF, temporal, and dispersion basis functions) were entered together into a second-level ANOVA. In this way, we could examine the HRF response for each group comparison without the assumption that the temporal and dispersion functions were constant (Henson et al., 2001; Henson and Penny, 2005). It is necessary to use an ANOVA to include the informed basis set. However, the resulting F-contrast does not provide information about the directionality of the effects. Thus, in order to determine the direction of the effects, the contrast estimates were examined and 95% confidence intervals for the HRF, for each significant cluster. Group-level F-contrasts were generated to examine: (1) brain activation time-locked to SW-coupled spindles, (2) brain activation time-locked to uncoupled spindles; (3) brain activation time-locked to uncoupled slow waves; (4) differences in activation between coupled and uncoupled spindles; (5) differences in activation between coupled and uncoupled slow waves. Finally, a conjunction analysis was also employed to test for similarities between uncoupled spindles and coupled SW-SP and also between uncoupled slow waves and coupled SW-SP. All group-level hypotheses were tested controlling for multiple comparisons using an uncorrected threshold of p < 0.001 and then family-wise error (FWE) corrected, at p < 0.05, at the cluster-level.

Results

Sleep architecture

All participants had more than a total of 14 min of sleep during the recording session in the MRI scanner, and on average, 46.10 min (SD = 23.42 min). The average sleep latency was 8.21 min (SD = 10.64 min), and the average time when participants fell asleep in the scanner was at 22:20 h (SD = 28 min). NREM data from 21 participants were included in the analyses (Table 1). The average duration of NREM sleep was 41.14 min (SD = 19.91 min).

TABLE 1

Table 1. Sleep architecture and spindle parameters for coupled and uncoupled spindles during EEG–fMRI recording.

Coupled SW-SP and uncoupled spindle parameters

Coupling status was determined using a previously reported and validated procedure developed by Mölle et al. (2002, 2011), Clemens et al. (2011), Staresina et al. (2015), Baena et al. (2020), Denis et al. (2021). Accordingly, spindles were marked as coupled if the spindle peak latency occurred within a 4 s window built around the slow wave negative peak. EEG trace example of coupled and uncoupled spindles are shown in Supplementary Figure 1. On average, 271 coupled spindles and 184 uncoupled spindles per participant were identified and included in the analyses (Table 1). In line with previous research (Baena et al., 2022; Malerba et al., 2022), spindle amplitude was significantly higher in coupled SW-SP in comparison to uncoupled spindles (t₍₂₀₎ = 2.65, p = 0.01). There were no differences between coupled SW-SP and uncoupled spindles in terms of the number (t₍₂₀₎ = 1.56, p = 0.13) duration (t₍₂₀₎ = −1.97, p = 0.06) or density (t₍₂₀₎ = 1.39, p = 0.18), hence suggesting they were evenly distributed.

Inspection of the temporal distribution of spindles co-occurring within the 4 sec window around the negative peak of the slow wave (Figure 2A) revealed that coupled spindles percentage during the downstate (up-to-down transition interval) ranging from −1 to 0 sec was significantly higher than the preceding interval ranging from −2 to −1 sec (t₍₂₀₎ = −5.45, p < 0.001). The percentage of coupled spindles during the upstate (down to up interval ranging from 0 to 1 sec), was also higher than the succeeding 1 to 2 sec interval (t₍₂₀₎ = 5.89, p < 0.001). Finally, the percentage of coupled spindles during the upstate was significantly higher than during the downstate (t₍₂₀₎ = −3.19, p = 0.005).

FIGURE 2

Figure 2. (A) Coupled spindle percentage histogram for all participants. Each bar represents the number of coupled spindles detected in an interval of 200 ms divided by the total number of spindles averaged across all participants. Average slow wave coupling window for all participants is superimposed In black. (B) Circular plot of preferred phase for each individual (slow wave phase at spindle amplitude maximum). Red dots denote an individual preferred phase (0° slow wave -upstate peak, ± 180° slow wave downstate peak). The direction of the line indicates the preferred direction of the grand average. Most individuals exhibit spindles adjacent to or immediately following the positive slow wave peak at 0°.

Next, to investigate the phase relationship between spindles and slow waves, the slow wave phase at the spindle peak for each SW-SP event was extracted and mean phase angle per individual was computed (Figure 2B). The coupling of spindle events within the slow wave cycle was maximal shortly before the upstate peak in 13 out of 21 participants [0° (slow wave positive peak); p < 0.001, V-test]. Further individual-level analyses revealed a non-uniform distribution (p < 0.05, Rayleigh test) of the preferred phases of SW-SP in 10 out of 21 participants, suggesting that spindles were coupled to slow waves preferentially adjacent to or immediately following the positive slow wave peak.

Cerebral activations time-locked to uncoupled (i.e., isolated) spindles

Activations time-locked to uncoupled spindles were observed in the bilateral thalamus, middle cingulate cortex (MCC), right posterior cingulate cortex (PCC), superior medial frontal cortex (sPFC), and bilateral occipital lobe (Figure 3A and Table 2).

FIGURE 3

Figure 3. Cerebral activations time-locked to: (A) SW-coupled spindles; (B) uncoupled spindles and (C) uncoupled slow waves. Statistical inferences were performed at p < 0.05, FWE corrected at the cluster level.

TABLE 2

Table 2. Cerebral activations time-locked to: (A) uncoupled spindles, (B) uncoupled slow waves, and (C) coupled SW-SP.

Cerebral activations time-locked to uncoupled (i.e., isolated) slow waves

Uncoupled slow waves showed a very similar pattern of activation (Figure 3C and Table 2) as compared to coupled SW-SP (see below, and Figure 3B), including the occipital lobe, bilateral postcentral gyrus, bilateral superior temporal gyrus (STG), bilateral IFG, bilateral thalamus, anterior cingulate cortex, PCC and bilateral hippocampus (Figure 3C and Table 2).

Cerebral activations time-locked to coupled SW-SP

Activations time-locked to SW-SP were observed in the occipital lobe, bilateral postcentral gyrus, bilateral STG, bilateral middle frontal gyrus (MFG), bilateral inferior frontal gyrus (IFG), bilateral thalamus, bilateral putamen, and hippocampus (Figure 3B and Table 2).

Similarities and differences between SW-SP and uncoupled (i.e., isolated) spindles

A conjunction analysis showed that both SW-coupled and uncoupled spindles presented robust activation of the bilateral thalamus and bilateral occipital lobe (Figure 4B). Further comparisons showed differences in fMRI activity between coupled SW-SP and uncoupled spindles were observed in the cerebellum, bilateral parahippocampus, anterior/middle/posterior cingulate cortex (ACC/MCC/PCC), precuneus and right postcentral gyrus (Figure 4A and Table 3). Inspection of the HRF contrasts verified that coupled SW-SP were associated with larger changes in fMRI signal than uncoupled spindles in the cerebellum and bilateral parahippocampus (Figure 4A). Uncoupled spindles on the other hand, showed greater activation in the cingulate cortex, precuneus (Pcu) and right postcentral gyrus compared to coupled SW-SP (Figure 4A). Thus suggesting that SW-SP coupled events differ from uncoupled spindles, notably in hippocampal regions, but share recruitment of areas critical for generating these events, and in sleep maintenance, such as the thalamus. BOLD-signal peristimulus plots for selected areas are included in Supplementary Figure 2.

FIGURE 4

Figure 4. Differences in cerebral activation between (A) coupled vs. uncoupled spindles, and, (B) conjunction between coupled and uncoupled spindles. Statistical inferences were performed at p < 0.05, FWE was corrected at the cluster level.

TABLE 3

Table 3. Differences in cerebral activations between: (A) coupled SW-SP vs. uncoupled spindles, and (B) coupled SW-SP and uncoupled slow waves.

Similarities and differences between SW-SP and uncoupled (i.e., isolated) slow waves

Conjunction analysis revealed the common areas in both coupled SW-SP and uncoupled slow waves, including the occipital lobe, bilateral postcentral gyrus, bilateral superior temporal gyrus (STG), bilateral IFG, bilateral thalamus and bilateral hippocampus (Figure 5B). Only the PCC was significantly greater in coupled SW-SP than that in response to uncoupled SW (Figure 5A and Table 3). Thus, suggesting that SW-SP coupled events are mostly indistinguishable from uncoupled slow waves, except for the PCC.

FIGURE 5

Figure 5. Differences in cerebral activation between (A) coupled vs. uncoupled slow waves, and (B) conjunction between coupled and uncoupled slow waves. Statistical inferences were performed at p < 0.05, FWE was corrected at the cluster level.

Discussion

In the present study, we investigated functional dissociations in cerebral activation that are specifically related to the occurrence of SW-SP coupling. Previously, converging evidence from multiple studies has implicated SW-SP coupling as an important mechanism supporting the consolidation of new information into long term memory (Novitskaya et al., 2016; Latchoumane et al., 2017; Bergmann and Born, 2018; Baena et al., 2020; Ngo et al., 2020). However, no study to date has specifically examined the patterns of functional brain activations of coupled SW-SP. The results of our study revealed a clear differentiation between regions recruited during coupled SW-SP versus uncoupled spindles as well as during, coupled SW-SP vs. uncoupled SWs, namely: (1) coupled SW-SP uniquely activated bilateral middle frontal gyrus and the bilateral putamen, (2) uncoupled spindles uniquely activated the MCC and sPFC, and, (3) uncoupled slow waves uniquely activated the ACC. Areas that were activated in all three contexts included the bilateral thalamus and bilateral occipital lobe. Furthermore, compared to uncoupled spindles, coupled SW-SP showed relatively greater activation of the cerebellum and bilateral parahippocampal cortex, and relatively lower activation of the ACC, MCC, PCC, precuneus, and right postcentral gyrus. Activation of the PCC was observed in both uncoupled spindles as well as uncoupled slow waves, while coupled slow waves showed greater deactivation in the PCC when compared to uncoupled slow waves. Activation observed in both SW-SP and uncoupled slow waves included bilateral postcentral gyrus, STG, IFG, and interestingly, the bilateral hippocampus. Taken together, these results support the view that the co-occurrence of spindles and slow waves acts as a unique neurophysiological mechanism, resulting in the activation of brain areas not observed when either spindles or slow waves occur in isolation.

Our results are consistent with previous studies using EEG-fMRI to measure cerebral activation and functional connectivity associated with sleep spindles (Schabus et al., 2007; Andrade et al., 2011; Caporro et al., 2012; Fogel et al., 2017b). The consensus from these studies suggests that the recruitment of the hippocampus associated with sleep spindles is a spindle-driven phenomenon. Our current results advance our understanding of this process and suggest that this relationship is more nuanced in important ways; in the absence of slow waves, spindles do not show hippocampal activation, and yet slow waves show hippocampal activation in the absence of spindles. Moreover, there were no specific brain structures that showed activation common to both coupled SW-SP and uncoupled spindles, except for those that were also activated during uncoupled slow waves and are involved in the generation of these events, and in sleep maintenance (i.e., thalamus). However, both coupled SW-SP and uncoupled slow waves showed hippocampal activation, whereas uncoupled spindles did not. In light of these findings, we postulate that hippocampal activation (and possibly, reactivation of brain areas) during sleep is not unique to spindles, but primarily driven by slow waves, and that the co-occurrence of spindles allows the recruitment of frontal areas, creating a unique functional activation pattern, which may be critical for the hippocampal-neocortical dialogue that preferentially occurs during sleep, and is known to support memory consolidation processes (Rasch and Born, 2013). However, the specific behavioral outcomes for reactivation of SW-SP coupling for memory consolidation, per se, remain to be directly investigated.

According to our results, the bilateral putamen were uniquely activated during coupled SW-SP. Putamen has been shown to play a central role in the consolidation of motor skills during sleep (Fogel et al., 2014; King et al., 2017; Vahdat et al., 2017) and also to be correlated with high-level cognitive functions beyond organization and execution of movements (Ketteler et al., 2008; Mestres-Missé et al., 2010). Based on related studies we speculate that activation of the putamen in this context may serve as a marker for coupling-related memory consolidation during sleep, perhaps via memory trace reactivation (Fogel et al., 2017b). Interestingly, while hippocampal activation did not differ between coupled and uncoupled spindles, coupled spindles did show greater activation of the parahippocampus and cerebellum. Thus, it may be that sleep spindles per se do not support hippocampal functions in isolation. Rather, this may preferentially occur when spindles and slow waves are coupled. This is consistent with the extant literature which suggests that SW-SP coupling supports the transfer of slow-wave-driven hippocampal signals into cortical and cerebellar structures, ultimately integrating memories into existing schema encoded in these areas. For example, slow wave-driven reactivation of a freshly encoded memory trace in the hippocampus may serve to maintain the integrity of the new memories’ neural code (Alger et al., 2012; Schlichting and Preston, 2015; Moyano et al., 2019). Coupling of spindles with slow waves then provides an opportunity for integration of the memory into the parahippocampal cortex to be output to other brain areas for long-term storage, where contextual information is processed and eventually integrated into a predominantly cortical, and more schematized memory trace (Schlichting and Preston, 2015; Kitamura et al., 2017; Noack et al., 2021). Moreover, the parahippocampal gyrus is of specific importance for memory of spatial configurations of objects, but not necessarily their identity (Bohbot et al., 2015); thus, suggesting that spindles may be more directly involved in memory processing for spatial aspects of memory processing (Albouy et al., 2013). While we cannot conclude specific mnemonic functions from the current study, this study lays the groundwork for future studies to directly investigate this at the behavioral level.

Uncoupled slow wave events showed an activation pattern very similar to the coupled SW-SP events with exception of the PCC region. This region has been previously associated with memory consolidation processes linking episodic and semantic information (Maguire et al., 1999). More recent studies using MRI confirm this integrating role, further detailing that PCC plays a crucial role in memory schema formation (Bird et al., 2015). Our results thus suggest that although hippocampal activation is shared between coupled SW-SP and uncoupled SWs, the lower level of PCC deactivation during coupled SW-SP could indicate that memory integration and schema formation may only be attributable to coupled SW-SP events. This, however, remains to be directly investigated.

There are several limitations worth mentioning in the present study. It was not possible to subdivide spindles into slow spindles (11–13 Hz) and fast spindles (13–16 Hz), nor was it possible to subdivide spindles from NREM2 and SWS due to the limited sleep duration and number of events recorded while participants were sleeping in the MRI scanner. In addition, we did not directly investigate any relationship between spindles, slow waves and SW-coupled spindles following new learning or associated with cognitive abilities. This remains to be directly investigated in future studies specifically designed to assess the reactivation of memory traces and the related behavioral outcomes. However, this study does represent a necessary first step in this line of investigation.

In summary, we investigated differences in cerebral activation linked to SW-coupled spindles and both their differences and similarities with uncoupled spindles and uncoupled SWs. The results demonstrate a clear differentiation in brain areas recruited time-locked to SW-coupled spindles and uncoupled spindles supporting that both events are functionally dissociable. By contrast, these results show largely similar activation patterns between SW-coupled spindles and uncoupled slow waves. For the first time, we provide direct evidence of hippocampal activation linked specifically to coupled SW-SP complexes in comparison to isolated spindles. This finding is consistent with other studies that show coordinated activity of these precisely timed complexes is specifically involved in memory consolidation via reactivation. Furthermore, and surprisingly, the present study revealed that this pattern of activation may be driven more so by slow waves than spindles. We propose that it is possible that hippocampal reactivation during sleep is not uniquely related to spindles. Rather, this pattern of activation may be primarily driven by slow waves. Our results also suggest that SW-spindle coupling may be critical for functional recruitment of the putamen. Importantly, the co-occurrence of spindles and slow waves may allow the recruitment of frontal areas, which may be critical for the hippocampal-neocortical dialogue that preferentially occurs during sleep. Future studies linking coupling-driven activation with memory performance are needed to further elucidate the dissociable contribution of isolated spindles, isolated slow waves and SW-SP complexes to memory consolidation and other functions.

Author’s note

Sleep is known to aid in memory consolidation. This involves transfer of information from the hippocampus to the neocortex and is facilitated by slow wave – spindle – hippocampal ripples complexes during sleep. However, no study to date has specifically examined if the patterns of functional brain activations of isolated spindles differ from slow wave coupled spindles. Here, we found a clear differentiation in brain areas recruited for each event. Coupled spindles recruited frontal areas and parahippocampus in comparison to uncoupled spindles. This suggests that coupled spindles may be critical for the hippocampal-neocortical dialogue that occurs during sleep. Our results lay the groundwork for future studies to differentiate between spindles according to coupling status.

Data availability statement

The datasets presented in this article are available by request to the corresponding author, following required ethical approval and inter-institutional data sharing agreements. Requests to access the datasets should be directed to SF, c2ZvZ2VsQHVvdHRhd2EuY2E=.

Ethics statement

The studies involving human participants were reviewed and approved by Western University Health Science Research Ethics Board. The patients/participants provided their written informed consent to participate in this study.

Author contributions

SF and AO designed the study and supervised the research. SF, ZF, DB, and LR carried out the research and collected the data and contributed to the data analyses. DB, SF, and LR prepared the figures and wrote the manuscript. AG, DS, JD, and SF revised the manuscript. All authors discussed the results and commented on the manuscript.

Funding

This research was funded by a Canada Excellence Research Chair Grant (#215063 to AO); the Natural Sciences and Engineering Research Council (NSERC) Discovery grant (RGPIN/2017–04328 to SF); and Ministry of Research and Innovation Early Researcher Award (ERA; ER17–13-023 to SF).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.1090045/full#supplementary-material

Footnotes

1. ^ https://github.com/jadesjardins/vised_marks
2. ^ https://github.com/stuartfogel/Period-Amplitude-Analysis
3. ^ https://github.com/stuartfogel/detect_spindles
4. ^ www.fil.ion.ucl.ac.uk/spm/software/spm12/

References

Albouy, G., Fogel, S. M., King, B. R., Laventure, S., Benali, H., Karni, A., et al. (2015). Maintaining vs. Enhancing motor sequence memories: Respective roles of striatal and hippocampal systems. NeuroImage 108, 423–434. doi: 10.1016/j.neuroimage.2014.12.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Albouy, G., Vandewalle, G., Sterpenich, V., Rauchs, G., Desseilles, M., Balteau, E., et al. (2013). Sleep stabilizes visuomotor adaptation memory: A functional magnetic resonance imaging study. J. Sleep Res. 22, 144–154. doi: 10.1111/j.1365-2869.2012.01059.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Alger, S. E., Lau, H., and Fishbein, W. (2012). Slow wave sleep during a daytime nap is necessary for protection from subsequent interference and long-term retention. Neurobiol. Learn. Mem. 98, 188–196. doi: 10.1016/j.nlm.2012.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, P. J., Josephs, O., and Turner, R. (2000). A method for removing imaging artifact from continuous EEG recorded during functional MRI. Neuroimage 12, 230–239. doi: 10.1006/nimg.2000.0599

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, P. J., Polizzi, G., Krakow, K., Fish, D. R., and Lemieux, L. (1998). Identification of EEG events in the MR scanner: The problem of pulse artifact and a method for its subtraction. NeuroImage 8, 229–239. doi: 10.1006/nimg.1998.0361

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrade, K. C., Spoormaker, V. I., Dresler, M., Wehrle, R., Holsboer, F., Sämann, P. G., et al. (2011). Sleep spindles and hippocampal functional connectivity in human NREM sleep. J. Neurosci. 31, 10331–10339. doi: 10.1523/JNEUROSCI.5660-10.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Baena, D., Cantero, J. L., Fuentemilla, L., and Atienza, M. (2020). Weakly encoded memories due to acute sleep restriction can be rescued after one night of recovery sleep. Sci. Rep. 10:1449. doi: 10.1038/s41598-020-58496-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Baena, D., Fang, Z., Ray, L. B., Owen, A. M., and Fogel, S. M. (2022). Brain activations time locked to slow wave-coupled sleep spindles correlates with intellectual abilities. Cereb. Cortex bhac428. doi: 10.1093/cercor/bhac428

PubMed Abstract | CrossRef Full Text | Google Scholar

Beck, A. T., Epstein, N., Brown, G., and Steer, R. A. (1988a). An inventory for measuring clinical anxiety: Psychometric properties. J. Consult. Clin. Psychol. 56, 893–897. doi: 10.1037//0022-006x.56.6.893

CrossRef Full Text | Google Scholar

Beck, A. T., Steer, R. A., and Carbin, M. G. (1988b). Psychometric properties of the Beck Depression Inventory: Twenty-five years of evaluation. Clin. Psychol. Rev. 8, 77–100. doi: 10.1016/0272-7358(88)90050-5

CrossRef Full Text | Google Scholar

Berens, P. (2009). CircStat: A MATLAB toolbox for circular statistics. J. Stat. Softw. 31, 1–21. doi: 10.18637/JSS.V031.I10

CrossRef Full Text | Google Scholar

Bergmann, T. O., and Born, J. (2018). Phase-amplitude coupling: A general mechanism for memory processing and synaptic plasticity? Neuron 97, 10–13. doi: 10.1016/J.NEURON.2017.12.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergmann, T. O., Mölle, M., Schmidt, M. A., Lindner, C., Marshall, L., Born, J., et al. (2012). EEG-guided transcranial magnetic stimulation reveals rapid shifts in motor cortical excitability during the human sleep slow oscillation. J. Neurosci. 32, 243–253. doi: 10.1523/JNEUROSCI.4792-11.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Bersagliere, A., and Achermann, P. (2010). Slow oscillations in human non-rapid eye movement sleep electroencephalogram: Effects of increased sleep pressure. J. Sleep Res. 19, 228–237. doi: 10.1111/j.1365-2869.2009.00775.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Betta, M., Handjaras, G., Leo, A., Federici, A., Farinelli, V., Ricciardi, E., et al. (2021). Cortical and subcortical hemodynamic changes during sleep slow waves in human light sleep. NeuroImage 236:118117. doi: 10.1016/J.NEUROIMAGE.2021.118117

PubMed Abstract | CrossRef Full Text | Google Scholar

Bird, C. M., Keidel, J. L., Ing, L. P., Horner, A. J., and Burgess, N. (2015). Consolidation of complex events via reinstatement in posterior cingulate cortex. J. Neurosci. 35, 14426–14434. doi: 10.1523/JNEUROSCI.1774-15.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Bohbot, V. D., Allen, J. J. B., Dagher, A., Dumoulin, S. O., Evans, A. C., Petrides, M., et al. (2015). Role of the parahippocampal cortex in memory for the configuration but not the identity of objects: Converging evidence from patients with selective thermal lesions and fMRI. Front. Hum. Neurosci. 9:431. doi: 10.3389/fnhum.2015.00431

PubMed Abstract | CrossRef Full Text | Google Scholar

Buzsáki, G. (2015). Hippocampal sharp wave-ripple: A cognitive biomarker for episodic memory and planning. Hippocampus 25, 1073–1188. doi: 10.1002/HIPO.22488

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporro, M., Haneef, Z., Yeh, H. J., Lenartowicz, A., Buttinelli, C., Parvizi, J., et al. (2012). Functional MRI of sleep spindles and K-complexes. Clin. Neurophysiol. 123, 303–309. doi: 10.1016/J.CLINPH.2011.06.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemens, Z., Mölle, M., Eross, L., Jakus, R., Rásonyi, G., Halász, P., et al. (2011). Fine-tuned coupling between human parahippocampal ripples and sleep spindles. Eur. J. Neurosci. 33, 511–520. doi: 10.1111/j.1460-9568.2010.07505.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang-Vu, T. T., Schabus, M., Desseilles, M., Albouy, G., Boly, M., Darsaud, A., et al. (2008). Spontaneous neural activity during human slow wave sleep. Proc. Natl. Acad. Sci. U.S.A. 105:15160. doi: 10.1073/pnas.0801819105

PubMed Abstract | CrossRef Full Text | Google Scholar

Delorme, A., and Makeig, S. (2004). EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21. doi: 10.1016/j.jneumeth.2003.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Denis, D., Mylonas, D., Poskanzer, C., Bursal, V., Payne, J. D., and Stickgold, R. (2021). Sleep spindles preferentially consolidate weakly encoded memories. J. Neurosci. 41, 4088–4099. doi: 10.1523/JNEUROSCI.0818-20.2021

PubMed Abstract | CrossRef Full Text | Google Scholar

Diekelmann, S., and Born, J. (2010). The memory function of sleep. Nat. Rev. Neurosci. 11, 114–126. doi: 10.1038/nrn2762

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglass, A. B., Bornstein, R. F., Nino-Murcia, G., Keenan, S., Miles, L., Zarcone, V. P., et al. (1994). The sleep disorders questionnaire. I: Creation and multivariate structure of SDQ. Sleep 17:160. doi: 10.1093/sleep/17.2.160

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudai, Y., Karni, A., and Born, J. (2015). The consolidation and transformation of memory. Neuron 88, 20–32. doi: 10.1016/J.NEURON.2015.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Eichenbaum, H. (2017). Memory: Organization and control. Annu. Rev. Psychol. 68, 19–45. doi: 10.1146/ANNUREV-PSYCH-010416-044131

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Z., Ray, L. B. B., Owen, A. M. M., and Fogel, S. M. (2019). Brain activation time-locked to sleep spindles associated with human cognitive abilities. Front. Neurosci. 13:46. doi: 10.3389/fnins.2019.00046

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, Z., Ray, L. B., Houldin, E., Smith, D., Owen, A. M., and Fogel, S. M. (2020). Sleep spindle-dependent functional connectivity correlates with cognitive abilities. J. Cogn. Neurosci. 32, 446–466. doi: 10.1162/jocn_a_01488

CrossRef Full Text | Google Scholar

Fogel, S., Vien, C., Karni, A., Benali, H., Carrier, J., and Doyon, J. (2017a). Sleep spindles: A physiological marker of age-related changes in gray matter in brain regions supporting motor skill memory consolidation. Neurobiol. Aging 49, 154–164. doi: 10.1016/J.NEUROBIOLAGING.2016.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Fogel, S., Albouy, G., King, B. R., Lungu, O., Vien, C., Boré, A., et al. (2017b). Reactivation or transformation? Motor memory consolidation associated with cerebral activation time-locked to sleep spindles. PLoS One 12:e0174755. doi: 10.1371/journal.pone.0174755

PubMed Abstract | CrossRef Full Text | Google Scholar

Fogel, S. M., Albouy, G., King, B. R., Vien, C., Karni, A., Benali, H., et al. (2014). Motor memory consolidation depends upon reactivation driven by the action of sleep spindles. J. Sleep Res. 23:47. doi: 10.1111/jsr.12213

CrossRef Full Text | Google Scholar

Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S. J., and Turner, R. (1996). Movement-related effects in fMRI time-series. Magn. Reson. Med. 35, 346–355. doi: 10.1002/MRM.1910350312

PubMed Abstract | CrossRef Full Text | Google Scholar

Friston, K., Josephs, O., Rees, G., and Turner, R. (1998). Nonlinear event-related responses in fMRI. Magn. Reson. Med. 39, 41–52.

Google Scholar

Helfrich, R. F., Lendner, J. D., Mander, B. A., Guillen, H., Paff, M., Mnatsakanyan, L., et al. (2019). Bidirectional prefrontal-hippocampal dynamics organize information transfer during sleep in humans. Nat. Commun. 10:3572. doi: 10.1038/s41467-019-11444-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Henson, R. N. A., and Penny, W. D. (2005). ANOVAs and SPM. London: Wellcome Department of Imaging Neuroscience.

Google Scholar

Henson, R. N. A., Price, C. J., Rugg, M. D., Turner, R., and Friston, K. J. (2002). Detecting latency differences in event-related BOLD responses: Application to words versus nonwords and initial versus repeated face presentations. NeuroImage 15, 83–97. doi: 10.1006/nimg.2001.0940

PubMed Abstract | CrossRef Full Text | Google Scholar

Henson, R., Rugg, M., and Friston, K. (2001). The choice of basis functions in event-related fMRI. Neuroimage 13:149. doi: 10.1016/S1053-8119(01)91492-2

CrossRef Full Text | Google Scholar

Iber, C., Ancoli-Israel, S., Chesson, A. L., and Quan, S. F. (2007). The AASM manual for the scoring of sleep and associated events: Rules, terminology and technical specifications, 1st Edn. Westchester, IL: American Academy of Sleep Medicine.

Google Scholar

Kaufmann, C., Wehrle, R., Wetter, T. C., Holsboer, F., Auer, D. P., Pollmächer, T., et al. (2006). Brain activation and hypothalamic functional connectivity during human non-rapid eye movement sleep: An EEG/fMRI study. Brain 129, 655–667. doi: 10.1093/BRAIN/AWH686

PubMed Abstract | CrossRef Full Text | Google Scholar

Ketteler, D., Kastrau, F., Vohn, R., and Huber, W. (2008). The subcortical role of language processing. High level linguistic features such as ambiguity-resolution and the human brain; an fMRI study. NeuroImage 39, 2002–2009. doi: 10.1016/j.neuroimage.2007.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

King, B. R., Saucier, P., Albouy, G., Fogel, S. M., Rumpf, J. J., Klann, J., et al. (2017). Cerebral activation during initial motor learning forecasts subsequent sleep-facilitated memory consolidation in older adults. Cereb. Cortex 27, 1588–1601. doi: 10.1093/CERCOR/BHV347

PubMed Abstract | CrossRef Full Text | Google Scholar

Kitamura, T., Ogawa, S. K., Roy, D. S., Okuyama, T., Morrissey, M. D., Smith, L. M., et al. (2017). Engrams and circuits crucial for systems consolidation of a memory. Science 356, 73–78. doi: 10.1126/SCIENCE.AAM6808/SUPPL_FILE/AAM6808-KITAMURA-SM.PDF

CrossRef Full Text | Google Scholar

Klinzing, J. G., Niethard, N., and Born, J. (2019). Mechanisms of systems memory consolidation during sleep. Nat. Neurosci. 22, 1598–1610. doi: 10.1038/s41593-019-0467-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Latchoumane, C.-F. F. V., Ngo, H.-V. V. V., Born, J., and Shin, H.-S. S. (2017). Thalamic spindles promote memory formation during sleep through triple phase-locking of cortical, thalamic, and hippocampal rhythms. Neuron 95, 424–435. doi: 10.1016/j.neuron.2017.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Laufs, H., Walker, M. C., and Lund, T. E. (2007). ‘Brain activation and hypothalamic functional connectivity during human non-rapid eye movement sleep: An EEG/fMRI study’-its limitations and an alternative approach. Brain 130:e75. doi: 10.1093/brain/awm084

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, P. A., and Durrant, S. J. (2011). Overlapping memory replay during sleep builds cognitive schemata. Trends Cogn. Sci. 15, 343–351. doi: 10.1016/J.TICS.2011.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Maguire, E. A., Frith, C. D., and Morris, R. G. M. (1999). The functional neuroanatomy of comprehension and memory: The importance of prior knowledge. Brain 122, 1839–1850. doi: 10.1093/BRAIN/122.10.1839

PubMed Abstract | CrossRef Full Text | Google Scholar

Malerba, P., Whitehurst, L., and Mednick, S. C. (2022). The space-time profiles of sleep spindles and their coordination with slow oscillations on the electrode manifold. Sleep 45:zsac132. doi: 10.1093/sleep/zsac132

PubMed Abstract | CrossRef Full Text | Google Scholar

Mantini, D., Perrucci, M. G., Cugini, S., Ferretti, A., Romani, G. L., and Gratta, C. D. (2007). Complete artifact removal for EEG recorded during continuous fMRI using independent component analysis. Neuroimage 34, 598–607. doi: 10.1016/j.neuroimage.2006.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Mestres-Missé, A., Rodriguez-Fornells, A., and Münte, T. F. (2010). Neural differences in the mapping of verb and noun concepts onto novel words. NeuroImage 49, 2826–2835. doi: 10.1016/j.neuroimage.2009.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mölle, M., Bergmann, T. O., Marshall, L., and Born, J. (2011). Fast and slow spindles during the sleep slow oscillation: Disparate coalescence and engagement in memory processing. Sleep 34, 1411–1421. doi: 10.5665/SLEEP.1290

PubMed Abstract | CrossRef Full Text | Google Scholar

Mölle, M., Marshall, L., Gais, S., and Born, J. (2002). Grouping of spindle activity during slow oscillations in human non-rapid eye movement sleep. J. Neurosci. 22, 10941–10947.

Google Scholar

Moyano, M. D., Diekelmann, S., Pedreira, M. E., and Forcato, C. (2019). Sleep accelerates re-stabilization of human declarative memories. Neurobiol. Learn. Mem. 162, 1–8. doi: 10.1016/J.NLM.2019.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulert, C., and Lemieux, L. (2010). EEG-fMRI: Physiological basis, technique and applications (illustrate ed.). Berlin: Springer.

Google Scholar

Mullinger, K. J., Yan, W. X., and Bowtell, R. (2011). Reducing the gradient artefact in simultaneous EEG-fMRI by adjusting the subject’s axial position. NeuroImage 54, 1942–1950. doi: 10.1016/j.neuroimage.2010.09.079

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevian, T., and Sakmann, B. (2006). Spine Ca2⁺ signaling in spike-timing-dependent plasticity. J. Neurosci. 26, 11001–11013. doi: 10.1523/JNEUROSCI.1749-06.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ngo, H. V., Fell, J., and Staresina, B. (2020). Sleep spindles mediate hippocampal-neocortical coupling during long-duration ripples. ELife 9: e57011. doi: 10.7554/ELIFE.57011

PubMed Abstract | CrossRef Full Text | Google Scholar

Niethard, N., Ngo, H.-V. V. V., Ehrlich, I., and Born, J. (2018). Cortical circuit activity underlying sleep slow oscillations and spindles. Proc. Natl. Acad. Sci. U.S.A. 115: E9220–E9229. doi: 10.1073/pnas.1805517115

PubMed Abstract | CrossRef Full Text | Google Scholar

Noack, H., Doeller, C. F., and Born, J. (2021). Sleep strengthens integration of spatial memory systems. Learn. Mem. 28, 162–170. doi: 10.1101/LM.053249.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Novitskaya, Y., Sara, S. J., Logothetis, N. K., and Eschenko, O. (2016). Ripple-triggered stimulation of the locus coeruleus during post-learning sleep disrupts ripple/spindle coupling and impairs memory consolidation. Learn. Mem. 23, 238–248. doi: 10.1101/LM.040923.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Penny, W. D., and Holmes, A. J. (2007). “Random effects analysis,” in Statistical parametric mapping: The analysis of functional brain images, (Amsterdam: Elsevier).

Google Scholar

Rasch, B., and Born, J. (2013). About sleep’s role in memory. Physiol. Rev. 93, 681–766. doi: 10.1152/physrev.00032.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ray, L. B., Sockeel, S., Soon, M., Bore, A., Myhr, A., Stojanoski, B., et al. (2015). Expert and crowd-sourced validation of an individualized sleep spindle detection method employing complex demodulation and individualized normalization. Front. Hum. Neurosci. 9:507. doi: 10.3389/fnhum.2015.00507

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadowski, J. H. L. P., Jones, M. W., and Mellor, J. R. (2016). Sharp-wave ripples orchestrate the induction of synaptic plasticity during reactivation of place cell firing patterns in the hippocampus. Cell Rep. 14, 1916–1929. doi: 10.1016/J.CELREP.2016.01.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Schabus, M., Dang-Vu, T. T., Albouy, G., Balteau, E., Boly, M., Carrier, J., et al. (2007). Hemodynamic cerebral correlates of sleep spindles during human non-rapid eye movement sleep. Proc. Natl. Acad. Sci. U.S.A. 104, 13164–13169. doi: 10.1073/PNAS.0703084104

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlichting, M. L., and Preston, A. R. (2015). Memory integration: Neural mechanisms and implications for behavior. Curr. Opin. Behav. Sci. 1, 1–8. doi: 10.1016/j.cobeha.2014.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Schreiner, T., Petzka, M., Staudigl, T., and Staresina, B. P. (2021). Endogenous memory reactivation during sleep in humans is clocked by slow oscillation-spindle complexes. Nat. Commun. 12:3112. doi: 10.1038/s41467-021-23520-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sejnowski, T. J., and Destexhe, A. (2000). Why do we sleep? Brain Res. 886, 208–223. doi: 10.1016/S0006-8993(00)03007-9

CrossRef Full Text | Google Scholar

Smith, D., Fang, Z., Thompson, K., and Fogel, S. M. (2020). Sleep and individual differences in intellectual abilities. Curr. Opin. Behav. Sci. 33, 126–131. doi: 10.1016/j.cobeha.2020.02.011

CrossRef Full Text | Google Scholar

Srivastava, G., Crottaz-Herbette, S., Lau, K. M., Glover, G. H., and Menon, V. (2005). ICA-based procedures for removing ballistocardiogram artifacts from EEG data acquired in the MRI scanner. NeuroImage 24, 50–60. doi: 10.1016/j.neuroimage.2004.09.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Staresina, B. P., Bergmann, T. O., Bonnefond, M., Meij, R. V. D., Jensen, O., Deuker, L., et al. (2015). Hierarchical nesting of slow oscillations, spindles and ripples in the human hippocampus during sleep. Nat. Neurosci. 18, 1679–1686. doi: 10.1038/nn.4119

PubMed Abstract | CrossRef Full Text | Google Scholar

Takehara-Nishiuchi, K. (2020). Neurobiology of systems memory consolidation. Eur. J. Neurosci. 54, 6850–6863. doi: 10.1111/EJN.14694

PubMed Abstract | CrossRef Full Text | Google Scholar

Timofeev, I., Grenier, F., Bazhenov, M., Houweling, A. R., Sejnowski, T. J., and Steriade, M. (2002). Short-and medium-term plasticity associated with augmenting responses in cortical slabs and spindles in intact cortex of cats in vivo. J. Physiol. 542, 583–598. doi: 10.1113/jphysiol.2001.013479

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyvaert, L., LeVan, P., Grova, C., Dubeau, F., and Gotman, J. (2008). Effects of fluctuating physiological rhythms during prolonged EEG-fMRI studies. Clin. Neurophysiol. 119, 2762–2774. doi: 10.1016/j.clinph.2008.07.284

PubMed Abstract | CrossRef Full Text | Google Scholar

Urbain, C., Schmitz, R., Schmidt, C., Cleeremans, A., Bogaert, P. V., Maquet, P., et al. (2013). Sleep-dependent neurophysiological processes in implicit sequence learning. J. Cogn. Neurosci. 25, 2003–2014. doi: 10.1162/jocn_a_00439

CrossRef Full Text | Google Scholar

Vahdat, S., Fogel, S., Benali, H., and Doyon, J. (2017). Network-wide reorganization of procedural memory during NREM sleep revealed by fMRI. ELife 6:e24987. doi: 10.7554/ELIFE.24987

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S.-H., and Morris, R. G. M. (2009). Hippocampal-neocortical interactions in memory formation, consolidation, and reconsolidation. Annu. Rev.Psychol. 61, 49–79. doi: 10.1146/ANNUREV.PSYCH.093008.100523

PubMed Abstract | CrossRef Full Text | Google Scholar

Witkowski, S., Schechtman, E., and Paller, K. A. (2020). Examining sleep’s role in memory generalization and specificity through the lens of targeted memory reactivation. Curr. Opin. Behav. Sci. 33, 86–91.

Google Scholar

Yordanova, J., Kirov, R., Verleger, R., and Kolev, V. (2017). Dynamic coupling between slow waves and sleep spindles during slow wave sleep in humans is modulated by functional pre-sleep activation. Sci. Rep. 7:14496. doi: 10.1038/S41598-017-15195-X

PubMed Abstract | CrossRef Full Text | Google Scholar