Interactions between membrane conductances underlying thalamocortical slow-wave oscillations

Mechanism of an Intrinsic Oscillation in Rat Geniculate Interneurons

Depolarizing current injections produced a rhythmic bursting of action potentials - a bursting oscillation - in a set of local interneurons in the lateral geniculate nucleus (LGN) of rats. The current dynamics underlying this firing pattern have not been determined, though this cell type constitutes an important cellular component of thalamocortical circuitry, and contributes to both pathologic and non-pathologic brain states. We thus investigated the source of the bursting oscillation using pharmacological manipulations in LGN slices in vitro and in silico. 1. Selective blockade of calcium channel subtypes revealed that high-threshold calcium currentsI L andI P contributed strongly to the oscillation. 2. Increased extracellular K+ concentration (decreased K+currents) eliminated the oscillation. 3. Selective blockade of K+ channel subtypes demonstrated that the calcium-sensitive potassium current (I A H P ) was of primary importance. A morphologically simplified, multicompartment model of the thalamic interneuron characterized the oscillation as follows: 1. The low-threshold calcium currentI T provided the strong initial burst characteristic of the oscillation. 2. Alternating fluxes through high-threshold calcium channels andI A H P then provided the continuing oscillation's burst and interburst periods respectively. This interplay betweenI L andI A H P contrasts with the current dynamics underlying oscillations in thalamocortical and reticularis neurons, which primarily involveI T andI H , orI T andI A H P respectively. These findings thus point to a novel electrophysiological mechanism for generating intrinsic oscillations in a major thalamic cell type. Because local interneurons can sculpt the behavior of thalamocortical circuits, these results suggest new targets for the manipulation of ascending thalamocortical network activity.

Two contrasting mediodorsal thalamic circuits target the mouse medial prefrontal cortex

At the heart of the prefrontal network is the mediodorsal (MD) thalamus. Despite the importance of MD in a broad range of behaviors and neuropsychiatric disorders, little is known about the physiology of neurons in MD. We injected the retrograde tracer cholera toxin subunit B (CTB) into the medial prefrontal cortex (mPFC) of adult wild-type mice. We prepared acute brain slices and used current clamp electrophysiology to measure and compare the intrinsic properties of the neurons in MD that project to mPFC (MD→mPFC neurons). We show that MD→mPFC neurons are located predominantly in the medial (MD-M) and lateral (MD-L) subnuclei of MD. MD-L→mPFC neurons had shorter membrane time constants and lower membrane resistance than MD-M→mPFC neurons. Relatively increased hyperpolarization-activated cyclic nucleotide-gated (HCN) channel activity in MD-L neurons accounted for the difference in membrane resistance. MD-L neurons had a higher rheobase that resulted in less readily generated action potentials compared with MD-M→mPFC neurons. In both cell types, HCN channels supported generation of burst spiking. Increased HCN channel activity in MD-L neurons results in larger after-hyperpolarization potentials compared with MD-M neurons. These data demonstrate that the two populations of MD→mPFC neurons have divergent physiologies and support a differential role in thalamocortical information processing and potentially behavior.NEW & NOTEWORTHY To realize the potential of circuit-based therapies for psychiatric disorders that localize to the prefrontal network, we need to understand the properties of the populations of neurons that make up this network. The mediodorsal (MD) thalamus has garnered attention for its roles in executive functioning and social/emotional behaviors mediated, at least in part, by its projections to the medial prefrontal cortex (mPFC). Here, we identify and compare the physiology of the projection neurons in the two MD subnuclei that provide ascending inputs to mPFC in mice. Differences in intrinsic excitability between the two populations of neurons suggest that neuromodulation strategies targeting the prefrontal thalamocortical network will have differential effects on these two streams of thalamic input to mPFC.

The role of the nucleus basalis of Meynert in neuromodulation therapy: a systematic review from the perspective of neural network oscillations

As a crucial component of the cerebral cholinergic system and the Papez circuit in the basal forebrain, dysfunction of the nucleus basalis of Meynert (NBM) is associated with various neurodegenerative disorders. However, no drugs, including existing cholinesterase inhibitors, have been shown to reverse this dysfunction. Due to advancements in neuromodulation technology, researchers are exploring the use of deep brain stimulation (DBS) therapy targeting the NBM (NBM-DBS) to treat mental and neurological disorders as well as the related mechanisms. Herein, we provided an update on the research progress on cognition-related neural network oscillations and complex anatomical and projective relationships between the NBM and other cognitive structures and circuits. Furthermore, we reviewed previous animal studies of NBM lesions, NBM-DBS models, and clinical case studies to summarize the important functions of the NBM in neuromodulation. In addition to elucidating the mechanism of the NBM neural network, future research should focus on to other types of neurons in the NBM, despite the fact that cholinergic neurons are still the key target for cell type-specific activation by DBS.

Canonical circuit computations for computer vision

Advanced computer vision mechanisms have been inspired by neuroscientific findings. However, with the focus on improving benchmark achievements, technical solutions have been shaped by application and engineering constraints. This includes the training of neural networks which led to the development of feature detectors optimally suited to the application domain. However, the limitations of such approaches motivate the need to identify computational principles, or motifs, in biological vision that can enable further foundational advances in machine vision. We propose to utilize structural and functional principles of neural systems that have been largely overlooked. They potentially provide new inspirations for computer vision mechanisms and models. Recurrent feedforward, lateral, and feedback interactions characterize general principles underlying processing in mammals. We derive a formal specification of core computational motifs that utilize these principles. These are combined to define model mechanisms for visual shape and motion processing. We demonstrate how such a framework can be adopted to run on neuromorphic brain-inspired hardware platforms and can be extended to automatically adapt to environment statistics. We argue that the identified principles and their formalization inspires sophisticated computational mechanisms with improved explanatory scope. These and other elaborated, biologically inspired models can be employed to design computer vision solutions for different tasks and they can be used to advance neural network architectures of learning.

High gamma coherence between task-responsive sensory-motor cortical regions in a motor reaction-time task

Electrical activity at high gamma frequencies (70-170 Hz) is thought to reflect the activity of small cortical ensembles. For example, high gamma activity (often quantified by spectral power) can increase in sensory-motor cortex in response to sensory stimuli or movement. On the other hand, synchrony of neural activity between cortical areas (often quantified by coherence) has been hypothesized as an important mechanism for inter-areal communication, thereby serving functional roles in cognition and behavior. Currently, high gamma activity has primarily been studied as a local amplitude phenomenon. We investigated the synchronization of high gamma activity within sensory-motor cortex and the extent to which underlying high gamma activity can explain coherence during motor tasks. We characterized high gamma coherence in sensory-motor networks and the relationship between coherence and power by analyzing electrocorticography (ECoG) data from human subjects as they performed a motor response to sensory cues. We found greatly increased high gamma coherence during the motor response compared with the sensory cue. High gamma power poorly predicted high gamma coherence, but the two shared a similar time course. However, high gamma coherence persisted longer than high gamma power. The results of this study suggest that high gamma coherence is a physiologically distinct phenomenon during a sensory-motor task, the emergence of which may require active task participation.NEW & NOTEWORTHY Motor action after auditory stimulus elicits high gamma responses in sensory-motor and auditory cortex, respectively. We show that high gamma coherence reliably and greatly increased during motor response, but not after auditory stimulus. Underlying high gamma power could not explain high gamma coherence. Our results indicate that high gamma coherence is a physiologically distinct sensory-motor phenomenon that may serve as an indicator of increased synaptic communication on short timescales (∼1 s).

Rapid thalamocortical network switching mediated by cortical synchronization underlies propofol-induced EEG signatures: a biophysical model

Propofol-mediated unconsciousness elicits strong alpha/low-beta and slow oscillations in the electroencephalogram (EEG) of patients. As anesthetic dose increases, the EEG signal changes in ways that give clues to the level of unconsciousness; the network mechanisms of these changes are only partially understood. Here, we construct a biophysical thalamocortical network involving brain stem influences that reproduces transitions in dynamics seen in the EEG involving the evolution of the power and frequency of alpha/low-beta and slow rhythm, as well as their interactions. Our model suggests that propofol engages thalamic spindle and cortical sleep mechanisms to elicit persistent alpha/low-beta and slow rhythms, respectively. The thalamocortical network fluctuates between two mutually exclusive states on the timescale of seconds. One state is characterized by continuous alpha/low-beta-frequency spiking in thalamus (C-state), whereas in the other, thalamic alpha spiking is interrupted by periods of co-occurring thalamic and cortical silence (I-state). In the I-state, alpha colocalizes to the peak of the slow oscillation; in the C-state, there is a variable relationship between an alpha/beta rhythm and the slow oscillation. The C-state predominates near loss of consciousness; with increasing dose, the proportion of time spent in the I-state increases, recapitulating EEG phenomenology. Cortical synchrony drives the switch to the I-state by changing the nature of the thalamocortical feedback. Brain stem influence on the strength of thalamocortical feedback mediates the amount of cortical synchrony. Our model implicates loss of low-beta, cortical synchrony, and coordinated thalamocortical silent periods as contributing to the unconscious state.NEW & NOTEWORTHY GABAergic anesthetics induce alpha/low-beta and slow oscillations in the EEG, which interact in dose-dependent ways. We constructed a thalamocortical model to investigate how these interdependent oscillations change with propofol dose. We find two dynamic states of thalamocortical coordination, which change on the timescale of seconds and dose-dependently mirror known changes in EEG. Thalamocortical feedback determines the oscillatory coupling and power seen in each state, and this is primarily driven by cortical synchrony and brain stem neuromodulation.

Responses in fast-spiking interneuron firing rates to parameter variations associated with degradation of perineuronal nets

The perineuronal nets (PNNs) are sugar coated protein structures that encapsulate certain neurons in the brain, such as parvalbumin positive (PV) inhibitory neurons. As PNNs are theorized to act as a barrier to ion transport, they may effectively increase the membrane charge-separation distance, thereby affecting the membrane capacitance. Tewari et al. (2018) found that degradation of PNNs induced a 25%-50% increase in membrane capacitance [Formula: see text] and a reduction in the firing rates of PV-cells. In the current work, we explore how changes in [Formula: see text] affects the firing rate in a selection of computational neuron models, ranging in complexity from a single compartment Hodgkin-Huxley model to morphologically detailed PV-neuron models. In all models, an increased [Formula: see text] lead to reduced firing, but the experimentally reported increase in [Formula: see text] was not alone sufficient to explain the experimentally reported reduction in firing rate. We therefore hypothesized that PNN degradation in the experiments affected not only [Formula: see text], but also ionic reversal potentials and ion channel conductances. In simulations, we explored how various model parameters affected the firing rate of the model neurons, and identified which parameter variations in addition to [Formula: see text] that are most likely candidates for explaining the experimentally reported reduction in firing rate.

Sedative Properties of Dexmedetomidine Are Mediated Independently from Native Thalamic Hyperpolarization-Activated Cyclic Nucleotide-Gated Channel Function at Clinically Relevant Concentrations

Dexmedetomidine is a selective α₂-adrenoceptor agonist and appears to disinhibit endogenous sleep-promoting pathways, as well as to attenuate noradrenergic excitation. Recent evidence suggests that dexmedetomidine might also directly inhibit hyperpolarization-activated cyclic-nucleotide gated (HCN) channels. We analyzed the effects of dexmedetomidine on native HCN channel function in thalamocortical relay neurons of the ventrobasal complex of the thalamus from mice, performing whole-cell patch-clamp recordings. Over a clinically relevant range of concentrations (1-10 µM), the effects of dexmedetomidine were modest. At a concentration of 10 µM, dexmedetomidine significantly reduced maximal I_h amplitude (relative reduction: 0.86 [0.78-0.91], n = 10, and p = 0.021), yet changes to the half-maximal activation potential V_1/2 occurred exclusively in the presence of the very high concentration of 100 µM (-4,7 [-7.5--4.0] mV, n = 10, and p = 0.009). Coincidentally, only the very high concentration of 100 µM induced a significant deceleration of the fast component of the HCN activation time course (τ_fast: +135.1 [+64.7-+151.3] ms, n = 10, and p = 0.002). With the exception of significantly increasing the membrane input resistance (starting at 10 µM), dexmedetomidine did not affect biophysical membrane properties and HCN channel-mediated parameters of neuronal excitability. Hence, the sedative qualities of dexmedetomidine and its effect on the thalamocortical network are not decisively shaped by direct inhibition of HCN channel function.

Enhancement of Hippocampal-Thalamocortical Temporal Coordination during Slow-Frequency Long-Duration Anterior Thalamic Spindles

Temporal nesting of cortical slow oscillations, thalamic spindles, and hippocampal ripples indicates multiregional neuronal interactions required for memory consolidation. However, how the thalamic activity during spindles organizes hippocampal dynamics remains largely undetermined. We analyzed simultaneous recordings of anterodorsal thalamus and CA1 in male mice to determine the contribution of thalamic spindles in cross-regional synchronization. Our results indicated that temporal hippocampo-thalamocortical coupling was more enhanced during slower and longer thalamic spindles. Additionally, spindles occurring closer to slow oscillation trough were more strongly coupled to ripples. We found that the temporal association between CA1 spiking/ripples and thalamic spindles was stronger following spatial exploration compared with baseline sleep. We further developed a hippocampal-thalamocortical model to explain the mechanism underlying the duration and frequency-dependent coupling of thalamic spindles to hippocampal activity. Our findings shed light on our understanding of the functional role of thalamic activity during spindles on multiregional information transfer.SIGNIFICANCE STATEMENT The contribution of thalamic spindles with differential properties to cross-regional synchronization and information transfer still remains poorly understood. Using simultaneous anterodorsal thalamic and hippocampal recordings from naturally sleeping mice before and after exploration, we found strong coupling of CA1 units to anterodorsal thalamic spindles and increase of this coupling following spatial experience. We further showed that the temporal coupling of CA1 units and hippocampal ripples with thalamic spindles and the spindle-associated modulation of CA1 units with ripples were stronger for spindles with slower frequency of oscillations. Our experimental as well as computational findings using a hippocampal-thalamocortical model provide the first demonstration that spindle frequency and duration can provide valuable information about the underlying multiregional interactions essential for memory consolidation computations.

Human Spindle Variability

In humans, sleep spindles are 10- to 16-Hz oscillations lasting approximately 0.5-2 s. Spindles, along with cortical slow oscillations, may facilitate memory consolidation by enabling synaptic plasticity. Early recordings of spindles at the scalp found anterior channels had overall slower frequency than central-posterior channels. This robust, topographical finding led to dichotomizing spindles as "slow" versus "fast," modeled as two distinct spindle generators in frontal versus posterior cortex. Using a large dataset of intracranial stereoelectroencephalographic (sEEG) recordings from 20 patients (13 female, 7 male) and 365 bipolar recordings, we show that the difference in spindle frequency between frontal and parietal channels is comparable to the variability in spindle frequency within the course of individual spindles, across different spindles recorded by a given site, and across sites within a given region. Thus, fast and slow spindles only capture average differences that obscure a much larger underlying overlap in frequency. Furthermore, differences in mean frequency are only one of several ways that spindles differ. For example, compared with parietal, frontal spindles are smaller, tend to occur after parietal when both are engaged, and show a larger decrease in frequency within-spindles. However, frontal and parietal spindles are similar in being longer, less variable, and more widespread than occipital, temporal, and Rolandic spindles. These characteristics are accentuated in spindles which are highly phase-locked to posterior hippocampal spindles. We propose that rather than a strict parietal-fast/frontal-slow dichotomy, spindles differ continuously and quasi-independently in multiple dimensions, with variability due about equally to within-spindle, within-region, and between-region factors.SIGNIFICANCE STATEMENT Sleep spindles are 10- to 16-Hz neural oscillations generated by cortico-thalamic circuits that promote memory consolidation. Spindles are often dichotomized into slow-anterior and fast-posterior categories for cognitive and clinical studies. Here, we show that the anterior-posterior difference in spindle frequency is comparable to that observed between different cycles of individual spindles, between spindles from a given site, or from different sites within a region. Further, we show that spindles vary on other dimensions such as duration, amplitude, spread, primacy and consistency, and that these multiple dimensions vary continuously and largely independently across cortical regions. These findings suggest that multiple continuous variables rather than a strict frequency dichotomy may be more useful biomarkers for memory consolidation or psychiatric disorders.

Cross-Frequency Slow Oscillation–Spindle Coupling in a Biophysically Realistic Thalamocortical Neural Mass Model

Sleep manifests itself by the spontaneous emergence of characteristic oscillatory rhythms, which often time-lock and are implicated in memory formation. Here, we analyze a neural mass model of the thalamocortical loop in which the cortical node can generate slow oscillations (approximately 1 Hz) while its thalamic component can generate fast sleep spindles of σ-band activity (12–15 Hz). We study the dynamics for different coupling strengths between the thalamic and cortical nodes, for different conductance values of the thalamic node's potassium leak and hyperpolarization-activated cation-nonselective currents, and for different parameter regimes of the cortical node. The latter are listed as follows: (1) a low activity (DOWN) state with noise-induced, transient excursions into a high activity (UP) state, (2) an adaptation induced slow oscillation limit cycle with alternating UP and DOWN states, and (3) a high activity (UP) state with noise-induced, transient excursions into the low activity (DOWN) state. During UP states, thalamic spindling is abolished or reduced. During DOWN states, the thalamic node generates sleep spindles, which in turn can cause DOWN to UP transitions in the cortical node. Consequently, this leads to spindle-induced UP state transitions in parameter regime (1), thalamic spindles induced in some but not all DOWN states in regime (2), and thalamic spindles following UP to DOWN transitions in regime (3). The spindle-induced σ-band activity in the cortical node, however, is typically the strongest during the UP state, which follows a DOWN state “window of opportunity” for spindling. When the cortical node is parametrized in regime (3), the model well explains the interactions between slow oscillations and sleep spindles observed experimentally during Non-Rapid Eye Movement sleep. The model is computationally efficient and can be integrated into large-scale modeling frameworks to study spatial aspects like sleep wave propagation.

Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia

In the arousal process, the brain restores its integrative activity from the synchronized state of slow wave activity (SWA). The mechanisms underpinning this state transition remain, however, to be elucidated. Here we simultaneously probed neuronal assemblies throughout the whole cortex with micro-electrocorticographic recordings in mice. We investigated the progressive shaping of propagating SWA at different levels of isoflurane. We found a form of memory of the wavefront shapes at deep anesthesia, tightly alternating posterior-anterior-posterior patterns. At low isoflurane, metastable patterns propagated in more directions, reflecting an increased complexity. The wandering across these mesostates progressively increased its randomness, as predicted by simulations of a network of spiking neurons, and confirmed in our experimental data. The complexity increase is explained by the elevated excitability of local assemblies with no modifications of the network connectivity. These results shed new light on the functional reorganization of the cortical network as anesthesia fades out.

•

Complexity of isoflurane-induced slow waves reliably determines anesthesia level

•

In deep anesthesia, the propagation strictly alternates between front-back-front patterns

•

In light anesthesia, there is a continuum of directions and faster propagation

•

Local excitability underpins the cortical reorganization in fading anesthesia

Can Corticomuscular Coherence Differentiate between REM Sleep Behavior Disorder with or without Parkinsonism?

REM sleep behavior disorder (RBD) could be a predictor of Parkinsonism even before development of typical motor symptoms. This study aims to characterize clinical features and corticomuscular and corticocortical coherence (CMC and CCC, respectively) during sleep in RBD patients with or without Parkinsonism. We enrolled a total of 105 subjects, including 20 controls, 54 iRBD, and 31 RBD+P patients, patients who were diagnosed as idiopathic RBD (iRBD) and RBD with Parkinsonism (RBD+P) in our neurology department. We analyzed muscle atonia index (MAI) and CMC between EEG and chin/limb muscle electromyography (EMG) and CCC during different sleep stages. Although differences in the CMC of iRBD group were observed only during REM sleep, MAI differences between groups were noted during both REM and NREM N2 stage sleep. During REM sleep, CMC was higher and MAI was reduced in iRBD patients compared to controls (p = 0.001, p < 0.001, respectively). Interestingly, MAI was more reduced in RBD+P compared to iRBD patients. In comparison, CCC was higher in iRBD patients compared to controls whereas CCC was lower in RBD+P groups compared to control and iRBD groups in various frequency bands during both NREM N2 and REM sleep stages. Among them, increased CMC during REM sleep revealed correlation between clinical severities of RBD symptoms. Our findings indicate that MAI, CMC, and CCC showed distinctive features in iRBD and RBD+P patients compared to controls, suggesting potential usefulness to understand possible links between these diseases.

Impact of In Utero Exposure to Antiepileptic Drugs on Neonatal Brain Function

In utero brain development underpins brain health across the lifespan but is vulnerable to physiological and pharmacological perturbation. Here, we show that antiepileptic medication during pregnancy impacts on cortical activity during neonatal sleep, a potent indicator of newborn brain health. These effects are evident in frequency-specific functional brain networks and carry prognostic information for later neurodevelopment. Notably, such effects differ between different antiepileptic drugs that suggest neurodevelopmental adversity from exposure to antiepileptic drugs and not maternal epilepsy per se. This work provides translatable bedside metrics of brain health that are sensitive to the effects of antiepileptic drugs on postnatal neurodevelopment and carry direct prognostic value.

The essential role of hippocampo-cortical connections in temporal coordination of spindles and ripples

The predominant activity of slow wave sleep is cortical slow oscillations (SOs), thalamic spindles and hippocampal sharp wave ripples. While the precise temporal nesting of these rhythms was shown to be essential for memory consolidation, the coordination mechanism is poorly understood. Here we develop a minimal hippocampo-cortico-thalamic network that can explain the mechanism underlying the SO-spindle-ripple coupling indicating of the succession of regional neuronal interactions. Further we verify the model predictions experimentally in naturally sleeping rodents showing our simple model provides a quantitative match to several experimental observations including the nesting of ripples in the spindle troughs and larger duration but lower amplitude of the ripples co-occurring with spindles or SOs compared to the isolated ripples. The model also predicts that the coupling of ripples to SOs and spindles monotonically enhances by increasing the strength of hippocampo-cortical connections while it is stronger at intermediate values of the cortico-hippocampal projections.

The impact of a closed-loop thalamocortical model on the spatiotemporal dynamics of cortical and thalamic traveling waves

Propagation of activity in spatially structured neuronal networks has been observed in awake, anesthetized, and sleeping brains. How these wave patterns emerge and organize across brain structures, and how network connectivity affects spatiotemporal neural activity remains unclear. Here, we develop a computational model of a two-dimensional thalamocortical network, which gives rise to emergent traveling waves similar to those observed experimentally. We illustrate how spontaneous and evoked oscillatory activity in space and time emerge using a closed-loop thalamocortical architecture, sustaining smooth waves in the cortex and staggered waves in the thalamus. We further show that intracortical and thalamocortical network connectivity, cortical excitation/inhibition balance, and thalamocortical or corticothalamic delay can independently or jointly change the spatiotemporal patterns (radial, planar and rotating waves) and characteristics (speed, direction, and frequency) of cortical and thalamic traveling waves. Computer simulations predict that increased thalamic inhibition induces slower cortical frequencies and that enhanced cortical excitation increases traveling wave speed and frequency. Overall, our results provide insight into the genesis and sustainability of thalamocortical spatiotemporal patterns, showing how simple synaptic alterations cause varied spontaneous and evoked wave patterns. Our model and simulations highlight the need for spatially spread neural recordings to uncover critical circuit mechanisms for brain functions.

A novel density-based neural mass model for simulating neuronal network dynamics with conductance-based synapses and membrane current adaptation

Despite its success in understanding brain rhythms, the neural mass model, as a low-dimensional mean-field network model, is phenomenological in nature, so that it cannot replicate some of rich repertoire of responses seen in real neuronal tissues. Here, using a colored-synapse population density method, we derived a novel neural mass model, termed density-based neural mass model (dNMM), as the mean-field description of network dynamics of adaptive exponential integrate-and-fire (aEIF) neurons, in which two critical neuronal features, i.e., voltage-dependent conductance-based synaptic interactions and adaptation of firing rate responses, were included. Our results showed that the dNMM was capable of correctly estimating firing rate responses of a neuronal population of aEIF neurons receiving stationary or time-varying excitatory and inhibitory inputs. Finally, it was also able to quantitatively describe the effect of spike-frequency adaptation in the generation of asynchronous irregular activity of excitatory-inhibitory cortical networks. We conclude that in terms of its biological reality and calculation efficiency, the dNMM is a suitable candidate to build significantly large-scale network models involving multiple brain areas, where the neuronal population is the smallest dynamic unit.

Travelling spindles create necessary conditions for spike-timing-dependent plasticity in humans

Sleep spindles facilitate memory consolidation in the cortex during mammalian non-rapid eye movement sleep. In rodents, phase-locked firing during spindles may facilitate spike-timing-dependent plasticity by grouping pre-then-post-synaptic cell firing within ~25 ms. Currently, microphysiological evidence in humans for conditions conducive for spike-timing-dependent plasticity during spindles is absent. Here, we analyze field potentials and unit firing from middle/upper layers during spindles from 10 × 10 microelectrode arrays at 400 μm pitch in humans. We report strong tonic and phase-locked increases in firing and co-firing within 25 ms during spindles, especially those co-occurring with down-to-upstate transitions. Co-firing, spindle co-occurrence, and spindle coherence are greatest within ~2 mm, and high co-firing of units on different contacts depends on high spindle coherence between those contacts. Spindles propagate at ~0.28 m/s in distinct patterns, with correlated cell co-firing sequences. Spindles hence organize spatiotemporal patterns of neuronal co-firing in ways that may provide pre-conditions for plasticity during non-rapid eye movement sleep.

Sleep spindles during non-rapid eye movement are important for memory consolidation and require specific neuronal firing conditions in non-human mammals. Here, the authors show these conditions are present in humans, potentially facilitating spike-timing-dependent plasticity.

Attenuation of Native Hyperpolarization-Activated, Cyclic Nucleotide-Gated Channel Function by the Volatile Anesthetic Sevoflurane in Mouse Thalamocortical Relay Neurons

As thalamocortical relay neurons are ascribed a crucial role in signal propagation and information processing, they have attracted considerable attention as potential targets for anesthetic modulation. In this study, we analyzed the effects of different concentrations of sevoflurane on the excitability of thalamocortical relay neurons and hyperpolarization-activated, cyclic-nucleotide gated (HCN) channels, which play a decisive role in regulating membrane properties and rhythmic oscillatory activity. The effects of sevoflurane on single-cell excitability and native HCN channels were investigated in acutely prepared brain slices from adult wild-type mice with the whole-cell patch-clamp technique, using voltage-clamp and current-clamp protocols. Sevoflurane dose-dependently depressed membrane biophysics and HCN-mediated parameters of neuronal excitability. Respective half-maximal inhibitory and effective concentrations ranged between 0.30 (95% CI, 0.18–0.50) mM and 0.88 (95% CI, 0.40–2.20) mM. We witnessed a pronounced reduction of HCN dependent I_h current amplitude starting at a concentration of 0.45 mM [relative change at −133 mV; 0.45 mM sevoflurane: 0.85 (interquartile range, 0.79–0.92), n = 12, p = 0.011; 1.47 mM sevoflurane: 0.37 (interquartile range, 0.34–0.62), n = 5, p < 0.001] with a half-maximal inhibitory concentration of 0.88 (95% CI, 0.40–2.20) mM. In contrast, effects on voltage-dependent channel gating were modest with significant changes only occurring at 1.47 mM [absolute change of half-maximal activation potential; 1.47 mM: −7.2 (interquartile range, −10.3 to −5.8) mV, n = 5, p = 0.020]. In this study, we demonstrate that sevoflurane inhibits the excitability of thalamocortical relay neurons in a concentration-dependent manner within a clinically relevant range. Especially concerning its effects on native HCN channel function, our findings indicate substance-specific differences in comparison to other anesthetic agents. Considering the importance of HCN channels, the observed effects might mechanistically contribute to the hypnotic properties of sevoflurane.

Synchronous Infra-Slow Oscillations Organize Ensembles of Accessory Olfactory Bulb Projection Neurons into Distinct Microcircuits

The accessory olfactory system controls social and sexual behavior. In the mouse accessory olfactory bulb, the first central stage of information processing along the accessory olfactory pathway, projection neurons (mitral cells) display infra-slow oscillatory discharge with remarkable periodicity. The physiological mechanisms that underlie this default output state, however, remain controversial. Moreover, whether such rhythmic infra-slow activity patterns exist in awake behaving mice and whether such activity reflects the functional organization of the accessory olfactory bulb circuitry remain unclear. Here, we hypothesize that mitral cell ensembles form synchronized microcircuits that subdivide the accessory olfactory bulb into segregated functional clusters. We use a miniature microscope to image the Ca²⁺ dynamics within the apical dendritic compartments of large mitral cell ensembles in vivo. We show that infra-slow periodic patterns of concerted neural activity, indeed, reflect the idle state of accessory olfactory bulb output in awake male and female mice. Ca²⁺ activity profiles are distinct and glomerulus-specific. Confocal time-lapse imaging in acute slices reveals that groups of mitral cells assemble into microcircuits that exhibit correlated Ca²⁺ signals. Moreover, electrophysiological profiling of synaptic connectivity indicates functional coupling between mitral cells. Our results suggest that both intrinsically rhythmogenic neurons and neurons entrained by fast synaptic drive are key elements in organizing the accessory olfactory bulb into functional microcircuits, each characterized by a distinct default pattern of infra-slow rhythmicity.

SIGNIFICANCE STATEMENT Information processing in the accessory olfactory bulb (AOB) plays a central role in conspecific chemosensory communication. Surprisingly, many basic physiological principles that underlie neuronal signaling in the AOB remain elusive. Here, we show that AOB projection neurons (mitral cells) form parallel synchronized ensembles both in vitro and in vivo. Infra-slow synchronous oscillatory activity within AOB microcircuits thus adds a new dimension to chemosensory coding along the accessory olfactory pathway.

Signal Propagation via Open-Loop Intrathalamic Architectures: A Computational Model

Propagation of signals across the cerebral cortex is a core component of many cognitive processes and is generally thought to be mediated by direct intracortical connectivity. The thalamus, by contrast, is considered to be devoid of internal connections and organized as a collection of parallel inputs to the cortex. Here, we provide evidence that "open-loop" intrathalamic pathways involving the thalamic reticular nucleus (TRN) can support propagation of oscillatory activity across the cortex. Recent studies support the existence of open-loop thalamo-reticulo-thalamic (TC-TRN-TC) synaptic motifs in addition to traditional closed-loop architectures. We hypothesized that open-loop structural modules, when connected in series, might underlie thalamic and, therefore cortical, signal propagation. Using a supercomputing platform to simulate thousands of permutations of a thalamocortical network based on physiological data collected in mice, rats, ferrets, and cats and in which select synapses were allowed to vary both by class and individually, we evaluated the relative capacities of closed-loop and open-loop TC-TRN-TC synaptic configurations to support both propagation and oscillation. We observed that (1) signal propagation was best supported in networks possessing strong open-loop TC-TRN-TC connectivity; (2) intrareticular synapses were neither primary substrates of propagation nor oscillation; and (3) heterogeneous synaptic networks supported more robust propagation of oscillation than their homogeneous counterparts. These findings suggest that open-loop, heterogeneous intrathalamic architectures might complement direct intracortical connectivity to facilitate cortical signal propagation.

A Computational Modeling Reveals That Strength of Inhibitory Input, E/I Balance, and Distance of Excitatory Input Modulate Thalamocortical Bursting Properties

The thalamus is a brain structure known to modulate sensory information before relaying to the cortex. The unique ability of a thalamocortical (TC) neuron to switch between the high frequency burst firing and single spike tonic firing has been implicated to have a key role in sensory modulation including pain. Of the two firing modes, burst firing, especially maintaining certain burst firing properties, was suggested to be critical in controlling nociceptive behaviors. Therefore, understanding the factors that influence burst firing properties would offer important insight into understanding sensory modulation. Using computational modeling, we investigated how the balance of excitatory and inhibitory inputs into a TC neuron influence TC bursting properties. We found that intensity of inhibitory inputs and the timing of excitatory input delivery control the dynamics of bursting properties. Then, to reflect a more realistic model, excitatory inputs delivered at different dendritic locations—proximal, intermediate, or distal—of a TC neuron were also investigated. Interestingly, excitatory input delivered into a distal dendrite, despite the furthest distance, had the strongest influence in shaping burst firing properties, suggesting that not all inputs equally contribute to modulating TC bursting properties. Overall, the results provide computational insights in understanding the detailed mechanism of the factors influencing temporal pattern of thalamic bursts.

Sleep Spindle-dependent Functional Connectivity Correlates with Cognitive Abilities

EEG studies have shown that interindividual differences in the electrophysiological properties of sleep spindles (e.g., density, amplitude, duration) are highly correlated with trait-like "reasoning" abilities (i.e., "fluid intelligence"; problem-solving skills; the ability to employ logic or identify complex patterns), but not interindividual differences in STM or "verbal" intellectual abilities. Previous simultaneous EEG-fMRI studies revealed brain activations time-locked to spindles. Our group has recently demonstrated that the extent of activation in a subset of these regions was related to interindividual differences in reasoning intellectual abilities, specifically. However, spindles reflect communication between spatially distant and functionally distinct brain areas. The functional communication among brain regions related to spindles and their relationship to reasoning abilities have yet to be investigated. Using simultaneous EEG-fMRI sleep recordings and psychophysiological interaction analysis, we identified spindle-related functional communication among brain regions in the thalamo-cortical-BG system, the salience network, and the default mode network. Furthermore, the extent of the functional connectivity of the cortical-striatal circuitry and the thalamo-cortical circuitry was specifically related to reasoning abilities but was unrelated to STM or verbal abilities, thus suggesting that individuals with higher fluid intelligence have stronger functional coupling among these brain areas during spontaneous spindle events. This may serve as a first step in further understanding the function of sleep spindles and the brain network functional communication, which support the capacity for fluid intelligence.

Sleep Modelling across Physiological Levels

Sleep and circadian rhythms are regulated across multiple functional, spatial and temporal levels: from genes to networks of coupled neurons and glial cells, to large scale brain dynamics and behaviour. The dynamics at each of these levels are complex and the interaction between the levels is even more so, so research have mostly focused on interactions within the levels to understand the underlying mechanisms—the so-called reductionist approach. Mathematical models were developed to test theories of sleep regulation and guide new experiments at each of these levels and have become an integral part of the field. The advantage of modelling, however, is that it allows us to simulate and test the dynamics of complex biological systems and thus provides a tool to investigate the connections between the different levels and study the system as a whole. In this paper I review key models of sleep developed at different physiological levels and discuss the potential for an integrated systems biology approach for sleep regulation across these levels. I also highlight the necessity of building mechanistic connections between models of sleep and circadian rhythms across these levels.

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.

Thalamocortical control of propofol phase-amplitude coupling

The anesthetic propofol elicits many different spectral properties on the EEG, including alpha oscillations (8-12 Hz), Slow Wave Oscillations (SWO, 0.1-1.5 Hz), and dose-dependent phase-amplitude coupling (PAC) between alpha and SWO. Propofol is known to increase GABAA inhibition and decrease H-current strength, but how it generates these rhythms and their interactions is still unknown. To investigate both generation of the alpha rhythm and its PAC to SWO, we simulate a Hodgkin-Huxley network model of a hyperpolarized thalamus and corticothalamic inputs. We find, for the first time, that the model thalamic network is capable of independently generating the sustained alpha seen in propofol, which may then be relayed to cortex and expressed on the EEG. This dose-dependent sustained alpha critically relies on propofol GABAA potentiation to alter the intrinsic spindling mechanisms of the thalamus. Furthermore, the H-current conductance and background excitation of these thalamic cells must be within specific ranges to exhibit any intrinsic oscillations, including sustained alpha. We also find that, under corticothalamic SWO UP and DOWN states, thalamocortical output can exhibit maximum alpha power at either the peak or trough of this SWO; this implies the thalamus may be the source of propofol-induced PAC. Hyperpolarization level is the main determinant of whether the thalamus exhibits trough-max PAC, which is associated with lower propofol dose, or peak-max PAC, associated with higher dose. These findings suggest: the thalamus generates a novel rhythm under GABAA potentiation such as under propofol, its hyperpolarization may determine whether a patient experiences trough-max or peak-max PAC, and the thalamus is a critical component of propofol-induced cortical spectral phenomena. Changes to the thalamus may be a critical part of how propofol accomplishes its effects, including unconsciousness.

Is more always better? How different 'doses' of exercise after incomplete spinal cord injury affects the membrane properties of deep dorsal horn interneurons

Interneurons in the deep dorsal horn (DDH) of the spinal cord process somatosensory input, and form an important link between upper and lower motoneurons to subsequently shape motor output. Exercise training after SCI is known to improve functional motor recovery, but little is known about the mechanisms within spinal cord neurons that underlie these improvements. Here we investigate how the properties of DDH interneurons are affected by spinal cord injury (SCI) alone, and SCI in combination with different 'doses' of treadmill exercise training (3, 6, and 9wks). In an adult mouse hemisection model of SCI we used whole-cell patch-clamp electrophysiology to record intrinsic, AP firing and gain modulation properties from DDH interneurons in a horizontal spinal cord slice preparation. We find that neurons within two segments of the injury, both ipsi- and contralateral to the hemisection, are similarly affected by SCI and SCI plus exercise. The passive intrinsic membrane properties input resistance (R_in) and rheobase are sensitive to the effects of recovery time and exercise training after SCI thus altering DDH interneuron excitability. Conversely, select active membrane properties are largely unaffected by either SCI or exercise training. SCI itself causes a mismatch in the expression of voltage-gated subthreshold currents and AP discharge firing type. Over time after SCI, and especially with exercise training (9wks), this mismatched expression is exacerbated. Lastly, amplification properties (i.e. gain of frequency-current relationship) of DDH interneurons are altered by SCI alone and recover spontaneously with no clear effect of exercise training. These results suggest a larger 'dose' of exercise training (9wks) has a strong and selective effect on specific membrane properties, and on the output of interneurons in the vicinity of a SCI. These electrophysiological data provide new insights into the plasticity of DDH interneurons and the mechanisms by which exercise therapy after SCI can improve recovery.

Using Oscillating Sounds to Manipulate Sleep Spindles

Introduction

EEG oscillations known as sleep spindles have been linked with various aspects of cognition, but the specific functions they signal remain controversial. Two types of EEG sleep spindles have been distinguished: slow spindles at 11-13.5 Hz and fast spindles at 13.5-16 Hz. Slow spindles exhibit a frontal scalp topography, whereas fast spindles exhibit a posterior scalp topography and have been preferentially linked with memory consolidation during sleep. To advance understanding beyond that provided from correlative studies of spindles, we aimed to develop a new method to systematically manipulate spindles.

Aims and Methods

We presented repeating bursts of oscillating white noise to people during a 90-min afternoon nap. During stage 2 and slow-wave sleep, oscillations were embedded within contiguous 10-s stimulation intervals, each comprising 2 s of white noise amplitude modulated at 12 Hz (targeting slow spindles), 15 Hz (targeting fast spindles), or 50 Hz followed by 8 s of constant white noise.

Results

During oscillating stimulation compared to constant stimulation, parietal EEG recordings showed more slow spindles in the 12-Hz condition, more fast spindles in the 15-Hz condition, and no change in the 50-Hz control condition. These effects were topographically selective, and were absent in frontopolar EEG recordings, where slow spindle density was highest. Spindles during stimulation were similar to spontaneous spindles in standard physiological features, including duration and scalp distribution.

Conclusions

These results define a new method to selectively and noninvasively manipulate spindles through acoustic resonance, while also providing new evidence for functional distinctions between the 2 types of EEG spindles.

Hippocampal gamma-slow oscillation coupling in macaques during sedation and sleep

Behavioral and neurophysiological evidence suggests that the slow (≤1 Hz) oscillation (SO) during sleep plays a role in consolidating hippocampal (HIPP)-dependent memories. The effects of the SO on HIPP activity have been studied in rodents and cats both during natural sleep and during anesthetic administration titrated to mimic sleep-like slow rhythms. In this study, we sought to document these effects in primates. First, HIPP field potentials were recorded during ketamine-dexmedetomidine sedation and during natural sleep in three rhesus macaques. Sedation produced regionally-specific slow and gamma (∼40 Hz) oscillations with strong coupling between the SO phase and gamma amplitude. These same features were seen in slow-wave sleep (SWS), but the coupling was weaker and the coupled gamma oscillation had a higher frequency (∼70 Hz) during SWS. Second, electrical stimuli were delivered to HIPP afferents in the parahippocampal gyrus (PHG) during sedation to assess the effects of sleep-like SO on excitability. Gamma bursts after the peak of SO cycles corresponded to periods of increased gain of monosynaptic connections between the PHG and HIPP. However, the two PHG-HIPP connectivity gains during sedation were both substantially lower than when the animal was awake. We conclude that the SO is correlated with rhythmic excitation and inhibition of the PHG-HIPP network, modulating connectivity and gamma generators intrinsic to this network. Ketamine-dexmedetomidine sedation produces a similar effect, but with a decreased contribution of the PHG to HIPP activity and gamma generation.

Increased Amplitude of Thalamocortical Low-Frequency Oscillations in Patients with Migraine

For many years, neurobiological theories have emphasized the importance of neuronal oscillations in the emergence of brain function. At the same time, clinical studies have shown that disturbances or irregularities in brain rhythms may relate to various common neurological conditions, including migraine. Increasing evidence suggests that the CNS plays a fundamental role in the predisposition to develop different forms of headache. Here, we present human imaging data that strongly support the presence of abnormal low-frequency oscillations (LFOs) in thalamocortical networks of patients in the interictal phase of migraine. Our results show that the main source of arrhythmic activity was localized to the higher-order thalamic relays of the medial dorsal nucleus. In addition, spontaneous LFOs in the thalamus were selectively associated with the headache attack frequency, meaning that the varying amplitude of dysrhythmia could predispose patients to recurrent attacks. Rhythmic cortical feedback to the thalamus is a major factor in the amplification of thalamocortical oscillations, making it a strong candidate for influencing neuronal excitability. We further speculate that the intrinsic dynamics of thalamocortical network oscillations are crucial for early sensory processing and therefore could underlie important pathophysiological processes involved in multisensory integration.

SIGNIFICANCE STATEMENT

In many cases, migraine attacks are thought to begin centrally. A major obstacle to studying intrinsic brain activity has been the identification of the precise anatomical structures and functional networks that are involved in migraine. Here, we present imaging data that strongly support the presence of abnormal low-frequency oscillations in thalamocortical networks of patients in the interictal phase of migraine. This arrhythmic activity was localized to the higher-order thalamic relays of the medial dorsal nucleus and was selectively associated with headache attack frequency. Rhythmic cortical feedback to the thalamus is a major factor in the amplification of thalamocortical oscillations, making it a strong candidate for influencing neuronal excitability and higher-level processes involved in multisensory integration.

Visually Evoked 3-5 Hz Membrane Potential Oscillations Reduce the Responsiveness of Visual Cortex Neurons in Awake Behaving Mice

Low-frequency membrane potential (V_m) oscillations were once thought to only occur in sleeping and anesthetized states. Recently, low-frequency V_m oscillations have been described in inactive awake animals, but it is unclear whether they shape sensory processing in neurons and whether they occur during active awake behavioral states. To answer these questions, we performed two-photon guided whole-cell V_m recordings from primary visual cortex layer 2/3 excitatory and inhibitory neurons in awake mice during passive visual stimulation and performance of visual and auditory discrimination tasks. We recorded stereotyped 3-5 Hz V_m oscillations where the V_m baseline hyperpolarized as the V_m underwent high amplitude rhythmic fluctuations lasting 1-2 s in duration. When 3-5 Hz V_m oscillations coincided with visual cues, excitatory neuron responses to preferred cues were significantly reduced. Despite this disruption to sensory processing, visual cues were critical for evoking 3-5 Hz V_m oscillations when animals performed discrimination tasks and passively viewed drifting grating stimuli. Using pupillometry and animal locomotive speed as indicators of arousal, we found that 3-5 Hz oscillations were not restricted to unaroused states and that they occurred equally in aroused and unaroused states. Therefore, low-frequency V_m oscillations play a role in shaping sensory processing in visual cortical neurons, even during active wakefulness and decision making.SIGNIFICANCE STATEMENT A neuron's membrane potential (V_m) strongly shapes how information is processed in sensory cortices of awake animals. Yet, very little is known about how low-frequency V_m oscillations influence sensory processing and whether they occur in aroused awake animals. By performing two-photon guided whole-cell recordings from layer 2/3 excitatory and inhibitory neurons in the visual cortex of awake behaving animals, we found visually evoked stereotyped 3-5 Hz V_m oscillations that disrupt excitatory responsiveness to visual stimuli. Moreover, these oscillations occurred when animals were in high and low arousal states as measured by animal speed and pupillometry. These findings show, for the first time, that low-frequency V_m oscillations can significantly modulate sensory signal processing, even in awake active animals.

Effects of Small-World Rewiring Probability and Noisy Synaptic Conductivity on Slow Waves: Cortical Network

Physiological rhythms play a critical role in the functional development of living beings. Many biological functions are executed with an interaction of rhythms produced by internal characteristics of scores of cells. While synchronized oscillations may be associated with normal brain functions, anomalies in these oscillations may cause or relate the emergence of some neurological or neuropsychological pathologies. This study was designed to investigate the effects of topological structure and synaptic conductivity noise on the spatial synchronization and temporal rhythmicity of the waves generated by cells in the network. Because of holding the ability of clustering and randomizing with change of parameters, small-world (SW) network topology was chosen. The oscillatory activity of network was tried out by manipulating an insulated SW, cortical network model whose morphology is very close to real world. According to the obtained results, it was observed that at the optimal probabilistic rates of conductivity noise and rewiring of SW, powerful synchronized oscillatory small waves are generated in relation to the internal dynamics of cells, which are in line with the network's input. These two parameters were observed to be quite effective on the excitation-inhibition balance of the network. Accordingly, it may be suggested that the topological dynamics of SW and noisy synaptic conductivity may be associated with the normal and abnormal development of neurobiological structure.

The Role of Neuromodulators in Cortical Plasticity. A Computational Perspective

Neuromodulators play a ubiquitous role across the brain in regulating plasticity. With recent advances in experimental techniques, it is possible to study the effects of diverse neuromodulatory states in specific brain regions. Neuromodulators are thought to impact plasticity predominantly through two mechanisms: the gating of plasticity and the upregulation of neuronal activity. However, the consequences of these mechanisms are poorly understood and there is a need for both experimental and theoretical exploration. Here we illustrate how neuromodulatory state affects cortical plasticity through these two mechanisms. First, we explore the ability of neuromodulators to gate plasticity by reshaping the learning window for spike-timing-dependent plasticity. Using a simple computational model, we implement four different learning rules and demonstrate their effects on receptive field plasticity. We then compare the neuromodulatory effects of upregulating learning rate versus the effects of upregulating neuronal activity. We find that these seemingly similar mechanisms do not yield the same outcome: upregulating neuronal activity can lead to either a broadening or a sharpening of receptive field tuning, whereas upregulating learning rate only intensifies the sharpening of receptive field tuning. This simple model demonstrates the need for further exploration of the rich landscape of neuromodulator-mediated plasticity. Future experiments, coupled with biologically detailed computational models, will elucidate the diversity of mechanisms by which neuromodulatory state regulates cortical plasticity.

Microbial Proteins as Novel Industrial Biotechnology Hosts to Treat Epilepsy

Epilepsy is characterized by the hyperexcitability of various neuronal circuits that results due to the imbalance between glutamate-mediated excitation of voltage-gated cation channels and γ-amino butyric acid (GABA)-mediated inhibition of anion channels leading to aberrant, sporadic oscillations or fluctuations in neuronal electrical activity. Epilepsy with a risk of mortality and around 65 million sufferers of all ages all over the world is limited therapeutically with high rates of adverse reactions, lack of complete seizure control, and over 30% patients with refractory epilepsy. The only alternative to medicines is to identify and surgically remove the seizure foci in the brain or to abort the seizures just as they begin using an implanted cerebral electrode. However, these alternatives are unable to precisely aim aberrant neuronal circuits while leaving others unaltered. Epilepsy animal models also constitute the identical constraint. Thus, a better target-specific approach is needed to study and treat epilepsy. Unicellular green algae Chlamydomonas reinhardtii expresses a channelrhodopsin-2 (ChR2) sodium ion channel protein that controls the phototaxis movement of algae in response to blue light. Similarly, archaeon Natronomonas pharaonis (NpHR) expresses a monovalent Cl^- channel protein halorhodopsin that responds to yellow light. These features of ChR2 and NpHR proteins can be used in optogenetic techniques to manipulate the bi-directional firing pattern of neuronal circuits in an attempt to better understand the pathophysiology of epileptic seizures as well as to discover novel potential drugs to treat epilepsy.

Long-range synchrony and emergence of neural reentry

Neural synchronization across long distances is a functionally important phenomenon in health and disease. In order to access the basis of different modes of long-range synchrony, we monitor spiking activities over centimetre scale in cortical networks and show that the mode of synchrony depends upon a length scale, λ, which is the minimal path that activity should propagate through to find its point of origin ready for reactivation. When λ is larger than the physical dimension of the network, distant neuronal populations operate synchronously, giving rise to irregularly occurring network-wide events that last hundreds of milliseconds to several seconds. In contrast, when λ approaches the dimension of the network, a continuous self-sustained reentry propagation emerges, a regular seizure-like mode that is marked by precise spatiotemporal patterns ('synfire chains') and may last many minutes. Termination of a reentry phase is preceded by a decrease of propagation speed to a halt. Stimulation decreases both propagation speed and λ values, which modifies the synchrony mode respectively. The results contribute to the understanding of the origin and termination of different modes of neural synchrony as well as their long-range spatial patterns, while hopefully catering to manipulation of the phenomena in pathological conditions.

Spatiotemporal properties of microsaccades: Model predictions and experimental tests

Microsaccades are involuntary and very small eye movements during fixation. Recently, the microsaccade-related neural dynamics have been extensively investigated both in experiments and by constructing neural network models. Experimentally, microsaccades also exhibit many behavioral properties. It's well known that the behavior properties imply the underlying neural dynamical mechanisms, and so are determined by neural dynamics. The behavioral properties resulted from neural responses to microsaccades, however, are not yet understood and are rarely studied theoretically. Linking neural dynamics to behavior is one of the central goals of neuroscience. In this paper, we provide behavior predictions on spatiotemporal properties of microsaccades according to microsaccade-induced neural dynamics in a cascading network model, which includes both retinal adaptation and short-term depression (STD) at thalamocortical synapses. We also successfully give experimental tests in the statistical sense. Our results provide the first behavior description of microsaccades based on neural dynamics induced by behaving activity, and so firstly link neural dynamics to behavior of microsaccades. These results indicate strongly that the cascading adaptations play an important role in the study of microsaccades. Our work may be useful for further investigations of the microsaccadic behavioral properties and of the underlying neural dynamical mechanisms responsible for the behavioral properties.

Callosal dysfunction explains injury sequelae in a computational network model of axonal injury

Mild traumatic brain injury (mTBI) often results in neurobehavioral aberrations such as impaired attention and increased reaction time. Diffusion imaging and postmortem analysis studies suggest that mTBI primarily affects myelinated axons in white matter tracts. In particular, corpus callosum, mediating interhemispheric information exchange, has been shown to be affected in mTBI. Yet little is known about the mechanisms linking the injury of myelinated callosal axons to the neurobehavioral sequelae of mTBI. To address this issue, we devised and studied a large, biologically plausible neuronal network model of cortical tissue. Importantly, the model architecture incorporated intra- and interhemispheric organization, including myelinated callosal axons and distance-dependent axonal conduction delays. In the resting state, the intact model network exhibited several salient features, including alpha-band (8-12 Hz) collective activity with low-frequency irregular spiking of individual neurons. The network model of callosal injury captured several clinical observations, including 1) "slowing down" of the network rhythms, manifested as an increased resting-state theta-to-alpha power ratio, 2) reduced response to attention-like network stimulation, manifested as a reduced spectral power of collective activity, and 3) increased population response time in response to stimulation. Importantly, these changes were positively correlated with injury severity, supporting proposals to use neurobehavioral indices as biomarkers for determining the severity of injury. Our modeling effort helps to understand the role played by the injury of callosal myelinated axons in defining the neurobehavioral sequelae of mTBI.

Transition between Functional Regimes in an Integrate-And-Fire Network Model of the Thalamus

The thalamus is a key brain element in the processing of sensory information. During the sleep and awake states, this brain area is characterized by the presence of two distinct dynamical regimes: in the sleep state activity is dominated by spindle oscillations (7 − 15 Hz) weakly affected by external stimuli, while in the awake state the activity is primarily driven by external stimuli. Here we develop a simple and computationally efficient model of the thalamus that exhibits two dynamical regimes with different information-processing capabilities, and study the transition between them. The network model includes glutamatergic thalamocortical (TC) relay neurons and GABAergic reticular (RE) neurons described by adaptive integrate-and-fire models in which spikes are induced by either depolarization or hyperpolarization rebound. We found a range of connectivity conditions under which the thalamic network composed by these neurons displays the two aforementioned dynamical regimes. Our results show that TC-RE loops generate spindle-like oscillations and that a minimum level of clustering (i.e. local connectivity density) in the RE-RE connections is necessary for the coexistence of the two regimes. We also observe that the transition between the two regimes occurs when the external excitatory input on TC neurons (mimicking sensory stimulation) is large enough to cause a significant fraction of them to switch from hyperpolarization-rebound-driven firing to depolarization-driven firing. Overall, our model gives a novel and clear description of the role that the two types of neurons and their connectivity play in the dynamical regimes observed in the thalamus, and in the transition between them. These results pave the way for the development of efficient models of the transmission of sensory information from periphery to cortex.

A Thalamocortical Neural Mass Model of the EEG during NREM Sleep and Its Response to Auditory Stimulation

Few models exist that accurately reproduce the complex rhythms of the thalamocortical system that are apparent in measured scalp EEG and at the same time, are suitable for large-scale simulations of brain activity. Here, we present a neural mass model of the thalamocortical system during natural non-REM sleep, which is able to generate fast sleep spindles (12–15 Hz), slow oscillations (<1 Hz) and K-complexes, as well as their distinct temporal relations, and response to auditory stimuli. We show that with the inclusion of detailed calcium currents, the thalamic neural mass model is able to generate different firing modes, and validate the model with EEG-data from a recent sleep study in humans, where closed-loop auditory stimulation was applied. The model output relates directly to the EEG, which makes it a useful basis to develop new stimulation protocols.

Sleep plays a pivotal role for the consolidation of memory. While REM sleep had originally been the focus of research due to its similarity with wakefulness, more recent studies suggest that different sleep stages are responsible for the consolidation of different types of memory. To better understand the changes in neuronal dynamics between the different sleep stages, neural mass models are a valuable tool, as they relate directly to the large-scale dynamics measured by an EEG. Here, we present a model of the sleeping thalamocortical system. We first show that the isolated thalamic submodule is able to generate different oscillatory behavior found in vivo. Furthermore, the full thalamocortical model reproduces the EEG of sleep stages N2 and N3 and preserves the temporal relationship between cortical K-complexes/slow oscillations and thalamic sleep spindles. A comparison with event related potential data from a recent sleep study in humans demonstrates its possible application in predicting the outcome of external stimulation on EEG rhythms. Our study shows, that a neural mass model incorporating few key mechanisms is sufficient to reproduce the complex brain dynamics observed during sleep.

Timing matters: open-loop stimulation does not improve overnight consolidation of word pairs in humans

The application of auditory clicks during non‐rapid eye movement (NREM) sleep phase‐locked to the up state of the slow oscillation (closed‐loop stimulation) has previously been shown to enhance the consolidation of declarative memories. We designed and applied sequences of three clicks during deep NREM sleep to achieve a quasi‐phase‐dependent open‐loop stimulation. This stimulation was successful in eliciting slow oscillation power in the stimulation period. Although fast and slow spindle power were markedly decreased during the stimulation period, memory consolidation did not differ from control. During putative up states fast spindle power remained, however, at control levels. We conclude that concurrence of slow oscillations and fast spindles suffices to maintain memory consolidation at control levels despite an overall decreased spindle activity.

Cellular and System Biology of Memory: Timing, Molecules, and Beyond

The storage of information in the mammalian nervous systems is dependent on a delicate balance between change and stability of neuronal networks. The induction and maintenance of processes that lead to changes in synaptic strength to a multistep process which can lead to long-lasting changes, which starts and ends with a highly choreographed and perfectly timed dance of molecules in different cell types of the central nervous system. This is accompanied by synchronization of specific networks, resulting in the generation of characteristic "macroscopic" rhythmic electrical fields, whose characteristic frequencies correspond to certain activity and information-processing states of the brain. Molecular events and macroscopic fields influence each other reciprocally. We review here cellular processes of synaptic plasticity, particularly functional and structural changes, and focus on timing events that are important for the initial memory acquisition, as well as mechanisms of short- and long-term memory storage. Then, we cover the importance of epigenetic events on the long-time range. Furthermore, we consider how brain rhythms at the network level participate in processes of information storage and by what means they participating in it. Finally, we examine memory consolidation at the system level during processes of sleep.

Characteristics and plasticity of electrical synaptic transmission

Electrical synapses are an omnipresent feature of nervous systems, from the simple nerve nets of cnidarians to complex brains of mammals. Formed by gap junction channels between neurons, electrical synapses allow direct transmission of voltage signals between coupled cells. The relative simplicity of this arrangement belies the sophistication of these synapses. Coupling via electrical synapses can be regulated by a variety of mechanisms on times scales ranging from milliseconds to days, and active properties of the coupled neurons can impart emergent properties such as signal amplification, phase shifts and frequency-selective transmission. This article reviews the biophysical characteristics of electrical synapses and some of the core mechanisms that control their plasticity in the vertebrate central nervous system.

Dynamic Analysis of the Conditional Oscillator Underlying Slow Waves in Thalamocortical Neurons

During non-REM sleep the EEG shows characteristics waves that are generated by the dynamic interactions between cortical and thalamic oscillators. In thalamic neurons, low-threshold T-type Ca(2+) channels play a pivotal role in almost every type of neuronal oscillations, including slow (< 1 Hz) waves, sleep spindles and delta waves. The transient opening of T channels gives rise to the low threshold spikes (LTSs), and associated high frequency bursts of action potentials, that are characteristically present during sleep spindles and delta waves, whereas the persistent opening of a small fraction of T channels, (i.e., ITwindow) is responsible for the membrane potential bistability underlying sleep slow oscillations. Surprisingly thalamocortical (TC) neurons express a very high density of T channels that largely exceed the amount required to generate LTSs and therefore, to support certain, if not all, sleep oscillations. Here, to clarify the relationship between T current density and sleep oscillations, we systematically investigated the impact of the T conductance level on the intrinsic rhythmic activities generated in TC neurons, combining in vitro experiments and TC neuron simulation. Using bifurcation analysis, we provide insights into the dynamical processes taking place at the transition between slow and delta oscillations. Our results show that although stable delta oscillations can be evoked with minimal T conductance, the full range of slow oscillation patterns, including groups of delta oscillations separated by Up states ("grouped-delta slow waves") requires a high density of T channels. Moreover, high levels of T conductance ensure the robustness of different types of slow oscillations.

Large-scale network-level processes during entrainment

Visual rhythmic stimulation evokes a robust power increase exactly at the stimulation frequency, the so-called steady-state response (SSR). Localization of visual SSRs normally shows a very focal modulation of power in visual cortex and led to the treatment and interpretation of SSRs as a local phenomenon. Given the brain network dynamics, we hypothesized that SSRs have additional large-scale effects on the brain functional network that can be revealed by means of graph theory. We used rhythmic visual stimulation at a range of frequencies (4–30 Hz), recorded MEG and investigated source level connectivity across the whole brain. Using graph theoretical measures we observed a frequency-unspecific reduction of global density in the alpha band “disconnecting” visual cortex from the rest of the network. Also, a frequency-specific increase of connectivity between occipital cortex and precuneus was found at the stimulation frequency that exhibited the highest resonance (30 Hz). In conclusion, we showed that SSRs dynamically re-organized the brain functional network. These large-scale effects should be taken into account not only when attempting to explain the nature of SSRs, but also when used in various experimental designs.

•

Visual entrainment is considered mostly to modulate cortical power locally.

•

Instead, we hypothesized large-scale effects in the brain functional network.

•

Graph theoretical analysis combined with MEG source localization.

•

Visual entrainment indeed yielded network-level effects.

Perturbation of Brain Oscillations after Ischemic Stroke: A Potential Biomarker for Post-Stroke Function and Therapy

Brain waves resonate from the generators of electrical current and propagate across brain regions with oscillation frequencies ranging from 0.05 to 500 Hz. The commonly observed oscillatory waves recorded by an electroencephalogram (EEG) in normal adult humans can be grouped into five main categories according to the frequency and amplitude, namely δ (1-4 Hz, 20-200 μV), θ (4-8 Hz, 10 μV), α (8-12 Hz, 20-200 μV), β (12-30 Hz, 5-10 μV), and γ (30-80 Hz, low amplitude). Emerging evidence from experimental and human studies suggests that groups of function and behavior seem to be specifically associated with the presence of each oscillation band, although the complex relationship between oscillation frequency and function, as well as the interaction between brain oscillations, are far from clear. Changes of brain oscillation patterns have long been implicated in the diseases of the central nervous system including ischemic stroke, in which the reduction of cerebral blood flow as well as the progression of tissue damage have direct spatiotemporal effects on the power of several oscillatory bands and their interactions. This review summarizes the current knowledge in behavior and function associated with each brain oscillation, and also in the specific changes in brain electrical activities that correspond to the molecular events and functional alterations observed after experimental and human stroke. We provide the basis of the generations of brain oscillations and potential cellular and molecular mechanisms underlying stroke-induced perturbation. We will also discuss the implications of using brain oscillation patterns as biomarkers for the prediction of stroke outcome and therapeutic efficacy.

Analysis of the role of the low threshold currents IT and Ih in intrinsic delta oscillations of thalamocortical neurons

Thalamocortical neurons are involved in the generation and maintenance of brain rhythms associated with global functional states. The repetitive burst firing of TC neurons at delta frequencies (1–4 Hz) has been linked to the oscillations recorded during deep sleep and during episodes of absence seizures. To get insight into the biophysical properties that are the basis for intrinsic delta oscillations in these neurons, we performed a bifurcation analysis of a minimal conductance-based thalamocortical neuron model including only the I_T channel and the sodium and potassium leak channels. This analysis unveils the dynamics of repetitive burst firing of TC neurons, and describes how the interplay between the amplifying variable m_T and the recovering variable h_T of the calcium channel I_T is sufficient to generate low threshold oscillations in the delta band. We also explored the role of the hyperpolarization activated cationic current I_h in this reduced model and determine that, albeit not required, I_h amplifies and stabilizes the oscillation.