Hz Oscillations Research Articles

Modulations of Cortical Power and Connectivity in Alpha and Beta Bands during the Preparation of Reaching Movements

Planning goal-directed movements towards different targets is at the basis of common daily activities (e.g., reaching), involving visual, visuomotor, and sensorimotor brain areas. Alpha (8-13 Hz) and beta (13-30 Hz) oscillations are modulated during movement preparation and are implicated in correct motor functioning. However, how brain regions activate and interact during reaching tasks and how brain rhythms are functionally involved in these interactions is still limitedly explored. Here, alpha and beta brain activity and connectivity during reaching preparation are investigated at EEG-source level, considering a network of task-related cortical areas. Sixty-channel EEG was recorded from 20 healthy participants during a delayed center-out reaching task and projected to the cortex to extract the activity of 8 cortical regions per hemisphere (2 occipital, 2 parietal, 3 peri-central, 1 frontal). Then, we analyzed event-related spectral perturbations and directed connectivity, computed via spectral Granger causality and summarized using graph theory centrality indices (in degree, out degree). Results suggest that alpha and beta oscillations are functionally involved in the preparation of reaching in different ways, with the former mediating the inhibition of the ipsilateral sensorimotor areas and disinhibition of visual areas, and the latter coordinating disinhibition of the contralateral sensorimotor and visuomotor areas.

The cerebellum regulates fear extinction through thalamo-prefrontal cortex interactions in male mice

Fear extinction is a form of inhibitory learning that suppresses the expression of aversive memories and plays a key role in the recovery of anxiety and trauma-related disorders. Here, using male mice, we identify a cerebello-thalamo-cortical pathway regulating fear extinction. The cerebellar fastigial nucleus (FN) projects to the lateral subregion of the mediodorsal thalamic nucleus (MD), which is reciprocally connected with the dorsomedial prefrontal cortex (dmPFC). The inhibition of FN inputs to MD in male mice impairs fear extinction in animals with high fear responses and increases the bursting of MD neurons, a firing pattern known to prevent extinction learning. Indeed, this MD bursting is followed by high levels of the dmPFC 4 Hz oscillations causally associated with fear responses during fear extinction, and the inhibition of FN-MD neurons increases the coherence of MD bursts and oscillations with dmPFC 4 Hz oscillations. Overall, these findings reveal a regulation of fear-related thalamo-cortical dynamics by the cerebellum and its contribution to fear extinction.

Random Tactile Noise Stimulation Reveals Beta-Rhythmic Impulse Response Function of the Somatosensory System.

Both passive tactile stimulation and motor actions result in dynamic changes in beta band (15-30 Hz Hz) oscillations over somatosensory cortex. Similar to alpha band (8-12 Hz) power decrease in the visual system, beta band power also decreases following stimulation of the somatosensory system. This relative suppression of α and β oscillations is generally interpreted as an increase in cortical excitability. Here, next to traditional single-pulse stimuli, we employed a random intensity continuous right index finger tactile stimulation (white noise), which enabled us to uncover an impulse response function (IRF) of the somatosensory system. Contrary to previous findings, we demonstrate a burst-like initial increase rather than decrease of beta activity following white noise stimulation (human participants, N = 18, 8 female). These β bursts, on average, lasted for three cycles and their frequency was correlated with resonant frequency of somatosensory cortex, as measured by a multi-frequency steady-state somatosensory evoked potential (SSSEP) paradigm. Furthermore, beta band bursts shared spectro-temporal characteristics with evoked and resting-state β oscillations. Taken together, our findings not only reveal a novel oscillatory signature of somatosensory processing that mimics the previously reported visual IRFs, but also point to a common oscillatory generator underlying spontaneous β bursts in the absence of tactile stimulation and phase-locked β bursts following stimulation, the frequency of which is determined by the resonance properties of the somatosensory system.SIGNIFICANCE STATEMENTThe investigation of the transient nature of oscillations has gained great popularity in recent years. The findings of bursting activity rather than sustained oscillations in the beta band has provided important insights into its role in movement planning, working memory, inhibition and reactivation of neural ensembles. In this study, we show that also in response to tactile stimulation the somatosensory system responds with ∼3 cycle oscillatory beta band bursts, whose spectro-temporal characteristics are shared with evoked and resting-state beta band oscillatory signatures of the somatosensory system. As similar bursts have been observed in the visual domain, these oscillatory signatures might reflect an important supramodal mechanism in sensory processing.

Theta-gamma phase-amplitude coupling in auditory cortex is modulated by language proficiency.

The coordination between the theta phase (3-7 Hz) and gamma power (25-35 Hz) oscillations (namely theta-gamma phase-amplitude coupling, PAC) in the auditory cortex has been proposed as an essential neural mechanism involved in speech processing. However, it has not been established how this mechanism is related to the efficiency with which a listener processes speech. Speech processing in a non-native language offers a useful opportunity to evaluate if theta-gamma PAC is modulated by the challenges imposed by the reception of speech input in a non-native language. The present study investigates how auditory theta-gamma PAC (recorded with magnetoencephalography) is modulated in both native and non-native speech reception. Participants were Spanish native (L1) speakers studying Basque (L2) at three different levels: beginner (Grade 1), intermediate (Grade 2), and advanced (Grade 3). We found that during L2 speech processing (i) theta-gamma PAC was more highly coordinated for intelligible compared to unintelligible speech; (ii) this coupling was modulated by proficiency in Basque being lower for beginners, higher for intermediate, and highest for advanced speakers (no difference observed in Spanish); (iii) gamma power did not differ between languages and groups. These findings highlight how the coordinated theta-gamma oscillatory activity is tightly related to speech comprehension: the stronger this coordination is, the more the comprehension system will proficiently parse the incoming speech input.

Detection of neuronal OFF periods as low amplitude neural activity segments

BackgroundDuring non-rapid eye movement sleep (NREM), alternating periods of synchronised high (ON period) and low (OFF period) neuronal activity are associated with high amplitude delta band (0.5–4 Hz) oscillations in neocortical electrophysiological signals termed slow waves. As this oscillation is dependent crucially on hyperpolarisation of cortical cells, there is an interest in understanding how neuronal silencing during OFF periods leads to the generation of slow waves and whether this relationship changes between cortical layers. A formal, widely adopted definition of OFF periods is absent, complicating their detection. Here, we grouped segments of high frequency neural activity containing spikes, recorded as multiunit activity from the neocortex of freely behaving mice, on the basis of amplitude and asked whether the population of low amplitude (LA) segments displayed the expected characteristics of OFF periods.ResultsAverage LA segment length was comparable to previous reports for OFF periods but varied considerably, from as short as 8 ms to > 1 s. LA segments were longer and occurred more frequently in NREM but shorter LA segments also occurred in half of rapid eye movement sleep (REM) epochs and occasionally during wakefulness. LA segments in all states were associated with a local field potential (LFP) slow wave that increased in amplitude with LA segment duration. We found that LA segments > 50 ms displayed a homeostatic rebound in incidence following sleep deprivation whereas short LA segments (< 50 ms) did not. The temporal organisation of LA segments was more coherent between channels located at a similar cortical depth.ConclusionWe corroborate previous studies showing neural activity signals contain uniquely identifiable periods of low amplitude with distinct characteristics from the surrounding signal known as OFF periods and attribute the new characteristics of vigilance-state-dependent duration and duration-dependent homeostatic response to this phenomenon. This suggests that ON/OFF periods are currently underdefined and that their appearance is less binary than previously considered, instead representing a continuum.

Peripheral blood flow estimated by laser doppler flowmetry provides additional information about sleep state beyond that provided by pulse rate variability.

Pulse rate variability (PRV), derived from Laser Doppler flowmetry (LDF) or photoplethysmography, has recently become widely used for sleep state assessment, although it cannot identify all the sleep stages. Peripheral blood flow (BF), also estimated by LDF, may be modulated by sleep stages; however, few studies have explored its potential for assessing sleep state. Thus, we aimed to investigate whether peripheral BF could provide information about sleep stages, and thus improve sleep state assessment. We performed electrocardiography and simultaneously recorded BF signals by LDF from the right-index finger and ear concha of 45 healthy participants (13 women; mean age, 22.5 ± 3.4years) during one night of polysomnographic recording. Time- and frequency-domain parameters of peripheral BF, and time-domain, frequency-domain, and non-linear indices of PRV and heart rate variability (HRV) were calculated. Finger-BF parameters in the time and frequency domains provided information about different sleep stages, some of which (such as the difference between N1 and rapid eye movement sleep) were not revealed by finger-PRV. In addition, finger-PRV patterns and HRV patterns were similar for most parameters. Further, both finger- and ear-BF results showed 0.2-0.3Hz oscillations that varied with sleep stages, with a significant increase in N3, suggesting a modulation of respiration within this frequency band. These results showed that peripheral BF could provide information for different sleep stages, some of which was complementary to the information provided by PRV. Furthermore, the combination of peripheral BF and PRV may be more advantageous than HRV alone in assessing sleep states and related autonomic nervous activity.

Optogenetic frequency scrambling of hippocampal theta oscillations dissociates working memory retrieval from hippocampal spatiotemporal codes

The precise temporal coordination of activity in the brain is thought to be fundamental for memory function. Inhibitory neurons in the medial septum provide a prominent source of innervation to the hippocampus and play a major role in controlling hippocampal theta (~8 Hz) oscillations. While pharmacological inhibition of medial septal neurons is known to disrupt memory, the exact role of septal inhibitory neurons in regulating hippocampal representations and memory is not fully understood. Here, we dissociate the role of theta rhythms in spatiotemporal coding and memory using an all-optical interrogation and recording approach. We find that optogenetic frequency scrambling stimulations abolish theta oscillations and modulate a portion of neurons in the hippocampus. Such stimulation decreased episodic and working memory retrieval while leaving hippocampal spatiotemporal codes intact. Our study suggests that theta rhythms play an essential role in memory but may not be necessary for hippocampal spatiotemporal codes.

Stimulus-evoked and resting-state alpha oscillations show a linked dependence on patterned visual experience for development

Persistent visual impairments after congenital blindness due to dense bilateral cataracts have been attributed to altered visual cortex development within a sensitive period. Occipital alpha (8–14 Hz) oscillations were found to be reduced after congenital cataract reversal, while participants performed visual motion tasks. However, it has been unclear whether reduced alpha oscillations were task-specific, or linked to impaired visual behavior in cataract-reversed individuals. Here, we compared resting-state and stimulus-evoked alpha activity between individuals who had been treated for dense bilateral congenital cataracts (CC, n = 13, mean duration of blindness = 11.0 years) and age-matched, normally sighted individuals (SC, n = 13). We employed the visual impulse response function, adapted from VanRullen and MacDonald (2012), to test for the characteristic alpha response to visual white noise. Participants observed white noise stimuli changing in luminance with equal power at frequencies between 0 and 30 Hz. Compared to SC individuals, CC individuals demonstrated a reduced likelihood of exhibiting an evoked alpha response. Moreover, stimulus-evoked alpha power was reduced and correlated with a corresponding reduction of resting-state alpha power in the same CC individuals. Finally, CC individuals with an above-threshold evoked alpha peak had better visual acuity than CC individual without an evoked alpha peak. Since alpha oscillations have been linked to feedback communication, we suggest that the concurrent impairment in resting-state and stimulus-evoked alpha oscillations indicates an altered interaction of top-down and bottom-up processing in the visual hierarchy, which likely contributes to incomplete behavioral recovery in individuals who experienced transient congenital blindness.

The Complexity of Electroencephalographic Signal Decreases during the Social Stress.

Social stress affects brain function. Trier social stress test (TSST) is a standard test to assess it. The study aimed to analyze the electroencephalographic (EEG) recording during and after TSST in healthy subjects. The EEG signals of 44 healthy men participating in the study were recorded in the control condition, during and after TSST and after 30 min of recovery. Salivary cortisol (SC) and the Emotional Visual Analog Scale (EVAS) score were measured in the control condition, after TSST, and after the recovery period. The false discovery rate correction was used to control the false positive of significance in EEG. In the comparison control condition, the SC and EVAS levels significantly increased after TSST. The relative Delta band frequency significantly increased during TSST. On the other hand, the Beta bands and, in less amount, the Theta and Gamma 1 (30-40 Hz) oscillations decreased, especially in the frontal region. The nonlinear features such as, approximate and spectral entropy, Katz fractal dimension behaved like Beta band oscillation. All changes returned to baseline after TSST except the increase of Katz in the F3 channel after the recovery period. Thus, stress on EEG increased low frequency (1-4 Hz), decreased high frequency (13-40 Hz), and complexity indices during TSST.

Mathematical model of the photoplethysmogram for testing methods of biological signals analysis

The purpose of this study was to develop a mathematical model of the photoplethysmogram, which can be used to test methods that introduce the instantaneous phases of the modulating signals. The model must reproduce statistical and spectral characteristics of the real photoplethysmogram, and explicitly incorporate the instantaneous phases of the modulating signals, so they can be used as a reference during testing. Methods. Anacrotic and catacrotic phases of the photoplethysmogram pulse wave were modeled as a sum of two density distributions for the skew normal distribution. The modulating signals were introduced as harmonic functions taken from the experimental instantaneous phases of the VLF (0.015...0.04 Hz), LF (0.04...0.15 Hz) and HF (0.15...0.4 Hz) oscillations in the real photoplethysmogram. The spectral power in the VLF, LF, and HF frequency ranges was calculated to compare the model and experimental data. Results. The model qualitatively reproduces the shape of the experimental photoplethysmogram pulse wave and shows less than 1% error when simulating the spectral properties of the signal. Conclusion. The proposed mathematical model can be used to test the methods for introduction of the instantaneous phases of the modulating signals in photoplethysmogram time-series.

Model-driven tACS in Alzheimer’s disease with personalized hybrid brain models

Abstract Animal models suggest that 40 Hz induced gamma oscillations reduce Aβ and p-Tau deposition, which is considered central to Alzheimer’s disease (AD) pathophysiology. Transcranial alternating current stimulation (tACS) has been proposed as a feasible option to induce gamma in AD patients and thus ease the cognitive burden of the disease. Here we investigate the potential of personalized whole-brain models to drive the stimulation in AD patients. We build a hybrid brain model that combines the (i) anatomical information, i.e., a parcellated brain, with subcortical parcels and their connections inferred from dMRI, (ii) functional data, i.e. fMRI-based functional connectivity and (iii) physiology based on Laminar Neural Mass Models (NMM). The NMM, personalized further with TMS-EEG data, can represent fast (gamma, 40-80 Hz) and slow (alpha/beta, 4-22 Hz) oscillations. Pathological disruptions have been incorporated in the most affected parcels based on FDG PET SUVR, which result in a decrease of gamma power in said parcels, consistent with findings in EEG and TMS-EEG. We have explored stimulation strategies by targeting small networks of parcels and specific regions, such as the prefrontal cortex (PFC). We show that the interplay of resonance phenomena and network effects can facilitate gamma entrainment and recover part of the gamma-power loss. This study is the first step into building model-driven optimization strategies for non-invasive stimulation in AD. Research Category and Technology and Methods Translational Research: 8. Transcranial Alternating Current Stimulation (tACS) Keywords: tACS, alzheimer's disease, computational neuroscience, resonance phenomena

Features of respiratory impedance in patients with respiratory tuberculosis with bronchial patency disorders

SummaryAim. To establish the features of respiratory im-pedance deviations detected by the method of forced oscillations in patients with tuberculosis of the respiratory organs and bronchial patency disorders. Маterials and methods. The parameters of spirometry and the method of forced oscillations were studied in 310 patients and in healthy subjects (group 3, n=22). At 8, 12, 16 Hz (oscillation frequencies), the respiratory impedance (RI) was measured: RI with calm breathing at inspiration — Rin; RI with calm breathing at exhalation — Rex; RI with respiratory volume — Rfo. In group 1 (n=229), the following subgroups were identified: 1A subgroup — 62 patients with initial obstructive disorders, 1B — 108 people with moderate, and 1C — 59 patients with significant obstruction. Patients without obstructive disorders comprised group 2 (n=81). Results. In group 1: the parameters of RI were higher than in groups 2 and 3, the frequency dependence of the RI was established. In subgroup 3B values of the parameters of RI, index of the frequency dependence of the RI, the number of RI indicators with pathological changes were higher than in 3A subgroup. Conclusion. In patients with tuberculosis of the respiratory organs with obstructive disorders an increase in respiratory impedance and the frequency dependence of respiratory impedance were observed. With the aggravation of obstructive disorders, the respiratory impedance and frequency dependence increase, pathologically altered impedance parameters are more often recorded

Time-Frequency Analysis of Somatosensory Evoked High-Frequency (600 Hz) Oscillations as an Early Indicator of Arousal Recovery after Hypoxic-Ischemic Brain Injury.

Cardiac arrest (CA) remains the leading cause of coma, and early arousal recovery indicators are needed to allocate critical care resources properly. High-frequency oscillations (HFOs) of somatosensory evoked potentials (SSEPs) have been shown to indicate responsive wakefulness days following CA. Nonetheless, their potential in the acute recovery phase, where the injury is reversible, has not been tested. We hypothesize that time-frequency (TF) analysis of HFOs can determine arousal recovery in the acute recovery phase. To test our hypothesis, eleven adult male Wistar rats were subjected to asphyxial CA (five with 3-min mild and six with 7-min moderate to severe CA) and SSEPs were recorded for 60 min post-resuscitation. Arousal level was quantified by the neurological deficit scale (NDS) at 4 h. Our results demonstrated that continuous wavelet transform (CWT) of SSEPs localizes HFOs in the TF domain under baseline conditions. The energy dispersed immediately after injury and gradually recovered. We proposed a novel TF-domain measure of HFO: the total power in the normal time-frequency space (NTFS) of HFO. We found that the NTFS power significantly separated the favorable and unfavorable outcome groups. We conclude that the NTFS power of HFOs provides earlier and objective determination of arousal recovery after CA.

Elucidating the oscillation instability of sessile drops triggered by surface acoustic waves

Surface acoustic waves (SAWs) are used in droplet microfluidics for operations on sessile droplets, such as internal mixing, directional motion, or nebulization. An unsolved problem has been to understand how high frequency (MegaHz) SAWs can counterintuitively cause low frequency (10-500 Hz) oscillations on the drop free surface, which is known to enhance drop mobility. Using simultaneous measurements of the free-surface dynamics and internal acoustic pressure, we find that the drop oscillation instability results from coupling between an intracavity acoustic mode excited by the SAW, and an inertia-capillary surface eigenmode, through amplitude modulation and delayed radiation pressure feedback.

Pulse Intensity Effects of Burst and Tonic Spinal Cord Stimulation on Neural Responses to Brushing in Patients With Neuropathic Pain

ObjectivesTonic spinal cord stimulation (SCS) is accompanied by paresthesia in affected body regions. Comparatively, the absence of paresthesia with burst SCS suggests different involvement of the dorsal column system conveying afferent impulses from low-threshold mechanoreceptors. This study evaluated cortical activation changes during gentle brushing of a pain-free leg during four SCS pulse intensities to assess the effect of intensity on recruitment of dorsal column system fibers during burst and tonic SCS. Materials and MethodsTwenty patients using SCS (11 burst, nine tonic) for neuropathic leg pain participated. Brushing was administered to a pain-free area of the leg during four SCS intensities: therapeutic (100%), medium (66%), low (33%), and no stimulation. Whole-brain electroencephalography was continuously recorded. Changes in spectral power during brushing were evaluated using the event-related desynchronization (ERD) method in theta (4–7 Hz), alpha (8–13 Hz), and beta (16–24 Hz) frequency bands. ResultsBrushing was accompanied by a suppression of cortical oscillations in the range 4–24 Hz. Stronger intensities of burst and tonic SCS led to less suppression of 4–7 Hz and 8–13 Hz bands in parietal electrodes, and in central electrodes in the 16–24 Hz band, with the strongest, statistically significant suppression at medium intensity. Tonic SCS showed a stronger reduction in 4–7 Hz oscillations over right sensorimotor electrodes, and over right frontal and left sensorimotor electrodes in the 8–13 Hz band, compared to burst SCS. ConclusionsResults suggest that burst and tonic SCS are mediated by both different and shared mechanisms. Attenuated brushing-related ERD with tonic SCS suggests a gating of cortical activation by afferent impulses in the dorsal column, whereas burst may engage different pathways. Diminished brushing-related ERD at medium and therapeutic intensities of burst and tonic SCS points towards a nonlinear effect of SCS on somatosensory processing.

Real-World 20-Hz IBR Subsynchronous Oscillations: Signatures and Mechanism Analysis

This article provides analyses of subsynchronous oscillations associated with solar photovoltaics (PVs) observed in real world. All have oscillation frequencies at approximately 20 Hz in voltage and current magnitudes. On the other hand, in phase voltage and current, the oscillations show very different signatures. Through constructing simplified linear models and conducting sensitivity analysis, we are able to associate dynamic phenomena with the underlying mechanism and speculate the possible causes of the 20-Hz oscillations. The analysis results have all been validated via simulation results from electromagnetic transient (EMT) test beds constructed in MATLAB/Simscape.

Modulation of cortical beta oscillations influences motor vigor: A rhythmic TMS-EEG study.

Previous electro- or magnetoencephalography (Electro/Magneto EncephaloGraphic; E/MEG) studies using a correlative approach have shown that β (13-30 Hz) oscillations emerging in the primary motor cortex (M1) are implicated in regulating motor response vigor and associated with an anti-kinetic role, that is, slowness of movement. However, the functional role of M1 β oscillations in regulation of motor responses remains unclear. To address this gap, we combined EEG with rhythmic TMS (rhTMS) delivered to M1 at the β (20 Hz) frequency shortly before subjects performed an isometric ramp-and-hold finger force production task at three force levels. rhTMS is a novel approach that can modulate rhythmic patterns of neural activity. β-rhTMS over M1 induced a modulation of neural oscillations to β frequency in the sensorimotor area and reduced peak force rate during the ramp-up period relative to sham and catch trials. Interestingly, this rhTMS effect occurred only in the large force production condition. To distinguish whether the effects of rhTMS on EEG and behavior stemmed from phase-resetting by each magnetic pulse or neural entrainment by the periodicity of rhTMS, we performed a control experiment using arrhythmic TMS (arTMS). arTMS did not induce changes in EEG oscillations nor peak force rate during the rump-up period. Our results provide novel evidence that β neural oscillations emerging the sensorimotor area influence the regulation of motor response vigor. Furthermore, our findings further demonstrate that rhTMS is a promising tool for tuning neural oscillations to the target frequency.

Corticomotor Plasticity Underlying Priming Effects of Motor Imagery on Force Performance.

The neurophysiological processes underlying the priming effects of motor imagery (MI) on force performance remain poorly understood. Here, we tested whether the priming effects of embedded MI practice involved short-term changes in corticomotor connectivity. In a within-subjects counterbalanced experimental design, participants (n = 20) underwent a series of experimental sessions consisting of successive maximal isometric contractions of elbow flexor muscles. During inter-trial rest periods, we administered MI, action observation (AO), and a control passive recovery condition. We collected electromyograms (EMG) from both agonists and antagonists of the force task, in addition to electroencephalographic (EEG) brain potentials during force trials. Force output was higher during MI compared to AO and control conditions (both p &lt; 0.01), although fatigability was similar across experimental conditions. We also found a weaker relationship between triceps brachii activation and force output during MI and AO compared to the control condition. Imaginary coherence topographies of alpha (8-12 Hz) oscillations revealed increased connectivity between EEG sensors from central scalp regions and EMG signals from agonists during MI, compared to AO and control. Present results suggest that the priming effects of MI on force performance are mediated by a more efficient cortical drive to motor units yielding reduced agonist/antagonist coactivation.

The psychotomimetic ketamine disrupts the transfer of late sensory information in the corticothalamic network.

In prodromal and early schizophrenia, disorders of attention and perception are associated with structural and chemical brain abnormalities and with dysfunctional corticothalamic networks exhibiting disturbed brain rhythms. The underlying mechanisms are elusive. The non-competitive NMDA receptor antagonist ketamine simulates the symptoms of prodromal and early schizophrenia, including disturbances in ongoing and task & sensory-related broadband beta-/gamma-frequency (17-29 Hz/30-80 Hz) oscillations in corticothalamic networks. In normal healthy subjects and rodents, complex integration processes, like sensory perception, induce transient, large-scale synchronised beta/gamma oscillations in a time window of a few hundred ms (200-700 ms) after the presentation of the object of attention (e.g., sensory stimulation). Our goal was to use an electrophysiological multisite network approach to investigate, in lightly anesthetised rats, the effects of a single psychotomimetic dose (2.5mg/kg, subcutaneous) of ketamine on sensory stimulus-induced oscillations. Ketamine transiently increased the power of baseline beta/gamma oscillations and decreased sensory-induced beta/gamma oscillations. In addition, it disrupted information transferability in both the somatosensory thalamus and the related cortex and decreased the sensory-induced thalamocortical connectivity in the broadband gamma range. The present findings support the hypothesis that NMDA receptor antagonism disrupts the transfer of perceptual information in the somatosensory cortico-thalamo-cortical system.

Large-scale biophysically detailed model of somatosensory thalamocortical circuits in NetPyNE.

The primary somatosensory cortex (S1) of mammals is critically important in the perception of touch and related sensorimotor behaviors. In 2015, the Blue Brain Project (BBP) developed a groundbreaking rat S1 microcircuit simulation with over 31,000 neurons with 207 morpho-electrical neuron types, and 37 million synapses, incorporating anatomical and physiological information from a wide range of experimental studies. We have implemented this highly detailed and complex S1 model in NetPyNE, using the data available in the Neocortical Microcircuit Collaboration Portal. NetPyNE provides a Python high-level interface to NEURON and allows defining complicated multiscale models using an intuitive declarative standardized language. It also facilitates running parallel simulations, automates the optimization and exploration of parameters using supercomputers, and provides a wide range of built-in analysis functions. This will make the S1 model more accessible and simpler to scale, modify and extend in order to explore research questions or interconnect to other existing models. Despite some implementation differences, the NetPyNE model preserved the original cell morphologies, electrophysiological responses and spatial distribution for all 207 cell types; and the connectivity properties of all 1941 pathways, including synaptic dynamics and short-term plasticity (STP). The NetPyNE S1 simulations produced reasonable physiological firing rates and activity patterns across all populations. When STP was included, the network generated a 1 Hz oscillation comparable to the original model in vitro-like state. By then reducing the extracellular calcium concentration, the model reproduced the original S1 in vivo-like states with asynchronous activity. These results validate the original study using a new modeling tool. Simulated local field potentials (LFPs) exhibited realistic oscillatory patterns and features, including distance- and frequency-dependent attenuation. The model was extended by adding thalamic circuits, including 6 distinct thalamic populations with intrathalamic, thalamocortical (TC) and corticothalamic connectivity derived from experimental data. The thalamic model reproduced single known cell and circuit-level dynamics, including burst and tonic firing modes and oscillatory patterns, providing a more realistic input to cortex and enabling study of TC interactions. Overall, our work provides a widely accessible, data-driven and biophysically-detailed model of the somatosensory TC circuits that can be employed as a community tool for researchers to study neural dynamics, function and disease.