- 1Department of Neurology, University of Colorado, Denver, Denver, CO, United States
- 2Pediatric Critical Care Medicine, University of Colorado, Denver, Denver, CO, United States
- 3Department of Neurology, University of Miami Miller School of Medicine, Miami, FL, United States
Slow wave activity (SWA) during sleep is associated with synaptic regulation and memory processing functions. Each cycle of non-rapid-eye-movement (NREM) sleep demonstrates a waxing and waning amount of SWA during the transitions between stages N2 and N3 sleep, and the deeper N3 sleep is associated with an increased density of SWA. Further, SWA is an amalgam of different types of slow waves, each identifiable by their temporal coupling to spindle subtypes with distinct physiological features. The objectives of this study were to better understand the neurobiological properties that distinguish different slow wave and spindle subtypes, and to examine the composition of SWA across cycles of NREM sleep. We further sought to explore changes in the composition of NREM cycles that occur among aging adults. To address these goals, we analyzed subsets of data from two well-characterized cohorts of healthy adults: (1) The DREAMS Subjects Database (n = 20), and (2) The Cleveland Family Study (n = 60). Our analyses indicate that slow wave/spindle coupled events can be characterized as frontal vs. central in their relative distribution between electroencephalography (EEG) channels. The frontal predominant slow waves are identifiable by their coupling to late-fast spindles and occur more frequently during stage N3 sleep. Conversely, the central-associated slow waves are identified by coupling to early-fast spindles and favor occurrence during stage N2 sleep. Together, both types of slow wave/spindle coupled events form the composite of SWA, and their relative contribution to the SWA rises and falls across cycles of NREM sleep in accordance with depth of sleep. Exploratory analyses indicated that older adults produce a different composition of SWA, with a shift toward the N3, frontal subtype, which becomes increasingly predominant during cycles of NREM sleep. Overall, these data demonstrate that subtypes of slow wave/spindle events have distinct cortical propagation patterns and differ in their distribution across lighter vs. deeper NREM sleep. Future efforts to understand how slow wave sleep and slow wave/spindle coupling impact memory performance and neurological disease may benefit from examining the composition of SWA to avoid potential confounds that may occur when comparing dissimilar neurophysiological events.
Introduction
The slow waves of non-rapid-eye-movement (NREM) sleep are hypothesized to coordinate synaptic remodeling and sleep-dependent memory processes (Walker, 2009; Born, 2010). Studies indicate that loss of normal slow wave physiology is related to cognitive impairment and pathophysiological changes of Alzheimer’s disease, suggesting that slow waves may have particular relevance to pathological aging processes (Backhaus et al., 2007; Scullin, 2013; Mander et al., 2015; Lee et al., 2020; Winer et al., 2020; Wunderlin et al., 2020). Efforts to understand the neurobiology of slow waves have noted that slow wave activity (SWA) is not uniform in nature and is instead a heterogeneous composition of different subtypes of slow waves. Within surface EEG, SWA has been separated into subcomponents by several slow wave characteristics, including frequency (Achermann and Borbely, 1997; Amzica and Steriade, 1998; Cox et al., 2019), local vs. global occurrence (Bernardi et al., 2018), and frontal vs. central localization (Murphy et al., 2009).
Slow waves can also be conceptualized as discrete neuronal communication events that mark time windows of hippocampal activity and reactivation of specific neuronal patterns. During these slow wave time windows, neuronal patterns recapitulate the original sequences of memory tracings that first formed during wake within the hippocampi and association cortexes (Wilson and McNaughton, 1994; Varela and Wilson, 2020; Schreiner et al., 2021). In this mechanistic framework, the EEG signals associated with slow waves, termed coupled sleep spindles, are additional components of a “memory replay mode” that occurs during NREM sleep (Latchoumane et al., 2017; Schreiner et al., 2021). Intracranial recordings further support a distinction between subtypes of slow waves based on their temporal coupling with anterior vs. posterior hippocampal activity, and the relative contribution of these distinct slow wave and coupled spindle components to SWA differs between stages N2 and N3 sleep (Jiang et al., 2019a,b). Within surface EEG, two distinct subtypes of slow waves can also be identified by their coupling to unique spindle subtypes, and congruent with intracranial observations, the contribution of each coupled slow wave/spindle subtype to SWA varies between stages N2 and N3 sleep (McConnell et al., 2021).
Experimental manipulation of SWA-associated memory functions in animal models further supports a physiological distinction between different subtypes of coupled slow wave/spindle events, with one subtype of coupled events promoting memory retention, while another is observed to facilitate weakening of memory representations (Kim et al., 2019). Notably, this observed combination of synaptic strengthening and weakening of memory content has been proposed as a possible mechanistic framework for sleep-dependent memory consolidation (Puentes-Mestril and Aton, 2017). Thus, the mixing of different subtypes of coupled slow wave/spindle events during cycles of NREM sleep may impact key and distinct memory consolidation processes underlying selective synaptic remodeling and synaptic homeostasis. Further, given that changes in the depth of sleep between N2 and N3 sleep are associated with shifting compositions of SWA subtypes, these transitions between N2 and N3 sleep across each NREM cycle may create ebbs and flows of different synaptic remodeling activities.
Analyzing the composition of NREM sleep cycles provides a tool to advance our fundamental understanding of sleep-dependent memory consolidation and the proposed neuroprotective aspects of slow wave sleep. Studies of sleep and memory have reported incongruencies regarding slow wave sleep’s memory promoting effects, with some data indicating that slow wave sleep strengthens overnight memory consolidation, while others fail to reproduce these observations (Scullin, 2013; Cordi and Rasch, 2021). Notably, these incongruencies remain poorly explained in a more limited, homogenous, conceptual framework of SWA that does not account for the differences in the composition of slow wave sleep.
Properties of SWA are also proposed to be a possible biomarker of neurodegenerative disease, given that the proposed neuroanatomical components of SWA-associated memory consolidation activities are directly impacted by Alzheimer’s disease processes including beta amyloid and tau accumulation (Mander et al., 2015; Winer and Mander, 2018; Winer et al., 2019, 2020). Disruptions in the specific timing relationships between slow waves and spindles correlate with worsening cognitive performance among aging adults (Helfrich et al., 2018; Djonlagic et al., 2021), and the presence of Alzheimer’s disease pathology may underly some these timing abnormalities (Winer et al., 2019). Thus far, these studies have been limited to analyzing SWA as a mixture of different slow wave and spindle subtypes with unknown composition, and given that the timing aspects of slow wave and spindle coupling are observed to differ between early-fast and late-fast spindles (McConnell et al., 2021), a more advanced approach to capturing the composition of SWA may provide tools and insights to further delineate the connections between sleep and Alzheimer’s disease processes.
Our study was designed to appraise whether structure and composition of the SWA may be fundamentally different depending on the depth of NREM sleep and time within a sleep session, with an overarching goal to examine how assessment of SWA may vary depending on experimental conditions. We leveraged two existing and well-characterized sleep recording datasets to examine the structure and composition of SWA across the night and within cycles of NREM sleep. We further examined how the composition of SWA through cycles of NREM sleep changes with age, building on prior observations that aging is associated with overall greater late-fast composition of NREM sleep (McConnell et al., 2021). Here we specifically expanded on this work to examine the temporal composition of SWA at discrete time intervals, and we measured the relative composition of SWA slow wave/spindle subtypes across the night for each younger and older research participant.
In light of studies suggesting that SWA and the coupling of slow waves to spindles may be an early biomarker of cognitive aging and neurodegenerative disease (Mander et al., 2015, 2016; Winer et al., 2019, 2020), efforts to advance these insights and to translate the science to clinical use will need to account for the fundamental changes in SWA composition across cycles of NREM. Here we characterize and illustrate the distinctions in the composition of SWA, and we demonstrate the remarkable variance in the composition of SWA between young adults and aging adults.
Materials and Methods
All data were collected with written informed consent under protocols approved by a local institutional review board at each institution. Participant data was accessed from the DREAMS database (Devuyst, 2005) and from the Cleveland Family Study (CFS) via the National Sleep Research Repository (Redline et al., 1995, 1999; Dean et al., 2016; Zhang et al., 2018). The DREAMS database was originally created to study automated sleep staging techniques among overnight sleep recordings among individuals with no identifiable sleep pathologies including sleep apnea. Recordings from DREAMS are characterized as containing very few EEG artifacts and were included in this database specifically for their recording clarity. We selected the DREAMS database for our sub-analysis based on the quality of the recordings and the availability of frontal and central recording electrodes (these recordings sites are not available within the CFS dataset). Overnight polysomnography was utilized from all of the available 20 normal control participants available within the DREAMS database.
The CFS was originally conducted to study sleep apnea among a well-characterized cohort of 2,284 individuals and we selected this study for our secondary analysis due to the availability of high-quality overnight polysomnography from participants in young adult and older adult age ranges (distinct age ranges are not available within the DREAMS database). Data quality was originally rated by sleep staging experts per standardized study protocol, and we selected recordings for analysis to contain at least 95% visually-scorable EEG data. Individuals from CFS were selected from ages 18–35 for young adults [n = 30; median age 24.19 (21.78–25.48 IQR)] and 60–70 for older adults [n = 30; median age 65 (62.09–68.13 IQR)], given prior data indicating that these age ranges differ in slow wave/spindle coupling (McConnell et al., 2021). To minimize possible effects of sleep apnea, we selected 60 participants within the CFS by lowest oxygen desaturation apneas and hypopneas ≥3% oxygen desaturation or arousal per hour of sleep (AHI) scores within each age desired range [median AHI for young adults = 0.21 (0.15–0.50 IQR); median AHI for older adults = 4.51 (2.80–9.15 IQR); Table 1].
EEG Acquisition
For Dream EEG files, recordings were acquired at the André Vésale hospital (Montigny-le-Tilleul, Belgium) using a digital 32-channel polygraph (BrainnetTM System of MEDATEC, Brussels, Belgium). Channels used for analysis include FP1 referenced to A2 and CZ referenced to A1, both recorded at a sampling rate of 200 Hz. Staging was performed by a single expert at the sleep laboratory following AASM criteria (Berry et al., 2012). Similarly for CFS EEG files, laboratory-based polysomnography was performed using a 14-channel Compumedics E-Series System (Abbotsford, Australia). Recordings were created using gold cup electrodes at a sampling rate of 128 Hz (filtered low pass 0.16 Hz and high pass 105 Hz). Electrodes at C3 and C4 were referenced to positions A2 and A1, respectively. Sleep staging was performed by trained technicians at a centralized reading center at Brigham and Women’s Hospital following study protocols and AASM criteria (Berry et al., 2012).
Slow Wave Detection
Slow waves were identified via automated zero-crossing detection as described in McConnell et al. (2021). Briefly, signals for slow wave detection were utilized from electrodes FP1, CZ, and C3 from stages N2 and N3 of slow wave sleep by AASM criteria (Berry et al., 2012). Channels FP1 and CZ were selected for their relatively close positions to slow wave sources described in localization experiments, which demonstrated a frontal source near the left insula and a central source near the cingulate gyrus (Murphy et al., 2009). Notably, while these sources were utilized to guide electrode selection, the neuroanatomical source of EEG signals within channels FP1 and CZ cannot be determined given the limitations of source localization inherent to single channel EEG measures. Electrode C3 was selected as having been well-characterized in our prior analysis of slow wave and spindle coupling among 582 individuals from CFS (McConnell et al., 2021).
Epochs with un-scorable data were excluded from analysis. Automated management of high amplitude artifacts was accomplished via exclusion of EEG segments exceeding 900 μV after detrending data with sliding window of 3 s across raw data. EEG data was subsequently detrended and band-pass filtered in a forward and backward direction using a 6th-order Butterworth filter between 0.16–4 Hz. Zero crossings were identified to detect negative and positive half-waves, and slow wave events were identified when the half-wave pairs approximated a frequency range of 0.4–4 Hz. Minimal and maximal half-wave amplitudes were measured, and slow waves with both positive and negative maximum amplitudes in the top 50% of all waves were selected for subsequent coupling analysis. An upper threshold of ± 200 μV for zero crossing pairs was utilized to reduce misidentification of non-slow wave events. A further reduction of false identifications was accomplished by rejecting all zero crossing pairs with peak/trough amplitudes exceeding four standard deviations from the mean min/max zero crossing pair values for each subject, and visual inspection was performed to ensure consistent performance of artifact rejection, as detailed in McConnell et al. (2021).
Time-Frequency Analysis
Time-frequency wavelet plots of coupled spindles were created via established methods (Mölle et al., 2011; Muehlroth et al., 2019; McConnell et al., 2021). Troughs of each slow wave were centered in 5-s intervals of EEG data and matched to 5-s baseline intervals immediately preceding slow wave events (excluding baseline segments containing slow wave events). Morlet-wavelet transformation (8 cycles) was applied to the unfiltered EEG for slow wave and baseline segments between 4 Hz and 20 Hz in steps of 0.25 Hz. The mean of baseline regions was used to normalize the amplitude of the mean Morlet-wavelet transformation of all 5-s peri-slow wave regions.
Slow Wave Separation by Co-identification of Spindle Events and Timing Alignment
Spindle event identification was consistent with parameters of prior reports (Mölle et al., 2011; Staresina et al., 2015; McConnell et al., 2021). EEG data was detrended and bandpass filtered in a forward and backward direction using a 3rd-order Butterworth filter between 10 and 13.5 Hz for late-fast spindles, and between 14.5 and 18 Hz for early-fast spindles. Note that early-fast and late-fast spindles share frequency ranges near 14 Hz (McConnell et al., 2021), therefore the range of 13.5–14.5 Hz was excluded from spindle identification to reduce cross-identification of spindle event types. The upper spindle envelopes were calculated and an amplitude threshold of 75% percentile of the root mean squared value with a length window of 0.5–3.0 s was used to define spindle events. An absolute threshold of 40 microvolts in range was used to eliminate artifacts and only spindle envelopes within 8 standard deviations from the mean amplitude values were used in further analysis. Individual slow wave events were sorted by their co-localization with each spindle type in the time domain. Spindles occurring within 0.2–1.5 s from the trough of each slow wave were classified as a coupled event and sorted as late-fast or early-fast coupled slow wave events. Note that the time window adjacent to the slow wave troughs was buffered by 0.2 s after the trough to avoid misidentification of slow spindles, which share frequency components with late-fast spindles and occur immediately prior to the slow wave trough (McConnell et al., 2021). These parameters we selected to reduce mis-matching of spindles that overlap in timing and frequency components, although given the restrictive quality of these parameters, a portion of slow wave events will necessarily remain unclassified.
Statistical Analysis
All statistical modeling was performed with SAS v9.4 (SAS Institute Inc., Cary, NC, United States).
Cross-channel comparison of SWA event count composition: Given that slow wave source localization experiments report both left frontal and central sources of slow waves (Murphy et al., 2009), we analyzed whether the occurrence of slow wave/spindle subtypes differs between left frontal vs. central recording locations. Specifically, to compare the occurrence of late-fast and early-fast slow wave/spindle events, counts of each event type across channels (CZ and FP1) were analyzed with negative binomial generalized estimating equations (GEE).
Cross-channel comparison of SWA ROI spectral composition: The normalized spectral power within frontal and central channels were further compared within the ROIs corresponding to early-fast and late-fast slow wave/spindle events. Here, repeated measures regression models were used to analyze ROI data across channels (CZ and FP1).
Comparison of sleep stage-specific SWA composition: The SWA composition of early-fast spindle/slow wave and late-fast spindle/slow wave event counts and percentage of late-fast spindle/slow wave were measured in stages N2 and N3 sleep for each participant to quantify the effect of shifting SWA composition between cycles of NREM sleep. We performed two-sample, two-sided T-tests of each channel configuration used to detect slow wave and spindle events to analyze differences in event counts and percentage of late-fast spindle slow wave events between sleep stages (given differences in the number of slow wave/spindle events detected, calculations were weighted by the number of detected events for each participant to improve statistical reliability).
Comparison of SWA event count composition between age groups: Prior work has demonstrated differences in early fast and late fast slow wave/spindle EEG power (McConnell et al., 2021), and here we sought to measure changes in the composition of SWA slow wave/spindle subtypes that may correspond with these age-related EEG power changes. To measure the composition of SWA, the ratio of early fast and late fast slow wave/spindle coupling types were compared between 18–35-year-old patients and 60–70 patients with a two-sample, two-sided T-test (similar to prior analyses, calculations were weighted by the number of detected events for each participant to improve statistical reliability).
Comparison of SWA event percentage composition between age groups: As an additional means to measure differences in the composition of SWA among younger and older individuals, we examined regions of the overnight recordings that were heavily composed of late-fast slow wave/spindle subtypes by counting the number of time bins where the SWA composite was 80% or more of the late-fast slow wave/spindle subtype. The count of bins greater than 80% were compared between 18–35-year-old patients and 60–70-year-old patients with a negative binomial count model.
Results
Coupling of slow waves to spindles demonstrates an interspersed composite of distinct slow wave/spindle events (DREAMS dataset). To aid in conceptualization of the composite structure of SWA, we identified each subtype of slow wave/spindle pair within the filtered EEG timeseries data of an individual subject (DREAMS Subject 3, see Table 1). Signal processing steps are illustrated via event-centered plots of each slow wave and spindle subtype (Figures 1A, B). Note that coupled spindles are not readily visible in slow wave-centered event average plots due to jitter in the phase lock between slow wave troughs and spindles (Figure 1B). Segments of EEG data, each with decrementing time durations, are illustrated alongside superimposed markers for each identified slow wave/spindle event to visually demonstrate the intermixed amalgam of SWA (Figure 1C). In contrast to event-averaged slow wave/spindle pairs, the coupled spindles of each individual slow wave/spindle event are discernable within the raw EEG timeseries data (Figure 1C, lower panel).
Figure 1. (A) Early-fast and late-fast spindle/slow wave couples are illustrated and color-coded as red and blue, respectively. (B) Original EEG signals for slow waves (bandpass filtered 0.16–4 Hz), early-fast spindles (bandpass filtered 14–18 Hz), late-fast spindles (bandpass filtered 10–13.5 Hz) and coupled slow wave + spindle events (bandpass filtered 0.16–18 Hz) are demonstrated for a single subject (DREAMS Subject 3). (C) EEG time-series data (bandpass filtered 0.16–18 Hz) are from a single subject (DREAMS Subject 3) are graphed at three levels of time duration to demonstrate the intermixing of early-fast spindle/slow waves (red markers) and late-fast spindle/slow waves (blue markers).
Subtypes of slow wave/spindle couples are differentially represented in frontal vs. central electrode locations (DREAMS dataset). We hypothesized that frontal slow waves may demonstrate distinct coupling patterns, compared to central slow waves, and we examined slow wave coupling from channels FP1 (left frontal) and CZ (central) to test our hypothesis. Total numbers of slow waves identified were not significantly different between frontal vs. central recording sites (p > 0.727), although late-fast spindle counts were ∼14.8% higher in the frontal channel [95% CI: (10.9, 18.9%), p < 0.0001] compared to the central channel, and early-fast spindle counts were ∼11.7% higher in the central channel compared to the frontal channel [95% CI: (7.7, 15.8%), p = 0.0004; Table 2]. Event counts for slow waves coupled to late-fast vs. early-fast spindles similarly demonstrated a frontal vs. central preference, with late-fast/slow wave (LF-SW) coupled events ∼18.3% higher in the frontal channel [95% CI: (11.9, 25.1%), p = 0.0003], and early-fast/slow wave (EF-SW) coupled events ∼20.8% higher in the central channel [95% CI: (13.7, 28.4%), p = 0.0003]. Together, these data indicate that although the number of slow waves identified in channels FP1 and CZ is similar, the mix of different slow wave/spindle couples differs between frontal vs. central recording locations (Table 2).
In addition to these differences in event counts, we also observed differences in the averaged EEG power of coupled spindles between frontal and central sites. Congruent with event count data, the normalized, averaged EEG power is ∼6.5% higher for coupled late-fast spindles in the frontal region [95% CI: (3.4, 9.5%), p = 0.0002], while the normalized, averaged early-fast coupled spindle power is ∼6.2% higher in the central region [95% CI: (1.8, 10.6%), p = 0.0079; Figure 2 and Table 2].
Figure 2. Frontal recording site FP1 demonstrates higher late-fast spindle EEG power (A), while central recording site Cz demonstrates higher early-fast spindle EEG power (B). Blue and red arrows indicate the direction of EEG power differential between electrode recording sites. (C) An electrode position diagram is provided for reference to the FP1 and Cz recording site locations [FP1 for panel (A) and Cz for panel (B)].
The composite of SWA is formed by a mixture of slow wave/spindle coupled events, and their relative contribution to the SWA rises and falls across cycles of NREM sleep (DREAMS dataset). In order to better characterize the changing mixture of distinct slow wave/spindle events, we measured the instances of each slow wave/spindle event subtype across cycles of NREM sleep. Congruent with measures of coupled spindle power, there were distinctions between stages N2 and N3 in the composition of slow wave/spindle subtypes. At the level of the individual, each cycle of NREM sleep is observed to shift in the composition of SWA, and a distinction between production of LF-SWs vs. EF-SWs moves through the night in accordance with the with the depth of sleep (illustrated for a single participant in Figure 3).
Figure 3. (A) An illustrative example of SWA composition for a single participant aged 24 years old (DREAMS Subject 3). A color density spectral array provides a graphic illustration of EEG power for an entire overnight sleep recording of a single subject. (B) The hypnogram marking sleep stages for this subject’s overnight sleep recording is illustrated. (C) The percentage of early-fast spindle/slow waves (red bars) and late-fast spindle/slow waves (blue bars) in each time bin of 2 min is graphed to illustrate the composite mix of SWA across a single night of sleep.
This undulation of SWA composite through NREM cycles is most clearly observed when slow wave/spindle subtypes are identified with respect to their location propensity, i.e., frontal for LF-SWs and central for EF-SWs. Identification of both subtypes from the same channel location yielded similar results, although the shifts between SWA composition were less pronounced (illustrated in Figure 4).
Figure 4. An illustrative example comparing SWA composition as measured from different channel pairs for a single participant (DREAMS Subject 3). (A) A hypnogram marking sleep stages for a single participant’s overnight sleep recording is illustrated. (B) The percentage of early-fast spindle/slow waves (red bars) detected in channel CZ and late-fast spindle/slow waves (blue bars) detected in channel FP1 are graphed in time bins of 2 min. (C) The percentage of early-fast spindle/slow waves (red bars) detected in channel FP1 and late-fast spindle/slow waves (blue bars) detected in channel FP1 are graphed in time bins of 2 min. (D) The percentage of early-fast spindle/slow waves (red bars) detected in channel CZ and late-fast spindle/slow waves (blue bars) detected in channel CZ are graphed in time bins of 2 min.
To further illustrate the shifts in SWA composition across individual subjects’ NREM sleep cycles, eight additional participants from the DREAMS database are presented with event histograms marking normalized percentage of total stage N2 and N3 sleep for each individual (Figure 5).
Figure 5. Individual differences in the shifting composition through stages N2 and N3 SWA. Event histograms marking the percentage of early-fast spindle/slow waves detected in channel CZ (red bars) and late-fast spindle/slow waves detected in channel FP1 (blue bars) are graphed in time bins of 2 min [(A–H); DREAMS Subjects 1, 2, 4–9].
Quantification of SWA composition between stages N2 and N3 sleep confirms stage-specific composition shifts underly undulations across cycles of NREM sleep (DREAMS dataset). Building on the observation that SWA composition shifts across cycles of NREM sleep, we next measured the number of each slow wave/spindle subtype in each stage of sleep from all 20 participant recordings in the DREAMS database to calculate the percentage of LF-SWs vs. EF-SWs for N2 and N3 sleep. Our visual representations of SWA composition suggest that the undulations of SWA shift toward greater LF-SW subtypes in N3 sleep, while EF-SW subtypes are more frequent in stage N2, congruent with prior observations that demonstrate a distinction in the late-fast spindle and early-fast spindle EEG power between stages N2 and N3 sleep (McConnell et al., 2021). Quantifying the percentage of LF-SW vs. EF-SW events from stages N2 and N3 sleep was congruent with this observation. Measuring LF-SW events from channel FP1 and EF-SW events from channel CZ, 59.93% [95% CI: (55.95, 63.92%)] of slow wave/spindle events during N3 sleep were LF-SW subtype, while only 49.25% [95% CI: (46.12, 52.39%)] of N2 sleep were the LF-SW subtype [FP1-CZ pair, estimated difference = 10.68 percent units, 95% CI: (6.78, 14.57), p < 0.0001; Table 3]. Measuring both LF-SW events and EF-SW events from channel CZ, 53.15% [95% CI: (48.36, 57.95%)] of slow wave/spindle events during N3 sleep were LF-SW subtype, while only 48.08% [95% CI: (44.70, 51.46%)] of N2 sleep were the LF-SW subtype [CZ-CZ pair, estimated difference = 5.08 percent units, 95% CI: (1.33, 8.82), p = 0.0114; Table 3]. Finally, measuring both LF-SW events and EF-SW events from channel FP1, 64.99% [95% CI: (60.88, 69.11%)] of slow wave/spindle events during N3 sleep were LF-SW subtype, while only 53.04% [95% CI: (50.04, 56.03%)] of N2 sleep were the LF-SW subtype [FP1-FP1 pair, estimated difference = 11.96 percent units, 95% CI: (7.66, 16.26), p < 0.0001; Table 3]. We further performed a comparative analysis of the differences in SWA composition between these channel pairs (FP1-CZ, CZ-CZ, and FP1-FP1). Additional statistical comparisons between channel pairs (FP1-CZ, CZ-CZ, and FP1-CZ) were performed for percentage of LF-SW events in both stages N2 and N3 sleep. Each permutation of channel comparison (i.e., FP1-CZ to CZ-CZ and FP1-FP1; CZ-CZ to FP1-CZ and FP1-FP1; FP1-FP1 to FP1-CZ and CZ-CZ) were statically significant (p < 0.001) for each comparison, except for LF-SW percentage in FP1-CZ vs. CZ-CZ (p = 0.06; Table 3).
Table 3. Stage N3- and N2-specific slow wave and spindle event counts and percentage of late-fast slow wave/spindle pairs from each channel pairing among the DREAMS subjects.
The composition of SWA among older adults favors a higher proportion of late-fast coupled slow wave/spindle events, and higher percentages of late-fast spindles during peak times (CFS dataset). To examine how the composition of SWA differs with age, we compared the slow wave/spindle subtype structure of NREM cycles from channel C3 between young adults and older adults from the CFS dataset. Overall, the mean proportion of LF-SW coupled events across the entire night is 57.1% [95% CI: (54.9, 59.3%)] in older adults, compared to 53.7% [95% CI: (52.0, 55.4%)] in younger adults [estimated difference = 3.43 percent units, 95% CI: (−0.7, −6.2%), p = 0.0146; Table 4]. Further, the peaks of high LF-SW event SWA composition frequently demonstrate a substantially higher amount of late-fast spindle coupling, when compared to the highest peaks of late-fast SWA among young adults (illustrated in Figure 6). To quantify this observation, a threshold was set to 80% normalized LF-SW events per 1% of the total normalized NREM time for each subject, and the number of bins exceeding this threshold was quantified among the young adults and older adults (Figure 6). Older adults averaged 10.9 [95% CI: (8.5, 14.0)] bins greater than or equal to 80% LF-SW events, while young adults averaged 2.9 [95% CI: (2.0, 4.2)] of these high LF-SW bins (older vs. younger ratio = 3.76, 95% CI: (2.41, 5.87), p < 0.0001; Table 4).
Table 4. Slow wave and spindle event counts, and mean number of bins exceeding a threshold of 80% composition of the late-fast subtype, for younger and older adults from the Cleveland family study.
Figure 6. A comparison of the SWA composition among younger and older adults. (A) A younger participant (CFS 800092), age 21-year-old, and (B) an older participant (CFS 800099), age 61-year-old, are presented as an illustrative example of SWA thresholding. The percentage of late-fast spindle/slow wave events (blue bars) and early-fast spindle/slow wave events (red bars) are assessed for bins with 80% or greater late-fast spindle/slow wave composition. (C) The number of bins at or exceeding 80% late-fast spindle/slow wave composition for each participant is quantified, and the younger adults (green-stripes; n = 30) are compared to older adults (yellow-checkered; n = 30) at the group level.
Discussion
Our analyses demonstrate that SWA can be decomposed into distinct subtypes based on coupling of slow wave and spindle pairs. Results from the study suggest that each subtype of slow wave/spindle event demonstrates unique neurophysiological properties that predict their rate of occurrence, including electrode location bias and relative depth of sleep. This work also progresses and expands our prior observations with regard to coupled spindle EEG power differences between stages N2 and N3 sleep (McConnell et al., 2021). Here we specifically demonstrate that the number of each distinct slow wave/spindle event type shifts the composition of SWA across cycles of NREM sleep via incorporation of whole-night event histograms. In addition, the temporal resolution provided by the event histograms suggests that age-related changes we previously reported in late-fast spindle coupling create windows of time during NREM sleep when the composition of SWA is strikingly biased toward the deep sleep-associated late-fast spindle/slow wave event subtype. Collectively our findings build upon and may provide increased clarity on prior work that demonstrates slow waves have anterior and central sources (Murphy et al., 2009), and that SWA differs in cortical propagation (Bernardi et al., 2018) and homeostatic properties (Werth et al., 1996) along the anterior-to-posterior axis.
Our whole-night illustration of SWA composition further demonstrate a striking rise and fall of slow wave/spindle events across cycles of NREM sleep, suggesting the possibility of competing or complementary processes that occur in a choreographed rise and fall of SWA composition states. In addition, our results demonstrate that the event composition of SWA is quantitatively distinct between stages N2 and N3 sleep, with deeper sleep favoring a higher percentage of late-fast slow wave/spindle events. This is consistent with our prior observations that demonstrate the normalized EEG power differs in the late fast vs. early fast ROI regions within slow wave/spindle events (McConnell et al., 2021), and suggests that the differences in EEG power may be, at least in part, attributable to differences in the number of each slow wave/spindle subtype that creates each sleep stage’s SWA composite.
Also congruent with our previous analyses of slow wave/spindle EEG power, (McConnell et al., 2021) we observed that aging is associated with a relative loss of early-fast slow wave/spindle events, as the late-fast coupled slow waves form a greater contribution to the mix of SWA and increasingly dominate the SWA event composition. Here we further demonstrate that extremes of SWA event composition are common among aging adults, and late-fast coupled slow wave/spindle pairs dominate these segments of skewed SWA composition. Notably, these results are seemingly counterintuitive within the context of a well-described reduction in stage N3 sleep observed among aging adults (Muehlroth, 2019); however, we and others have previously discussed the shortcomings of applying sleep staging rules to EEG data from older adults, which consistently results in misclassification of stage N3 sleep (Berry et al., 2012; Muehlroth, 2019; McConnell et al., 2021). Specifically, stage N3 sleep is increasingly misclassified as stage N2 sleep as age-related decreased in EEG amplitude progressively results in slow wave amplitudes that fall below criterion for slow wave classification according to AASM sleep staging rules (Berry et al., 2012). Thus, examining the composition of SWA may result in better characterization of NREM sleep state among aging adults, although additional studies will be required to determine how SWA composition may fit within existing sleep staging conceptual framework.
Understanding that SWA is not uniform, but instead shifts in composition across cycles of NREM, provides key insights to inform interpretations of prior studies and better design future experiments. One such scenario may be studies of slow wave coupling and memory performance, which have reported seemingly conflicting relationships (Scullin, 2013; Cordi and Rasch, 2021). Key differences in experimental design, particularly regarding the captured and analyzed depth of slow wave sleep, may be reflected in the composition of SWA among study participants and may account for discrepancies in outcome measures. Further, experimental manipulation of SWA may preferentially impact one type of slow wave/spindle event, creating an imbalance in the composite of SWA across portions of NREM sleep. Accounting for these changes in slow wave physiology may be critical to understanding the mechanisms underlying SWA manipulation.
While the function of early-fast vs. late-fast coupling with slow waves is yet to be determined, additional mechanistic insights may come from intracranial recording techniques. Human experiments among patients undergoing epilepsy surgery have described distinctions between anterior and posterior hippocampal slow wave/spindle coupling (Jiang et al., 2019a,b). Given that the anterior and posterior hippocampi are proposed to utilize distinct neuroanatomical connections for processing different forms of memory content (Poppenk et al., 2013; Strange et al., 2014), slow wave/spindle coupled memory processing events may also reflect different forms of memory processing in frontal vs. parietal cortical regions. Indeed, data from experimental manipulation of SWA in an animal model demonstrate two distinct forms of slow wave/spindle coupling, and their functions may be opposing to one another in synaptic regulation and memory processing functions (Kim et al., 2019).
Our analyses delineate key neurophysiological properties of late-fast and early-fast spindle coupling, although important limitations remain in characterizing the composition of SWA and the function of each distinct type of slow wave and spindle event. Signal processing experiments have utilized features of slow waves including the waveform characteristics and recording locations to describe multiple subtypes of slow waves as anterior vs. central (Murphy et al., 2009), global vs. local (Bernardi et al., 2018), and fast switching vs. slow switching (trough-to-peak transition frequency) (Maude et al., 2021). In addition, the coupling of slow waves with spindle and hippocampal ripples is reported to differ between anterior vs. posterior hippocampal locations, further suggesting that a division of coupled slow waves may follow an similar anterior vs. posterior pattern (Jiang et al., 2019a,b). Future experiments will be required to determine whether spindle coupling indeed separates slow wave events into these previously described slow wave subtypes. In addition, the use of coupling to identify slow waves has inherent limitations in the ability to label slow wave events, including an overlapping frequency component that precludes characterization of many spindles in this shared frequency range at approximately 14 Hz (McConnell et al., 2021). Further, many slow waves do not clearly demonstrate a coupled spindle on surface EEG, and there remain numerous undifferentiated slow waves that may represent additional subtypes of slow wave events. The incorporation of additional slow wave and spindle features may provide the next steps in better decomposing all SWA in future studies.
Additional limitations are inherent in the use of single channel surface EEG to describe neuroanatomical localization. While source localization from prior reports was used to guide electrode selection, neuroanatomical source localization of our signals is not possible due to the well-described limitations of the inverse problem (Fender, 1987). Further, while we previously described no significant differences in left vs. right laterality of slow wave/spindle coupling within channels C3 and C4 (McConnell et al., 2021), additional studies will be required to determine whether left/right laterality differences in SWA composition can be detected between hemispheres, and these analyses would be best suited for sleep recording data sets that can provide greater source localization.
Aging has been linked to changes in slow wave/spindle coupling (Helfrich et al., 2018; Muehlroth et al., 2019), and indeed, the continuum of the adult human lifespan demonstrates a shift in the composition of coupled spindle EEG power toward the late-fast subtype (McConnell et al., 2021). The results of this study support a conceptual model of a drift in the types of slow wave/spindle events that older adults produce, and a striking increase in the number of time periods that cycle disproportionately into late-fast spindle domination of SWA. Given the interest in spindle coupling as a biomarker of cognitive aging and risk of neurodegenerative disease (Mander et al., 2015, 2016; Winer et al., 2019, 2020), the potential impact of this late-fast spindle coupling bias among aging adults may provide mechanistic insights and facilitate the development of a predictive functional biomarker. Future studies may focus on the neuroanatomical basis of late-fast and early fast spindle to better understand the potential mechanisms that may produce a shift in composition of sleep.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://sleepdata.org/datasets/cfs and https://doi.org/10.5281/zenodo.2650142.
Ethics Statement
The studies involving human participants were reviewed and approved by André Vésale Hospital (Montigny-le-Tilleul, Belgium) for DREAM overnight EEG files, and Cleveland Family Study EEG files were acquired under Institutional Review Boards Multiple Hospitals (University Hospitals, the Louis Stokes Cleveland Department of Veterans Affairs Medical Center, and MetroHealth Medical Center). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
BM, EK, VK, SS, AR, RP, and BB designed the study. BM, SS, LM, and RP analyzed the data. BM and BB wrote the manuscript. EK, VK, and AR gave conceptual advice. All authors participated in editing and revision of the manuscript.
Funding
This work was supported in part by funding from the National Institutes of Health (R01AG058772). We are grateful for the open access to data provided by the DREAMS database (10.5281/zenodo.2650142) and the Cleveland Family Study (CFS) [supported by grants from the National Institutes of Health (HL46380, M01 RR00080- 39, T32-HL07567, and RO1-46380)], and the National Sleep Research Resource [supported by the National Heart, Lung, and Blood Institute (R24 HL114473 and RFP 75N92019R002)].
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.
References
Achermann, P., and Borbely, A. (1997). Low-frequency (< 1 Hz) oscillations in the human sleep electroencephalogram. Neuroscience 81, 213–222. doi: 10.1016/s0306-4522(97)00186-3
Amzica, F., and Steriade, M. (1998). Electrophysiological correlates of sleep delta waves. Electroencephal. Clin. Neurophysiol. 107, 69–83. doi: 10.1016/s0013-4694(98)00051-0
Backhaus, J., Born, J., Hoeckesfeld, R., Fokuhl, S., Hohagen, F., and Junghanns, K. (2007). Midlife decline in declarative memory consolidation is correlated with a decline in slow wave sleep. Learn. Memory 14, 336–341. doi: 10.1101/lm.470507
Bernardi, G., Siclari, F., Handjaras, G., Riedner, B. A., and Tononi, G. (2018). Local and widespread slow waves in stable NREM sleep: evidence for distinct regulation mechanisms. Front. Hum. Neurosci. 12:248. doi: 10.3389/fnhum.2018.00248
Berry, R. B., Brooks, R., Gamaldo, C. E., Harding, S. M., Marcus, C., and Vaughn, B. (2012). The AASM manual for the scoring of sleep and associated events. Rules, Terminology and Technical Specifications, Darien, Illinois. Darien: American Academy of Sleep Medicine.
Born, J. (2010). Slow-wave sleep and the consolidation of long-term memory. World J. Biol. Psychiatry 11, 16–21. doi: 10.3109/15622971003637637
Cordi, M. J., and Rasch, B. (2021). No evidence for intra-individual correlations between sleep-mediated declarative memory consolidation and slow-wave sleep. Sleep 44:zsab034. doi: 10.1093/sleep/zsab034
Cox, R., Rüber, T., Staresina, B. P., and Fell, J. (2019). Heterogeneous profiles of coupled sleep oscillations in human hippocampus. NeuroImage 202:116178. doi: 10.1016/j.neuroimage.2019.116178
Dean, D. A., Goldberger, A. L., Mueller, R., Kim, M., Rueschman, M., Mobley, D., et al. (2016). Scaling up scientific discovery in sleep medicine: the National Sleep Research Resource. Sleep 39, 1151–1164. doi: 10.5665/sleep.5774
Devuyst, S. (2005). The DREAMS Databases and Assessment Algorithm [Data set]. Zenodo. doi: 10.5281/zenodo.2650142
Djonlagic, I., Mariani, S., Fitzpatrick, A. L., Van Der Klei, V. M. G. T. H., Johnson, D. A., Wood, A. C., et al. (2021). Macro and micro sleep architecture and cognitive performance in older adults. Nat. Hum. Behav. 5, 123–145.
Fender, D. (1987). “Source localization of brain electrical activity,” in Methods of Analysis of Brain Electrical and Magnetic Signals: Handbook of Electroencephalography and Clinical Neurophysiology, eds A. Gevins and A. Remond, (Amsterdam: Elsevier), 355–403.
Helfrich, R. F., Mander, B. A., Jagust, W. J., Knight, R. T., and Walker, M. P. (2018). Old brains come uncoupled in sleep: slow wave-spindle synchrony, brain atrophy, and forgetting. Neuron 97, 221–230.e4. doi: 10.1016/j.neuron.2017.11.020
Jiang, X., Gonzalez-Martinez, J., and Halgren, E. (2019a). Coordination of Human Hippocampal Sharpwave Ripples during NREM Sleep with Cortical Theta Bursts, Spindles, Downstates, and Upstates. J. Neurosci. 39, 8744–8761. doi: 10.1523/JNEUROSCI.2857-18.2019
Jiang, X., Gonzalez-Martinez, J., and Halgren, E. (2019b). Posterior Hippocampal Spindle Ripples Co-occur with Neocortical Theta Bursts and Downstates-Upstates, and Phase-Lock with Parietal Spindles during NREM Sleep in Humans. J. Neurosci. 39, 8949–8968. doi: 10.1523/JNEUROSCI.2858-18.2019
Kim, J., Gulati, T., and Ganguly, K. (2019). Competing Roles of Slow Oscillations and Delta Waves in Memory Consolidation versus Forgetting. Cell 179, 514–526.e13. doi: 10.1016/j.cell.2019.08.040
Latchoumane, C.-F. V., Ngo, H.-V. V., Born, J., and Shin, H.-S. (2017). Thalamic spindles promote memory formation during sleep through triple phase-locking of cortical, thalamic, and hippocampal rhythms. Neuron 95, 424–435.e6. doi: 10.1016/j.neuron.2017.06.025
Lee, Y. F., Gerashchenko, D., Timofeev, I., Bacskai, B. J., and Kastanenka, K. V. (2020). Slow wave sleep is a promising intervention target for Alzheimer’s disease. Front. Neurosci. 14:705. doi: 10.3389/fnins.2020.00705
Mander, B. A., Marks, S. M., Vogel, J. W., Rao, V., Lu, B., Saletin, J. M., et al. (2015). β-amyloid disrupts human NREM slow waves and related hippocampus-dependent memory consolidation. Nat. Neurosci. 18, 1051–1057. doi: 10.1038/nn.4035
Mander, B. A., Winer, J. R., Jagust, W. J., and Walker, M. P. (2016). Sleep: a novel mechanistic pathway, biomarker, and treatment target in the pathology of Alzheimer’s disease? Trends Neurosci. 39, 552–566. doi: 10.1016/j.tins.2016.05.002
Maude, B., Jean-Marc, L., Pierre-Olivier, G., Lafrenière, A., Dubé, J., Gosselin, N., et al. (2021). Sleeping at the switch. eLife 10:e64337. doi: 10.7554/eLife.64337
McConnell, B. V., Kronberg, E., Teale, P. D., Sillau, S. H., Fishback, G. M., Kaplan, R. I., et al. (2021). The aging slow wave: a shifting amalgam of distinct slow wave and spindle coupling subtypes define slow wave sleep across the human lifespan. Sleep 44:zsab125. doi: 10.1093/sleep/zsab125
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
Muehlroth, B. E. (2019). Werkle-Bergner M. Understanding the interplay of sleep and aging: methodological challenges. Psychophysiology 0:e13523. doi: 10.1111/psyp.13523
Muehlroth, B. E., Sander, M. C., Fandakova, Y., Grandy, T. H., Rasch, B., Shing, Y. L., et al. (2019). Precise slow oscillation–spindle Coupling promotes Memory Consolidation in Younger and older Adults. Sci. Rep. 9:1940. doi: 10.1038/s41598-018-36557-z
Murphy, M., Riedner, B. A., Huber, R., Massimini, M., Ferrarelli, F., and Tononi, G. (2009). Source modeling sleep slow waves. Proc. Natl. Acad. Sci. 106, 1608–1613. doi: 10.1073/pnas.0807933106
Poppenk, J., Evensmoen, H. R., Moscovitch, M., and Nadel, L. (2013). Long-axis specialization of the human hippocampus. Trends Cogn. Sci. 17, 230–240.
Puentes-Mestril, C., and Aton, S. J. (2017). Linking network activity to synaptic plasticity during sleep: hypotheses and recent data. Front. Neural Circuits 11:61. doi: 10.3389/fncir.2017.00061
Redline, S., Tishler, P. V., Schluchter, M., Aylor, J., Clark, K., and Graham, G. (1999). Risk factors for sleep-disordered breathing in children: associations with obesity, race, and respiratory problems. Am. J. Respirat. Critical Care Med. 159, 1527–1532. doi: 10.1164/ajrccm.159.5.9809079
Redline, S., Tishler, P. V., Tosteson, T. D., Williamson, J., Kump, K., Browner, I., et al. (1995). The familial aggregation of obstructive sleep apnea. Am. J. Respirat. Critical Care Med. 151, 682–687.
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
Scullin, M. K. (2013). Sleep, memory, and aging: the link between slow-wave sleep and episodic memory changes from younger to older adults. Psychol. Aging 28, 105–114. doi: 10.1037/a0028830
Staresina, B. P., Bergmann, T. O., Bonnefond, M., van der Meij, R., 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
Strange, B. A., Witter, M. P., Lein, E. S., and Moser, E. I. (2014). Functional organization of the hippocampal longitudinal axis. Nat. Rev. Neurosci. 15, 655–669.
Varela, C., and Wilson, M. A. (2020). mPFC spindle cycles organize sparse thalamic activation and recently active CA1 cells during non-REM sleep. Elife 9:e48881. doi: 10.7554/eLife.48881
Walker, M. P. (2009). The role of slow wave sleep in memory processing. J. Clin. Sleep Med. 5, S20–S26.
Werth, E., Achermann, P., and Borbély, A. A. (1996). Brain topography of the human sleep EEG: antero-posterior shifts of spectral power. Neuroreport 8, 123–127. doi: 10.1097/00001756-199612200-00025
Wilson, M. A., and McNaughton, B. L. (1994). Reactivation of hippocampal ensemble memories during sleep. Science 265, 676–679. doi: 10.1126/science.8036517
Winer, J. R., and Mander, B. A. (2018). Waking up to the importance of sleep in the pathogenesis of Alzheimer disease. JAMA Neurol. 75, 654–656. doi: 10.1001/jamaneurol.2018.0005
Winer, J. R., Mander, B. A., Helfrich, R. F., Maass, A., Harrison, T. M., Baker, S. L., et al. (2019). Sleep as a potential biomarker of tau and β-amyloid burden in the human brain. J. Neurosci. 39, 6315–6324.
Winer, J. R., Mander, B. A., Kumar, S., Reed, M., Baker, S. L., Jagust, W. J., et al. (2020). Sleep disturbance forecasts β-amyloid accumulation across subsequent years. Curr. Biol. 30, 4291–4298.e3. doi: 10.1016/j.cub.2020.08.017
Wunderlin, M., Züst, M. A., Fehér, K. D., Klöppel, S., and Nissen, C. (2020). The role of slow wave sleep in the development of dementia and its potential for preventative interventions. Psychiatry Res. 306:111178. doi: 10.1016/j.pscychresns.2020.111178
Keywords: slow wave, sleep spindle, memory, coupling, biomarker, EEG
Citation: McConnell BV, Kronberg E, Medenblik LM, Kheyfets VO, Ramos AR, Sillau SH, Pulver RL and Bettcher BM (2022) The Rise and Fall of Slow Wave Tides: Vacillations in Coupled Slow Wave/Spindle Pairing Shift the Composition of Slow Wave Activity in Accordance With Depth of Sleep. Front. Neurosci. 16:915934. doi: 10.3389/fnins.2022.915934
Received: 08 April 2022; Accepted: 03 June 2022;
Published: 23 June 2022.
Edited by:
Ritchie Edward Brown, United States Department of Veterans Affairs, United StatesReviewed by:
Chiara Prosperetti, Neurocenter of Southern Switzerland, SwitzerlandFeng Han, The Pennsylvania State University (PSU), United States
Copyright © 2022 McConnell, Kronberg, Medenblik, Kheyfets, Ramos, Sillau, Pulver and Bettcher. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brice V. McConnell, brice.mcconnell@cuanschutz.edu