During sleep, the thalamus generates a characteristic pattern of transient, 11-15 Hz sleep spindle oscillations, which synchronize the cortex through large-scale thalamocortical loops. Spindles have been increasingly demonstrated to be critical for sleep-dependent consolidation of memory, but the specific neural mechanism for this process remains unclear. We show here that cortical spindles are spatiotemporally organized into circular wave-like patterns, organizing neuronal activity over tens of milliseconds, within the timescale for storing memories in large-scale networks across the cortex via spike-time dependent plasticity. These circular patterns repeat over hours of sleep with millisecond temporal precision, allowing reinforcement of the activity patterns through hundreds of reverberations. These results provide a novel mechanistic account for how global sleep oscillations and synaptic plasticity could strengthen networks distributed across the cortex to store coherent and integrated memories.

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.

Human intracranial electroencephalography (iEEG) recordings provide data with much greater spatiotemporal precision than is possible from data obtained using scalp EEG, magnetoencephalography (MEG), or functional MRI. Until recently, the fusion of anatomical data (MRI and computed tomography (CT) images) with electrophysiological data and their subsequent analysis have required the use of technologically and conceptually challenging combinations of software. Here, we describe a comprehensive protocol that enables complex raw human iEEG data to be converted into more readily comprehensible illustrative representations. The protocol uses an open-source toolbox for electrophysiological data analysis (FieldTrip). This allows iEEG researchers to build on a continuously growing body of scriptable and reproducible analysis methods that, over the past decade, have been developed and used by a large research community. In this protocol, we describe how to analyze complex iEEG datasets by providing an intuitive and rapid approach that can handle both neuroanatomical information and large electrophysiological datasets. We provide a worked example using an example dataset. We also explain how to automate the protocol and adjust the settings to enable analysis of iEEG datasets with other characteristics. The protocol can be implemented by a graduate student or postdoctoral fellow with minimal MATLAB experience and takes approximately an hour to execute, excluding the automated cortical surface extraction.