Individual Differences in White Matter Diffusion Affect Sleep Oscillations

The Journal of Neuroscience, January 2, 2013 • 33(1):227–233 • 227

Systems/Circuits

Individual Differences in White Matter Diffusion Affect Sleep Oscillations Giovanni Piantoni,1 Simon-Shlomo Poil,2 Klaus Linkenkaer-Hansen,2 Ilse M. Verweij,1 Jennifer R. Ramautar,1 Eus J. W. Van Someren,1,2* and Ysbrand D. Van Der Werf1,3* 1

Department of Sleep and Cognition, Netherlands Institute for Neuroscience, 1105 BA Amsterdam, The Netherlands, 2Departments of Integrative Neurophysiology and Medical Psychology, Center for Neurogenomics and Cognitive Research, and 3Department of Anatomy and Neurosciences, Neuroscience Campus Amsterdam, 1081 HV Amsterdam, The Netherlands

The characteristic oscillations of the sleeping brain, spindles and slow waves, show trait-like, within-subject stability and a remarkable interindividual variability that correlates with functionally relevant measures such as memory performance and intelligence. Yet, the mechanisms underlying these interindividual differences are largely unknown. Spindles and slow waves are affected by the recent history of learning and neuronal activation, indicating sensitivity to changes in synaptic strength and thus to the connectivity of the neuronal network. Because the structural backbone of this network is formed by white matter tracts, we hypothesized that individual differences in spindles and slow waves depend on the white matter microstructure across a distributed network. We recorded both diffusion-weighted magnetic resonance images and whole-night, high-density electroencephalography and investigated whether individual differences in sleep spindle and slow wave parameters were associated with diffusion tensor imaging metrics; white matter fractional anisotropy and axial diffusivity were quantified using tract-based spatial statistics. Individuals with higher spindle power had higher axial diffusivity in the forceps minor, the anterior corpus callosum, fascicles in the temporal lobe, and the tracts within and surrounding the thalamus. Individuals with a steeper rising slope of the slow wave had higher axial diffusivity in the temporal fascicle and frontally located white matter tracts (forceps minor, anterior corpus callosum). These results indicate that the profiles of sleep oscillations reflect not only the dynamics of the neuronal network at the synaptic level, but also the localized microstructural properties of its structural backbone, the white matter tracts.

Introduction People differ markedly with respect to their habitual sleep duration and timing. Electroencephalographic (EEG) recordings show large interindividual differences in sleep oscillation characteristics, such as slow waves and spindles, which, at the same time, have remarkable trait-like stability within an individual across nights (Kerkhof, 1991; Finelli et al., 2001; De Gennaro et al., 2005; Buckelmu¨ller et al., 2006; Gander et al., 2010; Botella-Soler et al., 2012). The topography of slow-wave expression over the scalp is unique to each person (Rusterholz and Achermann, 2011). Furthermore, subject-specific slow-wave expression profiles have been associated with functionally relevant attributes (Dijk, 2011) such as chronotype (Kerkhof, 1991), sleep maintenance (Dijk, 2006), and the risk of schizophrenia (Sarkar et al., 2010). SimiReceived April 27, 2012; revised Oct. 3, 2012; accepted Oct. 24, 2012. Author contributions: G.P., S.-S.P., K.L.-H., E.J.W.V.S., and Y.D.V.D.W. designed research; S.-S.P., K.L.-H., and Y.D.V.D.W. performed research; I.M.V. and J.R.R. contributed unpublished reagents/analytic tools; G.P. analyzed data; G.P., E.J.W.V.S., and Y.D.V.D.W. wrote the paper. This work was supported by Netherlands Organization for Scientific Research (NWO) Grants VIDI 016.095.359 (to Y.D.V.D.W.)andVICI453.07.001(toE.J.W.V.S.)andTopTalentGrant021.002.080(toS.-S.P.).WethankD.Stoffersfortechnical assistance and L. Sloot, D. Windhorst, and C. Zennig for their help with data acquisition. Part of the data were collected in collaboration with the Spinoza Centre for Neuroimaging, Amsterdam, The Netherlands, www.spinozacentre.nl. *E.J.W.V.S. and Y.D.V.D.W. are co-senior authors. Correspondence should be addressed to Giovanni Piantoni, Department of Sleep and Cognition, Netherlands Institute for Neuroscience, Meibergdreef 47, 1105BA Amsterdam, The Netherlands. E-mail: [email protected]. DOI:10.1523/JNEUROSCI.2030-12.2013 Copyright © 2013 the authors 0270-6474/13/330227-07$15.00/0

larly, sleep spindles have trait-like properties: the power profile in the spindle frequency range is remarkably constant over different nights within a person (Werth et al., 1997; De Gennaro et al., 2005, 2008) and indeed shows high heritability (De Gennaro et al., 2008). Spindle characteristics are likewise associated with meaningful measures, such as learning abilities (Schabus et al., 2008), cognitive performance (Bo´dizs et al., 2005), intelligence (Fogel and Smith, 2011), schizophrenia severity (Ferrarelli et al., 2007), and tolerance to environmental noise (Dang-Vu et al., 2010a). Despite this literature, little is known about the brain structural and functional mechanisms that cause the genetically specified interindividual differences and determine the properties of slow waves and spindles. Slow waves and spindles are affected by the recent history of induced neuronal activation and synaptic plasticity (for review, Tononi and Cirelli, 2006; Diekelmann and Born, 2010; Poe et al., 2010; Fogel and Smith, 2011). Because learning-induced plasticity mostly depends on synaptic changes, the profile of slow waves and spindles likely reflects synaptic strength (Esser et al., 2007; Vyazovskiy et al., 2007; Vyazovskiy et al., 2008; Vyazovskiy et al., 2009) and consequently connectivity dynamics of the neuronal network. Crucially, whereas synapses comprise the endpoints of the network, its structural backbone is formed by long-range myelinated white matter tracts. Therefore, we hypothesized that slow waves and spindles might not only depend on synaptic strength, but also on the structural properties of the white matter

228 • J. Neurosci., January 2, 2013 • 33(1):227–233

Piantoni et al. • White Matter Diffusion Affects Sleep Oscillations

Figure 1. Detection of slow-wave and spindle parameters. A, The black line shows a 20 s segment of sleep EEG referenced to the averaged mastoids. After 0.2– 4 Hz bandpass filtering (purple line), a slow wave was detected when the amplitude of the negative peak was less than ⫺60 ␮V (green line), the duration between two zero-crossings was between 0.2 and 1 s, and the negative peak-to-positive peak difference was ⬎75 ␮V. B, Three parameters were used to describe slow-wave activity: (1) slow-wave density, which is the number of slow waves in a minute; (2) the amplitude of the negative peak (blue line); and (3) the steepness of the rising slope, which is the amplitude of the negative peak divided by the interval between the negative peak and the following zero crossing (red line). C, To quantify spindles, the original signal shown in A was bandpass filtered between 11 and 16 Hz (blue line) and Hilbert-transformed into its envelope (purple line). A spindle was detected when the envelope of the band-passed signal was higher than the absolute threshold of 3 ␮V (green line) for 0.5–2 s. D, Two parameters were used to describe spindles: (1) spindle density, which is the number of spindles per minute; and (2) spindle power, which is the area outlined by the envelope (dark blue).

tracts. Indirect support for this hypothesis is provided by Buchmann et al. (2011), who applied magnetic resonance imaging (MRI) and voxel-based morphometry to demonstrate a correlation between corpus callosum volume and slow-wave power over 10 min of sleep. As indicated by Buchmann et al. (2011), more definite conclusions would require more sophisticated imaging methods, such as diffusion tensor imaging (DTI) (Jones, 2010). We, therefore, used DTI to investigate whether individual differences in the whole-brain white matter microstructure might contribute to the interindividual variability in slow wave and spindle characteristics as assessed by high-density sleep EEG. DTI allows for a detailed quantification of white matter properties, such as fractional anisotropy (FA), which is taken as a marker of white matter integrity (Beaulieu, 2010), and axial diffusivity (AD), which reflects the strength of diffusion along the axonal direction (Budde et al., 2007). Using this sensitive methodology, we here demonstrate, for the first time, partially overlapping and partially specific white matter networks associated with the subject-specific profiles of slow waves and spindles.

Materials and Methods Sleep recordings. Fifteen healthy males (mean age 25.4 ⫾ 4.8 SD) with no history of neurological disorders took part in this study after given written informed consent according to the Declaration of Helsinki. The protocol of this study was approved by the medical ethical review board of the VU University Medical Center (Amsterdam, The Netherlands). Participants followed a regular sleep schedule in the week preceding the recording, as verified from actigraphy and sleep logs, and were requested to abstain from alcohol and caffeine for 24 h before the recording session.

Whole-night 256-channel electroencephalography (Electrical Geodesics) was recorded referenced to the Cz electrode, high pass filtered at 0.01 Hz, low pass filtered at 200 Hz with a hardware Bessel filter, and digitized at 500 Hz. Using the FAST toolbox (Leclercq et al., 2011), two independent raters scored sleep according to standard criteria (Iber et al., 2007) in 30 s epochs. One participant was excluded because he did not reach sleep stage N3 (non-REM stage 3, i.e., deep sleep). Segments containing artifacts were visually detected and excluded from analysis. For each epoch, noisy channels were automatically excluded when the signals exceeded the absolute thresholds of 1000 ␮V or its variance was above 10000 ␮V 2. Slow wave and spindle detection was conducted after rereferencing to the averaged mastoid channels. Slow wave detection. Slow waves were automatically detected using methods described previously (Massimini et al., 2004; Dang-Vu et al., 2008) with custom code written in Matlab 7.13 (MathWorks). EEG potentials were averaged over four nonoverlapping areas of the scalp (on frontal, posterior, and lateral electrodes as in Massimini et al., 2004) and bandpass filtered between 0.2 and 4 Hz (Fig. 1A). The following criteria were applied independently to each local average to identify slow waves: (1) a negative zero crossing and a subsequent positive zero crossing separated by 0.2–1 s; (2) a negative peak between the two zero crossings with voltage less than ⫺60 ␮V; and (3) a peak-to-peak difference between the negative peak and the subsequent positive peak in the 1 s period after the negative peak of ⬎75 ␮V. The peak-to-peak threshold was based on the prescribed value for the detection of slow waves in the consensus sleep scoring manuals (Rechtschaffen and Kales, 1968; Iber et al., 2007). For each slow wave, we calculated (1) the maximum amplitude at the negative peak and (2) the slope between the negative peak and the following upward zero-crossing (Fig. 1B). Furthermore, we calculated slow wave density as the number of slow waves per minute. For


J. Neurosci., January 2, 2013 • 33(1):227–233 • 229

diffusion-weighted images at b ⫽ 1000 s/mm2 were acquired on a Philips Achieva 3T MRI scanner (Philips Medical Systems) with the following parameters: TR/TE ⫽ 8188/60 ms, flip angle ⫽ 90°, 224 ⫻ 224 ⫻ 120 mm FOV, 2 ⫻ 2 ⫻ 2 mm voxel size. The FMRIB Software Library (FSL) data analysis suite [FSL 4.1.9 (Smith et al., 2004); http://www.fmrib.ox.ac.uk/fsl/] was used for correction of eddy current distortions and motion artifacts and for the computation of FA and AD. Single-voxel statistical analysis of the FA and AD maps was carried out using tractbased spatial statistics (TBSS) (Smith et al., 2006). Briefly, individual maps are warped into the MNI152 template space using nonlinear registration. Individual maps are averaged across subjects, and voxels with FA ⬎ 0.2 are included in the analysis, resulting in the selection of only major fiber bundles. Individual DTI maps are then projected onto this skeleton of putative white matter tracts. Group statistics were computed separately for FA and AD on the slow-wave and spindle parameters. The design matrix included a column of ones to capture the intercept and a demeaned column for the sleep parameter, the latter representing the contrast of interest. In this way, the single-subject DTI values included in the white matter skeleton were correlated with the slow-wave and spindle parameters. Statistical tests were carried out using the FSL tool randomise, with a threshold-free cluster-enhanced correction for multiple comparisons [TFCE correction (Smith and Nichols, 2009)] using 10,000 permutations. The p values in the Tables refer to the most significant voxel in that particular contrast. The correlation maps with the slope of the slow wave in N3 and the Figure 2. White matter tracts where axial diffusivity was higher in individuals with a steeper slow-wave slope. The plot shows spindle power in N2 are plotted in Figures 2 the voxels with a positive significant correlation (thresholded at p ⱕ 0.05 after correction for multiple comparisons) superimposed and 3, respectively. These figures show the on the white matter skeleton (light gray) and on the study-specific MRI template (dark gray) following the neurological convention significant voxels thresholded at (corrected) (right is right). The average slope was computed from the slow waves detected in all the four regions of interest, while the color p ⱕ 0.05 after using tbss_fill overlaid on the coding indicates which region of interest was the most correlated with the AD at each voxel. The white matter tracts in which AD white matter skeleton in gray color and the was associated with the slope steepness are concentrated in the frontal white matter, including the forceps minor, the genu and study-specific MRI template. To identify anterior part of the corpus callosum (G CC and aCC, respectively). These frontal pathways connect to the superior longitudinal which electrodes were the most informative fascicle (sLF) and the temporal branch of the inferior fronto-occipital fascicle (iFOF). in determining the slow-wave slope, we calculated at each significant voxel which group of electrode had the highest correlation. The each participant, these values were then averaged separately for sleep four regions of interest are then color-coded in Figure 2. stages N2 and N3. The Matlab/FSL script for the complete analysis can be requested from Spindle detection. Spindles were automatically detected using an adapthe corresponding author. tation of the methods described previously (Mo¨lle et al., 2002; Ferrarelli et al., 2007). The raw EEG signal was averaged over the same four nonResults overlapping areas of the scalp as above, bandpass filtered between 11 and Association between properties of slow waves and white 16 Hz, and transformed into the envelope of the Hilbert transform (Fig. matter diffusion 1C,D). A spindle was detected when the signal envelope remained above Slow waves were automatically detected and their properties the threshold of 3 ␮V for a duration of 0.5–2 s. These two absolute criteria allow for a constant detection threshold that is consistent quantified by three parameters (Fig. 1 A, B): (1) the mean ampliacross subjects and simplifies the interpretation of between-subject tude of the negative peak, (2) the mean rising slope (Riedner et variability. al., 2007), and (3) the slow-wave density (the number of slow For each spindle we calculated the spindle power, which is defined as waves per minute). The values for these parameters were calcuthe area outlined by the envelope above the threshold (Fig. 1D). Furtherlated separately for sleep stages N2 and N3 (Iber et al., 2007) and more, we calculated spindle density as the number of spindles per minare shown in Table 1. ute. For each participant, these values were then averaged separately for The mean steepness of the rising slope of slow waves in stage sleep stages N2 and N3. N3 sleep was significantly positively correlated with AD, but not Diffusion tensor imaging. Diffusion properties in the white matter were FA, in several parts of the structural skeleton of major white assessed using DTI measures of fractional anisotropy and axial diffusivity matter bundles (corrected p value ⫽ 0.006 at the peak voxel, see from diffusion-weighted images. One image at b ⫽ 0 s/mm 2 and 64



tracts in the left middle temporal lobe (corrected p value ⫽ 0.04 at the peak voxel for both contrasts, details are reported in Table 2).

Table 1. Summary of the parameters for the slow waves, including mean and standard deviation Stage

Parameter

Mean

SD

PFA

PAD

N2

Slow wave density (SW per minute) Mean negative peak (␮V) Mean rising slope (␮V/s)

2.61 ⫺92.24 324.41

⫾0.56 ⫾5.96 ⫾31.04

0.07 0.26 0.59

0.11 0.68 0.13

N3

Slow wave density (SW per minute) Mean negative peak (␮V) Mean rising slope (␮V/s)

11.26 ⫺87.52 357.84

⫾3.77 ⫾4.23 ⫾28.05

0.27 0.56 0.14

0.11 0.93 0.006

The rightmost columns indicate the corrected p value for the most significant voxel of the correlation between the slow wave parameter (row) and the DTI parameter (column). A significant correlation between AD and the rising slope of slow waves during sleep stage N3 was observed, indicating that axial diffusivity was higher in individuals with a steeper slow wave slope (see Fig. 2).

Table 2. Summary of the parameters for the spindles, including mean and standard deviation Stage

Parameter

Mean

SD

PFA

PAD

N2

Spindle density (spindles per minute) Mean spindle power (␮V)

9.09 4.22

⫾2.84 ⫾0.55

0.33 0.50

0.04 0.02

N3

Spindle density (spindles per minute) Mean spindle power (␮V)

5.63 3.54

⫾3.29 ⫾0.50

0.31 0.58

0.04 0.11

The rightmost columns indicate the corrected p value for the most significant voxel of the correlation between the spindle parameter (row) and the DTI parameter (column). A significant correlation between AD and spindle power during sleep stage N2 was observed, indicating that axial diffusivity was higher in individuals with a stronger spindle power (Fig. 3). Furthermore, axial diffusivity was higher in individuals with a higher spindle density in sleep stages N2 and N3 in the temporal branch of the right inferior fronto-occipital fascicle (Fig. 3).

Fig. 2). A steeper rising slope was associated with higher AD in the forceps minor and neighboring white matter structures, including part of the uncinate fascicle, the genu and anterior part of the corpus callosum, the right superior longitudinal fascicle, and the temporal branch of the right inferior fronto-occipital fascicle. These areas comprise a mostly frontal interconnected network; major bundles in the parietal and occipital areas were below the significance threshold. To investigate which electrodes were the most informative in estimating the slow-wave slope, we calculated, at each significant voxel, the correlation between the slope and axial diffusivity for each of the four regions of interest and then selected the region of interest showing the highest correlation. The slow-wave slopes estimated from anterior electrodes showed the highest correlation with axial diffusivity, especially in the anterior white matter bundles (color-coded in Fig. 2). Neither AD nor FA correlated significantly with slow-wave density or the mean negative peak amplitude. No correlations were found for stage N2 sleep slow-wave parameters (details are reported in Table 1). Association between properties of spindles and white matter diffusion Spindles were automatically detected and their properties quantified by two parameters (Fig. 1 C,D): (1) the mean spindle power and (2) the mean spindle density (the mean number of spindles per minute). These parameters were calculated separately for sleep stages N2 and N3 and are shown in Table 2. Individuals with a higher spindle power during sleep stage N2 had a higher axial diffusivity in the forceps minor, the genu and anterior part of the corpus callosum, the temporal branch of the inferior longitudinal fascicles, and the areas surrounding and intrinsic to the thalamus, including parts of the internal and external capsule (corrected p value ⫽ 0.02 at the peak voxel; Fig. 3). Higher spindle density, calculated during N2 and N3, positively correlated with higher AD in specific parts of the white matter

Discussion Individual differences in the steepness of the rising slope of the slow wave during stage N3 sleep and in spindle power during stage N2 sleep are associated with axial diffusivity in partially overlapping major frontal white matter tracts that include the forceps minor, parts of the uncinate fascicle, and the anterior corpus callosum. Connected to this white matter backbone are the temporal branch of the inferior fronto-occipital fascicle and parts of the superior longitudinal fascicle, tracts in which AD was correlated with the profiles of both slow waves and spindles. In addition, participants with higher spindle power had higher axial diffusivity in white matter tracts within and surrounding the thalamus, including the internal and external capsule. Of note is that the slow-wave profile is mainly associated with frontal white matter, whereas the spindle profile is associated with white matter in a more extensive network that also includes subcortical regions. This difference is in line with the known anatomical and cellular substrate of both types of oscillatory events and with evidence from intracortical and scalp EEG, fMRI, and MRI studies (Schabus et al., 2007; Dang-Vu et al., 2008; Dang-Vu et al., 2010b). These results indicate that specific parameters of slow-wave and spindle activity do not only depend on the dynamics of the connectivity at the synaptic level (Vyazovskiy et al., 2007; Vyazovskiy et al., 2008; Vyazovskiy et al., 2009), but also reflect the more stable network properties that are determined by the structural backbone of white matter tracts. The correlation between the slow-wave slope and AD in the superior longitudinal fascicle and the temporal fascicle was more pronounced in the right than in the left hemisphere. It is unclear whether this lateralization is due to an underlying anatomical difference between hemispheres, such as the unilateral increase in diffusivity in the arcuate fasciculus for language-related skills and handedness (Glasser and Rilling, 2008; Lebel and Beaulieu, 2009; Propper et al., 2010), or to the predominance of slow-wave activity on the right hemisphere (Goldstein et al., 1972; Roth et al., 1999; Sekimoto et al., 2000). This latter fact or local magnetic field in homogeneities might reduce the signal-to-noise ratio, thereby weakening the correlation in the left hemisphere. In fact, a more liberal threshold reveals significant correlation in the homologous regions in the left hemisphere. Role of axial diffusivity in slow-wave synchronization and propagation Axial diffusivity is considered a marker of axonal integrity (Budde et al., 2007; Budde and Song, 2010). The present findings indicate that axonal integrity in frontal parts of the white matter bundles determines the steepness of the slow-wave slope. Slope steepness is thought to reflect the extent to which the neuronal network synchronizes activity during the transition from the downstate of hyperpolarization and neuronal silence to the upstate of depolarization and discharge of action potentials (Steriade et al., 1993; Vyazovskiy et al., 2007). Taken together, these observations suggest that higher AD in selected brain regions might contribute to the synchronization of the slow wave, resulting in a steeper slope. The effect of axonal integrity on slow-wave slope can be thought of in terms of known biophysical processes. We propose that while the intrinsic bistability of membrane potentials during slow-wave sleep eventually leads to some spontaneous isolated firing (Destexhe et al., 1999; Chauvette et al., 2010), white matter


J. Neurosci., January 2, 2013 • 33(1):227–233 • 231

travel over the cortex (Massimini et al., 2004; Botella-Soler et al., 2012). Although they can emerge in any part of the cortex, they are more likely to first appear in the anterior cingulate cortex and insula and then propagate along the anteroposterior axis over medial cortical areas (Murphy et al., 2009). Activity in the medial prefrontal cortex during slow waves (Dang-Vu et al., 2008) and spindles (Schabus et al., 2007) has been described using combined EEG-fMRI. Furthermore, the volume of the corpus callosum, measured with MRI, is associated with an individual’s power spectral density in slow-wave frequency band (Buchmann et al., 2011). By employing DTI, not only do we confirm the association between slow waves and frontal regions, such as the prefrontal and cingulate cortex, but we also demonstrate that the profiles of slow waves are specifically affected by the microstructural properties of the main interregional white matter bundles. Role of axial diffusivity in spindle profiles A similar mechanism might contribute to the enhanced spindle power in individuals with high AD in the areas surrounding and within the thalamus and the internal and external capsule, including the corona radiata. These areas are thought to be Figure 3. White matter tracts where axial diffusivity was higher in individuals with a stronger spindle power in sleep stage N2. involved in driving spindles by way of Some of the white matter tracts in which AD was associated with spindle power partially overlap with the tracts plotted in Figure projections to and from the cortex (Ste2, such as the forceps minor, the genu and anterior part of the corpus callosum (G CC and aCC, respectively), the superior longitudinal fascicle (sLF), and the inferior fronto-occipital fascicle (iFOF). The network of white matter tracts whose AD was correlated riade et al., 1985; Steriade, 2003a, Steriade, with spindle power was more extended in the temporal lobe than the network where AD was associated with the slow-wave slope 2003b). Corticothalamic and thalamocorand included the inferior longitudinal fascicle (iLF) as well. Other tracts more selectively associated with spindle power include the tical fibers run through the superior and white matter regions around the thalamus, such as the internal and external capsule (iC and eC, respectively) and the internal anterior thalamic peduncles that are lomedullary lamina (iML) of the thalamus. In addition, tracts in the middle temporal lobe, colored in blue, were associated with cated on the medial part of the internal spindle density in sleep stage N2. capsule and in the anterior limb of the internal capsule, respectively. Our findings tracts with higher axial diffusivity will promote a more synchrosuggest that a stronger AD in the thalamocortical loop that supnized switching into the upstate over multiple neurons, as reports the generation of spindles (Steriade et al., 1985; Steriade, flected in the steepness of the slow-wave slope. Crucially, 2003b) could enhance their power. Similarly, high AD in the according to our interpretation, strong axonal integrity of the frontal white matter network including the corpus callosum may bundles that connect prefrontal and limbic regions will catalyze facilitate interhemispheric synchronization (Achermann and the probability that sparse neuronal firing can ignite a widespread Borbe´ly, 1998), which is of functional relevance given its associavalanche of slow-wave activity across the neuronal network. ation with memory performance (Cantero et al., 2002; Mo¨lle et White matter tracts of particular importance for the expresal., 2004). sion of slow waves comprise, among others, the forceps minor, the genu of the corpus callosum, and the uncinate fascicle. These Temporal lobe bundles are ideal candidates to support the spreading of slow Expression profiles of slow waves and spindles consistently corwaves over widespread cortical areas across hemispheres (Fig. 2). related with axial diffusivity not only in frontal white matter, but Once a slow wave has been triggered in frontal regions, activity also in two major white matter tracts of the temporal lobe: the spreads toward the posterior part of the cortex along the cingulate temporal branch of the inferior fronto-occipital fascicle and the cortex and the longitudinal fascicles (Botella-Soler et al., 2012). inferior longitudinal fascicle. The association of AD in the temConduction over the corpus callosum might be involved in the poral lobe and spindle activity was particularly strong, as judged enhanced interhemispheric connectivity observed during sleep from the spindle power correlation in both hemispheres and the (Nir et al., 2008). spindle density correlation during both stages N2 and N3 sleep. This interpretation is supported by intracranial and highWhereas spindles are generated in thalamocortical networks density EEG experiments that have demonstrated that slow waves (Steriade, 2003b), their enhanced expression in individuals with are highly dynamic and are not confined to any brain region but higher AD in the inferior fronto-occipital fascicule and the infe-


rior longitudinal fascicle may be of particular functional relevance, because these fascicles connect the temporal lobes with widespread cortical areas. Several studies have shown that memory performance and intelligence are positively associated with spindle density (Gais et al., 2002; Clemens et al., 2005; Morin et al., 2008; Fogel and Smith, 2011) and also with white matter connectivity in the temporal lobe (Persson et al., 2006; Sasson et al., 2010; Kantarci et al., 2011). Individual differences in white matter connectivity may mediate the association of spindle expression with memory performance and intelligence. Thus, these different structural and neurophysiological traits might be related and result from a common genotype. An intriguing additional possible explanation is that the variability in AD in these white matter bundles may in part reflect plasticity-induced changes associated with prolonged activation. There is accumulating evidence of an increase in structural connectivity, assessed with DTI, in response to training (Keller and Just, 2009; Scholz et al., 2009; Landi et al., 2011). A possible mechanism underlying such strengthening might be the recently discovered changes in myelination in association with glutamate release onto oligodendrocytes along the axon during action potentials (Kukley et al., 2007; Araque and Navarrete, 2011; Wake et al., 2011). Based on these findings, a prolonged and repeated high level of spindle activity, for example triggered by extensive learning (Fogel and Smith, 2011) and measured in this study with high spindle power and high spindle density, will conceivably strengthen an individual’s temporal lobe white matter connectivity. Such increased connectivity might facilitate the synchronization between spindles and hippocampal ripples (Clemens et al., 2007) that underpins the communication between neocortex and hippocampus (Siapas and Wilson, 1998; Clemens et al., 2011) and bolsters memory consolidation during sleep (Diekelmann and Born, 2010). Conclusion Our findings indicate that interindividual variability in regionally specific white matter tracts contribute to the trait-like expression of slow waves and spindles, the most characteristic oscillations of the sleeping brain. The individual profiles of both these types of sleep oscillations correlate strongly with white matter axial diffusivity in frontal regions. In addition to these areas, spindle profiles show a marked association with axial diffusivity in subcortical white matter, including tracts within and around the thalamus. The regional distribution is in line with previous work on the neural sources of these oscillations. The findings complement earlier work which indicated that the expression profile of sleep oscillations reflect dynamic properties of network connectivity as a consequence of the relatively fast changes in synaptic strength (Tononi and Cirelli, 2006; Poe et al., 2010). White matter tracts, representing the structural backbone of neuronal network connectivity, appear equally important in determining the expression of sleep oscillations, reflecting the more stable and trait-like individual characteristics.

References Achermann P, Borbe´ly AA (1998) Temporal evolution of coherence and power in the human sleep electroencephalogram. J Sleep Res 7 [Suppl 1]:36 – 41. CrossRef Medline Araque A, Navarrete M (2011) Development: electrically driven insulation in the central nervous system. Science 333:1587–1588. CrossRef Medline Beaulieu C (2010) What makes diffusion anisotropic in the nervous system? In: Diffusion MRI: theory, methods, and application (Jones DK, ed), pp 92–109. Oxford: Oxford UP. Bo´dizs R, Kis T, La´za´r AS, Havra´n L, Rigo´ P, Clemens Z, Hala´sz P (2005)

Piantoni et al. • White Matter Diffusion Affects Sleep Oscillations Prediction of general mental ability based on neural oscillation measures of sleep. J Sleep Res 14:285–292. CrossRef Medline Botella-Soler V, Valderrama M, Cre´pon B, Navarro V, Le Van Quyen M (2012) Large-scale cortical dynamics of sleep slow waves. PLoS One 7:e30757. CrossRef Medline Buchmann A, Kurth S, Ringli M, Geiger A, Jenni OG, Huber R (2011) Anatomical markers of sleep slow wave activity derived from structural magnetic resonance images. J Sleep Res 20:506 –513. CrossRef Medline Buckelmu¨ller J, Landolt H-P, Stassen HH, Achermann P (2006) Trait-like individual differences in the human sleep electroencephalogram. Neuroscience 138:351–356. CrossRef Medline Budde MD, Song S-K (2010) Insights into diffusion tensor imaging from animal models of white matter pathology. In: Diffusion MRI: theory, methods, and application (Jones DK, ed), pp 670 –701. Oxford: Oxford UP. Budde MD, Kim JH, Liang H-F, Schmidt RE, Russell JH, Cross AH, Song S-K (2007) Toward accurate diagnosis of white matter pathology using diffusion tensor imaging. Magn Reson Med 57:688 – 695. CrossRef Medline Cantero JL, Atienza M, Salas RM, Dominguez-Marin E (2002) Effects of prolonged waking-auditory stimulation on electroencephalogram synchronization and cortical coherence during subsequent slow-wave sleep. J Neurosci 22:4702– 4708. Medline Chauvette S, Volgushev M, Timofeev I (2010) Origin of active states in local neocortical networks during slow sleep oscillation. Cereb Cortex 20: 2660 –2674. CrossRef Medline Clemens Z, Fabo´ D, Hala´sz P (2005) Overnight verbal memory retention correlates with the number of sleep spindles. Neuroscience 132:529 –535. CrossRef Medline Clemens Z, Mo¨lle M, Eross L, Barsi P, Hala´sz P, Born J (2007) Temporal coupling of parahippocampal ripples, sleep spindles and slow oscillations in humans. Brain 130:2868 –2878. CrossRef Medline Clemens Z, Mo¨lle M, Eross L, Jakus R, Ra´sonyi G, Hala´sz P, Born J (2011) Fine-tuned coupling between human parahippocampal ripples and sleep spindles. Eur J Neurosci 33:511–520. CrossRef Medline Dang-Vu TT, Schabus M, Desseilles M, Albouy G, Boly M, Darsaud A, Gais S, Rauchs G, Sterpenich V, Vandewalle G, Carrier J, Moonen G, Balteau E, Degueldre C, Luxen A, Phillips C, Maquet P (2008) Spontaneous neural activity during human slow wave sleep. Proc Natl Acad Sci U S A 105: 15160 –15165. CrossRef Medline Dang-Vu TT, McKinney SM, Buxton OM, Solet JM, Ellenbogen JM (2010a) Spontaneous brain rhythms predict sleep stability in the face of noise. Curr Biol 20:R626 –R627. CrossRef Medline Dang-Vu TT, Schabus M, Desseilles M, Sterpenich V, Bonjean M, Maquet P (2010b) Functional neuroimaging insights into the physiology of human sleep. Sleep 33:1589 –1603. Medline De Gennaro L, Ferrara M, Vecchio F, Curcio G, Bertini M (2005) An electroencephalographic fingerprint of human sleep. Neuroimage 26:114 – 122. CrossRef Medline De Gennaro L, Marzano C, Fratello F, Moroni F, Pellicciari MC, Ferlazzo F, Costa S, Couyoumdjian A, Curcio G, Sforza E, Malafosse A, Finelli LA, Pasqualetti P, Ferrara M, Bertini M, Rossini PM (2008) The electroencephalographic fingerprint of sleep is genetically determined: a twin study. Ann Neurol 64:455– 460. CrossRef Medline Destexhe A, Contreras D, Steriade M (1999) Spatiotemporal analysis of local field potentials and unit discharges in cat cerebral cortex during natural wake and sleep states. J Neurosci 19:4595– 4608. Medline Diekelmann S, Born J (2010) The memory function of sleep. Nat Rev Neurosci 11:114 –126. CrossRef Medline Dijk D-J (2006) Sleep of aging women and men: back to basics. Sleep 29:12–13. Medline Dijk D-J (2011) Phenotyping sleep. J Sleep Res 20:495. CrossRef Medline Esser SK, Hill SL, Tononi G (2007) Sleep homeostasis and cortical synchronization: I. Modeling the effects of synaptic strength on sleep slow waves. Sleep 30:1617–1630. Medline Ferrarelli F, Huber R, Peterson MJ, Massimini M, Murphy M, Riedner BA, Watson A, Bria P, Tononi G (2007) Reduced sleep spindle activity in schizophrenia patients. Am J Psychiatry 164:483– 492. CrossRef Medline Finelli LA, Achermann P, Borbe´ly AA (2001) Individual “fingerprints” in human sleep EEG topography. Neuropsychopharmacology 25:S57– 62. CrossRef Medline Fogel SM, Smith CT (2011) The function of the sleep spindle: a physiolog-

Piantoni et al. • White Matter Diffusion Affects Sleep Oscillations ical index of intelligence and a mechanism for sleep-dependent memory consolidation. Neurosci Biobehav Rev 35:1154 –1165. CrossRef Medline Gais S, Mo¨lle M, Helms K, Born J (2002) Learning-dependent increases in sleep spindle density. J Neurosci 22:6830 – 6834. Medline Gander P, Signal L, Van Dongen HP, Muller D, Van Den Berg M (2010) Stable inter-individual differences in slow-wave sleep during nocturnal sleep and naps. Sleep Biol Rhythms 8:239 –244. CrossRef Glasser MF, Rilling JK (2008) DTI tractography of the human brain’s language pathways. Cereb Cortex 18:2471–2482. CrossRef Medline Goldstein L, Stoltzfus NW, Gardocki JF (1972) Changes in interhemispheric amplitude relationships in the EEG during sleep. Physiol Behav 8:811– 815. CrossRef Medline Iber C, Ancoli-Israel S, Chesson AL, Quan SF (2007) The AASM manual for the scoring of sleep and associated events: rules, terminology, and technical specifications, Ed 1. Westchester, IL: American Academy of Sleep Medicine. Jones DK ed. (2010) Diffusion MRI: theory, methods, and application. Oxford, UK: Oxford UP. Kantarci K, Senjem ML, Avula R, Zhang B, Samikoglu AR, Weigand SD, Przybelski SA, Edmonson HA, Vemuri P, Knopman DS, Boeve BF, Ivnik RJ, Smith GE, Petersen RC, Jack CR Jr (2011) Diffusion tensor imaging and cognitive function in older adults with no dementia. Neurology 77: 26 –34. CrossRef Medline Keller TA, Just MA (2009) Altering cortical connectivity: remediationinduced changes in the white matter of poor readers. Neuron 64:624 – 631. CrossRef Medline Kerkhof GA (1991) Differences between morning-types and evening-types in the dynamics of EEG slow wave activity during night sleep. Electroencephalogr Clin Neurophysiol 78:197–202. CrossRef Medline Kukley M, Capetillo-Zarate E, Dietrich D (2007) Vesicular glutamate release from axons in white matter. Nat Neurosci 10:311–320. CrossRef Medline Landi SM, Baguear F, Della-Maggiore V (2011) One week of motor adaptation induces structural changes in primary motor cortex that predict long-term memory one year later. J Neurosci 31:11808 –11813. CrossRef Medline Lebel C, Beaulieu C (2009) Lateralization of the arcuate fasciculus from childhood to adulthood and its relation to cognitive abilities in children. Hum Brain Mapp 30:3563–3573. CrossRef Medline Leclercq Y, Schrouff J, Noirhomme Q, Maquet P, Phillips C (2011) fMRI artefact rejection and sleep scoring toolbox. Comput Intell Neurosci 2011:598206. CrossRef Medline Massimini M, Huber R, Ferrarelli F, Hill S, Tononi G (2004) The sleep slow oscillation as a traveling wave. J Neurosci 24:6862– 6870. CrossRef Medline Mo¨lle M, Marshall L, Gais S, Born J (2002) Grouping of spindle activity during slow oscillations in human non-rapid eye movement sleep. J Neurosci 22:10941–10947. Medline Mo¨lle M, Marshall L, Gais S, Born J (2004) Learning increases human electroencephalographic coherence during subsequent slow sleep oscillations. Proc Natl Acad Sci U S A 101:13963–13968. CrossRef Medline Morin A, Doyon J, Dostie V, Barakat M, Hadj Tahar A, Korman M, Benali H, Karni A, Ungerleider LG, Carrier J (2008) Motor sequence learning increases sleep spindles and fast frequencies in post-training sleep. Sleep 31:1149 –1156. Medline Murphy M, Riedner BA, Huber R, Massimini M, Ferrarelli F, Tononi G (2009) Source modeling sleep slow waves. Proc Natl Acad Sci U S A 106: 1608 –1613. CrossRef Medline Nir Y, Mukamel R, Dinstein I, Privman E, Harel M, Fisch L, Gelbard-Sagiv H, Kipervasser S, Andelman F, Neufeld MY, Kramer U, Arieli A, Fried I, Malach R (2008) Interhemispheric correlations of slow spontaneous neuronal fluctuations revealed in human sensory cortex. Nat Neurosci 11:1100 –1108. CrossRef Medline Persson J, Nyberg L, Lind J, Larsson A, Nilsson L-G, Ingvar M, Buckner RL (2006) Structure-function correlates of cognitive decline in aging. Cereb Cortex 16:907–915. CrossRef Medline Poe GR, Walsh CM, Bjorness TE (2010) Cognitive neuroscience of sleep. Prog Brain Res 185:1–19. CrossRef Medline Propper RE, O’Donnell LJ, Whalen S, Tie Y, Norton IH, Suarez RO, Zollei L, Radmanesh A, Golby AJ (2010) A combined fMRI and DTI examination of functional language lateralization and arcuate fasciculus structure:

J. Neurosci., January 2, 2013 • 33(1):227–233 • 233 Effects of degree versus direction of hand preference. Brain Cogn 73:85– 92. CrossRef Medline Rechtschaffen A, Kales A (1968) A manual of standardized terminology, techniques and scoring system for sleep stages of human subjects. Los Angeles: UCLA Brain Information Service/Brain Research Institute. Riedner BA, Vyazovskiy VV, Huber R, Massimini M, Esser S, Murphy M, Tononi G (2007) Sleep homeostasis and cortical synchronization: III. A high-density EEG study of sleep slow waves in humans. Sleep 30:1643–1657. Medline Roth C, Achermann P, Borbe´ly AA (1999) Frequency and state specific hemispheric asymmetries in the human sleep EEG. Neurosci Lett 271: 139 –142. CrossRef Medline Rusterholz T, Achermann P (2011) Topographical aspects in the dynamics of sleep homeostasis in young men: individual patterns. BMC Neurosci 12:84. CrossRef Medline Sarkar S, Katshu MZUH, Nizamie SH, Praharaj SK (2010) Slow wave sleep deficits as a trait marker in patients with schizophrenia. Schizophr Res 124:127–133. CrossRef Medline Sasson E, Doniger GM, Pasternak O, Assaf Y (2010) Structural correlates of memory performance with diffusion tensor imaging. Neuroimage 50:1231–1242. CrossRef Medline Schabus M, Dang-Vu TT, Albouy G, Balteau E, Boly M, Carrier J, Darsaud A, Degueldre C, Desseilles M, Gais S, Phillips C, Rauchs G, Schnakers C, Sterpenich V, Vandewalle G, Luxen A, Maquet P (2007) Hemodynamic cerebral correlates of sleep spindles during human non-rapid eye movement sleep. Proc Natl Acad Sci U S A 104:13164 –13169. CrossRef Medline Schabus M, Hoedlmoser K, Pecherstorfer T, Anderer P, Gruber G, Parapatics S, Sauter C, Kloesch G, Klimesch W, Saletu B, Zeitlhofer J (2008) Interindividual sleep spindle differences and their relation to learning-related enhancements. Brain Res 1191:127–135. CrossRef Medline Scholz J, Klein MC, Behrens TEJ, Johansen-Berg H (2009) Training induces changes in white matter architecture. Nat Neurosci 12:1370 –1371. CrossRef Medline Sekimoto M, Kato M, Kajimura N, Watanabe T, Takahashi K, Okuma T (2000) Asymmetric interhemispheric delta waves during all-night sleep in humans. Clin Neurophysiol 111:924 –928. CrossRef Medline Siapas AG, Wilson MA (1998) Coordinated interactions between hippocampal ripples and cortical spindles during slow-wave sleep. Neuron 21:1123–1128. CrossRef Medline Smith SM, Nichols TE (2009) Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage 44:83–98. CrossRef Medline Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TEJ, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage [23 Suppl] 1:S208 –S219. CrossRef Medline Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE, Watkins KE, Ciccarelli O, Cader MZ, Matthews PM, Behrens TEJ (2006) Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage 31:1487–1505. CrossRef Medline Steriade M (2003a) The corticothalamic system in sleep. Front Biosci 8:d878 –99. CrossRef Medline Steriade M (2003b) Neuronal substrates of sleep and epilepsy. Cambridge: Cambridge University UP. Steriade M, Descheˆnes M, Domich L, Mulle C (1985) Abolition of spindle oscillations in thalamic neurons disconnected from nucleus reticularis thalami. J Neurophysiol 54:1473–1497. Medline Steriade M, Nun˜ez A, Amzica F (1993) A novel slow (⬍1 Hz) oscillation of neocortical neurons in vivo: depolarizing and hyperpolarizing components. J Neurosci 13:3252–3265. 8340806 Tononi G, Cirelli C (2006) Sleep function and synaptic homeostasis. Sleep Med Rev 10:49 – 62. CrossRef Medline Vyazovskiy VV, Riedner BA, Cirelli C, Tononi G (2007) Sleep homeostasis and cortical synchronization: II. A local field potential study of sleep slow waves in the rat. Sleep 30:1631–1642. Medline Vyazovskiy VV, Cirelli C, Pfister-Genskow M, Faraguna U, Tononi G (2008) Molecular and electrophysiological evidence for net synaptic potentiation in wake and depression in sleep. Nat Neurosci 11:200 –208. CrossRef Medline

233a • J. Neurosci., January 2, 2013 • 33(1):227–233 Vyazovskiy VV, Olcese U, Lazimy YM, Faraguna U, Esser SK, Williams JC, Cirelli C, Tononi G (2009) Cortical firing and sleep homeostasis. Neuron 63:865– 878. CrossRef Medline Wake H, Lee PR, Fields RD (2011) Control of local protein synthesis and initial events in myelination by action potentials. Science 333:1647–1651. CrossRef Medline

Piantoni et al. • White Matter Diffusion Affects Sleep Oscillations Werth E, Achermann P, Dijk DJ, Borbe´ly AA (1997) Spindle frequency activity in the sleep EEG: individual differences and topographic distribution. Electroencephalogr Clin Neurophysiol 103:535–542. CrossRef Medline