Relationship between the mechanisms of gamma rhythm generation and the magnitude of the macroscopic phase response function in a population of excitatory and inhibitory modified quadratic integrate-and-fire neurons

Anatomically resolved oscillatory bursts reveal dynamic motifs of thalamocortical activity during naturalistic stimulus viewing

Natural vision requires circuit mechanisms which process complex spatiotemporal stimulus features in parallel. In the mammalian forebrain, one signature of circuit activation is fast oscillatory dynamics, reflected in the local field potential (LFP). Using data from the Allen Neuropixels Visual Coding project, we show that local visual features in naturalistic stimuli induce in mouse primary visual cortex (V1) retinotopically specific oscillations in various frequency bands and V1 layers. Specifically, layer 4 (L4) narrowband gamma was linked to luminance, low-gamma to optic flow, and L4/L5 epsilon oscillations to contrast. These feature-specific oscillations were associated with distinct translaminar spike-phase coupling patterns, which were conserved across a range of stimuli containing the relevant visual features, suggesting that they might constitute feature-specific circuit motifs. Our findings highlight visually induced fast oscillations as markers of dynamic circuit motifs, which may support differential and multiplexed coding of complex visual input and thalamocortical information propagation.

Distributed Synaptic Connection Strength Changes Dynamics in a Population Firing Rate Model in Response to Continuous External Stimuli

Neural network complexity allows for diverse neuronal population dynamics and realizes higherorder brain functions such as cognition and memory. Complexity is enhanced through chemical synapses with exponentially decaying conductance and greater variation in the neuronal connection strength due to synaptic plasticity. However, in the macroscopic neuronal population model, synaptic connections are often described by spike connections, and connection strengths within the population are assumed to be uniform. Thus, the effects of synaptic connections variation on network synchronization remain unclear. Based on recent advances in mean field theory for the quadratic integrate-and-fire neuronal network model, we introduce synaptic conductance and variation of connection strength into the excitatory and inhibitory neuronal population model and derive the macroscopic firing rate equations for faithful modeling. We then introduce a heuristic switching rule of the dynamic system with respect to the mean membrane potentials to avoid divergences in the computation caused by variations in the neuronal connection strength. We show that the switching rule agrees with the numerical computation of the microscopic level model. In the derived model, variations in synaptic conductance and connection strength strongly alter the stability of the solutions to the equations, which is related to the mechanism of synchronous firing. When we apply physiologically plausible values from layer 4 of the mammalian primary visual cortex to the derived model, we observe event-related desynchronization at the alpha and beta frequencies and event-related synchronization at the gamma frequency over a wide range of balanced external currents. Our results show that the introduction of complex synaptic connections and physiologically valid numerical values into the low-dimensional mean field equations reproduces dynamic changes such as eventrelated (de)synchronization, and provides a unique mathematical insight into the relationship between synaptic strength variation and oscillatory mechanism.

Entrainment by transcranial alternating current stimulation: Insights from models of cortical oscillations and dynamical systems theory

Signature of neuronal oscillations can be found in nearly every brain function. However, abnormal oscillatory activity is linked with several brain disorders. Transcranial alternating current stimulation (tACS) is a non-invasive brain stimulation technique that can potentially modulate neuronal oscillations and influence behavior both in health and disease. Yet, a complete understanding of how interacting networks of neurons are affected by tACS remains elusive. Entrainment effects by which tACS synchronizes neuronal oscillations is one of the main hypothesized mechanisms, as evidenced in animals and humans. Computational models of cortical oscillations may shed light on the entrainment effects of tACS, but current modeling studies lack specific guidelines to inform experimental investigations. This study addresses the existing gap in understanding the mechanisms of tACS effects on rhythmogenesis within the brain by providing a comprehensive overview of both theoretical and experimental perspectives. We explore the intricate interactions between oscillators and periodic stimulation through the lens of dynamical systems theory. Subsequently, we present a synthesis of experimental findings that demonstrate the effects of tACS on both individual neurons and collective oscillatory patterns in animal models and humans. Our review extends to computational investigations that elucidate the interplay between tACS and neuronal dynamics across diverse cortical network models. To illustrate these concepts, we conclude with a simple oscillatory neuron model, showcasing how fundamental theories of oscillatory behavior derived from dynamical systems, such as phase response of neurons to external perturbation, can account for the entrainment effects observed with tACS. Studies reviewed here render the necessity of integrated experimental and computational approaches for effective neuromodulation by tACS in health and disease.

Analyzing top-down visual attention in the context of gamma oscillations: a layer- dependent network-of- networks approach

Top-down visual attention is a fundamental cognitive process that allows individuals to selectively attend to salient visual stimuli in the environment. Recent empirical findings have revealed that gamma oscillations participate in the modulation of visual attention. However, computational studies face challenges when analyzing the attentional process in the context of gamma oscillation due to the unstable nature of gamma oscillations and the complexity induced by the layered fashion in the visual cortex. In this study, we propose a layer-dependent network-of-networks approach to analyze such attention with gamma oscillations. The model is validated by reproducing empirical findings on orientation preference and the enhancement of neuronal response due to top-down attention. We perform parameter plane analysis to classify neuronal responses into several patterns and find that the neuronal response to sensory and attention signals was modulated by the heterogeneity of the neuronal population. Furthermore, we revealed a counter-intuitive scenario that the excitatory populations in layer 2/3 and layer 5 exhibit opposite responses to the attentional input. By modification of the original model, we confirmed layer 6 plays an indispensable role in such cases. Our findings uncover the layer-dependent dynamics in the cortical processing of visual attention and open up new possibilities for further research on layer-dependent properties in the cerebral cortex.

Macroscopic Gamma Oscillation With Bursting Neuron Model Under Stochastic Fluctuation

Gamma oscillations are thought to play a role in information processing in the brain. Bursting neurons, which exhibit periodic clusters of spiking activity, are a type of neuron that are thought to contribute largely to gamma oscillations. However, little is known about how the properties of bursting neurons affect the emergence of gamma oscillation, its waveforms, and its synchronized characteristics, especially when subjected to stochastic fluctuations. In this study, we proposed a bursting neuron model that can analyze the bursting ratio and the phase response function. Then we theoretically analyzed the neuronal population dynamics composed of bursting excitatory neurons, mixed with inhibitory neurons. The bifurcation analysis of the equivalent Fokker-Planck equation exhibits three types of gamma oscillations of unimodal firing, bimodal firing in the inhibitory population, and bimodal firing in the excitatory population under different interaction strengths. The analyses of the macroscopic phase response function by the adjoint method of the Fokker-Planck equation revealed that the inhibitory doublet facilitates synchronization of the high-frequency oscillations. When we keep the strength of interactions constant, decreasing the bursting ratio of the individual neurons increases the relative high-gamma component of the populational phase-coupling functions. This also improves the ability of the neuronal population model to synchronize with faster oscillatory input. The analytical frameworks in this study provide insight into nontrivial dynamics of the population of bursting neurons, which further suggest that bursting neurons have an important role in rhythmic activities.

A framework for macroscopic phase-resetting curves for generalised spiking neural networks

Brain rhythms emerge from synchronization among interconnected spiking neurons. Key properties of such rhythms can be gleaned from the phase-resetting curve (PRC). Inferring the PRC and developing a systematic phase reduction theory for large-scale brain rhythms remains an outstanding challenge. Here we present a theoretical framework and methodology to compute the PRC of generic spiking networks with emergent collective oscillations. We adopt a renewal approach where neurons are described by the time since their last action potential, a description that can reproduce the dynamical feature of many cell types. For a sufficiently large number of neurons, the network dynamics are well captured by a continuity equation known as the refractory density equation. We develop an adjoint method for this equation giving a semi-analytical expression of the infinitesimal PRC. We confirm the validity of our framework for specific examples of neural networks. Our theoretical framework can link key biological properties at the individual neuron scale and the macroscopic oscillatory network properties. Beyond spiking networks, the approach is applicable to a broad class of systems that can be described by renewal processes.

The formation of oscillatory neuronal assemblies at the network level has been hypothesized to be fundamental to many cognitive and motor functions. One prominent tool to understand the dynamics of oscillatory activity response to stimuli, and hence the neural code for which it is a substrate, is a nonlinear measure called Phase-Resetting Curve (PRC). At the network scale, the PRC defines the measure of how a given synaptic input perturbs the timing of next upcoming volley of spike assemblies: either advancing or delaying this timing. As a further application, one can use PRCs to make unambiguous predictions about whether communicating networks of neurons will phase-lock as it is often observed across the cortical areas and what would be this stable phase-configuration: synchronous, asynchronous or with asymmetric phase-shifts. The latter configuration also implies a preferential flow of information form the leading network to the follower, thereby giving causal signatures of directed functional connectivity. Because of the key position of the PRC in studying synchrony, information flow and entrainment to external forcing, it is crucial to move toward a theory that allows to compute the PRCs of network-wide oscillations not only for a restricted class of models, as has been done in the past, but to network descriptions that are generalized and can reflect flexibly single cell properties. In this manuscript, we tackle this issue by showing how the PRC for network oscillations can be computed using the adjoint systems of partial differential equations that define the dynamics of the neural activity density.

Paradoxical phase response of gamma rhythms facilitates their entrainment in heterogeneous networks

The synchronization of different γ-rhythms arising in different brain areas has been implicated in various cognitive functions. Here, we focus on the effect of the ubiquitous neuronal heterogeneity on the synchronization of ING (interneuronal network gamma) and PING (pyramidal-interneuronal network gamma) rhythms. The synchronization properties of rhythms depends on the response of their collective phase to external input. We therefore determine the macroscopic phase-response curve for finite-amplitude perturbations (fmPRC) of ING- and PING-rhythms in all-to-all coupled networks comprised of linear (IF) or quadratic (QIF) integrate-and-fire neurons. For the QIF networks we complement the direct simulations with the adjoint method to determine the infinitesimal macroscopic PRC (imPRC) within the exact mean-field theory. We show that the intrinsic neuronal heterogeneity can qualitatively modify the fmPRC and the imPRC. Both PRCs can be biphasic and change sign (type II), even though the phase-response curve for the individual neurons is strictly non-negative (type I). Thus, for ING rhythms, say, external inhibition to the inhibitory cells can, in fact, advance the collective oscillation of the network, even though the same inhibition would lead to a delay when applied to uncoupled neurons. This paradoxical advance arises when the external inhibition modifies the internal dynamics of the network by reducing the number of spikes of inhibitory neurons; the advance resulting from this disinhibition outweighs the immediate delay caused by the external inhibition. These results explain how intrinsic heterogeneity allows ING- and PING-rhythms to become synchronized with a periodic forcing or another rhythm for a wider range in the mismatch of their frequencies. Our results identify a potential function of neuronal heterogeneity in the synchronization of coupled γ-rhythms, which may play a role in neural information transfer via communication through coherence.

The interaction of a large number of oscillating units can lead to the emergence of a collective, macroscopic oscillation in which many units oscillate in near-unison or near-synchrony. This has been exploited technologically, e.g., to combine many coherently interacting, individual lasers to form a single powerful laser. Collective oscillations are also important in biology. For instance, the circadian rhythm of animals is controlled by the near-synchronous dynamics of a large number of individually oscillating cells. In animals and humans brain rhythms reflect the coherent dynamics of a large number of neurons and are surmised to play an important role in the communication between different brain areas. To be functionally relevant, these rhythms have to respond to external inputs and have to be able to synchronize with each other. We show that the ubiquitous heterogeneity in the properties of the individual neurons in a network can contribute to that ability. It can allow the external inputs to modify the internal network dynamics such that the network can follow these inputs over a wider range of frequencies. Paradoxically, while an external perturbation may delay individual neurons, their ensuing within-network interaction can overcompensate this delay, leading to an overall advance of the rhythm.

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.

Theta-Nested Gamma Oscillations in Next Generation Neural Mass Models

Theta-nested gamma oscillations have been reported in many areas of the brain and are believed to represent a fundamental mechanism to transfer information across spatial and temporal scales. In a series of recent experiments in vitro it has been possible to replicate with an optogenetic theta frequency stimulation several features of cross-frequency coupling (CFC) among theta and gamma rhythms observed in behaving animals. In order to reproduce the main findings of these experiments we have considered a new class of neural mass models able to reproduce exactly the macroscopic dynamics of spiking neural networks. In this framework, we have examined two set-ups able to support collective gamma oscillations: namely, the pyramidal interneuronal network gamma (PING) and the interneuronal network gamma (ING). In both set-ups we observe the emergence of theta-nested gamma oscillations by driving the system with a sinusoidal theta-forcing in proximity of a Hopf bifurcation. These mixed rhythms always display phase amplitude coupling. However, two different types of nested oscillations can be identified: one characterized by a perfect phase locking between theta and gamma rhythms, corresponding to an overall periodic behavior; another one where the locking is imperfect and the dynamics is quasi-periodic or even chaotic. From our analysis it emerges that the locked states are more frequent in the ING set-up. In agreement with the experiments, we find theta-nested gamma oscillations for forcing frequencies in the range [1:10] Hz, whose amplitudes grow proportionally to the forcing intensity and which are clearly modulated by the theta phase. Furthermore, analogously to the experiments, the gamma power and the frequency of the gamma-power peak increase with the forcing amplitude. At variance with experimental findings, the gamma-power peak does not shift to higher frequencies by increasing the theta frequency. This effect can be obtained, in our model, only by incrementing, at the same time, also the stimulation power. An effect achieved by increasing the amplitude either of the noise or of the forcing term proportionally to the theta frequency. On the basis of our analysis both the PING and the ING mechanism give rise to theta-nested gamma oscillations with almost identical features.

Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits

Macroscopic oscillations of different brain regions show multiple phase relationships that are persistent across time and have been implicated in routing information. While multiple cellular mechanisms influence the network oscillatory dynamics and structure the macroscopic firing motifs, one of the key questions is to identify the biophysical neuronal and synaptic properties that permit such motifs to arise. A second important issue is how the different neural activity coherence states determine the communication between the neural circuits. Here we analyse the emergence of phase-locking within bidirectionally delayed-coupled spiking circuits in which global gamma band oscillations arise from synaptic coupling among largely excitable neurons. We consider both the interneuronal (ING) and the pyramidal-interneuronal (PING) population gamma rhythms and the inter coupling targeting the pyramidal or the inhibitory neurons. Using a mean-field approach together with an exact reduction method, we reduce each spiking network to a low dimensional nonlinear system and derive the macroscopic phase resetting-curves (mPRCs) that determine how the phase of the global oscillation responds to incoming perturbations. This is made possible by the use of the quadratic integrate-and-fire model together with a Lorentzian distribution of the bias current. Depending on the type of gamma (PING vs. ING), we show that incoming excitatory inputs can either speed up the macroscopic oscillation (phase advance; type I PRC) or induce both a phase advance and a delay (type II PRC). From there we determine the structure of macroscopic coherence states (phase-locking) of two weakly synaptically-coupled networks. To do so we derive a phase equation for the coupled system which links the synaptic mechanisms to the coherence states of the system. We show that a synaptic transmission delay is a necessary condition for symmetry breaking, i.e. a non-symmetric phase lag between the macroscopic oscillations. This potentially provides an explanation to the experimentally observed variety of gamma phase-locking modes. Our analysis further shows that symmetry-broken coherence states can lead to a preferred direction of signal transfer between the oscillatory networks where this directionality also depends on the timing of the signal. Hence we suggest a causal theory for oscillatory modulation of functional connectivity between cortical circuits.

Large scale brain oscillations emerge from synaptic interactions within neuronal circuits. Over the past years, such macroscopic rhythms have been suggested to play a crucial role in routing the flow of information across cortical regions, resulting in a functional connectome. The underlying mechanism is cortical oscillations that bind together following a well-known motif called phase-locking. While there is significant experimental support for multiple phase-locking modes in the brain, it is still unclear what is the underlying mechanism that permits macroscopic rhythms to phase lock. In the present paper we take up with this issue, and to show that, one can study the emergent macroscopic phase-locking within the mathematical framework of weakly coupled oscillators. We find that under synaptic delays, fully symmetrically coupled networks can display symmetry-broken states of activity, where one network starts to lead in phase the second (also sometimes known as stuttering states). When we analyse how incoming transient signals affect the coupled system, we find that in the symmetry-broken state, the effect depends strongly on which network is targeted (the leader or the follower) as well as the timing of the input. Hence we show how the dynamics of the emergent phase-locked activity imposes a functional directionality on how signals are processed. We thus offer clarification on the synaptic and circuit properties responsible for the emergence of multiple phase-locking patterns and provide support for its functional implication in information transfer.

Analogue closed-loop optogenetic modulation of hippocampal pyramidal cells dissociates gamma frequency and amplitude

Gamma-band oscillations are implicated in modulation of attention, integration of sensory information and flexible communication among anatomically connected brain areas. How networks become entrained is incompletely understood. Specifically, it is unclear how the spectral and temporal characteristics of network oscillations can be altered on rapid timescales needed for efficient communication. We use closed-loop optogenetic modulation of principal cell excitability in mouse hippocampal slices to interrogate the dynamical properties of hippocampal oscillations. Gamma frequency and amplitude can be modulated bi-directionally, and dissociated, by phase-advancing or delaying optogenetic feedback to pyramidal cells. Closed-loop modulation alters the synchrony rather than average frequency of action potentials, in principle avoiding disruption of population rate-coding of information. Modulation of phasic excitatory currents in principal neurons is sufficient to manipulate oscillations, suggesting that feed-forward excitation of pyramidal cells has an important role in determining oscillatory dynamics and the ability of networks to couple with one another.