Timely coupling of sleep spindles and slow waves linked to early amyloid-β burden and predicts memory decline
Abstract
Sleep alteration is a hallmark of ageing and emerges as a risk factor for Alzheimer’s disease (AD). While the fine-tuned coalescence of sleep microstructure elements may influence age-related cognitive trajectories, its association with AD processes is not fully established. Here, we investigated whether the coupling of spindles and slow waves (SW) is associated with early amyloid-β (Aβ) brain burden, a hallmark of AD neuropathology, and cognitive change over 2 years in 100 healthy individuals in late-midlife (50–70 years; 68 women). We found that, in contrast to other sleep metrics, earlier occurrence of spindles on slow-depolarisation SW is associated with higher medial prefrontal cortex Aβ burden (p=0.014, r²β*=0.06) and is predictive of greater longitudinal memory decline in a large subsample (p=0.032, r²β*=0.07, N=66). These findings unravel early links between sleep, AD-related processes, and cognition and suggest that altered coupling of sleep microstructure elements, key to its mnesic function, contributes to poorer brain and cognitive trajectories in ageing.
Editor's evaluation
This paper is of interest to neuroscientists studying sleep, memory, and neurodegeneration. The authors found that an altered pattern of brain wave during NREM sleep, changes in the coupling of spindles and slow waves, correlates with amyloid burden and predicts memory decline over time in healthy older individuals. The results suggest that sleep brain waves may be a useful tool in identifying older adults at risk for future cognitive impairment in the earliest stage.
https://doi.org/10.7554/eLife.78191.sa0Introduction
Alterations in sleep quality are typical of the ageing process with a more fragmented and less intense (or shallower) sleep detected as early as the fifth decade of life (Carrier et al., 2011). Beyond healthy ageing, alterations in sleep are predictive of the risk of developing Alzheimer’s disease (AD) over the next 5–10 years (Lim et al., 2013; Musiek and Holtzman, 2016). Similarly, sleep disorders such as insomnia and obstructive sleep apnoea syndrome are associated with increased odds for AD diagnosis (Baril et al., 2018; Elias et al., 2018). Brain burdens of aggregated amyloid-β (Aβ) and tau proteins, hallmarks of AD pathophysiology, have been linked with a reduced sleep intensity, as indexed by the overall production of slow waves (SW) during sleep, but also to worse objective sleep efficiency and subjective sleep quality, in healthy and cognitively normal older individuals aged >70 years (Mander et al., 2015; Ju et al., 2017; Lucey et al., 2019; Kang et al., 2009; Ju et al., 2013; Sprecher et al., 2017). Sleep alteration may, in turn, contribute to the aggregation of Aβ and tau proteins: experimental sleep deprivation and sleep fragmentation (disturbance in the production of SW during sleep) lead to increased concentration of Aβ in the CSF, both in animal models and in healthy human populations (Ju et al., 2017; Kang et al., 2009; Lucey et al., 2018; Ooms et al., 2014). Overall, a bidirectional detrimental relationship between sleep quality and the neuropathology of AD is emerging in the literature. Sleep may therefore constitute a modifiable risk factor which one could act upon to prevent or delay the neuropathological processes associated to AD and favour successful cognitive trajectories over the lifespan (Van Egroo et al., 2019b; Ju et al., 2014; Wang and Holtzman, 2020). Hitherto, however, sleep is not yet widely recognised as an independent risk factor for AD, and the mechanistic associations between sleep and early AD neuropathology are not yet fully established.
Sleep microstructure elements, such as sleep spindles and SW, are essential correlates of the cognitive functions of sleep, as higher density of both elements during post-learning sleep has been linked to a better overnight memory consolidation (Ulrich, 2016; Gais et al., 2002; Schmidt et al., 2006; Mednick et al., 2013; Miyamoto et al., 2017). Furthermore, the slower portion of SW activity (i.e. a power measure combining the density and amplitude of sleep SW over sleep cycles in the 0.6–1 Hz band) was reported to modulate the regression between the burden of Aβ over the prefrontal cortex and a lower overnight memory consolidation in cognitively normal older individuals (Mander et al., 2015). The fine-tuned coupling of spindles and SW has further been reported to be altered in ageing, with an earlier occurrence of the spindle relative to the SW depolarisation phase in the older compared to younger individuals, and to predict overnight memory retention (Helfrich et al., 2018). However, whether this change in the spindle-SW coupling in ageing is associated to AD-pathological processes and to individual cognitive trajectories is not fully established. A study in 31 individuals, aged around 75 years, found a link between the brain deposit of tau protein and the coupling of spindles and SW (Winer, 2019). By contrast, another research failed to find a link between this coupling and Aβ brain burden (Winer et al., 2020). Here, we argue that, on top of potential statistical power issues, the difficulty to detect this link may be due to the fact that the assessments were carried out late over the lifespan (i.e. >70 years), when subtle associations may be masked by concurrent brain alterations, and by the heterogeneity of sleep SW.
The low frequency oscillations of the electroencephalography (EEG) during non-Rapid Eye Movement (NREM) sleep have been divided in slow oscillations (≤1 Hz) and delta waves (1–4 Hz) for decades in humans, notably based on the theoretical framework of the generation of SWs (Achermann and Borbély, 1997; Lee et al., 2004; Hubbard et al., 2020). As one of the main features of the SW resides in their transition from down- to up-state, reflecting synchronised depolarisation, a recent work proposed the transition frequency of the down-to-up state as an objective way of distinguishing between slow and fast switcher SWs (Bouchard et al., 2021). Compared to young adults, on top of exhibiting a typical overall lower density of SW, older individuals reportedly exhibit higher probabilities of producing slow as compared to fast switcher SWs, providing important insights into age-related changes in sleep microstructure.
Investigating the coupling of spindles and SW, appropriately split between the slow and fast switcher SW, in late middle-aged healthy adults may be the best approach to gain insight into the biology underlying the early relationship between sleep and AD-related processes. In a longitudinal study, we therefore tested whether the coupling of spindles with the slow and the fast switcher SWs is differently associated with the early brain burden of Aβ in a large sample (N=100) of healthy and cognitively normal individuals in their late midlife (50–70 years). We further explored whether spindle coupling onto slow and fast switcher SWs would be associated with the change in memory performance at 2 years in a large subsample (N=66). We recorded habitual sleep in these individuals devoid of sleep disorders of both sexes (59.5±5 years; 68 women) under EEG and extracted the density and coupling of the spindles and fast and slow switcher SWs over frontal derivations. The burden of Aβ was assessed using PET tracers ([18F]Flutemetamol, N=97; [18F]Florbetapir, N=3) over the medial prefrontal areas, known to be an early site for Aβ deposits and the most important generator of SWs during sleep (Mander et al., 2015; Dang-Vu et al., 2005; Dang-Vu et al., 2010; Saletin et al., 2013). Performance to the Mnemonic Similarity Task (MST), a memory task highly sensitive to early signs of cognitive decline (Stark et al., 2013; Marks et al., 2017) was assessed in all participants concomitantly to EEG and PET measurements as well as at follow-up, 2 years later (N=66). We hypothesised that our large sample of individuals, positioned relatively early in the ageing process, would allow to detect subtle differential associations between the impaired fine-tuned coupling of spindles and slow and fast switcher SWs, and both the early Aβ burden and memory performance decline over 2 years.
Results
Spindle onset on slow switcher SWs is linked to prefrontal Aβ burden
Following SW and spindle detections, we segregated SWs into those showing a slow or a fast down-to-up state transition, i.e., slow and fast switcher SWs (Figure 1a). Of the 341,836 detected slow switcher SWs, 75,910 co-occurred with a spindle (22%); while of the 78,235 fast switcher SWs, 26,912 (34%) were found to co-occur with a spindle. Regarding spindles, 563,928 spindles were detected over all the recordings, of which 102,822 were coupled to a SW (18%), 75,910 to slow switcher SWs (13%), and 26,912 to fast switcher SWs (5%). Table 1 gathers average characteristics of SW types and spindles in addition the demographic characteristics of our sample.
Our primary analysis tested whether the phase of coupling of spindles onto SWs, as the dependent variable, was differently associated with the accumulation of Aβ protein over the medial prefrontal cortex (MPFC) (Figure 1b) for slow and fast switcher SWs, controlling for age, sex, and total sleep time (TST) as well as for SW type overall duration, as it significantly differed across SW types and could condition the time available for the co-occurrence of a spindle and an SW (Table 1). The generalised linear mixed model (GLMM) yielded a significant difference between SW types (main effect of SW type: F1179.6=34.97, p<0.0001; r²β*=0.16) and no direct association with the burden of Aβ in the MPFC (main effect of Aβ burden: F198.49=0.09, p=0.76), while, importantly, the interaction between the burden of Aβ in the MPFC and the type of SWs was significant (Aβ burden by SW type interaction: F196=7.05, p=0.009; r²β*=0.07), controlling for all other covariates (age: F197.3=0.02, p=0.9; sex: F199.1=.75, p=0.4; TST: F199.1=0.75, p=0.4, SW duration: F1192.6=8.73, p=0.0035; r²β*=0.04). Post-hoc tests indicated that the link between the coupling of the spindles onto the SWs and the MPFC Aβ burden was significant for the slow switcher type (t149.8=−2.00, p=0.047) and not for the fast switcher type (t149.8=0.56, p=0.57) (Figure 1c–d). This primary result is in line with our initial hypothesis. In the following lines, we report additional statistical tests that were computed for a full understanding of this outcome.
We assessed whether spindles showed a preferential phase of anchoring with both slow and fast switcher SWs. Qualitative appreciation of the distributions suggests that there is no preferential phase of anchoring of spindles onto the fast switcher SWs, while spindle initiation onto slow switcher SWs would show a clear preferred phase (Figure 1e). Watson U² tests indicate, however, that the phase of anchoring onto slow and fast SWs is non-uniformly distributed (see Materials and methods; slow switcher SWs: U²=904.29, p<0.001; fast switcher SWs: U²=136.76, p<0.001), i.e., they both show some phase preference. Importantly, further statistical analysis with Watson’s U² test showed that the distribution of spindles anchoring phase was significantly different between slow and fast switcher SWs (U²=71.143, p<0.001). This finding reinforces the idea that slow and fast switcher SWs constitute distinct realisations of NREM oscillations that are differently associated with brain aggregation of Aβ during the ageing process. This has likely contributed to previous failures to detect links between the coupling of spindles and SWs and the deposit of Aβ. In fact, when testing the association between the coupling of spindles and SWs, irrespective of the type of SWs, and PET Aβ burden over the MPFC, the statistical analysis only yields a weak negative association between the phase of the coupling of the spindles onto the SWs and Aβ burden (main effect of Aβ uptake: F196=3.96, p=0.049, r²β*=0.04; main effect of sex: F196=4.33, p=0.04, r²β*=0.04; main effect of age: F1,96=0.05, p=0.83), which could arguably go undetected in a smaller or different sample.
We further assessed the specificity of our main finding and found that the density of slow switcher SWs (main effect of Aβ PET uptake: F196=0.14, p=0.71) or of spindles (main effect of Aβ PET uptake: F196=0.09, p=0.76) was not associated with the MPFC Aβ burden (Figure 1f–g) after correcting for age (F195≤3.71, p≥0.06), sex (F195≤7.78, p≥0.006), and TST (F195≤0.36, p≥0.55), further reinforcing the idea that it is the coupling of sleep microstructure elements that matters rather than their individual occurrence. Likewise, unlike previous reports (Mander et al., 2015; Winer, 2019), we did not find any association between the MPFC Aβ burden and the density of SWs, irrespective of their type (F195=1.18, p=0.28) or slow wave energy (SWE) – the cumulated power generated in the delta band: (F195=0.84, p=0.36), after correcting for age (F195≥4.2, p≤0.04, r²β*≥0.04), sex (F195≥6.63, p≤0.01, r²β*≥0.07), and TST (F195≤0.36, p≥0.34). This is true also when attempting to reproduce a previous finding made in a smaller and older sample than in the present study (Mander et al., 2015): in our sample, power in the 0.5–1 Hz EEG frequency band as well as the ratio between the power in the 0.5–1 Hz and 1–4 Hz bands were not associated with the burden of Aβ (main effect of Aβ burden: F195<2.85, p≥0.10) after correcting for age (F195≥3.28 p=0.014, r²β*=0.03), sex (F195≤0.25, p≥0.12), and TST (F195≤2.7, p≥0.10) (Figure 1h).
The anchoring of spindles onto slow-switcher SWs is associated to memory change over 2 years
A large subsample (N=66) completed a follow-up visit at 2 years which included the MST consisting in a pattern separation task targeting the ability to distinguish between highly resembling memory events, a hippocampus-dependent task which is very sensitive to early cognitive decline (Stark et al., 2013; Marks et al., 2017). We first observed an overall decline in performance between the baseline and follow-up performance at the MST (t65=2.19, p=0.032). We then explored whether the coupling of spindles with the slow switcher SWs was associated with memory performance decline over 2 years (Figure 2a). Statistical analyses revealed a significant negative link between the relative change in memory performance and the phase of spindle anchoring onto slow switcher SWs, indicating that an earlier spindle onset is predictive of a memory worsening over 2 years (main effect of spindle-slow switcher SW coupling: F161=4.80, p=0.032, r²β*=0.07), after correcting for age (F161=0.25, p=0.62), sex (F161=0.20, p=0.66), and education (F161=0.25, p=0.62) (Figure 2d–e). No such association was detected when considering spindle coupling to fast switcher SWs (main effect of spindle-fast switcher SW coupling: F161=2.51, p=0.12; main effect of age: F161=0.33, p=0.57; main effect of sex: F161=0.18, p=0.68; and main effect of education: F161=1.11, p=0.30). Further statistical analyses show that the memory performance change over the 2 year was not significantly related to the MPFC Aβ burden (main effect of Aβ burden: F160=2.33, p=0.13; main effect of age: F160=1.27, p=0.26; main effect of sex: F160=0.03, p=0.87; and main effect of education: F160=0.41, p=0.53). When considering baseline performance to the MST across the entire sample (N=100; i.e. assessed at the same time as the sleep measures), we found no significant link between the coupling of the spindles onto both SW types and the performance on the task (main effect of spindle-SW coupling: F196≥0.2, p≤0.61; main effect of age: F196≤0.36, p≥0.55; main effect of sex: F196≤0.54, p≥0.46; and main effect of education: F196≤0.47, p≥0.50 (Figure 2b–c)).
Discussion
In order to unravel early associations between the microstructure of sleep and the burden of Aβ in the brain, and their cognitive implications, we collected polysomnography, PET, and behavioural data in a relatively large sample of individuals without cognitive impairments or sleep disorders. To this end, we recruited individuals in late middle age (50–70 years) that could in most cases only present limited age-related alterations in sleep and accumulation of Aβ protein in the brain (Grothe et al., 2017). We investigated whether the coupling of spindles onto SWs, showing a slower and a faster frequency of transition from the down to the up states (slow and fast switcher SWs), was associated to the accumulation of Aβ over the MPFC. We further explored whether the coupling of spindles onto SWs was associated with the performance to a sensitive memory test, assessed at the time of the sleep and PET recordings and, longitudinally, 2 years later. The coupling of spindles onto the slow, but not the fast, switcher SWs was significantly associated with the Aβ PET signal assessed over the MPFC. Moreover, this coupling between spindles and slow switcher SWs was significantly linked to the memory performance change detected 2 years after the initial assessment. Overall, our results provide compelling evidence that the link between sleep and the accumulation of Aβ over the MPFC, an early AD-related brain features, involves the precise and timely coupling of two key elements of NREM sleep, spindles and SWs, and that this coupling bears a predictive value for the subsequent decline in memory performance. Our study does not indicate, at least not in these healthy and relatively young older adults, that the number of spindles or SWs generated overnight is associated with the accumulation of Aβ over the MPFC.
Sleep SWs provide a readout of the homeostatic sleep pressure and are more prevalent at the beginning relative to the end of the night (Achermann et al., 1993). In addition, both the density of spindles and SWs have separately been related to overnight consolidation of memory (Miyamoto et al., 2017; Fernandez and Lüthi, 2020). They actively take part in information transfer from hippocampic to neocortical networks and in synaptic plasticity (Ulrich, 2016; Miyamoto et al., 2017). Recent research has put forward the importance of their precise phase coupling during sleep, and reported an age-related difference in that coupling (Helfrich et al., 2018). In the younger individuals, spindles tend to reach their maximum around the cortical up-state of the SWs, whereas in older individuals, spindles occur earlier on the depolarisation phase of the SWs. This earlier coupling between spindles and SWs is related to a poorer overnight memory retention (Helfrich et al., 2018), suggesting a suboptimal neuronal interplay for the exchange of information during sleep in older individuals. In line with a presumed suboptimal coupling between spindles and SWs, we find that, when controlling for age, individuals for which the spindles occur earlier during the transition phase of the slow switcher SWs show higher Aβ burden over the MPFC and a worse memory change over time. We did not find any significant effect of the age of participants on the phase of the coupling between spindles and SWs. This is likely due, in part, to the restricted age range of our participants, but also to the variability existing between the individuals in the changes they undergo in their sleep during ageing. Importantly, the relationship we observed between the phase of the spindle onset onto the SWs and the MPFC Aβ burden shows the same directionality (i.e. earlier spindles onto SWs) as the changes previously reported in older individuals.
Sleep SWs are classically divided into slow oscillations and delta waves based on whether their overall frequency lay between ~0.5 and 1 Hz or between 1 and 4 Hz, respectively (Steriade, 2003). Although meaningful, the definition of these frequency bands is arguably arbitrary. SWs were recently divided into two categories based on an objective functional feature consisting in the frequency of their transition from the down- to up-state, which reflects the relative synchronisation of the depolarisation of the neurons when generating a SW (Bouchard et al., 2021). Beyond the well-known decrease in the production of SWs in ageing, the slow switcher SWs were relatively preserved in the older individuals compared to the fast switcher SWs. This finding suggests that the two populations of SWs constitute distinct elements of sleep microstructure. Three present results confirm that the two types of SWs – slow and fast switchers – behave differentially. First, sleep spindles show a difference in their preferential coupling with the transition period from down-to-up state of the slow and fast switcher SWs. While spindles occurring concomitantly to slow switchers SWs show a clear preference to the late part of the depolarisation phase, spindles co-occurring with fast switcher SWs show more widespread phase of coupling (that is still not randomly/uniformly distributed). Furthermore, only the coupling of spindles onto slow switcher SWs was significantly associated to the early accumulation of Aβ in the brain. Finally, only the coupling of spindles and slow switcher SWs was predictive of the memory change after 2 years. Our results support that slow switcher SWs, the type previously reported to be relatively spared during ageing (Bouchard et al., 2021), are important to the development of AD-related pathological changes, at least in the form of Aβ protein accumulation, and for the subsequent development of subtle alteration in cognitive abilities, leastways over the memory domain.
Sleep spindles are considered as thalamic events. They are generated through the interplay between the inhibitory reticular nucleus and the excitatory thalamocortical neurons, which project to the cortical neurons that feedback in turns to the thalamus (Ulrich, 2016). In contrast, SWs are intrinsically cortical. They consist in the spontaneous alternations between down (hyperpolarised) and up (depolarised) neuronal states (Adamantidis et al., 2019). The SWs undergo, however, an influence from the thalamus reticular nucleus that contributes to the synchronisation of distant neuronal populations (Lee et al., 2004; Adamantidis et al., 2019). Hence, although the exact neurophysiological origin of the functional associations between spindles and SWs remain to be established, the thalamus reticular neurons could arguably be involved. Our results would therefore indicate that the early aggregation of Aβ protein in the MPFC could disturb the thalamocortical interplay, driving the coalescence of spindles and SWs. This chronic disturbance would then trigger the cognitive changes detectable after 2 years. The cross-sectional nature of our imaging data does not, however, allow to make inferences regarding the directionality of the relationship between the spindle-slow switcher SW coupling and the burden of Aβ protein. Evidence accumulates to show that if AD-neuropathological hallmarks can affect the quality of sleep, sleep can also impinge onto these hallmarks (Van Egroo et al., 2019b). Specifically relevant in the context of this study, the occurrence of the SWs has been associated with a transient increase in the glymphatic flow related to local variations in neuromodulator concentrations (Xie et al., 2013). It is therefore possible that the changes in the density of SWs occurring during ageing affect both their coupling with spindles and the early accumulation of Aβ protein.
We found that the density of either spindles or SWs is not related to the early accumulation of the Aβ protein in the MPFC. This suggests a specific role for the coupling of both elements of sleep microstructure. Moreover and as previously reported in an intermediate analysis of a subsample of the same study (Van Egroo et al., 2019a), we do not find significant associations between more macroscopic measures, which are typically used to characterise sleep, and the accumulation of Aβ protein. The cumulated power of the oscillations generated over the entire 0.5–4 Hz delta band (i.e. the SWE) or only over the slower 0.5–1 Hz portion of SWE (as in Mander et al., 2015) as well as the ratio between slower (0.5–1 Hz) and faster (1–4 Hz) SWE (as in Mander et al., 2015) was not associated with PET measures. This finding argues against the idea that these rougher measures of sleep disruption constitute the earliest manifestations of the association between sleep and AD-related processes and contrasts with previous reports in a smaller sample of individuals older (>70 years) than our sample. Altogether, these discrepancies reinforce the idea that the alteration in the microstructure of sleep, consisting in the coupling of the spindles onto a specific subpopulation of SWs, as reported here, but also in the occurrence of microarousals during sleep we previously reported based on the same sample (Chylinski et al., 2021), shows a prior association with AD-related processes compared with the amount of slow brain oscillations generated during overnight sleep. The latter may only be significantly associated at a later age, when the pathophysiological changes are already more substantial. In addition, the coupling of spindles onto slow switcher SWs is predictive of the future change in memory performance. Sleep microstructure could therefore constitute a promising early marker of future-cognitive and brain-ageing trajectory (Chylinski et al., 2021). We did not evaluate whether distinct links between the SW types and slow and fast spindles are observed. As some reports describe that fast spindles are rather coupled to the up-state of the SWs and slow spindles tend to occur on the waning depolarisation phase of the SWs (Mölle et al., 2011), we could hypothesise that the associations we observe are probably rather driven by fast spindles. Future investigations are, however, needed to confirm this hypothesis.
One should bear in mind the potential limitations of our study. First, although we collected data in a relatively large sample, we may have insufficient power to detect some associations with other sleep measures. One can nevertheless frame our findings in relative terms such that association between spindle-SWs coupling and the early accumulation of Aβ protein is at least stronger than the association with the coupling of spindles onto fast switchers SW, the density of SWs and spindles, the SWE, etc. Also, effect sizes of the significant associations we found were relatively small, in line with the complex and multifactorial nature of AD. In addition, the longitudinal aspect of our study is relatively short-termed and only concerned the performance to a sensitive mnesic task while it did not include sleep EEG and the PET assessments. Further studies should replicate our findings and evaluate the predictive value of such parameters on longer longitudinal protocols, and the evolution of the sleep EEG and the PET parameters as well as their generalisability over other precociously impacted cognitive abilities. Finally, given that our protocol does not include manipulation of the coupling of the spindles onto the SWs, it precludes any inference on the causality of one aspect onto the other. It may be that cognition and the coupling of spindles and SWs are sensitive to the same age-related or AD-related phenomenon (e.g. presence of amyloid oligomers, or non-amyloid processes, that would go mostly undetected using common PET scan Aβ radioligand; Yamin and Teplow, 2017).
Together, our findings reveal that the timely occurrence of spindles onto a specific type of SWs showing a relative preservation in ageing may play an important role in ageing trajectory, both at the cognitive level and with regards to structural brain integrity. These findings may help to unravel early links between sleep, AD-related pathophysiology, and cognitive trajectories in ageing and warrants future clinical trials attempting at manipulating sleep microstructure or Aβ protein accumulation.
Materials and methods
Study design and participants
Request a detailed protocol101 healthy participants aged from 50 to 70 years (68 women; mean ± SD = 59.4±5.3 years) were enrolled between 15 June 2016 and 2 October 2019 for a multimodal cross-sectional study taking place at the GIGA-Cyclotron Research Centre/In Vivo Imaging of the University of Liège (Cognitive fitness in ageing – COFITAGE – study) which has already led to several scientific publications (e.g. Chylinski et al., 2021; Rizzolo et al., 2021). One participant was excluded from analyses due to the lack of PET imaging data. The exclusion criteria were as follows: clinical symptoms of cognitive impairment (Mattis Dementia Rating Scale >130; Mini-Mental State Evaluation >27); recent psychiatric history, or severe brain trauma; self-reported or clinically diagnosed sleep disorder; ≤18 and ≥29; use of medication affecting sleep or the CNS; smoking; excessive alcohol (>14 units/week) or caffeine (>5 cups/day) consumption; shift work in the 6 months or transmeridian travel in the 2 months preceding the study. All participants gave their written informed consent prior to their participation. The study was registered with EudraCT 2016-001436-35. All procedures were approved by the Hospital-Faculty Ethics Committee of ULiège. All participants signed an informed consent prior to participating in the study.
Sleep assessment
Request a detailed protocolA first night of sleep was recorded at the laboratory under full polysomnography to avoid potential first night effects and exclude volunteers with sleep apnoea (Apnea Hypopnea Index ≥15 /hr). A second night of sleep was recorded with EEG, following 1 week of regular sleep-wake schedule based on each participant’s preferred bed and wake-up time (compliance was verified by actimetry and sleep diary – Actiwatch, Cambridge Neurotechnology, UK). Sleep was recorded with N7000 amplifiers (EMBLA, Natus, Planegg, Germany). The recording comprised 11 EEG derivations, placed according to the 10–20 system (F3, Fz, and F4; C3, Cz, and C4; P3, Pz, and P4; O1 and O2), two bipolar electrooculogram, and two bipolar submental electromyogram electrodes. Sampling was set at 200 Hz, and the signal was re-referenced to the mean of the two mastoids. Recordings were scored for sleep stages in 30 s windows using a validated automatic algorithm (ASEEGA, Physip, Paris, France) (Berthomier et al., 2007; Peter-Derex et al., 2021). Automatic arousal and artefact detection (Chylinski et al., 2020; ’t Wallant et al., 2016) were performed in order to remove EEG segments containing artefacts and arousals from further analysis.
Slow wave and spindle detections
Request a detailed protocolOnly the frontal electrodes were considered because the frontal cortex is an early site showing Aβ deposit and is the primary generator of the SWs during sleep (Mander et al., 2015; Dang-Vu et al., 2005; Dang-Vu et al., 2010; Saletin et al., 2013) as well as to facilitate interpretations of future large-scale studies using headband EEG restricted to frontal electrodes (Lucey et al., 2019). SWs were automatically detected during N2 and N3 30s-epochs of NREM sleep devoid of artefacts/arousals >5 s long, using a previously developed algorithm (Rosinvil et al., 2021). Data were first band-filtered between 0.3 and 4.0 Hz with a linear phase finite impulse response (FIR) filter. Following recent work, SW detection criteria was adapted for age and sex (Rosinvil et al., 2021): peak to peak amplitude ≥70 µV (resp. ≥60.5 µV) and negative amplitude ≤–37 µV (resp. ≤–32 µV) was used for women (resp. for men), instead of the standard ≥75 µV and ≤–40 µV. The duration of the negative deflection had to fit in the range 125–1500 ms, and the duration of the positive deflection could not exceed 1000 ms. The SWs were sorted according to their transition frequency (Bouchard et al., 2021) (inverse of the duration between the hyperpolarised and depolarised state) into either slow or fast switchers (the critical value for distinguishing between the two types being the intersection between two Gaussian, around 1.2 Hz) (Figure 3).
Sleep spindles were also automatically detected over the same N2 and N3 epochs with a previously published methods (Gaudreault et al., 2018; Lafortune et al., 2014; Martin et al., 2013). The EEG signal was bandpass filtered between 10 and 16 Hz with a linear phase FIR filter (–3 dB at 10 and 16 Hz). The envelope amplitude of the Hilbert transform of this band-limited signal was smoothed, and a threshold was set at the 75th percentile. All events of duration between 0.5 and 3 s were then selected as a spindle.
After detection of SWs and spindles, analysis of their coincidence was performed. A coincidence was defined as to occurrence of the ignition of a spindle within the time frame of a SW: SW ignition at zero µV=phase 0°, SW maximum hyperpolarisation=π/2, zero crossing=π, SW maximum depolarisation=1.5π, and SW termination at zero µV=2π. This criterion was used on slow and fast switchers. All metrics related to SWs and spindles and to their coupling were averaged over F3, F4, and Fz derivations.
MRI data
Request a detailed protocolQuantitative multi-parametric MRI acquisition was performed on a 3-Tesla MR scanner (Siemens MAGNETOM Prisma, Siemens Healthineers, Erlangen, Germany). Quantitative maps were obtained by combining the images using different parameters sensitive to distinct tissue properties. The multi-parameter mapping was based on multi-echo 3D fast low angle shot at 1 mm isotropic resolution (Weiskopf and Helms, 2008). This included three datasets with T1, proton density (PD), and magnetisation transfer (MT)-weighted contrasts imposed by the choice of the flip angle (FA = 6° for PD and MT, 21° for T1) and the application of an additional off-resonance Gaussian-shaped Radio Frequency (RF) pulse for the MT-weighted acquisition. MRI multi-parameter maps were processed with the hMRI toolbox (Tabelow et al., 2019) (http://hmri.info) and SPM12 (Welcome Trust Centre for Neuroimaging, London, UK) to obtain notably a quantitative MT map, which was segmented into grey matter, white matter, and CSF using unified segmentation (Ashburner and Friston, 2005). Flow-field deformation parameters obtained from DARTEL spatial normalisation of the individual MT maps were applied to the averaged co-registered PET images (Ashburner, 2007). The volumes of interest were determined using the automated anatomical labelling atlas (Tzourio-Mazoyer et al., 2002).
PET scan
Request a detailed protocolAβ PET imaging was performed using [18F] Flutemetamol, except for three volunteers for which [18F] Florbetapir was used. PET scans were performed on an ECAT EXACT +HR scanner (Siemens, Erlangen, Germany). Participants received a single dose of the radioligand in the antecubital vein (target dose 185±10% MBq); image acquisition started 85 min after the injection and consisted of four frames of 5 min, followed by a 10 min transmission scan using 68Ge line sources. Images were reconstructed using a filtered back-projection algorithm including corrections for the measured attenuation, dead time, random events, and scatter using standard software (Siemens ECAT – HR +V7.1, Siemens/CTI Knoxville, TN, USA). Individual PET average images were produced using all frames and were then manually reoriented according to MT-weighted structural MRI volumes and co-registered to the individual space structural MT map. Standardised uptake value ratio (SUVR) was computed using the whole cerebellum as reference region (Klunk et al., 2015). As images were acquired using two different radioligands, their SUVR values were converted into Centiloid units (Klunk et al., 2015) (the validation of the procedure in our sample was previously published Narbutas et al., 2021). The Aβ burden was averaged over a mask covering the MPFC, previously reported to undergo the earliest aggregation sites for Aβ pathology (Grothe et al., 2017).
Cognitive assessments
Request a detailed protocolAs part of an extensive neuropsychological assessment, participants were administered the MST (Stark et al., 2015), a visual recognition memory task. After an incidental encoding phase during which participants were randomly presented 128 common objects for a period of 2 s and were instructed to determine whether the object presented on the screen was rather an ‘indoor’ or ‘outdoor’ item, the recognition memory phase consisted in the presentation of 192 objects (64 old, presented previously – target items; 64 similar but not identical to the previously presented stimuli – lure; and 64 new objects – foil items). In this phase, participants were instructed to determine whether the presented object was new (foil), previously presented (old), or similar but not perfectly identical (lure). For statistical analyses, the recognition memory (RM) score was used, computed as the difference between the rate of calling a target item ‘old’ minus the rate of calling a foil item ‘old’ (P[‘old’|target]-P[‘old’|foil]) (Stark et al., 2013; Rizzolo et al., 2021).
The MST was administered at two timepoints: the first time, the day preceding the baseline night, during a cognitive evaluation performed ~6.5 hr before habitual bedtime. The second neuropsychological evaluation was carried out ~24 months after the first one (mean 767±54 days) between 4 and 10h after wake up time. The memory decline score was computed as the baseline performance minus the follow-up performance, divided by the baseline performance, so that a higher score indicates a higher decline over the 2 years.
Statistics
Our primary analysis tested in a single model whether the phase of coupling, as dependent variable, would differ between SW types and would be associated with Aβ burden (independent variables), controlling for sex, age, TST, and SW duration. The phase of spindle-SW coupling was set as the phase of the onset of the spindle on the SW converted to its cosine value, to deal with the circularity of the phase variable and perform linear statistics (analysis using the phase in degrees yielded the same outcomes). As we tested our hypothesis in a single test, it did not require correction for multiple testing, and significance was set at p<0.05. Our main exploratory analysis was to assess whether spindle-SW coupling to a specific SW type was associated with longitudinal change in the performance to the memory test, controlling for sex, age, education, and TST. The remaining analyses were aimed a better characterising, the primary analysis, or the main exploratory analyses.
Statistical analyses were performed using GLMMs in SAS 9.4 (SAS Institute, Cary, NC). The distribution of dependent variables was verified in MATLAB 2013a, and the GLMMs were adapted accordingly. Subject was treated as a random factor and each model was corrected for age and sex effects. Kenward-Roger’s correction was used to determine the degrees of freedom. Cook’s distance was used to assess the potential presence of outliers driving the associations, and as values ranged below 0.45, no data point was excluded from the analyses (a Cook’s distant >1 is typically considered to reflect outlier value). Semi-partial R² (R²β*) values were computed to estimate the effect sizes of significant fixed effects and statistical trends in all GLMMs (Jaeger et al., 2017). p-values in post-hoc contrasts (difference of least square means) were adjusted for multiple testing using Tukey’s procedure. Watson’s non-parametric two-sample U² test for circular-normal data was performed in MATLAB 2019 to assess the difference between the distribution of spindle onset on the phase of SWs for slow and fast switcher SWs. We further assessed whether the distribution of spindle onset on the phase of SWs per type was different from a uniform distribution. For each SW type, we generated series of uniformly distributed random values composed of the same number of values spanning the same ranges. Watson’s non-parametric two-sample U² test compared this random series to the actual values.
Optimal sensitivity and power analyses in GLMMs remain under investigation (e.g. Kain et al., 2015). We nevertheless computed a prior sensitivity analysis to get an indication of the minimum detectable effect size in our main analyses given our sample size. According to G*Power 3 (version 3.1.9.4) (Faul et al., 2009), taking into account a power of 0.8, an error rate α of 0.05, and a sample size of 100 allowed us to detect small effect sizes r>0.27 (two-sided; absolute values; CI: 0.07–0.44; R²>0.07, R² CI:0.005–0.19) within a linear multiple regression framework including one tested predictor (Aβ) and two covariates (age and sex).
Data availability
Request a detailed protocolThe data and analysis scripts supporting the results included in this manuscript are publicly available via the following open repository: https://gitlab.uliege.be/CyclotronResearchCentre/Public/fasst/slow-wave-spindle-coupling-and-amyloid; Vandewalle, 2022. We used Matlab script for MRI and PET data processing and to detect slow wave and spindles as well as their coupling, while we used SAS for statistical analyses. The raw data could be identified and linked to a single subject and represent a huge amount of data (>200 Gb). Researchers willing to access to the raw should send a request to the corresponding author (GV). Data sharing will require evaluation of the request by the local Research Ethics Board and the signature of a data transfer agreement (DTA).
Data availability
The data and analysis scripts supporting the results included in this manuscript are publicly available via the following open repository: https://gitlab.uliege.be/CyclotronResearchCentre/Public/fasst/slow-wave-spindle-coupling-and-amyloid (copy archived at swh:1:rev:a63699b5f06c98ade284689af5671e13c77751f2). We used Matlab script for MRI and PET data processing and to detect slow waves and spindles as well as their coupling, while we used SAS for statistical analyses. The raw data could be identified and linked to a single subject and represent a huge amount of data (> 200 Gb). Researchers willing to access the raw data should send a request to the corresponding author (GV). Data sharing will require evaluation of the request by the local Research Ethics Board of the Faculty of Medicine at the University of Liège, Belgium and the signature of a data transfer agreement (DTA).
References
-
Oscillating circuitries in the sleeping brainNature Reviews. Neuroscience 20:746–762.https://doi.org/10.1038/s41583-019-0223-4
-
Biomarkers of dementia in obstructive sleep apneaSleep Medicine Reviews 42:139–148.https://doi.org/10.1016/j.smrv.2018.08.001
-
Sleep slow wave changes during the middle years of lifeThe European Journal of Neuroscience 33:758–766.https://doi.org/10.1111/j.1460-9568.2010.07543.x
-
Risk of Alzheimer’s Disease in Obstructive Sleep Apnea Syndrome: Amyloid-β and Tau ImagingJournal of Alzheimer’s Disease 66:733–741.https://doi.org/10.3233/JAD-180640
-
Statistical power analyses using G*Power 3.1: tests for correlation and regression analysesBehavior Research Methods 41:1149–1160.https://doi.org/10.3758/BRM.41.4.1149
-
Sleep Spindles: Mechanisms and FunctionsPhysiological Reviews 100:805–868.https://doi.org/10.1152/physrev.00042.2018
-
Learning-dependent increases in sleep spindle densityThe Journal of Neuroscience 22:6830–6834.
-
An R2statistic for fixed effects in the generalized linear mixed modelJournal of Applied Statistics 44:1086–1105.https://doi.org/10.1080/02664763.2016.1193725
-
Sleep quality and preclinical Alzheimer diseaseJAMA Neurology 70:587–593.https://doi.org/10.1001/jamaneurol.2013.2334
-
Sleep and Alzheimer disease pathology--a bidirectional relationshipNature Reviews. Neurology 10:115–119.https://doi.org/10.1038/nrneurol.2013.269
-
Amyloid-beta dynamics are regulated by orexin and the sleep-wake cycleScience (New York, N.Y.) 326:1005–1007.https://doi.org/10.1126/science.1180962
-
Effect of sleep on overnight cerebrospinal fluid amyloid β kineticsAnnals of Neurology 83:197–204.https://doi.org/10.1002/ana.25117
-
Reduced non-rapid eye movement sleep is associated with tau pathology in early Alzheimer’s diseaseScience Translational Medicine 11:eaau6550.https://doi.org/10.1126/scitranslmed.aau6550
-
β-amyloid disrupts human NREM slow waves and related hippocampus-dependent memory consolidationNature Neuroscience 18:1051–1057.https://doi.org/10.1038/nn.4035
-
Topography of age-related changes in sleep spindlesNeurobiology of Aging 34:468–476.https://doi.org/10.1016/j.neurobiolaging.2012.05.020
-
The critical role of sleep spindles in hippocampal-dependent memory: a pharmacology studyThe Journal of Neuroscience 33:4494–4504.https://doi.org/10.1523/JNEUROSCI.3127-12.2013
-
The Roles of Cortical Slow Waves in Synaptic Plasticity and Memory ConsolidationFrontiers in Neural Circuits 11:1–8.https://doi.org/10.3389/fncir.2017.00092
-
Mechanisms linking circadian clocks, sleep, and neurodegenerationScience (New York, N.Y.) 354:1004–1008.https://doi.org/10.1126/science.aah4968
-
Automatic analysis of single-channel sleep EEG in a large spectrum of sleep disordersJournal of Clinical Sleep Medicine 17:393–402.https://doi.org/10.5664/jcsm.8864
-
Encoding difficulty promotes postlearning changes in sleep spindle activity during nappingThe Journal of Neuroscience 26:8976–8982.https://doi.org/10.1523/JNEUROSCI.2464-06.2006
-
Stability of age-related deficits in the mnemonic similarity task across task variationsBehavioral Neuroscience 129:257–268.https://doi.org/10.1037/bne0000055
-
The corticothalamic system in sleepFrontiers in Bioscience 8:d878–d899.https://doi.org/10.2741/1043
-
Sleep Spindles as Facilitators of Memory Formation and LearningNeural Plasticity 2016:1796715.https://doi.org/10.1155/2016/1796715
-
SoftwareSlow wave-spindle coupling and amyloid, version swh:1:rev:a63699b5f06c98ade284689af5671e13c77751f2Software Heritage.
-
Multi-parameter mapping of the human brain at 1mm resolution in less than 20 minutesProceeding Int. Soc. Magetic Reson. Med 16:2241.
-
BookSleep as a Potential Biomarker of Tau and β-Amyloid Burden in the Human Brain Abbreviated Title Authors Center for Human Sleep ScienceUniversit: Department of Psychology , University of California Berkeley , Department of Psychiatry and Human Behavior.
-
Sleep drives metabolite clearance from the adult brainScience (New York, N.Y.) 342:373–377.https://doi.org/10.1126/science.1241224
-
Pittsburgh Compound-B (PiB) binds amyloid β-protein protofibrilsJournal of Neurochemistry 140:210–215.https://doi.org/10.1111/jnc.13887
-
Automatic artifacts and arousals detection in whole-night sleep EEG recordingsJournal of Neuroscience Methods 258:124–133.https://doi.org/10.1016/j.jneumeth.2015.11.005
Article and author information
Author details
Funding
Fonds De La Recherche Scientifique - FNRS (FRSM 3.4516.11)
- Gilles Vandewalle
Fédération Wallonie-Bruxelles (ARC-SLEEPDEM 17/09)
- Daphne Chylinski
- Maxime Van Egroo
- Justinas Narbutas
- Christine Bastin
- Christophe Phillips
- Fabienne Collette
- Pierre Maquet
- Gilles Vandewalle
European Regional Development Fund (RAdiomed)
- Daphne Chylinski
- Maxime Van Egroo
- Justinas Narbutas
- Vincenzo Muto
- Mohamed Ali Bahri
- Eric Salmon
- Christine Bastin
- Christophe Phillips
- Fabienne Collette
- Pierre Maquet
- Gilles Vandewalle
Canadian Institutes of Health Research (grant number 190750)
- Julie Carrier
- Jean-Marc Lina
General Electric (ISS290)
- Eric Salmon
- Christine Bastin
- Christophe Phillips
- Fabienne Collette
- Pierre Maquet
- Gilles Vandewalle
Fonds De La Recherche Scientifique - FNRS
- Maxime Van Egroo
- Christine Bastin
- Christophe Phillips
- Fabienne Collette
- Gilles Vandewalle
Fondation Recherche Alzheimer
- Fabienne Collette
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank M Blanpain, P Cardone, M Cerasuolo, E Lambot, P Ghaemmaghami, C Hagelstein, S Laloux, E Balteau, A Claes, C Degueldre, B Herbillon, P Hawotte, B Lauricella, A Lesoine, A Luxen, X Pepin, E Tezel, D Marzoli, C Schmidt, and P Villar González for their help in different steps of the study. This work was supported by Fonds National de la Recherche Scientifique (FRS-FNRS, FRSM 3451611, Belgium), Actions de Recherche Concertées (ARC SLEEPDEM 17/27–09) of the Fédération Wallonie-Bruxelles, University of Liège (ULiège), Fondation Simone et Pierre Clerdent, European Regional Development Fund (ERDF, Radiomed Project), and the Canadian Institutes of Health Research (CIHR) (grant number 190750). [18 F]Flutemetamol doses were provided and cost covered by GE Healthcare Ltd (Little Chalfont, UK) as part of an investigator sponsored study (ISS290) agreement. This agreement had no influence on the protocol and results of the study reported here MVE, CBastin, FC, CP, and GV are/were supported by the FRS-FNRS Belgium.
Ethics
Human subjects: The study was registered with EudraCT 2016-001436-35. All procedures were approved by the Hospital-Faculty Ethics Committee of ULiège. All participants signed an informed consent prior to participating in the study.
Copyright
© 2022, Chylinski et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,664
- views
-
- 449
- downloads
-
- 27
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
Complex macro-scale patterns of brain activity that emerge during periods of wakeful rest provide insight into the organisation of neural function, how these differentiate individuals based on their traits, and the neural basis of different types of self-generated thoughts. Although brain activity during wakeful rest is valuable for understanding important features of human cognition, its unconstrained nature makes it difficult to disentangle neural features related to personality traits from those related to the thoughts occurring at rest. Our study builds on recent perspectives from work on ongoing conscious thought that highlight the interactions between three brain networks – ventral and dorsal attention networks, as well as the default mode network. We combined measures of personality with state-of-the-art indices of ongoing thoughts at rest and brain imaging analysis and explored whether this ‘tri-partite’ view can provide a framework within which to understand the contribution of states and traits to observed patterns of neural activity at rest. To capture macro-scale relationships between different brain systems, we calculated cortical gradients to describe brain organisation in a low-dimensional space. Our analysis established that for more introverted individuals, regions of the ventral attention network were functionally more aligned to regions of the somatomotor system and the default mode network. At the same time, a pattern of detailed self-generated thought was associated with a decoupling of regions of dorsal attention from regions in the default mode network. Our study, therefore, establishes that interactions between attention systems and the default mode network are important influences on ongoing thought at rest and highlights the value of integrating contemporary perspectives on conscious experience when understanding patterns of brain activity at rest.
-
- Neuroscience
Motivation depends on dopamine, but might be modulated by acetylcholine which influences dopamine release in the striatum, and amplifies motivation in animal studies. A corresponding effect in humans would be important clinically, since anticholinergic drugs are frequently used in Parkinson’s disease, a condition that can also disrupt motivation. Reward and dopamine make us more ready to respond, as indexed by reaction times (RT), and move faster, sometimes termed vigour. These effects may be controlled by preparatory processes that can be tracked using electroencephalography (EEG). We measured vigour in a placebo-controlled, double-blinded study of trihexyphenidyl (THP), a muscarinic antagonist, with an incentivised eye movement task and EEG. Participants responded faster and with greater vigour when incentives were high, but THP blunted these motivational effects, suggesting that muscarinic receptors facilitate invigoration by reward. Preparatory EEG build-up (contingent negative variation [CNV]) was strengthened by high incentives and by muscarinic blockade, although THP reduced the incentive effect. The amplitude of preparatory activity predicted both vigour and RT, although over distinct scalp regions; frontal activity predicted vigour, whereas a larger, earlier, central component predicted RT. The incentivisation of RT was partly mediated by the CNV, though vigour was not. Moreover, the CNV mediated the drug’s effect on dampening incentives, suggesting that muscarinic receptors underlie the motivational influence on this preparatory activity. Taken together, these findings show that a muscarinic blocker impairs motivated action in healthy people, and that medial frontal preparatory neural activity mediates this for RT.