What determines the frequency of fast network oscillations with irregular neural discharges? I. Synaptic dynamics and excitation-inhibition balance

Distinct effects of heterogeneity and noise on gamma oscillation in a model of neuronal network with different reversal potential

Gamma oscillation is crucial in brain functions such as attentional selection, and is inextricably linked to both heterogeneity and noise (or so-called stochastic fluctuation) in neuronal networks. However, under coexistence of these factors, it has not been clarified how the synaptic reversal potential modulates the entraining of gamma oscillation. Here we show distinct effects of heterogeneity and noise in a population of modified theta neurons randomly coupled via GABAergic synapses. By introducing the Fokker-Planck equation and circular cumulants, we derive a set of two-cumulant macroscopic equations. In bifurcation analyses, we find a stabilizing effect of heterogeneity and a nontrivial effect of noise that results in promoting, diminishing, and shifting the oscillatory region, and is largely dependent on the reversal potential of GABAergic synapses. These findings are verified by numerical simulations of a finite-size neuronal network. Our results reveal that slight changes in reversal potential and magnitude of stochastic fluctuations can lead to immediate control of gamma oscillation, which would results in complex spatio-temporal dynamics for attentional selection and recognition.

Differential auditory brain response abnormalities in two intellectual disability conditions: SYNGAP1 mutations and Down syndrome

OBJECTIVE

Altered sensory processing is common in intellectual disability (ID). Here, we study electroencephalographic responses to auditory stimulation in human subjects presenting a rare condition (mutations in SYNGAP1) which causes ID, epilepsy and autism.

METHODS

Auditory evoked potentials, time-frequency and inter-trial coherence analyses were used to compare subjects with SYNGAP1 mutations with Down syndrome (DS) and neurotypical (NT) participants (N = 61 ranging from three to 19 years of age).

RESULTS

Altered synchronization in the brain responses to sound were found in both ID groups. The SYNGAP1 mutations group showed less phase-locking in early time windows and lower frequency bands compared to NT, and in later time windows compared to NT and DS. Time-frequency analysis showed more power in beta-gamma in the SYNGAP1 group compared to NT participants.

CONCLUSIONS

This study indicated reduced synchronization as well as more high frequencies power in SYNGAP1 mutations, while maintained synchronization was found in the DS group. These results might reflect dysfunctional sensory information processing caused by excitation/inhibition imbalance, or an imperfect compensatory mechanism in SYNGAP1 mutations individuals.

SIGNIFICANCE

Our study is the first to reveal brain response abnormalities in auditory sensory processing in SYNGAP1 mutations individuals, that are distinct from DS, another ID condition.

Lévy walk dynamics explain gamma burst patterns in primate cerebral cortex

Lévy walks describe patterns of intermittent motion with variable step sizes. In complex biological systems, Lévy walks (non-Brownian, superdiffusive random walks) are associated with behaviors such as search patterns of animals foraging for food. Here we show that Lévy walks also describe patterns of oscillatory activity in primate cerebral cortex. We used a combination of empirical observation and modeling to investigate high-frequency (gamma band) local field potential activity in visual motion-processing cortical area MT of marmoset monkeys. We found that gamma activity is organized as localized burst patterns that propagate across the cortical surface with Lévy walk dynamics. Lévy walks are fundamentally different from either global synchronization, or regular propagating waves, because they include large steps that enable activity patterns to move rapidly over cortical modules. The presence of Lévy walk dynamics therefore represents a previously undiscovered mode of brain activity, and implies a novel way for the cortex to compute. We apply a biophysically realistic circuit model to explain that the Lévy walk dynamics arise from critical-state transitions between asynchronous and localized propagating wave states, and that these dynamics yield optimal spatial sampling of the cortical sheet. We hypothesise that Lévy walk dynamics could help the cortex to efficiently process variable inputs, and to find links in patterns of activity among sparsely spiking populations of neurons.

FNS allows efficient event-driven spiking neural network simulations based on a neuron model supporting spike latency

Neural modelling tools are increasingly employed to describe, explain, and predict the human brain's behavior. Among them, spiking neural networks (SNNs) make possible the simulation of neural activity at the level of single neurons, but their use is often threatened by the resources needed in terms of processing capabilities and memory. Emerging applications where a low energy burden is required (e.g. implanted neuroprostheses) motivate the exploration of new strategies able to capture the relevant principles of neuronal dynamics in reduced and efficient models. The recent Leaky Integrate-and-Fire with Latency (LIFL) spiking neuron model shows some realistic neuronal features and efficiency at the same time, a combination of characteristics that may result appealing for SNN-based brain modelling. In this paper we introduce FNS, the first LIFL-based SNN framework, which combines spiking/synaptic modelling with the event-driven approach, allowing us to define heterogeneous neuron groups and multi-scale connectivity, with delayed connections and plastic synapses. FNS allows multi-thread, precise simulations, integrating a novel parallelization strategy and a mechanism of periodic dumping. We evaluate the performance of FNS in terms of simulation time and used memory, and compare it with those obtained with neuronal models having a similar neurocomputational profile, implemented in NEST, showing that FNS performs better in both scenarios. FNS can be advantageously used to explore the interaction within and between populations of spiking neurons, even for long time-scales and with a limited hardware configuration.

Aspartame-True or False? Narrative Review of Safety Analysis of General Use in Products

Aspartame is a sweetener introduced to replace the commonly used sucrose. It was discovered by James M. Schlatter in 1965. Being 180-200 times sweeter than sucrose, its intake was expected to reduce obesity rates in developing countries and help those struggling with diabetes. It is mainly used as a sweetener for soft drinks, confectionery, and medicines. Despite its widespread use, its safety remains controversial. This narrative review investigates the existing literature on the use of aspartame and its possible effects on the human body to refine current knowledge. Taking to account that aspartame is a widely used artificial sweetener, it seems appropriate to continue research on safety. Studies mentioned in this article have produced very interesting results overall, the current review highlights the social problem of providing visible and detailed information about the presence of aspartame in products. The studies involving the impact of aspartame on obesity, diabetes mellitus, children and fetus, autism, neurodegeneration, phenylketonuria, allergies and skin problems, its cancer properties and its genotoxicity were analyzed. Further research should be conducted to ensure clear information about the impact of aspartame on health.

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.

Regulating synchronous oscillations of cerebellar granule cells by different types of inhibition

Synchronous oscillations in neural populations are considered being controlled by inhibitory neurons. In the granular layer of the cerebellum, two major types of cells are excitatory granular cells (GCs) and inhibitory Golgi cells (GoCs). GC spatiotemporal dynamics, as the output of the granular layer, is highly regulated by GoCs. However, there are various types of inhibition implemented by GoCs. With inputs from mossy fibers, GCs and GoCs are reciprocally connected to exhibit different network motifs of synaptic connections. From the view of GCs, feedforward inhibition is expressed as the direct input from GoCs excited by mossy fibers, whereas feedback inhibition is from GoCs via GCs themselves. In addition, there are abundant gap junctions between GoCs showing another form of inhibition. It remains unclear how these diverse copies of inhibition regulate neural population oscillation changes. Leveraging a computational model of the granular layer network, we addressed this question to examine the emergence and modulation of network oscillation using different types of inhibition. We show that at the network level, feedback inhibition is crucial to generate neural oscillation. When short-term plasticity was equipped on GoC-GC synapses, oscillations were largely diminished. Robust oscillations can only appear with additional gap junctions. Moreover, there was a substantial level of cross-frequency coupling in oscillation dynamics. Such a coupling was adjusted and strengthened by GoCs through feedback inhibition. Taken together, our results suggest that the cooperation of distinct types of GoC inhibition plays an essential role in regulating synchronous oscillations of the GC population. With GCs as the sole output of the granular network, their oscillation dynamics could potentially enhance the computational capability of downstream neurons.

A Complex-Valued Oscillatory Neural Network for Storage and Retrieval of Multidimensional Aperiodic Signals

Recurrent neural networks with associative memory properties are typically based on fixed-point dynamics, which is fundamentally distinct from the oscillatory dynamics of the brain. There have been proposals for oscillatory associative memories, but here too, in the majority of cases, only binary patterns are stored as oscillatory states in the network. Oscillatory neural network models typically operate at a single/common frequency. At multiple frequencies, even a pair of oscillators with real coupling exhibits rich dynamics of Arnold tongues, not easily harnessed to achieve reliable memory storage and retrieval. Since real brain dynamics comprises of a wide range of spectral components, there is a need for oscillatory neural network models that operate at multiple frequencies. We propose an oscillatory neural network that can model multiple time series simultaneously by performing a Fourier-like decomposition of the signals. We show that these enhanced properties of a network of Hopf oscillators become possible by operating in the complex-variable domain. In this model, the single neural oscillator is modeled as a Hopf oscillator, with adaptive frequency and dynamics described over the complex domain. We propose a novel form of coupling, dubbed "power coupling," between complex Hopf oscillators. With power coupling, expressed naturally only in the complex-variable domain, it is possible to achieve stable (normalized) phase relationships in a network of multifrequency oscillators. Network connections are trained either by Hebb-like learning or by delta rule, adapted to the complex domain. The network is capable of modeling N-channel electroencephalogram time series with high accuracy and shows the potential as an effective model of large-scale brain dynamics.

Differential Contributions of Inhibitory Subnetwork to Visual Cortical Modulations Identified via Computational Model of Working Memory

Extrastriate visual neurons show no firing rate change during a working memory (WM) task in the absence of sensory input, but both αβ oscillations and spike phase locking are enhanced, as is the gain of sensory responses. This lack of change in firing rate is at odds with many models of WM, or attentional modulation of sensory networks. In this article we devised a computational model in which this constellation of results can be accounted for via selective activation of inhibitory subnetworks by a top-down working memory signal. We confirmed the model prediction of selective inhibitory activation by segmenting cells in the experimental neural data into putative excitatory and inhibitory cells. We further found that this inhibitory activation plays a dual role in influencing excitatory cells: it both modulates the inhibitory tone of the network, which underlies the enhanced sensory gain, and also produces strong spike-phase entrainment to emergent network oscillations. Using a phase oscillator model we were able to show that inhibitory tone is principally modulated through inhibitory network gain saturation, while the phase-dependent efficacy of inhibitory currents drives the phase locking modulation. The dual contributions of the inhibitory subnetwork to oscillatory and non-oscillatory modulations of neural activity provides two distinct ways for WM to recruit sensory areas, and has relevance to theories of cortical communication.

Computational exploration of dynamic mechanisms of steady state visual evoked potentials at the whole brain level

Periodic visual stimulation can induce stable steady-state visual evoked potentials (SSVEPs) distributed in multiple brain regions and has potential applications in both neural engineering and cognitive neuroscience. However, the underlying dynamic mechanisms of SSVEPs at the whole-brain level are still not completely understood. Here, we addressed this issue by simulating the rich dynamics of SSVEPs with a large-scale brain model designed with constraints of neuroimaging data acquired from the human brain. By eliciting activity of the occipital areas using an external periodic stimulus, our model was capable of replicating both the spatial distributions and response features of SSVEPs that were observed in experiments. In particular, we confirmed that alpha-band (8-12 Hz) stimulation could evoke stronger SSVEP responses; this frequency sensitivity was due to nonlinear entrainment and resonance, and could be modulated by endogenous factors in the brain. Interestingly, the stimulus-evoked brain networks also exhibited significant superiority in topological properties near this frequency-sensitivity range, and stronger SSVEP responses were demonstrated to be supported by more efficient functional connectivity at the neural activity level. These findings not only provide insights into the mechanistic understanding of SSVEPs at the whole-brain level but also indicate a bright future for large-scale brain modeling in characterizing the complicated dynamics and functions of the brain.

Role of synaptic nonlinearity in persistent firing rate shifts caused by external periodic forcing

Information storage and processing in the brain largely relies on the neural population coding principle. In this framework, information is reflected in the population firing rate that reflects asynchronous irregular spiking of its constituent neurons. Periodic modulations of neural activity can lead to neural activity oscillations. Data indicate that such oscillations are ubiquitous in brain activity and are modulated, in frequency and amplitude, in a functionally meaningful manner. The relationship between oscillations and the population rate code remains an open issue. While ample works show how changes in the mean firing rate may alter neural oscillations, the reverse connection is unclear. One notable possibility is that oscillatory activity impinging on a neural population modulates its mean firing rate, thereby impacting information processing. We suggest that such modulation requires nonlinearities and propose nonlinear excitatory coupling via slow N-methyl-D-aspartate (NMDA) receptors as the prevalent mechanism. The aim of our paper is to theoretically explore to what extent the NMDA-related mechanism could account for oscillation-induced mean firing rate changes. We consider a mean-field model of a neural circuit containing an excitatory and an inhibitory population with linear transfer functions. Along with NMDA excitation, the model included fast recurrent excitatory and inhibitory connectivity. To explicitly study the effects of impinging oscillation on the rate dynamics, we subjected the circuit to a sinusoidal input signal imitating an input from distant brain regions or from a larger network into which the circuit is embedded. Using time-scale separation and time-averaging techniques, we developed a geometric method to determine the oscillation-induced mean firing rate shifts and validated it by numeric simulations of the model. Our results indicate that a large-amplitude stable firing rate shift requires nonlinear NMDA synapses on both the excitatory and the inhibitory populations. Our results delineate specific neural synaptic properties that enable neural oscillations to act as flexible modulators of the population rate code.

Individual Alpha Frequency Determines the Impact of Bottom-Up Drive on Visual Processing

Endogenous alpha oscillations propagate from higher-order to early visual cortical regions, consistent with the observed modulation of these oscillations by top-down factors. However, bottom-up manipulations also influence alpha oscillations, and little is known about how these top-down and bottom-up processes interact to impact behavior. To address this, participants performed a detection task while viewing a stimulus flickering at multiple alpha band frequencies. Bottom-up drive at a participant's endogenous alpha frequency either impaired or enhanced perception, depending on the frequency, but not amplitude, of their endogenous alpha oscillation. Fast alpha drive impaired perceptual performance in participants with faster endogenous alpha oscillations, while participants with slower oscillations displayed enhanced performance. This interaction was reflected in slower endogenous oscillatory dynamics in participants with fast alpha oscillations and more rapid dynamics in participants with slow endogenous oscillations when receiving high-frequency bottom-up drive. This central tendency may suggest that driving visual circuits at alpha band frequencies that are away from the peak alpha frequency improves perception through dynamical interactions with the endogenous oscillation. As such, studies that causally manipulate neural oscillations via exogenous stimulation should carefully consider interacting effects of bottom-up drive and endogenous oscillations on behavior.

Coherence resonance in neuronal populations: Mean-field versus network model

The counterintuitive phenomenon of coherence resonance describes a nonmonotonic behavior of the regularity of noise-induced oscillations in the excitable regime, leading to an optimal response in terms of regularity of the excited oscillations for an intermediate noise intensity. We study this phenomenon in populations of FitzHugh-Nagumo (FHN) neurons with different coupling architectures. For networks of FHN systems in an excitable regime, coherence resonance has been previously analyzed numerically. Here we focus on an analytical approach studying the mean-field limits of the globally and locally coupled populations. The mean-field limit refers to an averaged behavior of a complex network as the number of elements goes to infinity. We apply the mean-field approach to the globally coupled FHN network. Further, we derive a mean-field limit approximating the locally coupled FHN network with low noise intensities. We study the effects of the coupling strength and noise intensity on coherence resonance for both the network and the mean-field models. We compare the results of the mean-field and network frameworks and find good agreement in the globally coupled case, where the correspondence between the two approaches is sufficiently good to capture the emergence of coherence resonance, as well as of anticoherence resonance.

Transmission delays and frequency detuning can regulate information flow between brain regions

Brain networks exhibit very variable and dynamical functional connectivity and flexible configurations of information exchange despite their overall fixed structure. Brain oscillations are hypothesized to underlie time-dependent functional connectivity by periodically changing the excitability of neural populations. In this paper, we investigate the role of the connection delay and the detuning between the natural frequencies of neural populations in the transmission of signals. Based on numerical simulations and analytical arguments, we show that the amount of information transfer between two oscillating neural populations could be determined by their connection delay and the mismatch in their oscillation frequencies. Our results highlight the role of the collective phase response curve of the oscillating neural populations for the efficacy of signal transmission and the quality of the information transfer in brain networks.

Malleability of gamma rhythms enhances population-level correlations

An important problem in systems neuroscience is to understand how information is communicated among brain regions, and it has been proposed that communication is mediated by neuronal oscillations, such as rhythms in the gamma band. We sought to investigate this idea by using a network model with two components, a source (sending) and a target (receiving) component, both built to resemble local populations in the cerebral cortex. To measure the effectiveness of communication, we used population-level correlations in spike times between the source and target. We found that after correcting for a response time that is independent of initial conditions, spike-time correlations between the source and target are significant, due in large measure to the alignment of firing events in their gamma rhythms. But, we also found that regular oscillations cannot produce the results observed in our model simulations of cortical neurons. Surprisingly, it is the irregularity of gamma rhythms, the absence of internal clocks, together with the malleability of these rhythms and their tendency to align with external pulses - features that are known to be present in gamma rhythms in the real cortex - that produced the results observed. These findings and the mechanistic explanations we offered are our primary results. Our secondary result is a mathematical relationship between correlations and the sizes of the samples used for their calculation. As improving technology enables recording simultaneously from increasing numbers of neurons, this relationship could be useful for interpreting results from experimental recordings.

Computation of the electroencephalogram (EEG) from network models of point neurons

The electroencephalogram (EEG) is a major tool for non-invasively studying brain function and dysfunction. Comparing experimentally recorded EEGs with neural network models is important to better interpret EEGs in terms of neural mechanisms. Most current neural network models use networks of simple point neurons. They capture important properties of cortical dynamics, and are numerically or analytically tractable. However, point neurons cannot generate an EEG, as EEG generation requires spatially separated transmembrane currents. Here, we explored how to compute an accurate approximation of a rodent's EEG with quantities defined in point-neuron network models. We constructed different approximations (or proxies) of the EEG signal that can be computed from networks of leaky integrate-and-fire (LIF) point neurons, such as firing rates, membrane potentials, and combinations of synaptic currents. We then evaluated how well each proxy reconstructed a ground-truth EEG obtained when the synaptic currents of the LIF model network were fed into a three-dimensional network model of multicompartmental neurons with realistic morphologies. Proxies based on linear combinations of AMPA and GABA currents performed better than proxies based on firing rates or membrane potentials. A new class of proxies, based on an optimized linear combination of time-shifted AMPA and GABA currents, provided the most accurate estimate of the EEG over a wide range of network states. The new linear proxies explained 85-95% of the variance of the ground-truth EEG for a wide range of network configurations including different cell morphologies, distributions of presynaptic inputs, positions of the recording electrode, and spatial extensions of the network. Non-linear EEG proxies using a convolutional neural network (CNN) on synaptic currents increased proxy performance by a further 2-8%. Our proxies can be used to easily calculate a biologically realistic EEG signal directly from point-neuron simulations thus facilitating a quantitative comparison between computational models and experimental EEG recordings.

Modulation of Coordinated Activity across Cortical Layers by Plasticity of Inhibitory Synapses

In the neocortex, synaptic inhibition shapes all forms of spontaneous and sensory evoked activity. Importantly, inhibitory transmission is highly plastic, but the functional role of inhibitory synaptic plasticity is unknown. In the mouse barrel cortex, activation of layer (L) 2/3 pyramidal neurons (PNs) elicits strong feedforward inhibition (FFI) onto L5 PNs. We find that FFI involving parvalbumin (PV)-expressing cells is strongly potentiated by postsynaptic PN burst firing. FFI plasticity modifies the PN excitation-to-inhibition (E/I) ratio, strongly modulates PN gain, and alters information transfer across cortical layers. Moreover, our LTPi-inducing protocol modifies firing of L5 PNs and alters the temporal association of PN spikes to γ-oscillations both in vitro and in vivo. All of these effects are captured by unbalancing the E/I ratio in a feedforward inhibition circuit model. Altogether, our results indicate that activity-dependent modulation of perisomatic inhibitory strength effectively influences the participation of single principal cortical neurons to cognition-relevant network activity.

•

Feedforward inhibition (FFI) of layer 5 pyramidal neurons (PNs) is highly plastic

•

Long-term potentiation of FFI modulates spiking activity of layer 5 PNs

•

LTPi affects information transfer across cortical layers

•

The strength of LTPi determines the phase locking of PN firing to γ-oscillations

Spectral signatures of L-DOPA-induced dyskinesia depend on L-DOPA dose and are suppressed by ketamine

L-DOPA-induced dyskinesias (LID) are debilitating motor symptoms of dopamine-replacement therapy for Parkinson's disease (PD) that emerge after years of L-DOPA treatment. While there is an abundance of research into the cellular and synaptic origins of LID, less is known about how LID impacts systems-level circuits and neural synchrony, how synchrony is affected by the dose and duration of L-DOPA exposure, or how potential novel treatments for LID, such as sub-anesthetic ketamine, alter this activity. Sub-anesthetic ketamine treatments have recently been shown to reduce LID, and ketamine is known to affect neural synchrony. To investigate these questions, we measured movement and local-field potential (LFP) activity from the motor cortex (M1) and the striatum of preclinical rodent models of PD and LID. In the first experiment, we investigated the effect of the LID priming procedures and L-DOPA dose on neural signatures of LID. Two common priming procedures were compared: a high-dose procedure that exposed unilateral 6-hydroxydopamine-lesioned rats to 12 mg/kg L-DOPA for 7 days, and a low-dose procedure that exposed rats to 7 mg/kg L-DOPA for 21 days. Consistent with reports from other groups, 12 mg/kg L-DOPA triggered LID and 80-Hz oscillations; however, these 80-Hz oscillations were not observed after 7 mg/kg administration despite clear evidence of LID, indicating that 80-Hz oscillations are not an exclusive signature of LID. We also found that weeks-long low-dose priming resulted in the emergence of non-oscillatory broadband gamma activity (> 30 Hz) in the striatum and theta-to-high-gamma cross-frequency coupling (CFC) in M1. In a second set of experiments, we investigated how ketamine exposure affects spectral signatures of low-dose L-DOPA priming. During each neural recording session, ketamine was delivered through 5 injections (20 mg/kg, i.p.) administered every 2 h. We found that ketamine exposure suppressed striatal broadband gamma associated with LID but enhanced M1 broadband activity. We also found that M1 theta-to-high-gamma CFC associated with the LID on-state was suppressed by ketamine. These results suggest that ketamine's therapeutic effects are region specific. Our findings also have clinical implications, as we are the first to report novel oscillatory signatures of the common low-dose LID priming procedure that more closely models dopamine replacement therapy in individuals with PD. We also identify neural correlates of the anti-dyskinetic activity of sub-anesthetic ketamine treatment.

Neuronal heterogeneity modulates phase synchronization between unidirectionally coupled populations with excitation-inhibition balance

Several experiments and models have highlighted the importance of neuronal heterogeneity in brain dynamics and function. However, how such a cell-to-cell diversity can affect cortical computation, synchronization, and neuronal communication is still under debate. Previous studies have focused on the effect of neuronal heterogeneity in one neuronal population. Here we are specifically interested in the effect of neuronal variability on the phase relations between two populations, which can be related to different cortical communication hypotheses. It has been recently shown that two spiking neuron populations unidirectionally connected in a sender-receiver configuration can exhibit anticipated synchronization (AS), which is characterized by a negative phase lag. This phenomenon has been reported in electrophysiological data of nonhuman primates and human EEG during a visual discrimination cognitive task. In experiments, the unidirectional coupling could be accessed by Granger causality and can be accompanied by either positive or negative phase difference between cortical areas. Here we propose a model of two coupled populations in which the neuronal heterogeneity can determine the dynamical relation between the sender and the receiver and can reproduce phase relations reported in experiments. Depending on the distribution of parameters characterizing the neuronal firing patterns, the system can exhibit both AS and the usual delayed synchronization regime (DS, with positive phase) as well as a zero-lag synchronization regime and phase bistability between AS and DS. Furthermore, we show that our network can present diversity in their phase relations maintaining the excitation-inhibition balance.

Selective postnatal excitation of neocortical pyramidal neurons results in distinctive behavioral and circuit deficits in adulthood

In genetic and pharmacological models of neurodevelopmental disorders, and human data, neural activity is altered within the developing neocortical network. This commonality begs the question of whether early enhancement in excitation might be a common driver, across etiologies, of characteristic behaviors. We tested this concept by chemogenetically driving cortical pyramidal neurons during postnatal days 4–14. Hyperexcitation of Emx1-, but not dopamine transporter-, parvalbumin-, or Dlx5/6-expressing neurons, led to decreased social interaction and increased grooming activity in adult animals. In vivo optogenetic interrogation in adults revealed decreased baseline but increased stimulus-evoked firing rates of pyramidal neurons and impaired recruitment of inhibitory neurons. Slice recordings in adults from prefrontal cortex layer 5 pyramidal neurons revealed decreased intrinsic excitability and increased synaptic E/I ratio. Together these results support the prediction that enhanced pyramidal firing during development, in otherwise normal cortex, can selectively drive altered adult circuit function and maladaptive changes in behavior.

•

BL-OG allows chemogenetic activation and optogenetic interrogation in the same animal

•

Developmental hyperexcitation in normal mice leads to neurodevelopmental disorders

•

In these mice adult neurons show reduced baseline activity and increased excitability

•

Reduced activity-triggered coherence and altered oscillations in cortex and striatum

Reduced coupling of somatosensory gating and gamma oscillation in panic disorder

Previous studies have reported that patients with panic disorder (PD) exhibited an aberrant level of GABA concentration, an inhibitory neurotransmitter in the human brain. However, it remains substantially unclear whether the inhibitory function regarding the neurophysiological characteristics is altered in this disease. Sensory gating (SG) is considered as an automatic inhibitory function in the sensory cortex. In addition, brain's gamma oscillation within the sensory cortex is another index to reflect inhibitory function. Here we aimed to investigate whether the patients with PD showed altered inhibitory function in the somatosensory system, including the primary (SI) and secondary (SII) somatosensory cortices. A total of 20 healthy controls and 21 patients with PD underwent magnetoencephalographic recordings. Paired-pulse and single-pulse paradigms were used to study SG and gamma oscillations, respectively. There were no significant between-group differences in the SG function in the SI and SII. However, patients with PD demonstrated a reduced gamma power in the SI. Among the healthy individuals, strong associations between SG ratios and gamma frequency values were observed in the SI. However, such a functional relationship disappeared among the patients with PD. We suggested the reduced coupling of SG and gamma oscillation as one of the neural signatures in PD.

Neuronal oscillations and the mouse prefrontal cortex

The mouse prefrontal cortex (PFC) encompasses a collection of agranual brain regions in the rostral neocortex and is considered to be critically involved in the neuronal computations underlying intentional behaviors. Flexible behavioral responses demand coordinated integration of sensory inputs with state, goal and memory information in brain-wide neuronal networks. Neuronal oscillations are proposed to provide a temporal scaffold for coordination of neuronal network activity and routing of information. In the present book chapter, we review findings on the role neuronal oscillations in prefrontal functioning, with a specific focus on research in mice. We discuss discoveries pertaining to local prefrontal processing, as well to interactions with other brain regions. We also discuss how the recent discovery of brain-wide respiration-entrained rhythms (RR) warrant re-evaluation of certain findings on slow oscillations (<10Hz) in prefrontal functioning.

Gamma Oscillations Facilitate Effective Learning in Excitatory-Inhibitory Balanced Neural Circuits

Gamma oscillation in neural circuits is believed to associate with effective learning in the brain, while the underlying mechanism is unclear. This paper aims to study how spike-timing-dependent plasticity (STDP), a typical mechanism of learning, with its interaction with gamma oscillation in neural circuits, shapes the network dynamics properties and the network structure formation. We study an excitatory-inhibitory (E-I) integrate-and-fire neuronal network with triplet STDP, heterosynaptic plasticity, and a transmitter-induced plasticity. Our results show that the performance of plasticity is diverse in different synchronization levels. We find that gamma oscillation is beneficial to synaptic potentiation among stimulated neurons by forming a special network structure where the sum of excitatory input synaptic strength is correlated with the sum of inhibitory input synaptic strength. The circuit can maintain E-I balanced input on average, whereas the balance is temporal broken during the learning-induced oscillations. Our study reveals a potential mechanism about the benefits of gamma oscillation on learning in biological neural circuits.

Empirically constrained network models for contrast-dependent modulation of gamma rhythm in V1

Gamma oscillations are thought to play a key role in neuronal network function and neuronal communication, yet the underlying generating mechanisms have not been fully elucidated to date. At least partly, this may be due to the fact that even in simple network models of interconnected inhibitory (I) and excitatory (E) neurons, many parameters remain unknown and are set based on practical considerations or by convention. Here, we mitigate this problem by requiring PING (Pyramidal Interneuron Network Gamma) models to simultaneously satisfy a broad set of criteria for realistic behaviour based on empirical data spanning both the single unit (spikes) and local population (LFP) levels while unknown parameters are varied. By doing so, we were able to constrain the parameter ranges and select empirically valid models. The derived model constraints implied weak rather than strong PING as the generating mechanism for gamma, connectivity between E and I neurons within specific bounds, and variations of the external input to E but not I neurons. Constrained models showed valid behaviours, including gamma frequency increases with contrast and power saturation or decay at high contrasts. Using an empirically-validated model we studied the route to gamma instability at high contrasts. This involved increased heterogeneity of E neurons with increasing input triggering a breakdown of I neuron pacemaker function. Further, we illustrate the model's capacity to resolve disputes in the literature concerning gamma oscillation properties and GABA conductance proxies. We propose that the models derived in our study will be useful for other modelling studies, and that our approach to the empirical constraining of PING models can be expanded when richer empirical datasets become available. As local gamma networks are the building blocks of larger networks that aim to understand complex cognition through their interactions, there is considerable value in improving our models of these building blocks.

Transcranial direct current stimulation modulates brain functional connectivity in autism

Autism spectrum disorder (ASD) is characterized by deficits in social interactions, impairments in language and communication, and highly restricted behavioral interests. Transcranial direct current stimulation (tDCS) is a widely used form of noninvasive stimulation and may have therapeutic potential for ASD. So far, despite the widespread use of this technique in the neuroscience field, its effects on network-level neural activity and the underlying mechanisms of any effects are still unclear. In the present study, we used electroencephalography (EEG) to investigate tDCS induced brain network changes in children with ASD before and after active and sham stimulation. We recorded 5 min of resting state EEG before and after a single session of tDCS (of approximately 20 min) over dorsolateral prefrontal cortex (DLPFC). Two network-based methods were applied to investigate tDCS modulation on brain networks: 1) temporal network dynamics were analyzed by comparing "flexibility" changes before vs after stimulation, and 2) frequency specific network changes were identified using non-negative matrix factorization (NMF). We found 1) an increase in network flexibility following tDCS (rapid network configuration of dynamic network communities), 2) specific increase in interhemispheric connectivity within the alpha frequency band following tDCS. Together, these results demonstrate that tDCS could help modify both local and global brain network dynamics, and highlight stimulation-induced differences in the manifestation of network reconfiguration. Meanwhile, frequency-specific subnetworks, as a way to index local and global information processing, highlight the core modulatory effects of tDCS on the modular architecture of the functional connectivity patterns within higher frequency bands.

On the Validity of Neural Mass Models

Modeling the dynamics of neural masses is a common approach in the study of neural populations. Various models have been proven useful to describe a plenitude of empirical observations including self-sustained local oscillations and patterns of distant synchronization. We discuss the extent to which mass models really resemble the mean dynamics of a neural population. In particular, we question the validity of neural mass models if the population under study comprises a mixture of excitatory and inhibitory neurons that are densely (inter-)connected. Starting from a network of noisy leaky integrate-and-fire neurons, we formulated two different population dynamics that both fall into the category of seminal Freeman neural mass models. The derivations contained several mean-field assumptions and time scale separation(s) between membrane and synapse dynamics. Our comparison of these neural mass models with the averaged dynamics of the population reveals bounds in the fraction of excitatory/inhibitory neuron as well as overall network degree for a mass model to provide adequate estimates. For substantial parameter ranges, our models fail to mimic the neural network's dynamics proper, be that in de-synchronized or in (high-frequency) synchronized states. Only around the onset of low-frequency synchronization our models provide proper estimates of the mean potential dynamics. While this shows their potential for, e.g., studying resting state dynamics obtained by encephalography with focus on the transition region, we must accept that predicting the more general dynamic outcome of a neural network via its mass dynamics requires great care.

A Connectome-Based, Corticothalamic Model of State- and Stimulation-Dependent Modulation of Rhythmic Neural Activity and Connectivity

Rhythmic activity in the brain fluctuates with behaviour and cognitive state, through a combination of coexisting and interacting frequencies. At large spatial scales such as those studied in human M/EEG, measured oscillatory dynamics are believed to arise primarily from a combination of cortical (intracolumnar) and corticothalamic rhythmogenic mechanisms. Whilst considerable progress has been made in characterizing these two types of neural circuit separately, relatively little work has been done that attempts to unify them into a single consistent picture. This is the aim of the present paper. We present and examine a whole-brain, connectome-based neural mass model with detailed long-range cortico-cortical connectivity and strong, recurrent corticothalamic circuitry. This system reproduces a variety of known features of human M/EEG recordings, including spectral peaks at canonical frequencies, and functional connectivity structure that is shaped by the underlying anatomical connectivity. Importantly, our model is able to capture state- (e.g., idling/active) dependent fluctuations in oscillatory activity and the coexistence of multiple oscillatory phenomena, as well as frequency-specific modulation of functional connectivity. We find that increasing the level of sensory drive to the thalamus triggers a suppression of the dominant low frequency rhythms generated by corticothalamic loops, and subsequent disinhibition of higher frequency endogenous rhythmic behaviour of intracolumnar microcircuits. These combine to yield simultaneous decreases in lower frequency and increases in higher frequency components of the M/EEG power spectrum during states of high sensory or cognitive drive. Building on this, we also explored the effect of pulsatile brain stimulation on ongoing oscillatory activity, and evaluated the impact of coexistent frequencies and state-dependent fluctuations on the response of cortical networks. Our results provide new insight into the role played by cortical and corticothalamic circuits in shaping intrinsic brain rhythms, and suggest new directions for brain stimulation therapies aimed at state-and frequency-specific control of oscillatory brain activity.

Effect of diverse recoding of granule cells on optokinetic response in a cerebellar ring network with synaptic plasticity

We consider a cerebellar ring network for the optokinetic response (OKR), and investigate the effect of diverse recoding of granule (GR) cells on OKR by varying the connection probability p_c from Golgi to GR cells. For an optimal value of p_c^∗(=0.06), individual GR cells exhibit diverse spiking patterns which are in-phase, anti-phase, or complex out-of-phase with respect to their population-averaged firing activity. Then, these diversely-recoded signals via parallel fibers (PFs) from GR cells are effectively depressed by the error-teaching signals via climbing fibers from the inferior olive which are also in-phase ones. Synaptic weights at in-phase PF-Purkinje cell (PC) synapses of active GR cells are strongly depressed via strong long-term depression (LTD), while those at anti-phase and complex out-of-phase PF-PC synapses are weakly depressed through weak LTD. This kind of "effective" depression (i.e., strong/weak LTD) at the PF-PC synapses causes a big modulation in firings of PCs, which then exert effective inhibitory coordination on the vestibular nucleus (VN) neuron (which evokes OKR). For the firing of the VN neuron, the learning gain degree L_g, corresponding to the modulation gain ratio, increases with increasing the learning cycle, and it saturates at about the 300th cycle. By varying p_c from p_c^∗, we find that a plot of saturated learning gain degree L_g^∗ versus p_c forms a bell-shaped curve with a peak at p_c^∗ (where the diversity degree in spiking patterns of GR cells is also maximum). Consequently, the more diverse in recoding of GR cells, the more effective in motor learning for the OKR adaptation.

Hopf Bifurcation in Mean Field Explains Critical Avalanches in Excitation-Inhibition Balanced Neuronal Networks: A Mechanism for Multiscale Variability

Cortical neural circuits display highly irregular spiking in individual neurons but variably sized collective firing, oscillations and critical avalanches at the population level, all of which have functional importance for information processing. Theoretically, the balance of excitation and inhibition inputs is thought to account for spiking irregularity and critical avalanches may originate from an underlying phase transition. However, the theoretical reconciliation of these multilevel dynamic aspects in neural circuits remains an open question. Herein, we study excitation-inhibition (E-I) balanced neuronal network with biologically realistic synaptic kinetics. It can maintain irregular spiking dynamics with different levels of synchrony and critical avalanches emerge near the synchronous transition point. We propose a novel semi-analytical mean-field theory to derive the field equations governing the network macroscopic dynamics. It reveals that the E-I balanced state of the network manifesting irregular individual spiking is characterized by a macroscopic stable state, which can be either a fixed point or a periodic motion and the transition is predicted by a Hopf bifurcation in the macroscopic field. Furthermore, by analyzing public data, we find the coexistence of irregular spiking and critical avalanches in the spontaneous spiking activities of mouse cortical slice in vitro, indicating the universality of the observed phenomena. Our theory unveils the mechanism that permits complex neural activities in different spatiotemporal scales to coexist and elucidates a possible origin of the criticality of neural systems. It also provides a novel tool for analyzing the macroscopic dynamics of E-I balanced networks and its relationship to the microscopic counterparts, which can be useful for large-scale modeling and computation of cortical dynamics.

Differences in visually induced MEG oscillations reflect differences in deep cortical layer activity

Neural activity is organized at multiple scales, ranging from the cellular to the whole brain level. Connecting neural dynamics at different scales is important for understanding brain pathology. Neurological diseases and disorders arise from interactions between factors that are expressed in multiple scales. Here, we suggest a new way to link microscopic and macroscopic dynamics through combinations of computational models. This exploits results from statistical decision theory and Bayesian inference. To validate our approach, we used two independent MEG datasets. In both, we found that variability in visually induced oscillations recorded from different people in simple visual perception tasks resulted from differences in the level of inhibition specific to deep cortical layers. This suggests differences in feedback to sensory areas and each subject's hypotheses about sensations due to differences in their prior experience. Our approach provides a new link between non-invasive brain imaging data, laminar dynamics and top-down control.

Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture

Complex cognitive functions such as working memory and decision-making require information maintenance over seconds to years, from transient sensory stimuli to long-term contextual cues. While theoretical accounts predict the emergence of a corresponding hierarchy of neuronal timescales, direct electrophysiological evidence across the human cortex is lacking. Here, we infer neuronal timescales from invasive intracranial recordings. Timescales increase along the principal sensorimotor-to-association axis across the entire human cortex, and scale with single-unit timescales within macaques. Cortex-wide transcriptomic analysis shows direct alignment between timescales and expression of excitation- and inhibition-related genes, as well as genes specific to voltage-gated transmembrane ion transporters. Finally, neuronal timescales are functionally dynamic: prefrontal cortex timescales expand during working memory maintenance and predict individual performance, while cortex-wide timescales compress with aging. Thus, neuronal timescales follow cytoarchitectonic gradients across the human cortex and are relevant for cognition in both short and long terms, bridging microcircuit physiology with macroscale dynamics and behavior.

Characterizing Inscapes and resting-state in MEG: Effects in typical and atypical development

Examining the brain at rest is a powerful approach used to understand the intrinsic properties of typical and disordered human brain function, yet task-free paradigms are associated with greater head motion, particularly in young and/or clinical populations such as autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD). Inscapes, a non-social and non-verbal movie paradigm, has been introduced to increase attention, thus mitigating head motion, while reducing the task-induced activations found during typical movie watching. Inscapes has not yet been validated for use in magnetoencephalography (MEG), and it has yet to be shown whether its effects are stable in clinical populations. Across typically developing (N = 32) children and adolescents and those with ASD (N = 46) and ADHD (N = 42), we demonstrate that head motion is reduced during Inscapes. Due to the task state evoked by movie paradigms, we also expectedly observed concomitant modulations in local neural activity (oscillatory power) and functional connectivity (phase and envelope coupling) in intrinsic resting-state networks and across the frequency spectra compared to a fixation cross resting-state. Increases in local activity were accompanied by decreases in low-frequency connectivity within and between resting-state networks, primarily the visual network, suggesting that task-state evoked by Inscapes moderates ongoing and spontaneous cortical inhibition that forms the idling intrinsic networks found during a fixation cross resting-state. Importantly, these effects were similar in ASD and ADHD, making Inscapes a well-suited advancement for investigations of resting brain function in young and clinical populations.

Generative modelling of the thalamo-cortical circuit mechanisms underlying the neurophysiological effects of ketamine

Cortical recordings of task-induced oscillations following subanaesthetic ketamine administration demonstrate alterations in amplitude, including increases at high-frequencies (gamma) and reductions at low frequencies (theta, alpha). To investigate the population-level interactions underlying these changes, we implemented a thalamo-cortical model (TCM) capable of recapitulating broadband spectral responses. Compared with an existing cortex-only 4-population model, Bayesian Model Selection preferred the TCM. The model was able to accurately and significantly recapitulate ketamine-induced reductions in alpha amplitude and increases in gamma amplitude. Parameter analysis revealed no change in receptor time-constants but significant increases in select synaptic connectivity with ketamine. Significantly increased connections included both AMPA and NMDA mediated connections from layer 2/3 superficial pyramidal cells to inhibitory interneurons and both GABAA and NMDA mediated within-population gain control of layer 5 pyramidal cells. These results support the use of extended generative models for explaining oscillatory data and provide in silico support for ketamine's ability to alter local coupling mediated by NMDA, AMPA and GABA-A.

Synchronization, Stochasticity, and Phase Waves in Neuronal Networks With Spatially-Structured Connectivity

Oscillations in the beta/low gamma range (10–45 Hz) are recorded in diverse neural structures. They have successfully been modeled as sparsely synchronized oscillations arising from reciprocal interactions between randomly connected excitatory (E) pyramidal cells and local interneurons (I). The synchronization of spatially distant oscillatory spiking E–I modules has been well-studied in the rate model framework but less so for modules of spiking neurons. Here, we first show that previously proposed modifications of rate models provide a quantitative description of spiking E–I modules of Exponential Integrate-and-Fire (EIF) neurons. This allows us to analyze the dynamical regimes of sparsely synchronized oscillatory E–I modules connected by long-range excitatory interactions, for two modules, as well as for a chain of such modules. For modules with a large number of neurons (> 10⁵), we obtain results similar to previously obtained ones based on the classic deterministic Wilson-Cowan rate model, with the added bonus that the results quantitatively describe simulations of spiking EIF neurons. However, for modules with a moderate (~ 10⁴) number of neurons, stochastic variations in the spike emission of neurons are important and need to be taken into account. On the one hand, they modify the oscillations in a way that tends to promote synchronization between different modules. On the other hand, independent fluctuations on different modules tend to disrupt synchronization. The correlations between distant oscillatory modules can be described by stochastic equations for the oscillator phases that have been intensely studied in other contexts. On shorter distances, we develop a description that also takes into account amplitude modes and that quantitatively accounts for our simulation data. Stochastic dephasing of neighboring modules produces transient phase gradients and the transient appearance of phase waves. We propose that these stochastically-induced phase waves provide an explanative framework for the observations of traveling waves in the cortex during beta oscillations.

Neurostimulation stabilizes spiking neural networks by disrupting seizure-like oscillatory transitions

An improved understanding of the mechanisms underlying neuromodulatory approaches to mitigate seizure onset is needed to identify clinical targets for the treatment of epilepsy. Using a Wilson–Cowan-motivated network of inhibitory and excitatory populations, we examined the role played by intrinsic and extrinsic stimuli on the network’s predisposition to sudden transitions into oscillatory dynamics, similar to the transition to the seizure state. Our joint computational and mathematical analyses revealed that such stimuli, be they noisy or periodic in nature, exert a stabilizing influence on network responses, disrupting the development of such oscillations. Based on a combination of numerical simulations and mean-field analyses, our results suggest that high variance and/or high frequency stimulation waveforms can prevent multi-stability, a mathematical harbinger of sudden changes in network dynamics. By tuning the neurons’ responses to input, stimuli stabilize network dynamics away from these transitions. Furthermore, our research shows that such stabilization of neural activity occurs through a selective recruitment of inhibitory cells, providing a theoretical undergird for the known key role these cells play in both the healthy and diseased brain. Taken together, these findings provide new vistas on neuromodulatory approaches to stabilize neural microcircuit activity.

Dynamics of multiple interacting excitatory and inhibitory populations with delays

A network consisting of excitatory and inhibitory (EI) neurons is a canonical model for understanding local cortical network activity. In this study, we extended the local circuit model and investigated how its dynamical landscape can be enriched when it interacts with another excitatory (E) population with long transmission delays. Through analysis of a rate model and numerical simulations of a corresponding network of spiking neurons, we studied the transition from stationary to oscillatory states by analyzing the Hopf bifurcation structure in terms of two network parameters: (1) transmission delay between the EI subnetwork and the E population and (2) inhibitory couplings that induced oscillatory activity in the EI subnetwork. We found that the critical coupling strength can strongly modulate as a function of transmission delay, and consequently the stationary state can be interwoven intricately with the oscillatory state. Such a dynamical landscape gave rise to an isolated stationary state surrounded by multiple oscillatory states that generated different frequency modes, and cross-frequency coupling developed naturally at the bifurcation points. We identified the network motifs with short- and long-range inhibitory connections that underlie the emergence of oscillatory states with multiple frequencies. Thus, we provided a mechanistic explanation of how the transmission delay to and from the additional E population altered the dynamical landscape. In summary, our results demonstrated the potential role of long-range connections in shaping the network activity of local cortical circuits.

Generation of Sharp Wave-Ripple Events by Disinhibition

Sharp wave-ripple complexes (SWRs) are hippocampal network phenomena involved in memory consolidation. To date, the mechanisms underlying their occurrence remain obscure. Here, we show how the interactions between pyramidal cells, parvalbumin-positive (PV⁺) basket cells, and an unidentified class of anti-SWR interneurons can contribute to the initiation and termination of SWRs. Using a biophysically constrained model of a network of spiking neurons and a rate-model approximation, we demonstrate that SWRs emerge as a result of the competition between two interneuron populations and the resulting disinhibition of pyramidal cells. Our models explain how the activation of pyramidal cells or PV⁺ cells can trigger SWRs, as shown in vitro, and suggests that PV⁺ cell-mediated short-term synaptic depression influences the experimentally reported dynamics of SWR events. Furthermore, we predict that the silencing of anti-SWR interneurons can trigger SWRs. These results broaden our understanding of the microcircuits supporting the generation of memory-related network dynamics.

SIGNIFICANCE STATEMENT The hippocampus is a part of the mammalian brain that is crucial for episodic memories. During periods of sleep and inactive waking, the extracellular activity of the hippocampus is dominated by sharp wave-ripple events (SWRs), which have been shown to be important for memory consolidation. The mechanisms regulating the emergence of these events are still unclear. We developed a computational model to study the emergence of SWRs and to explain the roles of different cell types in regulating them. The model accounts for several previously unexplained features of SWRs and thus advances the understanding of memory-related dynamics.

Firing rate-dependent phase responses of Purkinje cells support transient oscillations

Both spike rate and timing can transmit information in the brain. Phase response curves (PRCs) quantify how a neuron transforms input to output by spike timing. PRCs exhibit strong firing-rate adaptation, but its mechanism and relevance for network output are poorly understood. Using our Purkinje cell (PC) model, we demonstrate that the rate adaptation is caused by rate-dependent subthreshold membrane potentials efficiently regulating the activation of Na+ channels. Then, we use a realistic PC network model to examine how rate-dependent responses synchronize spikes in the scenario of reciprocal inhibition-caused high-frequency oscillations. The changes in PRC cause oscillations and spike correlations only at high firing rates. The causal role of the PRC is confirmed using a simpler coupled oscillator network model. This mechanism enables transient oscillations between fast-spiking neurons that thereby form PC assemblies. Our work demonstrates that rate adaptation of PRCs can spatio-temporally organize the PC input to cerebellar nuclei.

A recurrent circuit implements normalization, simulating the dynamics of V1 activity

The normalization model has been applied to explain neural activity in diverse neural systems including primary visual cortex (V1). The model's defining characteristic is that the response of each neuron is divided by a factor that includes a weighted sum of activity of a pool of neurons. Despite the success of the normalization model, there are three unresolved issues. 1) Experimental evidence supports the hypothesis that normalization in V1 operates via recurrent amplification, i.e., amplifying weak inputs more than strong inputs. It is unknown how normalization arises from recurrent amplification. 2) Experiments have demonstrated that normalization is weighted such that each weight specifies how one neuron contributes to another's normalization pool. It is unknown how weighted normalization arises from a recurrent circuit. 3) Neural activity in V1 exhibits complex dynamics, including gamma oscillations, linked to normalization. It is unknown how these dynamics emerge from normalization. Here, a family of recurrent circuit models is reported, each of which comprises coupled neural integrators to implement normalization via recurrent amplification with arbitrary normalization weights, some of which can recapitulate key experimental observations of the dynamics of neural activity in V1.

Reconfigurations within resonating communities of brain regions following TMS reveal different scales of processing

An overarching goal of neuroscience research is to understand how heterogeneous neuronal ensembles cohere into networks of coordinated activity to support cognition. To investigate how local activity harmonizes with global signals, we measured electroencephalography (EEG) while single pulses of transcranial magnetic stimulation (TMS) perturbed occipital and parietal cortices. We estimate the rapid network reconfigurations in dynamic network communities within specific frequency bands of the EEG, and characterize two distinct features of network reconfiguration, flexibility and allegiance, among spatially distributed neural sources following TMS. Using distance from the stimulation site to infer local and global effects, we find that alpha activity (8–12 Hz) reflects concurrent local and global effects on network dynamics. Pairwise allegiance of brain regions to communities on average increased near the stimulation site, whereas TMS-induced changes to flexibility were generally invariant to distance and stimulation site. In contrast, communities within the beta (13–20 Hz) band demonstrated a high level of spatial specificity, particularly within a cluster comprising paracentral areas. Together, these results suggest that focal magnetic neurostimulation to distinct cortical sites can help identify both local and global effects on brain network dynamics, and highlight fundamental differences in the manifestation of network reconfigurations within alpha and beta frequency bands.

TMS may be used to probe the causal link between local regional activity and global brain dynamics. Using simultaneous TMS-EEG and dynamic community detection, we introduce what we call “resonating communities” or frequency band-specific clusters in the brain, as a way to index local and global processing. These resonating communities within the alpha and beta bands display both global (or integrating) behavior and local specificity, highlighting fundamental differences in the manifestation of network reconfigurations.

Portrait of visual cortical circuits for generating neural oscillation dynamics

The mouse primary visual cortex (V1) has emerged as a classical system to study neural circuit mechanisms underlying visual function and plasticity. A variety of efferent-afferent neuronal connections exists within the V1 and between the V1 and higher visual cortical areas or thalamic nuclei, indicating that the V1 system is more than a mere receiver in information processing. Sensory representations in the V1 are dynamically correlated with neural activity oscillations that are distributed across different cortical layers in an input-dependent manner. Circuits consisting of excitatory pyramidal cells (PCs) and inhibitory interneurons (INs) are the basis for generating neural oscillations. In general, INs are clustered with their adjacent PCs to form specific microcircuits that gate or filter the neural information. The interaction between these two cell populations has to be coordinated within a local circuit in order to preserve neural coding schemes and maintain excitation-inhibition (E-I) balance. Phasic alternations of the E-I balance can dynamically regulate temporal rhythms of neural oscillation. Accumulating experimental evidence suggests that the two major sub-types of INs, parvalbumin-expressing (PV⁺) cells and somatostatin-expressing (SOM⁺) INs, are active in controlling slow and fast oscillations, respectively, in the mouse V1. The review summarizes recent experimental findings on elucidating cellular or circuitry mechanisms for the generation of neural oscillations with distinct rhythms in either developing or matured mouse V1, mainly focusing on visual relaying circuits and distinct local inhibitory circuits.

Facilitating the propagation of spiking activity in feedforward networks by including feedback

Transient oscillations in network activity upon sensory stimulation have been reported in different sensory areas of the brain. These evoked oscillations are the generic response of networks of excitatory and inhibitory neurons (EI-networks) to a transient external input. Recently, it has been shown that this resonance property of EI-networks can be exploited for communication in modular neuronal networks by enabling the transmission of sequences of synchronous spike volleys (’pulse packets’), despite the sparse and weak connectivity between the modules. The condition for successful transmission is that the pulse packet (PP) intervals match the period of the modules’ resonance frequency. Hence, the mechanism was termed communication through resonance (CTR). This mechanism has three severe constraints, though. First, it needs periodic trains of PPs, whereas single PPs fail to propagate. Second, the inter-PP interval needs to match the network resonance. Third, transmission is very slow, because in each module, the network resonance needs to build up over multiple oscillation cycles. Here, we show that, by adding appropriate feedback connections to the network, the CTR mechanism can be improved and the aforementioned constraints relaxed. Specifically, we show that adding feedback connections between two upstream modules, called the resonance pair, in an otherwise feedforward modular network can support successful propagation of a single PP throughout the entire network. The key condition for successful transmission is that the sum of the forward and backward delays in the resonance pair matches the resonance frequency of the network modules. The transmission is much faster, by more than a factor of two, than in the original CTR mechanism. Moreover, it distinctly lowers the threshold for successful communication by synchronous spiking in modular networks of weakly coupled networks. Thus, our results suggest a new functional role of bidirectional connectivity for the communication in cortical area networks.

The cortex is organized as a modular system, with the modules (cortical areas) communicating via weak long-range connections. It has been suggested that the intrinsic resonance properties of population activities in these areas might contribute to enabling successful communication. A module’s intrinsic resonance appears in the damped oscillatory response to an incoming spike volley, enabling successful communication during the peaks of the oscillation. Such communication can be exploited in feedforward networks, provided the participating networks have similar resonance frequencies. This, however, is not necessarily true for cortical networks. Moreover, the communication is slow, as it takes several oscillation cycles to build up the response in the downstream network. Also, only periodic trains of spikes volleys (and not single volleys) with matching intervals can propagate. Here, we present a novel mechanism that alleviates these shortcomings and enables propagation of synchronous spiking across weakly connected networks with not necessarily identical resonance frequencies. In this framework, an individual spike volley can propagate by local amplification through reverberation in a loop between two successive networks, connected by feedforward and feedback connections: the resonance pair. This overcomes the need for activity build-up in downstream networks, causing the volley to propagate distinctly faster and more reliably.

Intrinsic excitation-inhibition imbalance affects medial prefrontal cortex differently in autistic men versus women

Excitation-inhibition (E:I) imbalance is theorized as an important pathophysiological mechanism in autism. Autism affects males more frequently than females and sex-related mechanisms (e.g., X-linked genes, androgen hormones) can influence E:I balance. This suggests that E:I imbalance may affect autism differently in males versus females. With a combination of in-silico modeling and in-vivo chemogenetic manipulations in mice, we first show that a time-series metric estimated from fMRI BOLD signal, the Hurst exponent (H), can be an index for underlying change in the synaptic E:I ratio. In autism we find that H is reduced, indicating increased excitation, in the medial prefrontal cortex (MPFC) of autistic males but not females. Increasingly intact MPFC H is also associated with heightened ability to behaviorally camouflage social-communicative difficulties, but only in autistic females. This work suggests that H in BOLD can index synaptic E:I ratio and that E:I imbalance affects autistic males and females differently.

Dependency analysis of frequency and strength of gamma oscillations on input difference between excitatory and inhibitory neurons

It has been found that gamma oscillations and the oscillation frequencies are regulated by the properties of external stimuli in many biology experimental researches. To unveil the underlying mechanism, firstly, we reproduced the experimental observations in an excitatory/inhibitory (E/I) neuronal network that the oscillation became stronger and moved to a higher frequency band (gamma band) with the increasing of the input difference between E/I neurons. Secondly, we found that gamma oscillation was induced by the unbalance between positive and negative synaptic currents, which was caused by the input difference between E/I neurons. When this input difference became greater, there would be a stronger gamma oscillation (i.e., a higher peak power in the power spectrum of the population activity of neurons). Further investigation revealed that the frequency dependency of gamma oscillation on the input difference between E/I neurons could be explained by the well-known mechanisms of inter-neuron-gamma (ING) and pyramidal-interneuron-gamma (PING). Finally, we derived mathematical analysis to verify the mechanism of frequency regulations and the results were consistent with the simulation results. The results of this paper provide a possible mechanism for the external stimuli-regulated gamma oscillations.

Parallel Neural Multiprocessing with Gamma Frequency Latencies

The Poisson variability in cortical neural responses has been typically modeled using spike averaging techniques, such as trial averaging and rate coding, since such methods can produce reliable correlates of behavior. However, mechanisms that rely on counting spikes could be slow and inefficient and thus might not be useful in the brain for computations at timescales in the 10 millisecond range. This issue has motivated a search for alternative spike codes that take advantage of spike timing and has resulted in many studies that use synchronized neural networks for communication. Here we focus on recent studies that suggest that the gamma frequency may provide a reference that allows local spike phase representations that could result in much faster information transmission. We have developed a unified model (gamma spike multiplexing) that takes advantage of a single cycle of a cell's somatic gamma frequency to modulate the generation of its action potentials. An important consequence of this coding mechanism is that it allows multiple independent neural processes to run in parallel, thereby greatly increasing the processing capability of the cortex. System-level simulations and preliminary analysis of mouse cortical cell data are presented as support for the proposed theoretical model.

A kernel-based method to calculate local field potentials from networks of spiking neurons

BACKGROUND

The local field potential (LFP) is usually calculated from current sources arising from transmembrane currents, in particular in asymmetric cellular morphologies such as pyramidal neurons.

NEW METHOD

Here, we adopt a different point of view and relate the spiking of neurons to the LFP through efferent synaptic connections and provide a method to calculate LFPs.

RESULTS

We show that the so-called unitary LFPs (uLFP) provide the key to such a calculation. We show experimental measurements and simulations of uLFPs in neocortex and hippocampus, for both excitatory and inhibitory neurons. We fit a "kernel" function to measurements of uLFPs, and we estimate its spatial and temporal spread by using simulations of morphologically detailed reconstructions of hippocampal pyramidal neurons. Assuming that LFPs are the sum of uLFPs generated by every neuron in the network, the LFP generated by excitatory and inhibitory neurons can be calculated by convolving the trains of action potentials with the kernels estimated from uLFPs. This provides a method to calculate the LFP from networks of spiking neurons, even for point neurons for which the LFP is not easily defined. We show examples of LFPs calculated from networks of point neurons and compare to the LFP calculated from synaptic currents.

CONCLUSIONS

The kernel-based method provides a practical way to calculate LFPs from networks of point neurons.

Effects of Transcranial Magnetic Stimulation Therapy on Evoked and Induced Gamma Oscillations in Children with Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a behaviorally diagnosed neurodevelopmental condition of unknown pathology. Research suggests that abnormalities of elecltroencephalogram (EEG) gamma oscillations may provide a biomarker of the condition. In this study, envelope analysis of demodulated waveforms for evoked and induced gamma oscillations in response to Kanizsa figures in an oddball task were analyzed and compared in 19 ASD and 19 age/gender-matched neurotypical children. The ASD group was treated with low frequency transcranial magnetic stimulation (TMS), (1.0 Hz, 90% motor threshold, 18 weekly sessions) targeting the dorsolateral prefrontal cortex. In ASD subjects, as compared to neurotypicals, significant differences in evoked and induced gamma oscillations were evident in higher magnitude of gamma oscillations pre-TMS, especially in response to non-target cues. Recordings post-TMS treatment in ASD revealed a significant reduction of gamma responses to task-irrelevant stimuli. Participants committed fewer errors post-TMS. Behavioral questionnaires showed a decrease in irritability, hyperactivity, and repetitive behavior scores. The use of a novel metric for gamma oscillations. i.e., envelope analysis using wavelet transformation allowed for characterization of the impedance of the originating neuronal circuit. The results suggest that gamma oscillations may provide a biomarker reflective of the excitatory/inhibitory balance of the cortex and a putative outcome measure for interventions in autism.