Unified thalamic model generates multiple distinct oscillations with state-dependent entrainment by stimulation

Visual light flicker stimulation: enhancing alertness in sleep-deprived rats

Introduction

In the evolving field of neurophysiological research, visual light flicker stimulation is recognized as a promising non-invasive intervention for cognitive enhancement, particularly in sleep-deprived conditions.

Methods

This study explored the effects of specific flicker frequencies (40 Hz and 20-30 Hz random flicker) on alertness recovery in sleep-deprived rats. We employed a multidisciplinary approach that included behavioral assessments with the Y-maze, in vivo electrophysiological recordings, and molecular analyses such as c-FOS immunohistochemistry and hormone level measurements.

Results

Both 40 Hz and 20-30 Hz flicker significantly enhanced behavioral performance in the Y-maze test, suggesting an improvement in alertness. Neurophysiological data indicated activation of neural circuits in key brain areas like the thalamus and hippocampus. Additionally, flicker exposure normalized cortisol and serotonin levels, essential for stress response and mood regulation. Notably, increased c-FOS expression in brain regions related to alertness and cognitive functions suggested heightened neural activity.

Discussion

These findings underscore the potential of light flicker stimulation not only to mitigate the effects of sleep deprivation but also to enhance cognitive functions. The results pave the way for future translational research into light-based therapies in human subjects, with possible implications for occupational health and cognitive ergonomics.

Thalamic control of sensory processing and spindles in a biophysical somatosensory thalamoreticular circuit model of wakefulness and sleep

Thalamoreticular circuitry plays a key role in arousal, attention, cognition, and sleep spindles, and is linked to several brain disorders. A detailed computational model of mouse somatosensory thalamus and thalamic reticular nucleus has been developed to capture the properties of over 14,000 neurons connected by 6 million synapses. The model recreates the biological connectivity of these neurons, and simulations of the model reproduce multiple experimental findings in different brain states. The model shows that inhibitory rebound produces frequency-selective enhancement of thalamic responses during wakefulness. We find that thalamic interactions are responsible for the characteristic waxing and waning of spindle oscillations. In addition, we find that changes in thalamic excitability control spindle frequency and their incidence. The model is made openly available to provide a new tool for studying the function and dysfunction of the thalamoreticular circuitry in various brain states.

State-dependent modulation of thalamocortical oscillations by gamma light flicker with different frequencies, intensities, and duty cycles

Rhythmic light flickers have emerged as useful tools to modulate cognition and rescue pathological oscillations related to neurological disorders by entrainment. However, a mechanistic understanding of the entrainment for different brain oscillatory states and light flicker parameters is lacking. To address this issue, we proposed a biophysical neural network model for thalamocortical oscillations (TCOs) and explored the stimulation effects depending on the thalamocortical oscillatory states and stimulation parameters (frequency, intensity, and duty cycle) using the proposed model and electrophysiology experiments. The proposed model generated alpha, beta, and gamma oscillatory states (with main oscillation frequences at 9, 25, and 35 Hz, respectively), which were successfully transmitted from the thalamus to the cortex. By applying light flicker stimulation, we found that the entrainment was state-dependent and it was more prone to induce entrainment if the flicker perturbation frequency was closer to the endogenous oscillatory frequency. In addition, endogenous oscillation would be accelerated, whereas low-frequency oscillatory power would be suppressed by gamma (30–50 Hz) flickers. Notably, the effects of intensity and duty cycle on entrainment were complex; a high intensity of light flicker did not mean high entrainment possibility, and duty cycles below 50% could induce entrainment easier than those above 50%. Further, we observed entrainment discontinuity during gamma flicker stimulations with different frequencies, attributable to the non-linear characteristics of the network oscillations. These results provide support for the experimental design and clinical applications of the modulation of TCOs by gamma (30–50 Hz) light flicker.

Role of the locus coeruleus and basal forebrain in arousal and attention

Experimental evidence has implicated multiple neurotransmitter systems in either the direct or indirect modulation of cortical arousal and attention circuitry. In this review, we selectively focus on three such systems: 1) norepinephrine (NE)-containing neurons of the locus coeruleus (LC), 2) acetylcholine (ACh)-containing neurons of the basal forebrain (BF), and 3) parvalbumin (PV)-containing gamma-aminobutyric acid neurons of the BF. Whereas BF-PV neurons serve as a rapid and transient arousal system, LC-NE and BF-ACh neuromodulation are typically activated on slower but longer-lasting timescales. Recent findings suggest that the BF-PV system serves to rapidly respond to even subtle sensory stimuli with a microarousal. We posit that salient sensory stimuli, such as those that are threatening or predict the need for a response, will quickly activate the BF-PV system and subsequently activate both the BF-ACh and LC-NE systems if the circumstances require longer periods of arousal and vigilance. We suggest that NE and ACh have overlapping psychological functions with the main difference being the precise internal/environmental sensory situations/contexts that recruit each neurotransmitter system - a goal for future research to determine. Implications of dysfunction of each of these three attentional systems for our understanding of neuropsychiatric conditions are considered. Finally, the contemporary availability of research tools to selectively manipulate and measure the activity of these distinctive neuronal populations promises to answer longstanding questions, such as how various arousal systems influence downstream decision-making and motor responding.

State-dependent effects of neural stimulation on brain function and cognition

Invasive and non-invasive brain stimulation methods are widely used in neuroscience to establish causal relationships between distinct brain regions and the sensory, cognitive and motor functions they subserve. When combined with concurrent brain imaging, such stimulation methods can reveal patterns of neuronal activity responsible for regulating simple and complex behaviours at the level of local circuits and across widespread networks. Understanding how fluctuations in physiological states and task demands might influence the effects of brain stimulation on neural activity and behaviour is at the heart of how we use these tools to understand cognition. Here we review the concept of such 'state-dependent' changes in brain activity in response to neural stimulation, and consider examples from research on altered states of consciousness (for example, sleep and anaesthesia) and from task-based manipulations of selective attention and working memory. We relate relevant findings from non-invasive methods used in humans to those obtained from direct electrical and optogenetic stimulation of neuronal ensembles in animal models. Given the widespread use of brain stimulation as a research tool in the laboratory and as a means of augmenting or restoring brain function, consideration of the influence of changing physiological and cognitive states is crucial for increasing the reliability of these interventions.

Multiscale neural modeling of resting-state fMRI reveals executive-limbic malfunction as a core mechanism in major depressive disorder

Major depressive disorder (MDD) represents a grand challenge to human health and society, but the underlying pathophysiological mechanisms remain elusive. Previous neuroimaging studies have suggested that MDD is associated with abnormal interactions and dynamics in two major neural systems including the default mode - salience (DMN-SAL) network and the executive - limbic (EXE-LIM) network, but it is not clear which network plays a central role and which network plays a subordinate role in MDD pathophysiology. To address this question, we refined a newly developed Multiscale Neural Model Inversion (MNMI) framework and applied it to test whether MDD is more affected by impaired circuit interactions in the DMN-SAL network or the EXE-LIM network. The model estimates the directed connection strengths between different neural populations both within and between brain regions based on resting-state fMRI data collected from normal healthy subjects and patients with MDD. Results show that MDD is primarily characterized by abnormal circuit interactions in the EXE-LIM network rather than the DMN-SAL network. Specifically, we observe reduced frontoparietal effective connectivity that potentially contributes to hypoactivity in the dorsolateral prefrontal cortex (dlPFC), and decreased intrinsic inhibition combined with increased excitation from the superior parietal cortex (SPC) that potentially lead to amygdala hyperactivity, together resulting in activation imbalance in the PFC-amygdala circuit that pervades in MDD. Moreover, the model reveals reduced PFC-to-hippocampus excitation but decreased SPC-to-thalamus inhibition in MDD population that potentially lead to hypoactivity in the hippocampus and hyperactivity in the thalamus, consistent with previous experimental data. Overall, our findings provide strong support for the long-standing limbic-cortical dysregulation model in major depression but also offer novel insights into the multiscale pathophysiology of this debilitating disease.

Impact of GABAA and GABAB Inhibition on Cortical Dynamics and Perturbational Complexity during Synchronous and Desynchronized States

Quantitative estimations of spatiotemporal complexity of cortical activity patterns are used in the clinic as a measure of consciousness levels, but the cortical mechanisms involved are not fully understood. We used a version of the perturbational complexity index (PCI) adapted to multisite recordings from the ferret (either sex) cerebral cortex in vitro (sPCI) to investigate the role of GABAergic inhibition in cortical complexity. We studied two dynamical states: slow-wave activity (synchronous state) and desynchronized activity, that express low and high causal complexity respectively. Progressive blockade of GABAergic inhibition during both regimes revealed its impact on the emergent cortical activity and on sPCI. Gradual GABAA receptor blockade resulted in higher synchronization, being able to drive the network from a desynchronized to a synchronous state, with a progressive decrease of complexity (sPCI). Blocking GABAB receptors also resulted in a reduced sPCI, in particular when in a synchronous, slow wave state. Our findings demonstrate that physiological levels of inhibition contribute to the generation of dynamical richness and spatiotemporal complexity. However, if inhibition is diminished or enhanced, cortical complexity decreases. Using a computational model, we explored a larger parameter space in this relationship and demonstrate a link between excitatory/inhibitory balance and the complexity expressed by the cortical network.SIGNIFICANCE STATEMENT The spatiotemporal complexity of the activity expressed by the cerebral cortex is a highly revealing feature of the underlying network's state. Complexity varies with physiological brain states: it is higher during awake than during sleep states. But it also informs about pathologic states: in disorders of consciousness, complexity is lower in an unresponsive wakefulness syndrome than in a minimally conscious state. What are the network parameters that modulate complexity? Here we investigate how inhibition, mediated by either GABAA or GABAA receptors, influences cortical complexity. And we do this departing from two extreme functional states: a highly synchronous, slow-wave state, and a desynchronized one that mimics wakefulness. We find that there is an optimal level of inhibition in which complexity is highest.

Thalamocortical Spectral Transmission Relies on Balanced Input Strengths

The thalamus is a key element of sensory transmission in the brain, as it gates and selects sensory streams through a modulation of its internal activity. A preponderant role in these functions is played by its internal activity in the alpha range ([8–14] Hz), but the mechanism underlying this process is not completely understood. In particular, how do thalamocortical connections convey stimulus driven information selectively over the back-ground of thalamic internally generated activity? Here we investigate this issue with a spiking network model of feedforward connectivity between thalamus and primary sensory cortex reproducing the local field potential of both areas. We found that in a feedforward network, thalamic oscillations in the alpha range do not entrain cortical activity for two reasons: (i) alpha range oscillations are weaker in neurons projecting to the cortex, (ii) the gamma resonance dynamics of cortical networks hampers oscillations over the 10–20 Hz range thus weakening alpha range oscillations. This latter mechanism depends on the balance of the strength of thalamocortical connections toward excitatory and inhibitory neurons in the cortex. Our results highlight the relevance of corticothalamic feedback to sustain alpha range oscillations and pave the way toward an integrated understanding of the sensory streams traveling between the periphery and the cortex.

Transcranial alternating current stimulation entrains alpha oscillations by preferential phase synchronization of fast-spiking cortical neurons to stimulation waveform

Computational modeling and human studies suggest that transcranial alternating current stimulation (tACS) modulates alpha oscillations by entrainment. Yet, a direct examination of how tACS interacts with neuronal spiking activity that gives rise to the alpha oscillation in the thalamo-cortical system has been lacking. Here, we demonstrate how tACS entrains endogenous alpha oscillations in head-fixed awake ferrets. We first show that endogenous alpha oscillations in the posterior parietal cortex drive the primary visual cortex and the higher-order visual thalamus. Spike-field coherence is largest for the alpha frequency band, and presumed fast-spiking inhibitory interneurons exhibit strongest coupling to this oscillation. We then apply alpha-tACS that results in a field strength comparable to what is commonly used in humans (<0.5 mV/mm). Both in these ferret experiments and in a computational model of the thalamo-cortical system, tACS entrains alpha oscillations by following the theoretically predicted Arnold tongue. Intriguingly, the fast-spiking inhibitory interneurons exhibit a stronger entrainment response to tACS in both the ferret experiments and the computational model, likely due to their stronger endogenous coupling to the alpha oscillation. Our findings demonstrate the in vivo mechanism of action for the modulation of the alpha oscillation by tACS.

Evidence of state-dependence in the effectiveness of responsive neurostimulation for seizure modulation

Background:

An implanted device for brain-responsive neurostimulation (RNS^® System) is approved as an effective treatment to reduce seizures in adults with medically-refractory focal epilepsy. Clinical trials of the RNS System demonstrate population-level reduction in average seizure frequency, but therapeutic response is highly variable.

Hypothesis:

Recent evidence links seizures to cyclical fluctuations in underlying risk. We tested the hypothesis that effectiveness of responsive neurostimulation varies based on current state within cyclical risk fluctuations.

Methods:

We analyzed retrospective data from 25 adults with medically-refractory focal epilepsy implanted with the RNS System. Chronic electrocorticography was used to record electrographic seizures, and hidden Markov models decoded seizures into fluctuations in underlying risk. State-dependent associations of RNS System stimulation parameters with changes in risk were estimated.

Results:

Higher charge density was associated with improved outcomes, both for remaining in a low seizure risk state and for transitioning from a high to a low seizure risk state. The effect of stimulation frequency depended on initial seizure risk state: when starting in a low risk state, higher stimulation frequencies were associated with remaining in a low risk state, but when starting in a high risk state, lower stimulation frequencies were associated with transition to a low risk state. Findings were consistent across bipolar and monopolar stimulation configurations.

Conclusion:

The impact of RNS on seizure frequency exhibits state-dependence, such that stimulation parameters which are effective in one seizure risk state may not be effective in another. These findings represent conceptual advances in understanding the therapeutic mechanism of RNS, and directly inform current practices of RNS tuning and the development of next-generation neurostimulation systems.

Revealing the Computational Meaning of Neocortical Interarea Signals

To understand the function of the neocortex, which is a hierarchical distributed network, it is useful giving meaning to the signals transmitted between these areas from the computational viewpoint. The overall anatomical structure or organs related to this network, including the neocortex, thalamus, and basal ganglia, has been roughly revealed, and much physiological knowledge, though often fragmentary, is being accumulated. The computational theories involving the neocortex have also been developed considerably. By introducing the assumption “The signals transmitted by interarea axonal projections of pyramidal cells in the neocortex carry different meanings for each cell type, common to all areas,” derived from its nature as a distributed network in the neocortex, allows us to specify the computational meanings of interarea signals. In this paper, first, the types of signals exchanged between neocortical areas are investigated, taking into account biological constraints, and employing theories such as predictive coding, reinforcement learning, representation emulation theory, and BDI logic as theoretical starting points, two types of feedforward signals (observation and deviation) and three types of feedback signals (prediction, plan, and intention) are identified. Next, based on the anatomical knowledge of the neocortex and thalamus, the pathways connecting the areas are organized and summarized as three corticocortical pathways and two thalamocortical pathways. Using this summation as preparation, this paper proposes a hypothesis that gives meaning to each type of signals transmitted in the different pathways in the neocortex, from the viewpoint of their functions. This hypothesis reckons that the feedforward corticocortical pathway transmits observation signals, the feedback corticocortical pathway transmits prediction signals, and the corticothalamic pathway mediated by core relay cells transmits deviation signals. The thalamocortical pathway, which is mediated by matrix relay cells, would be responsible for transmitting the signals that activate a part of prediction signals as intentions, due to the reason that the nature of the other available feedback pathways are not sufficient for conveying plans and intentions as signals. The corticocortical pathway, which is projected from various IT cells to the first layer, would be responsible for transmitting signals that activate a part of prediction signals as plans.

Relations between large-scale brain connectivity and effects of regional stimulation depend on collective dynamical state

At the macroscale, the brain operates as a network of interconnected neuronal populations, which display coordinated rhythmic dynamics that support interareal communication. Understanding how stimulation of different brain areas impacts such activity is important for gaining basic insights into brain function and for further developing therapeutic neurmodulation. However, the complexity of brain structure and dynamics hinders predictions regarding the downstream effects of focal stimulation. More specifically, little is known about how the collective oscillatory regime of brain network activity—in concert with network structure—affects the outcomes of perturbations. Here, we combine human connectome data and biophysical modeling to begin filling these gaps. By tuning parameters that control collective system dynamics, we identify distinct states of simulated brain activity and investigate how the distributed effects of stimulation manifest at different dynamical working points. When baseline oscillations are weak, the stimulated area exhibits enhanced power and frequency, and due to network interactions, activity in this excited frequency band propagates to nearby regions. Notably, beyond these linear effects, we further find that focal stimulation causes more distributed modifications to interareal coherence in a band containing regions’ baseline oscillation frequencies. Importantly, depending on the dynamical state of the system, these broadband effects can be better predicted by functional rather than structural connectivity, emphasizing a complex interplay between anatomical organization, dynamics, and response to perturbation. In contrast, when the network operates in a regime of strong regional oscillations, stimulation causes only slight shifts in power and frequency, and structural connectivity becomes most predictive of stimulation-induced changes in network activity patterns. In sum, this work builds upon and extends previous computational studies investigating the impacts of stimulation, and underscores the fact that both the stimulation site, and, crucially, the regime of brain network dynamics, can influence the network-wide responses to local perturbations.

Stimulation can be used to alter brain activity and is a therapeutic option for certain neurological conditions. However, predicting the distributed effects of local perturbations is difficult. Previous studies show that responses to stimulation depend on anatomical (or structural) coupling. In addition to structure, here we consider how stimulation effects also depend on the brain’s collective dynamical (or functional) state, arising from the coordination of rhythmic activity across large-scale networks. In a whole-brain computational model, we show that global responses to regional stimulation can indeed be contingent upon and differ across various dynamical working points. Notably, depending on the network’s oscillatory regime, stimulation can accelerate the activity of the stimulated site, and lead to widespread effects at both the new, excited frequency, as well as in a much broader frequency range including areas’ baseline frequencies. While structural connectivity is a good predictor of “excited band” changes, in some states “baseline band” effects can be better predicted by functional connectivity, which depends upon the system’s oscillatory regime. By integrating and extending past efforts, our results thus indicate that dynamical—in additional to structural—brain organization plays a role in governing how focal stimulation modulates interactions between distributed network elements.

An Attempt at a Unified Theory of the Neocortical Microcircuit in Sensory Cortex

The neocortex performs a wide range of functions, including working memory, sensory perception, and motor planning. Despite this diversity in function, evidence suggests that the neocortex is made up of repeating subunits ("macrocolumns"), each of which is largely identical in circuitry. As such, the specific computations performed by these macrocolumns are of great interest to neuroscientists and AI researchers. Leading theories of this microcircuit include models of predictive coding, hierarchical temporal memory (HTM), and Adaptive Resonance Theory (ART). However, these models have not yet explained: (1) how microcircuits learn sequences input with delay (i.e., working memory); (2) how networks of columns coordinate processing on precise timescales; or (3) how top-down attention modulates sensory processing. I provide a theory of the neocortical microcircuit that extends prior models in all three ways. Additionally, this theory provides a novel working memory circuit that extends prior models to support simultaneous multi-item storage without disrupting ongoing sensory processing. I then use this theory to explain the functional origin of a diverse set of experimental findings, such as cortical oscillations.

Converging intracortical signatures of two separated processing timescales in human early auditory cortex

Neural oscillations in auditory cortex are argued to support parsing and representing speech constituents at their corresponding temporal scales. Yet, how incoming sensory information interacts with ongoing spontaneous brain activity, what features of the neuronal microcircuitry underlie spontaneous and stimulus-evoked spectral fingerprints, and what these fingerprints entail for stimulus encoding, remain largely open questions. We used a combination of human invasive electrophysiology, computational modeling and decoding techniques to assess the information encoding properties of brain activity and to relate them to a plausible underlying neuronal microarchitecture. We analyzed intracortical auditory EEG activity from 10 patients while they were listening to short sentences. Pre-stimulus neural activity in early auditory cortical regions often exhibited power spectra with a shoulder in the delta range and a small bump in the beta range. Speech decreased power in the beta range, and increased power in the delta-theta and gamma ranges. Using multivariate machine learning techniques, we assessed the spectral profile of information content for two aspects of speech processing: detection and discrimination. We obtained better phase than power information decoding, and a bimodal spectral profile of information content with better decoding at low (delta-theta) and high (gamma) frequencies than at intermediate (beta) frequencies. These experimental data were reproduced by a simple rate model made of two subnetworks with different timescales, each composed of coupled excitatory and inhibitory units, and connected via a negative feedback loop. Modeling and experimental results were similar in terms of pre-stimulus spectral profile (except for the iEEG beta bump), spectral modulations with speech, and spectral profile of information content. Altogether, we provide converging evidence from both univariate spectral analysis and decoding approaches for a dual timescale processing infrastructure in human auditory cortex, and show that it is consistent with the dynamics of a simple rate model.

Effects of Cholinergic Neuromodulation on Thalamocortical Rhythms During NREM Sleep: A Model Study

It has been suggested that cholinergic neurons shape the oscillatory activity of the thalamocortical (TC) network in behavioral and electrophysiological experiments. However, theoretical modeling demonstrating how cholinergic neuromodulation of thalamocortical rhythms during non-rapid eye movement (NREM) sleep might occur has been lacking. In this paper, we first develop a novel computational model (TC-ACH) by incorporating a cholinergic neuron population (CH) into the classical thalamo-cortical circuitry, where connections between populations are modeled in accordance with existing knowledge. The neurotransmitter acetylcholine (ACH) released by neurons in CH, which is able to change the discharge activity of thalamocortical neurons, is the primary focus of our work. Simulation results with our TC-ACH model reveal that the cholinergic projection activity is a key factor in modulating oscillation patterns in three ways: (1) transitions between different patterns of thalamocortical oscillations are dramatically modulated through diverse projection pathways; (2) the model expresses a stable spindle oscillation state with certain parameter settings for the cholinergic projection from CH to thalamus, and more spindles appear when the strength of cholinergic input from CH to thalamocortical neurons increases; (3) the duration of oscillation patterns during NREM sleep including K-complexes, spindles, and slow oscillations is longer when cholinergic input from CH to thalamocortical neurons becomes stronger. Our modeling results provide insights into the mechanisms by which the sleep state is controlled, and provide a theoretical basis for future experimental and clinical studies.

Primate somatosensory cortical neurons are entrained to both spontaneous and peripherally evoked spindle oscillations

Recurrent thalamocortical circuits produce a number of rhythms critical to brain function. In slow-wave sleep, spindles (7-16 Hz) are a prominent spontaneous oscillation generated by thalamic circuits and triggered by cortical slow waves. In wakefulness and under anesthesia, brief peripheral sensory stimuli can evoke 10-Hz reverberations due potentially to similar thalamic mechanisms. Functionally, sleep spindles and peripherally evoked spindles may play a role in memory consolidation and perception, respectively. Yet, rarely have the circuits involved in these two rhythms been compared in the same animals and never in primates. Here, we investigated the entrainment of primary somatosensory cortex (S1) neurons to both rhythms in ketamine-sedated macaques. First, we compared spontaneous spindles in sedation and natural sleep to validate the model. Then, we quantified entrainment with spike-field coherence and phase-locking statistics. We found that S1 neurons entrained to spontaneous sleep spindles were also entrained to the evoked spindles, although entrainment strength and phase systematically differed. Our results indicate that the spindle oscillations triggered by top-down spontaneous cortical activity and bottom-up peripheral input share a common cortical substrate.NEW & NOTEWORTHY Brief sensory stimuli evoke 10-Hz oscillations in thalamocortical neuronal activity and in perceptual thresholds. The mechanisms underlying this evoked rhythm are not well understood but are thought to be similar to those generating sleep spindles. We directly compared the entrainment of cortical neurons to both spontaneous spindles and peripherally evoked oscillations in sedated monkeys. We found that the entrainment strengths to each rhythm were positively correlated, although with differing entrainment phases, implying involvement of similar networks.

Identification of Abnormal Circuit Dynamics in Major Depressive Disorder via Multiscale Neural Modeling of Resting-State fMRI

Resting-state functional magnetic resonance imaging (rs-fMRI) studies have focused primarily on characterizing functional or effective connectivity of discrete brain regions. A major drawback of this approach is that it does not provide a mechanistic understanding of brain cognitive function or dysfunction at cellular and circuit levels. To overcome this limitation, we combined the methods of computational neuroscience with traditional macroscale connectomic analysis and developed a Multiscale Neural Model Inversion (MNMI) framework that links microscale circuit interaction with macroscale network dynamics and estimates both local coupling and inter-regional connections via stochastic optimization based on blood oxygen-level dependent (BOLD) rs-fMRI. We applied this method to the rs-fMRI data of 66 normal healthy subjects and 66 individuals with major depressive disorder (MDD) to identify potential biomarkers at both local circuit and global network level. Results suggest that the recurrent excitation and inhibition within the dorsal lateral prefrontal cortex (dlPFC) might be disrupted in MDD, consistent with the commonly accepted hypothetical model of MDD. In addition, recurrent excitation in the thalamus was found to be abnormally elevated, which may be responsible to abnormal thalamocortical oscillations often observed in MDD. Overall, our modeling approach holds the promise to overcome the limitation of traditional large-scale connectome modeling by providing hidden mechanistic insights into neuroanatomy, circuit dynamics and pathophysiology.

Spectral organization of focal seizures within the thalamotemporal network

Objective

To investigate dynamic changes in neural activity between the anterior nucleus of the thalamus (ANT) and the seizure onset zone (SOZ) in patients with drug‐resistant temporal lobe epilepsy (TLE) based on anatomic location, seizure subtype, and state of vigilance (SOV).

Methods

Eleven patients undergoing stereoelectroencephalography for seizure localization were recruited prospectively for local field potential (LFP) recording directly from the ANT. The SOZ was identified using line length and epileptogenicity index. Changes in power spectral density (PSD) were compared between the two anatomic sites as seizures (N = 53) transitioned from interictal baseline to the posttermination stage.

Results

At baseline, the thalamic LFPs were significantly lower and distinct from the SOZ with the presence of higher power in the fast ripple band (P < 0.001). Temporal changes in ictal power of neural activity within ANT mimic those of the SOZ, are increased significantly at seizure onset (P < 0.05), and are distinct for seizures that impaired awareness or that secondarily generalized (P < 0.05). The onset of seizure was preceded by a decrease in the mean power spectral density (PSD) in ANT and SOZ (P < 0.05). Neural activity correlated with different states of vigilance at seizure onset within the ANT but not in the SOZ (P = 0.005).

Interpretation

The ANT can be recruited at the onset of mesial temporal lobe seizures, and the recruitment pattern differs with seizure subtypes. Furthermore, changes in neural dynamics precede seizure onset and are widespread to involve temporo‐thalamic regions, thereby providing an opportunity to intervene early with closed‐loop DBS.

Low-frequency direct cortical stimulation of left superior frontal gyrus enhances working memory performance

The neural substrates of working memory are spread across prefrontal, parietal and cingulate cortices and are thought to be coordinated through low frequency cortical oscillations in the theta (3-8 Hz) and alpha (8-12 Hz) frequency bands. While the functional role of many subregions have been elucidated using neuroimaging studies, the role of superior frontal gyrus (SFG) is not yet clear. Here, we combined electrocorticography and direct cortical stimulation in three patients implanted with subdural electrodes to assess if superior frontal gyrus is indeed involved in working memory. We found left SFG exhibited task-related modulation of oscillations in the theta and alpha frequency bands specifically during the encoding epoch. Stimulation at the frequency matched to the endogenous oscillations resulted in reduced reaction times in all three participants. Our results provide evidence for SFG playing a functional role in working memory and suggest that SFG may coordinate working memory through low-frequency oscillations thus bolstering the feasibility of using intracranial electric stimulation for restoring cognitive function.

Neuromodulation-dependent effect of gated high-frequency, LFMS-like electric field stimulation in mouse cortical slices

Low-field magnetic stimulation (LFMS) is a gated high-frequency non-invasive brain stimulation method (500 Hz gated at 2 Hz) with a proposed antidepressant effect. However, it has remained unknown how such stimulation paradigms modulate neuronal network activity and how the induced changes depend on network state. Here we examined the immediate and outlasting effects of the gated high-frequency electric field associated with LFMS on the cortical activity as a function of neuromodulatory tone that defines network state. We used a sham-controlled study design to investigate effects of stimulation (20 min of 0.5 s trains of 500 Hz charge-balanced pulse stimulation patterned at 0.5 Hz) on neural activity in mouse medial prefrontal cortex in vitro. Bath application of cholinergic and noradrenergic agents enabled us to examine the stimulation effects as a function of neuromodulatory tone. The stimulation attenuated the increase in firing rate of layer V cortical neurons during the post-stimulation period in the presence of cholinergic activation. The same stimulation had no significant immediate or outlasting effect in the absence of exogenous neuromodulators or in the presence of noradrenergic activation. These results provide electrophysiological insights into the neuromodulatory-dependent effects of gated high-frequency stimulation. More broadly, our results are the first to provide a mechanistic demonstration of how behavioral states and arousal levels may modify the effects of non-invasive brain stimulation.

Rhythmic modulation of thalamic oscillations depends on intrinsic cellular dynamics

OBJECTIVE

Rhythmic brain stimulation has emerged as a powerful tool to modulate cognition and to target pathological oscillations related to neurological and psychiatric disorders. However, we lack a systematic understanding of how periodic stimulation interacts with endogenous neural activity as a function of the brain state and target.

APPROACH

To address this critical issue, we applied periodic stimulation to a unified biophysical thalamic network model that generates multiple distinct oscillations, and examined thoroughly the impact of rhythmic stimulation on different oscillatory states.

MAIN RESULTS

We found that rhythmic perturbation induces four basic response mechanisms: entrainment, acceleration, resonance and suppression. Importantly, the appearance and expression of these mechanisms depend highly on the intrinsic cellular dynamics in each state. Specifically, the low-threshold bursting of thalamocortical cells (TCs) in delta (δ) oscillation renders the network relatively insensitive to entrainment; the high-threshold bursting of TCs in alpha (α) oscillation leads to widespread oscillation suppression while the tonic spiking of TC cells in gamma (γ) oscillation results in prominent entrainment and resonance. In addition, we observed entrainment discontinuity during α oscillation that is mediated by firing pattern switching of high-threshold bursting TC cells. Furthermore, we demonstrate that direct excitatory stimulation of the lateral geniculate nucleus (LGN) entrains thalamic oscillations via an asymmetric Arnold tongue that favors higher frequency entrainment and resonance, while stimulation of the inhibitory circuit, the reticular nucleus, induces much weaker and more symmetric entrainment and resonance. These results support the notion that rhythmic stimulation engages brain oscillations in a state- and target-dependent manner.

SIGNIFICANCE

Overall, our study provides, for the first time, insights into how the biophysics of thalamic oscillations guide the emergence of complex, state-dependent mechanisms of target engagement, which can be leveraged for the future rational design of novel therapeutic stimulation modalities.

Spiking patterns and synchronization of thalamic neurons along the sleep-wake cycle

Spiking patterns and synchronization dynamics of thalamic neurons along the sleep-wake cycle are studied in a minimal model of four coupled conductance-based neurons. The model simulates two thalamic neurons coupled via a gap junction and driven by a synaptic input from a two-neuron model of sleep regulation by the hypothalamus. In accord with experimental data, the model shows that during sleep, when hypothalamic wake-active neurons are silent, the thalamic neurons discharge bursts of spikes. During wake, the excitatory synaptic input from the hypothalamus drives the coupled thalamic neurons to a state of tonic firing (single spikes). In the deterministic case, the thalamic neurons synchronize in-phase in the bursting regime but demonstrate multi-stability of out-of-phase, in-phase, and asynchronous states in the tonic firing. However, along the sleep-wake cycle, once the neurons synchronize in-phase during sleep (bursting), they stay synchronized in wake (tonic firing). It is thus found that noise is needed to reproduce the experimentally observed transitions between synchronized bursting during sleep and asynchronous tonic firing during wake. Overall, synchronization of bursting is found to be more robust to noise than synchronization of tonic firing, where a small disturbance is sufficient to desynchronize the thalamic neurons. The model predicts that the transitions between sleep and wake happen via chaos because a single thalamic neuron exhibits chaos between regular bursting and tonic activity. The results of this study suggest that the sleep- and wake-related dynamics in the thalamus may be generated at a level of gap junction-coupled clusters of thalamic neurons driven from the hypothalamus which would then propagate throughout the thalamus and cortex via axonal long-range connections.

Stochastic resonance mediates the state-dependent effect of periodic stimulation on cortical alpha oscillations

Brain stimulation can be used to engage and modulate rhythmic activity in brain networks. However, the outcomes of brain stimulation are shaped by behavioral states and endogenous fluctuations in brain activity. To better understand how this intrinsic oscillatory activity controls the susceptibility of the brain to stimulation, we analyzed a computational model of the thalamo-cortical system in two distinct states (rest and task-engaged) to identify the mechanisms by which endogenous alpha oscillations (8Hz–12Hz) are modulated by periodic stimulation. Our analysis shows that the different responses to stimulation observed experimentally in these brain states can be explained by a passage through a bifurcation combined with stochastic resonance — a mechanism by which irregular fluctuations amplify the response of a nonlinear system to weak periodic signals. Indeed, our findings suggest that modulation of brain oscillations is best achieved in states of low endogenous rhythmic activity, and that irregular state-dependent fluctuations in thalamic inputs shape the susceptibility of cortical population to periodic stimulation.