A stochastic model of input effectiveness during irregular gamma rhythms

The stabilized supralinear network accounts for the contrast dependence of visual cortical gamma oscillations

When stimulated, neural populations in the visual cortex exhibit fast rhythmic activity with frequencies in the gamma band (30-80 Hz). The gamma rhythm manifests as a broad resonance peak in the power-spectrum of recorded local field potentials, which exhibits various stimulus dependencies. In particular, in macaque primary visual cortex (V1), the gamma peak frequency increases with increasing stimulus contrast. Moreover, this contrast dependence is local: when contrast varies smoothly over visual space, the gamma peak frequency in each cortical column is controlled by the local contrast in that column's receptive field. No parsimonious mechanistic explanation for these contrast dependencies of V1 gamma oscillations has been proposed. The stabilized supralinear network (SSN) is a mechanistic model of cortical circuits that has accounted for a range of visual cortical response nonlinearities and contextual modulations, as well as their contrast dependence. Here, we begin by showing that a reduced SSN model without retinotopy robustly captures the contrast dependence of gamma peak frequency, and provides a mechanistic explanation for this effect based on the observed non-saturating and supralinear input-output function of V1 neurons. Given this result, the local dependence on contrast can trivially be captured in a retinotopic SSN which however lacks horizontal synaptic connections between its cortical columns. However, long-range horizontal connections in V1 are in fact strong, and underlie contextual modulation effects such as surround suppression. We thus explored whether a retinotopically organized SSN model of V1 with strong excitatory horizontal connections can exhibit both surround suppression and the local contrast dependence of gamma peak frequency. We found that retinotopic SSNs can account for both effects, but only when the horizontal excitatory projections are composed of two components with different patterns of spatial fall-off with distance: a short-range component that only targets the source column, combined with a long-range component that targets columns neighboring the source column. We thus make a specific qualitative prediction for the spatial structure of horizontal connections in macaque V1, consistent with the columnar structure of cortex.

The Thalamocortical Mechanism Underlying the Generation and Regulation of the Auditory Steady-State Responses in Awake Mice

The auditory steady-state response (ASSR) is a cortical oscillation induced by trains of 40 Hz acoustic stimuli. While the ASSR has been widely used in clinic measurement, the underlying neural mechanism remains poorly understood. In this study, we investigated the contribution of different stages of auditory thalamocortical pathway-medial geniculate body (MGB), thalamic reticular nucleus (TRN), and auditory cortex (AC)-to the generation and regulation of 40 Hz ASSR in C57BL/6 mice of both sexes. We found that the neural response synchronizing to 40 Hz sound stimuli was most prominent in the GABAergic neurons in the granular layer of AC and the ventral division of MGB (MGBv), which were regulated by optogenetic manipulation of TRN neurons. Behavioral experiments confirmed that disrupting TRN activity has a detrimental effect on the ability of mice to discriminate 40 Hz sounds. These findings revealed a thalamocortical mechanism helpful to interpret the results of clinical ASSR examinations.Significance Statement Our study contributes to clarifying the thalamocortical mechanisms underlying the generation and regulation of the auditory steady-state response (ASSR), which is commonly used in both clinical and neuroscience research to assess the integrity of auditory function. Combining a series of electrophysiological and optogenetic experiments, we demonstrate that the generation of cortical ASSR is dependent on the lemniscal thalamocortical projections originating from the ventral division of medial geniculate body to the GABAergic interneurons in the granule layer of the auditory cortex. Furthermore, the thalamocortical process for ASSR is strictly regulated by the activity of thalamic reticular nucleus (TRN) neurons. Behavioral experiments confirmed that dysfunction of TRN would cause a disruption of mice's behavioral performance in the auditory discrimination task.

Mutual information resonances in delay-coupled limit cycle and quasi-cycle brain rhythms

We elucidate how coupling delays and noise impact phase and mutual information relationships between two stochastic brain rhythms. This impact depends on the dynamical regime of each PING-based rhythm, as well as on network heterogeneity and coupling asymmetry. The number of peaks at positive and negative time lags in the delayed mutual information between the two bi-directionally communicating rhythms defines our measure of flexibility of information sharing and reflects the number of ways in which the two networks can alternately lead one another. We identify two distinct mechanisms for the appearance of qualitatively similar flexible information sharing. The flexibility in the quasi-cycle regime arises from the coupling delay-induced bimodality of the phase difference distribution, and the related bimodal mutual information. It persists in the presence of asymmetric coupling and heterogeneity but is limited to two routes of information sharing. The second mechanism in noisy limit cycle regime is not induced by the delay. However, delay-coupling and heterogeneity enable communication routes at multiple time lags. Noise disrupts the shared compromise frequency, allowing the expression of individual network frequencies which leads to a slow beating pattern. Simulations of an envelope-phase description for delay-coupled quasi-cycles yield qualitatively similar properties as for the full system. Near the bifurcation from in-phase to out-of-phase behaviour, a single preferred phase difference can coexist with two information sharing routes; further, the phase laggard can be the mutual information leader, or vice versa. Overall, the coupling delay endows a two-rhythm system with an array of lead-lag relationships and mutual information resonances that exist in spite of the noise and across the Hopf bifurcation. These beg to be mapped out experimentally with the help of our predictions.

Determinants of Brain Rhythm Burst Statistics

Brain rhythms recorded in vivo, such as gamma oscillations, are notoriously variable both in amplitude and frequency. They are characterized by transient epochs of higher amplitude known as bursts. It has been suggested that, despite their short-life and random occurrence, bursts in gamma and other rhythms can efficiently contribute to working memory or communication tasks. Abnormalities in bursts have also been associated with e.g. motor and psychiatric disorders. It is thus crucial to understand how single cell and connectivity parameters influence burst statistics and the corresponding brain states. To address this problem, we consider a generic stochastic recurrent network of Pyramidal Interneuron Network Gamma (PING) type. Using the stochastic averaging method, we derive dynamics for the phase and envelope of the amplitude process, and find that they depend on only two meta-parameters that combine all the model parameters. This allows us to identify an optimal parameter regime of healthy variability with similar statistics to those seen in vivo; in this regime, oscillations and bursts are supported by synaptic noise. The probability density for the rhythm’s envelope as well as the mean burst duration are then derived using first passage time analysis. Our analysis enables us to link burst attributes, such as duration and frequency content, to system parameters. Our general approach can be extended to different frequency bands, network topologies and extra populations. It provides the much needed insight into the biophysical determinants of rhythm burst statistics, and into what needs to be changed to correct rhythms with pathological statistics.

Non-monotonic accumulation of spike time variance during membrane potential oscillations

A spike-phase neural code has been proposed as a mechanism to encode stimuli based on the precise timing of spikes relative to the phase of membrane potential oscillations. This form of coding has been reported in both in vivo and in vitro experiments across several regions of the brain, yet there are concerns that such precise timing may be compromised by an effect referred to as variance accumulation, wherein spike timing variance increases over the phase of an oscillation. Here, we provide a straightforward explanation of this effect based on the theoretical spike time variance. The proposed theory is consistent with recordings of mitral neurons. It shows that spike time variance can increase in a nonlinear fashion with spike number, in a way that is dependent upon the frequency and amplitude of the oscillation. Further, non-monotonic accumulation of variance can arise from different combinations of oscillation parameters. Nonlinear accumulation sometimes leads to lower variance than that of a mean rate-matched homogeneous Poisson process, particularly for spikes that occur in later phases of oscillation. However, such an advantage is limited to a narrow range of oscillation amplitudes and frequencies. These results suggest fundamental constraints on spike-phase coding, and reveal how certain spikes in a sequence may exhibit increased firing time precision relative to their neighbors.