The overall goal of the proposed research is to determine how brain events at the biophysical level in the human brain influence the macroscopic recordings outside the skull. The specific focus is on sleep spindles, the best studied sleep rhythm, for which we will obtain simultaneous EEG, MEG and depth electrode recordings. These empirical data will be analyzed to determine the laminar sources of the currents within the cortex associated with spindles and these results will be matched to predictions from computational models. We will develop three types of related network models to cover a wide range of spatial scales. The working hypotheses are, first, that spindles are generated by multiple loosely-coupled cortical regions through rhythmic activation of thalamocortical afferents, onto superficial cortical layers from the matrix system, and onto middle layers from the core system, and second, that spindle oscillations in restricted thalamocortical domains in the core system are spread via the matrix system. Empirical Specific Aims and Hypotheses 1) Record MEG and EEG simultaneously during sleep spindles. Hypothesis: MEG and EEG will record very different activity patterns at the sensor level and in their inferred sources. Poor correlation will be apparent as inconstant phase and amplitude relations within spindle discharges, as well as only loose correlations as to when spindles occur. 2) Record from intracranial EEG macro-electrodes during sleep spindles, with simultaneous scalp EEG and MEG recordings. Hypothesis: Recordings from different cortical generators will be only loosely coupled with each other, or with scalp EEG, thus supporting the general view suggested by MEG. 3) Record from intracranial microelectrode arrays during sleep spindles. Hypothesis: Neuronal generator currents vary during each spindle between superficial and middle layers. 4) Reconstruct supragranular and infragranular pyramidal cells from association cortex in humans, determine their laminar distribution, and estimate the likely termination zones of core and matrix thalamocortical projections. Hypothesis: Significant anatomical differences will be found between human association cortex and the rodent sensory cortex. Modeling Specific Aims and Hypotheses 1) Construct accurate, realistic EEG/MEG forward solutions based on cortical reconstruction from structural MRI. Hypothesis: Basic parameters of extracranial EEG/MEG can be replicated by modeling multiple thalamocortical domains with varying synchrony between domains. The values of the fit parameters to EEG/MEG will match those inferred from intracranial macroelectrode recordings. 2) Analyze the recorded cortical Current Source Density (CSD) using Principal Components Analysis (PCA). Model the CSD patterns expected from the matrix and core thalamocortical afferents using reconstructions of supragranular and infragranular pyramidal cells, their population distribution, and terminations of matrix and core afferents. Hypothesis: The main spatiotemporal CSD components contributing to the spindle identified with PCA will correspond to the CSD components modeled to result from activation of matrix and core thalamocortical afferents. 3) Construct neuronal models based on Hodgkin-Huxley ionic currents that include cortical cells, thalamic reticular nuclear cells, matrix and core thalamic relay cells, and their interconnections in a minicolumn. Hypothesis: The spatiotemporal patterns of currents predicted in different cortical layers will match those obtained using the methods described in modeling aim 2. 4) Scale up the cortical minicolumn network model using simplified neural models to a spatially accurate cortical model that can generate EEG/MEG patterns. Hypothesis: The predicted EEG and MEG from the model will match those observed in recordings, with the involvement of matrix vs core thalamocortical systems corresponding to the parameters derived in modeling aim 1. 5) Develop a statistical model that matches the properties of the scaled up model and analyze it with methods from statistical physics. Hypothesis: The matrix system controls the effective cortical connectivity and the core system controls the local correlation length. Collaborative research Recordings from humans will be performed at MGH (Cash), with analysis and modeling at UCSD (Halgren and Sejnowski). These 3 teams of PIs, students and postdoctoral fellows will interact on a daily basis during the research and will meet formally at least once a year to assess progress and plan new experiments.

National Institute of Health (NIH)
National Institute of Biomedical Imaging and Bioengineering (NIBIB)
Research Project (R01)
Project #
Application #
Study Section
Special Emphasis Panel (ZRG1-IFCN-B (50))
Program Officer
Peng, Grace
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of California San Diego
Biostatistics & Other Math Sci
Schools of Arts and Sciences
La Jolla
United States
Zip Code
Mak-McCully, Rachel A; Rolland, Matthieu; Sargsyan, Anna et al. (2017) Coordination of cortical and thalamic activity during non-REM sleep in humans. Nat Commun 8:15499
Tele?czuk, Bartosz; Dehghani, Nima; Le Van Quyen, Michel et al. (2017) Local field potentials primarily reflect inhibitory neuron activity in human and monkey cortex. Sci Rep 7:40211
Lainscsek, Claudia; Muller, Lyle E; Sampson, Aaron L et al. (2017) Analytical Derivation of Nonlinear Spectral Effects and 1/f Scaling Artifact in Signal Processing of Real-World Data. Neural Comput 29:2004-2020
Das, Anup; Sampson, Aaron L; Lainscsek, Claudia et al. (2017) Interpretation of the Precision Matrix and Its Application in Estimating Sparse Brain Connectivity during Sleep Spindles from Human Electrocorticography Recordings. Neural Comput 29:603-642
Piantoni, Giovanni; Halgren, Eric; Cash, Sydney S (2017) Spatiotemporal characteristics of sleep spindles depend on cortical location. Neuroimage 146:236-245
Piantoni, Giovanni; Halgren, Eric; Cash, Sydney S (2016) The Contribution of Thalamocortical Core and Matrix Pathways to Sleep Spindles. Neural Plast 2016:3024342
Dehghani, Nima; Peyrache, Adrien; Telenczuk, Bartosz et al. (2016) Dynamic Balance of Excitation and Inhibition in Human and Monkey Neocortex. Sci Rep 6:23176
Wei, Yina; Krishnan, Giri P; Bazhenov, Maxim (2016) Synaptic Mechanisms of Memory Consolidation during Sleep Slow Oscillations. J Neurosci 36:4231-47
Muller, Lyle; Piantoni, Giovanni; Koller, Dominik et al. (2016) Rotating waves during human sleep spindles organize global patterns of activity that repeat precisely through the night. Elife 5:
Krishnan, Giri P; Chauvette, Sylvain; Shamie, Isaac et al. (2016) Cellular and neurochemical basis of sleep stages in the thalamocortical network. Elife 5:

Showing the most recent 10 out of 34 publications