Temporal characteristics of gamma rhythm constrain properties of noise in an inhibition-stabilized network model

Gamma rhythm refers to oscillatory neural activity between 30 and 80 Hz, induced in visual cortex by stimuli such as iso-luminant hues or gratings. The power and peak frequency of gamma depend on the properties of the stimulus such as size and contrast. Gamma waveform is typically arch-shaped, with narrow troughs and broad peaks, and can be replicated in a self-oscillating Wilson-Cowan (WC) model operating in an appropriate regime. However, oscillations in this model are infinitely long, unlike physiological gamma that occurs in short bursts. Further, unlike the model, gamma is faster after stimulus onset and slows down over time. Here, we first characterized gamma burst duration in local field potential data recorded from two monkeys as they viewed full screen iso-luminant hues. We then added different types of noise in the inputs to the WC model and tested how that affected duration and temporal dynamics of gamma. While the model failed with the often-used Poisson noise, Ornstein-Uhlenbeck noise applied to both the excitatory and the inhibitory populations replicated the duration and slowing of gamma and replicated the shape and stimulus dependencies. Thus, the temporal dynamics of gamma oscillations put constraints on the type and properties of underlying neural noise.

Hierarchical Wilson-Cowan Models and Connection Matrices

This work aims to study the interplay between the Wilson-Cowan model and connection matrices. These matrices describe cortical neural wiring, while Wilson-Cowan equations provide a dynamical description of neural interaction. We formulate Wilson-Cowan equations on locally compact Abelian groups. We show that the Cauchy problem is well posed. We then select a type of group that allows us to incorporate the experimental information provided by the connection matrices. We argue that the classical Wilson-Cowan model is incompatible with the small-world property. A necessary condition to have this property is that the Wilson-Cowan equations be formulated on a compact group. We propose a p-adic version of the Wilson-Cowan model, a hierarchical version in which the neurons are organized into an infinite rooted tree. We present several numerical simulations showing that the p-adic version matches the predictions of the classical version in relevant experiments. The p-adic version allows the incorporation of the connection matrices into the Wilson-Cowan model. We present several numerical simulations using a neural network model that incorporates a p-adic approximation of the connection matrix of the cat cortex.

Beyond Wilson-Cowan dynamics: oscillations and chaos without inhibition

Fifty years ago, Wilson and Cowan developed a mathematical model to describe the activity of neural populations. In this seminal work, they divided the cells in three groups: active, sensitive and refractory, and obtained a dynamical system to describe the evolution of the average firing rates of the populations. In the present work, we investigate the impact of the often neglected refractory state and show that taking it into account can introduce new dynamics. Starting from a continuous-time Markov chain, we perform a rigorous derivation of a mean-field model that includes the refractory fractions of populations as dynamical variables. Then, we perform bifurcation analysis to explain the occurrence of periodic solutions in cases where the classical Wilson-Cowan does not predict oscillations. We also show that our mean-field model is able to predict chaotic behavior in the dynamics of networks with as little as two populations.

Dynamical Mechanism Analysis of Three Neuroregulatory Strategies on the Modulation of Seizures

This paper attempts to explore and compare the regulatory mechanisms of optogenetic stimulation (OS), deep brain stimulation (DBS) and electromagnetic induction on epilepsy. Based on the Wilson-Cowan model, we first demonstrate that the external input received by excitatory and inhibitory neural populations can induce rich dynamic bifurcation behaviors such as Hopf bifurcation, and make the system exhibit epileptic and normal states. Then, both OS and DBS are shown to be effective in controlling the epileptic state to a normal low-level state, and the stimulus parameters have a broad effective range. However, electromagnetic induction cannot directly control epilepsy to this desired state, even if it can significantly reduce the oscillation frequency of neural populations. One main difference worth noting is that the high spatiotemporal specificity of OS allows it to target inhibitory neuronal populations, whereas DBS and electromagnetic induction can only stimulate excitatory as well as inhibitory neuronal populations together. Next, the propagation behavior of epilepsy is explored under a typical three-node feedback loop structure. An increase in coupling strength accelerates and exacerbates epileptic activity in other brain regions. Finally, OS and DBS applied to the epileptic focus play similar positive roles in controlling the behavior of the area of seizure propagation, while electromagnetic induction still only achieves unsatisfactory effects. It is hoped that these dynamical results can provide insights into the treatment of epilepsy as well as other neurological disorders.

Phase response approaches to neural activity models with distributed delay

In weakly coupled neural oscillator networks describing brain dynamics, the coupling delay is often distributed. We present a theoretical framework to calculate the phase response curve of distributed-delay induced limit cycles with infinite-dimensional phase space. Extending previous works, in which non-delayed or discrete-delay systems were investigated, we develop analytical results for phase response curves of oscillatory systems with distributed delay using Gaussian and log-normal delay distributions. We determine the scalar product and normalization condition for the linearized adjoint of the system required for the calculation of the phase response curve. As a paradigmatic example, we apply our technique to the Wilson-Cowan oscillator model of excitatory and inhibitory neuronal populations under the two delay distributions. We calculate and compare the phase response curves for the Gaussian and log-normal delay distributions. The phase response curves obtained from our adjoint calculations match those compiled by the direct perturbation method, thereby proving that the theory of weakly coupled oscillators can be applied successfully for distributed-delay-induced limit cycles. We further use the obtained phase response curves to derive phase interaction functions and determine the possible phase locked states of multiple inter-coupled populations to illuminate different synchronization scenarios. In numerical simulations, we show that the coupling delay distribution can impact the stability of the synchronization between inter-coupled gamma-oscillatory networks.

Divisive gain modulation enables flexible and rapid entrainment in a neocortical microcircuit model

Neocortical circuits exhibit a rich dynamic repertoire, and their ability to achieve entrainment (adjustment of their frequency to match the input frequency) is thought to support many cognitive functions and indicate functional flexibility. Although previous studies have explored the influence of various circuit properties on this phenomenon, the role of divisive gain modulation (or divisive inhibition) is unknown. This gain control mechanism is thought to be delivered mainly by the soma-targeting interneurons in neocortical microcircuits. In this study, we use a neural mass model of the neocortical microcircuit (extended Wilson-Cowan model) featuring both soma-targeting and dendrite-targeting interneuronal subpopulations to investigate the role of divisive gain modulation in entrainment. Our results demonstrate that the presence of divisive inhibition in the microcircuit, as delivered by the soma-targeting interneurons, enables its entrainment to a wider range of input frequencies. Divisive inhibition also promotes a faster entrainment, with the microcircuit needing less time to converge to the fully entrained state. We suggest that divisive inhibition, working alongside subtractive inhibition, allows for more adaptive oscillatory responses in neocortical circuits and, thus, supports healthy brain functioning.NEW & NOTEWORTHY We introduce a computational neocortical microcircuit model that features two inhibitory neural populations, with one providing subtractive and the other divisive inhibition to the excitatory population. We demonstrate that divisive inhibition widens the range of input frequencies to which the microcircuit can become entrained and diminishes the time needed to reach full entrainment. We suggest that divisive inhibition enables more adaptive oscillatory activity, with important implications for both normal and pathological brain function.

Noise-induced precursors of state transitions in the stochastic Wilson-cowan model

The Wilson-Cowan neural field equations describe the dynamical behavior of a 1-D continuum of excitatory and inhibitory cortical neural aggregates, using a pair of coupled integro-differential equations. Here we use bifurcation theory and small-noise linear stochastics to study the range of a phase transitions-sudden qualitative changes in the state of a dynamical system emerging from a bifurcation-accessible to the Wilson-Cowan network. Specifically, we examine saddle-node, Hopf, Turing, and Turing-Hopf instabilities. We introduce stochasticity by adding small-amplitude spatio-temporal white noise, and analyze the resulting subthreshold fluctuations using an Ornstein-Uhlenbeck linearization. This analysis predicts divergent changes in correlation and spectral characteristics of neural activity during close approach to bifurcation from below. We validate these theoretical predictions using numerical simulations. The results demonstrate the role of noise in the emergence of critically slowed precursors in both space and time, and suggest that these early-warning signals are a universal feature of a neural system close to bifurcation. In particular, these precursor signals are likely to have neurobiological significance as early warnings of impending state change in the cortex. We support this claim with an analysis of the in vitro local field potentials recorded from slices of mouse-brain tissue. We show that in the period leading up to emergence of spontaneous seizure-like events, the mouse field potentials show a characteristic spectral focusing toward lower frequencies concomitant with a growth in fluctuation variance, consistent with critical slowing near a bifurcation point. This observation of biological criticality has clear implications regarding the feasibility of seizure prediction.

On the Influence of Amplitude on the Connectivity between Phases

In recent studies, functional connectivities have been reported to display characteristics of complex networks that have been suggested to concur with those of the underlying structural, i.e., anatomical, networks. Do functional networks always agree with structural ones? In all generality, this question can be answered with "no": for instance, a fully synchronized state would imply isotropic homogeneous functional connections irrespective of the "real" underlying structure. A proper inference of structure from function and vice versa requires more than a sole focus on phase synchronization. We show that functional connectivity critically depends on amplitude variations, which implies that, in general, phase patterns should be analyzed in conjunction with the corresponding amplitude. We discuss this issue by comparing the phase synchronization patterns of interconnected Wilson-Cowan models vis-à-vis Kuramoto networks of phase oscillators. For the interconnected Wilson-Cowan models we derive analytically how connectivity between phases explicitly depends on the generating oscillators' amplitudes. In consequence, the link between neurophysiological studies and computational models always requires the incorporation of the amplitude dynamics. Supplementing synchronization characteristics by amplitude patterns, as captured by, e.g., spectral power in M/EEG recordings, will certainly aid our understanding of the relation between structural and functional organizations in neural networks at large.