Impact of neuronal heterogeneity on correlated colored noise-induced synchronization

Coherent olfactory bulb gamma oscillations arise from coupling independent columnar oscillators

Spike timing-based representations of sensory information depend on embedded dynamical frameworks within neuronal networks that establish the rules of local computation and interareal communication. Here, we investigated the dynamical properties of olfactory bulb circuitry in mice of both sexes using microelectrode array recordings from slice and in vivo preparations. Neurochemical activation or optogenetic stimulation of sensory afferents evoked persistent gamma oscillations in the local field potential. These oscillations arose from slower, GABA(A) receptor-independent intracolumnar oscillators coupled by GABA(A)-ergic synapses into a faster, broadly coherent network oscillation. Consistent with the theoretical properties of coupled-oscillator networks, the spatial extent of zero-phase coherence was bounded in slices by the reduced density of lateral interactions. The intact in vivo network, however, exhibited long-range lateral interactions that suffice in simulation to enable zero-phase gamma coherence across the olfactory bulb. The timing of action potentials in a subset of principal neurons was phase-constrained with respect to evoked gamma oscillations. Coupled-oscillator dynamics in olfactory bulb thereby enable a common clock, robust to biological heterogeneities, that is capable of supporting gamma-band spike synchronization and phase coding across the ensemble of activated principal neurons.NEW & NOTEWORTHY Odor stimulation evokes rhythmic gamma oscillations in the field potential of the olfactory bulb, but the dynamical mechanisms governing these oscillations have remained unclear. Establishing these mechanisms is important as they determine the biophysical capacities of the bulbar circuit to, for example, maintain zero-phase coherence across a spatially extended network, or coordinate the timing of action potentials in principal neurons. These properties in turn constrain and suggest hypotheses of sensory coding.

Induced synchronization by endogenous noise modulation in finite-size random neural networks: A stochastic mean-field study

Event-related synchronization and desynchronization (ERS/ERD) are well-known features found experimentally in brain signals during cognitive tasks. Their understanding promises to have much better insights into neural information processes in cognition. Under the hypothesis that neural information affects the endogenous neural noise level in populations, we propose to employ a stochastic mean-field model to explain ERS/ERD in the γ-frequency range. The work extends previous mean-field studies by deriving novel effects from finite network size. Moreover, numerical simulations of ERS/ERD and their analytical explanation by the mean-field model suggest several endogenous noise modulation schemes, which may modulate the system's synchronization.

The steady state and response to a periodic stimulation of the firing rate for a theta neuron with correlated noise

The stochastic activity of neurons is caused by various sources of correlated fluctuations and can be described in terms of simplified, yet biophysically grounded, integrate-and-fire models. One paradigmatic model is the quadratic integrate-and-fire model and its equivalent phase description by the theta neuron. Here we study the theta neuron model driven by a correlated Ornstein-Uhlenbeck noise and by periodic stimuli. We apply the matrix-continued-fraction method to the associated Fokker-Planck equation to develop an efficient numerical scheme to determine the stationary firing rate as well as the stimulus-induced modulation of the instantaneous firing rate. For the stationary case, we identify the conditions under which the firing rate decreases or increases by the effect of the colored noise and compare our results to existing analytical approximations for limit cases. For an additional periodic signal we demonstrate how the linear and nonlinear response terms can be computed and report resonant behavior for some of them. We extend the method to the case of two periodic signals, generally with incommensurable frequencies, and present a particular case for which a strong mixed response to both signals is observed, i.e. where the response to the sum of signals differs significantly from the sum of responses to the single signals. We provide Python code for our computational method: https://github.com/jannikfranzen/theta_neuron .

The Effects of Background Noise on a Biophysical Model of Olfactory Bulb Mitral Cells

The spiking activity of mitral cells (MC) in the olfactory bulb is a key attribute in olfactory sensory information processing to downstream cortical areas. A more detailed understanding of the modulation of MC spike statistics could shed light on mechanistic studies of olfactory bulb circuits and olfactory coding. We study the spike response of a recently developed single-compartment biophysical MC model containing seven known ionic currents and calcium dynamics subject to constant current input with background white noise. We observe rich spiking dynamics even with constant current input, including multimodal peaks in the interspike interval distribution (ISI). Although weak-to-moderate background noise for a fixed current input does not change the firing rate much, the spike dynamics can change dramatically, exhibiting non-monotonic spike variability not commonly observed in standard neuron models. We explain these dynamics with a phenomenological model of the ISI probability density function. Our study clarifies some of the complexities of MC spiking dynamics.

Deep brain stimulation in the subthalamic nucleus for Parkinson's disease can restore dynamics of striatal networks

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is highly effective in alleviating movement disability in patients with Parkinson’s disease (PD). However, its therapeutic mechanism of action is unknown. The healthy striatum exhibits rich dynamics resulting from an interaction of beta, gamma, and theta oscillations. These rhythms are essential to selection and execution of motor programs, and their loss or exaggeration due to dopamine (DA) depletion in PD is a major source of behavioral deficits. Restoring the natural rhythms may then be instrumental in the therapeutic action of DBS. We develop a biophysical networked model of a BG pathway to study how abnormal beta oscillations can emerge throughout the BG in PD and how DBS can restore normal beta, gamma, and theta striatal rhythms. Our model incorporates STN projections to the striatum, long known but understudied, found to preferentially target fast-spiking interneurons (FSI). We find that DBS in STN can normalize striatal medium spiny neuron activity by recruiting FSI dynamics and restoring the inhibitory potency of FSIs observed in normal conditions. We also find that DBS allows the reexpression of gamma and theta rhythms, thought to be dependent on high DA levels and thus lost in PD, through cortical noise control. Our study highlights that DBS effects can go beyond regularizing BG output dynamics to restoring normal internal BG dynamics and the ability to regulate them. It also suggests how gamma and theta oscillations can be leveraged to supplement DBS treatment and enhance its effectiveness.

Quantitative comparison of the mean-return-time phase and the stochastic asymptotic phase for noisy oscillators

Seminal work by A. Winfree and J. Guckenheimer showed that a deterministic phase variable can be defined either in terms of Poincaré sections or in terms of the asymptotic (long-time) behaviour of trajectories approaching a stable limit cycle. However, this equivalence between the deterministic notions of phase is broken in the presence of noise. Different notions of phase reduction for a stochastic oscillator can be defined either in terms of mean-return-time sections or as the argument of the slowest decaying complex eigenfunction of the Kolmogorov backwards operator. Although both notions of phase enjoy a solid theoretical foundation, their relationship remains unexplored. Here, we quantitatively compare both notions of stochastic phase. We derive an expression relating both notions of phase and use it to discuss differences (and similarities) between both definitions of stochastic phase for (i) a spiral sink motivated by stochastic models for electroencephalograms, (ii) noisy limit-cycle systems-neuroscience models, and (iii) a stochastic heteroclinic oscillator inspired by a simple motor-control system.

Afterhyperpolarization Promotes the Firing of Mitral Cells through a Voltage-Dependent Modification of Action Potential Threshold

In the olfactory bulb, mitral cells (MCs) display a spontaneous firing that is characterized by bursts of action potentials (APs) intermixed with silent periods. Intraburst firing frequency and duration are heterogeneous among MCs and increase with membrane depolarization. By using patch-clamp recording on rat slices, we dissected out the intrinsic properties responsible for this bursting activity. We showed that the threshold of AP generation dynamically changes as a function of the preceding trajectory of the membrane potential. In fact, the AP threshold became more negative when the membrane was hyperpolarized and had a recovery rate inversely proportional to the membrane repolarization rate. Such variations appeared to be produced by changes in the inactivation state of voltage-dependent Na⁺ channels. Thus, AP initiation was favored by hyperpolarizing events, such as negative membrane oscillations or inhibitory synaptic input. After the first AP, the following fast afterhyperpolarization (AHP) brought the threshold to more negative values and then promoted the emission of the following AP. This phenomenon was repeated for each AP of the burst making the fast AHP a regenerative mechanism that sustained the firing, AHP with larger amplitudes and faster repolarizations being associated with larger and higher-frequency bursts. Burst termination was found to be because of the development of a slow repolarization component of the AHP (slow AHP). Overall, the AHP characteristics appeared as a major determinant of the bursting properties.

Synchronization Through Uncorrelated Noise in Excitatory-Inhibitory Networks

Gamma rhythms play a major role in many different processes in the brain, such as attention, working memory, and sensory processing. While typically considered detrimental, counterintuitively noise can sometimes have beneficial effects on communication and information transfer. Recently, Meng and Riecke showed that synchronization of interacting networks of inhibitory neurons in the gamma band (i.e., gamma generated through an ING mechanism) increases while synchronization within these networks decreases when neurons are subject to uncorrelated noise. However, experimental and modeling studies point towardz an important role of the pyramidal-interneuronal network gamma (PING) mechanism in the cortex. Therefore, we investigated the effect of uncorrelated noise on the communication between excitatory-inhibitory networks producing gamma oscillations via a PING mechanism. Our results suggest that, at least in a certain range of noise strengths and natural frequency differences between the regions, synaptic noise can have a supporting role in facilitating inter-regional communication, similar to the ING case for a slightly larger parameter range. Furthermore, the noise-induced synchronization between networks is generated via a different mechanism than when synchronization is mediated by strong synaptic coupling. Noise-induced synchronization is achieved by lowering synchronization within networks which allows the respective other network to impose its own gamma rhythm resulting in synchronization between networks.

Analytical approach to the mean-return-time phase of isotropic stochastic oscillators

One notion of phase for stochastic oscillators is based on the mean return-time (MRT): a set of points represents a certain phase if the mean time to return from any point in this set to this set after one rotation is equal to the mean rotation period of the oscillator (irrespective of the starting point). For this so far only algorithmically defined phase, we derive here analytical expressions for the important class of isotropic stochastic oscillators. This allows us to evaluate cases from the literature explicitly and to study the behavior of the MRT phase in the limits of strong noise. We also use the same formalism to show that lines of constant return time variance (instead of constant mean return time) can be defined, and that they in general differ from the MRT isochrons.

Heterogeneity in GnRH and kisspeptin neurons and their significance in vertebrate reproductive biology

Vertebrate reproduction is essentially controlled by the hypothalamus-pituitary-gonadal (HPG) axis, which is a central dogma of reproductive biology. Two major hypothalamic neuroendocrine cell groups containing gonadotropin-releasing hormone (GnRH) and kisspeptin are crucial for control of the HPG axis in vertebrates. GnRH and kisspeptin neurons exhibit high levels of heterogeneity including their cellular morphology, biochemistry, neurophysiology and functions. However, the molecular foundation underlying heterogeneities in GnRH and kisspeptin neurons remains unknown. More importantly, the biological and physiological significance of their heterogeneity in reproductive biology is poorly understood. In this review, we first describe the recent advances in the neuroendocrine functions of kisspeptin-GnRH pathways. We then view the recent emerging progress in the heterogeneity of GnRH and kisspeptin neurons using morphological and single-cell transcriptomic analyses. Finally, we discuss our views on the significance of functional heterogeneity of reproductive endocrine cells and their potential relevance to reproductive health.

Cell and circuit origins of fast network oscillations in the mammalian main olfactory bulb

Neural synchrony generates fast network oscillations throughout the brain, including the main olfactory bulb (MOB), the first processing station of the olfactory system. Identifying the mechanisms synchronizing neurons in the MOB will be key to understanding how network oscillations support the coding of a high-dimensional sensory space. Here, using paired recordings and optogenetic activation of glomerular sensory inputs in MOB slices, we uncovered profound differences in principal mitral cell (MC) vs. tufted cell (TC) spike-time synchrony: TCs robustly synchronized across fast- and slow-gamma frequencies, while MC synchrony was weaker and concentrated in slow-gamma frequencies. Synchrony among both cell types was enhanced by shared glomerular input but was independent of intraglomerular lateral excitation. Cell-type differences in synchrony could also not be traced to any difference in the synchronization of synaptic inhibition. Instead, greater TC than MC synchrony paralleled the more periodic firing among resonant TCs than MCs and emerged in patterns consistent with densely synchronous network oscillations. Collectively, our results thus reveal a mechanism for parallel processing of sensory information in the MOB via differential TC vs. MC synchrony, and further contrast mechanisms driving fast network oscillations in the MOB from those driving the sparse synchronization of irregularly firing principal cells throughout cortex.

Interspike interval correlations in neuron models with adaptation and correlated noise

The generation of neural action potentials (spikes) is random but nevertheless may result in a rich statistical structure of the spike sequence. In particular, contrary to the popular renewal assumption of theoreticians, the intervals between adjacent spikes are often correlated. Experimentally, different patterns of interspike-interval correlations have been observed and computational studies have identified spike-frequency adaptation and correlated noise as the two main mechanisms that can lead to such correlations. Analytical studies have focused on the single cases of either correlated (colored) noise or adaptation currents in combination with uncorrelated (white) noise. For low-pass filtered noise or adaptation, the serial correlation coefficient can be approximated as a single geometric sequence of the lag between the intervals, providing an explanation for some of the experimentally observed patterns. Here we address the problem of interval correlations for a widely used class of models, multidimensional integrate-and-fire neurons subject to a combination of colored and white noise sources and a spike-triggered adaptation current. Assuming weak noise, we derive a simple formula for the serial correlation coefficient, a sum of two geometric sequences, which accounts for a large class of correlation patterns. The theory is confirmed by means of numerical simulations in a number of special cases including the leaky, quadratic, and generalized integrate-and-fire models with colored noise and spike-frequency adaptation. Furthermore we study the case in which the adaptation current and the colored noise share the same time scale, corresponding to a slow stochastic population of adaptation channels; we demonstrate that our theory can account for a nonmonotonic dependence of the correlation coefficient on the channel's time scale. Another application of the theory is a neuron driven by network-noise-like fluctuations (green noise). We also discuss the range of validity of our weak-noise theory and show that by changing the relative strength of white and colored noise sources, we can change the sign of the correlation coefficient. Finally, we apply our theory to a conductance-based model which demonstrates its broad applicability.

Noise effect on the temporal patterns of neural synchrony

Neural synchrony in the brain is often present in an intermittent fashion, i.e., there are intervals of synchronized activity interspersed with intervals of desynchronized activity. A series of experimental studies showed that this kind of temporal patterning of neural synchronization may be very specific and may be correlated with behaviour (even if the average synchrony strength is not changed). Prior studies showed that a network with many short desynchronized intervals may be functionally different from a network with few long desynchronized intervals as it may be more sensitive to synchronizing input signals. In this study, we investigated the effect of channel noise on the temporal patterns of neural synchronization. We employed a small network of conductance-based model neurons that were mutually connected via excitatory synapses. The resulting dynamics of the network was studied using the same time-series analysis methods as used in prior experimental and computational studies. While it is well known that synchrony strength generally degrades with noise, we found that noise also affects the temporal patterning of synchrony. Noise, at a sufficient intensity (yet too weak to substantially affect synchrony strength), promotes dynamics with predominantly short (although potentially very numerous) desynchronizations. Thus, channel noise may be one of the mechanisms contributing to the short desynchronization dynamics observed in multiple experimental studies.

Heterogeneities in intrinsic excitability and frequency-dependent response properties of granule cells across the blades of the rat dentate gyrus

The dentate gyrus (DG), the input gate to the hippocampus proper, is anatomically segregated into three different sectors, namely, the suprapyramidal blade, the crest region, and the infrapyramidal blade. Although there are well-established differences between these sectors in terms of neuronal morphology, connectivity patterns, and activity levels, differences in electrophysiological properties of granule cells within these sectors have remained unexplored. Here, employing somatic whole cell patch-clamp recordings from the rat DG, we demonstrate that granule cells in these sectors manifest considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, these neurons showed positive temporal summation of their responses to inputs mimicking excitatory postsynaptic currents and showed little to no sag in their voltage responses to pulse currents. Consistently, the impedance amplitude profile manifested low-pass characteristics and the impedance phase profile lacked positive phase values at all measured frequencies and voltages and for all sectors. Granule cells in all sectors exhibited class I excitability, with broadly linear firing rate profiles, and granule cells in the crest region fired significantly fewer action potentials compared with those in the infrapyramidal blade. Finally, we found weak pairwise correlations across the 18 different measurements obtained individually from each of the three sectors, providing evidence that these measurements are indeed reporting distinct aspects of neuronal physiology. Together, our analyses show that granule cells act as integrators of afferent information and emphasize the need to account for the considerable physiological heterogeneities in assessing their roles in information encoding and processing.NEW & NOTEWORTHY We employed whole cell patch-clamp recordings from granule cells in the three subregions of the rat dentate gyrus to demonstrate considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, granule cells did not express membrane potential resonance, and their impedance profiles lacked inductive phase leads at all measured frequencies. Our analyses also show that granule cells manifest class I excitability characteristics, categorizing them as integrators of afferent information.

Analysis of heterogeneous cardiac pacemaker tissue models and traveling wave dynamics

The sinoatrial-node (SAN) is a complex heterogeneous tissue that generates a stable rhythm in healthy hearts, yet a general mechanistic explanation for when and how this tissue remains stable is lacking. Although computational and theoretical analyses could elucidate these phenomena, such methods have rarely been used in realistic (large-dimensional) gap-junction coupled heterogeneous pacemaker tissue models. In this study, we adapt a recent model of pacemaker cells (Severi et al., 2012), incorporating biophysical representations of ion channel and intracellular calcium dynamics, to capture physiological features of a heterogeneous population of pacemaker cells, in particular "center" and "peripheral" cells with distinct intrinsic frequencies and action potential morphology. Large-scale simulations of the SAN tissue, represented by a heterogeneous tissue structure of pacemaker cells, exhibit a rich repertoire of behaviors, including complete synchrony, traveling waves of activity originating from periphery to center, and transient traveling waves originating from the center. We use phase reduction methods that do not require fully simulating the large-scale model to capture these observations. Moreover, the phase reduced models accurately predict key properties of the tissue electrical dynamics, including wave frequencies when synchronization occurs, and wave propagation direction in a variety of tissue models. With the reduced phase models, we analyze the relationship between cell distributions and coupling strengths and the resulting transient dynamics. Further, the reduced phase model predicts parameter regimes of irregular electrical dynamics. Thus, we demonstrate that phase reduced oscillator models applied to realistic pacemaker tissue is a useful tool for investigating the spatial-temporal dynamics of cardiac pacemaker activity.

Degeneracy in the robust expression of spectral selectivity, subthreshold oscillations, and intrinsic excitability of entorhinal stellate cells

Biological heterogeneities are ubiquitous and play critical roles in the emergence of physiology at multiple scales. Although neurons in layer II (LII) of the medial entorhinal cortex (MEC) express heterogeneities in channel properties, the impact of such heterogeneities on the robustness of their cellular-scale physiology has not been assessed. Here, we performed a 55-parameter stochastic search spanning nine voltage- or calcium-activated channels to assess the impact of channel heterogeneities on the concomitant emergence of 10 in vitro electrophysiological characteristics of LII stellate cells (SCs). We generated 150,000 models and found a heterogeneous subpopulation of 449 valid models to robustly match all electrophysiological signatures. We employed this heterogeneous population to demonstrate the emergence of cellular-scale degeneracy in SCs, whereby disparate parametric combinations expressing weak pairwise correlations resulted in similar models. We then assessed the impact of virtually knocking out each channel from all valid models and demonstrate that the mapping between channels and measurements was many-to-many, a critical requirement for the expression of degeneracy. Finally, we quantitatively predict that the spike-triggered average of SCs should be endowed with theta-frequency spectral selectivity and coincidence detection capabilities in the fast gamma-band. We postulate this fast gamma-band coincidence detection as an instance of cellular-scale-efficient coding, whereby SC response characteristics match the dominant oscillatory signals in LII MEC. The heterogeneous population of valid SC models built here unveils the robust emergence of cellular-scale physiology despite significant channel heterogeneities, and forms an efficacious substrate for evaluating the impact of biological heterogeneities on entorhinal network function.

NEW & NOTEWORTHY We assessed the impact of heterogeneities in channel properties on the robustness of cellular-scale physiology of medial entorhinal cortical stellate neurons. We demonstrate that neuronal models with disparate channel combinations were endowed with similar physiological characteristics, as a consequence of the many-to-many mapping between channel properties and the physiological characteristics that they modulate. We predict that the spike-triggered average of stellate cells should be endowed with theta-frequency spectral selectivity and fast gamma-band coincidence detection capabilities.

A quantitative theory of gamma synchronization in macaque V1

Gamma-band synchronization coordinates brief periods of excitability in oscillating neuronal populations to optimize information transmission during sensation and cognition. Commonly, a stable, shared frequency over time is considered a condition for functional neural synchronization. Here, we demonstrate the opposite: instantaneous frequency modulations are critical to regulate phase relations and synchronization. In monkey visual area V1, nearby local populations driven by different visual stimulation showed different gamma frequencies. When similar enough, these frequencies continually attracted and repulsed each other, which enabled preferred phase relations to be maintained in periods of minimized frequency difference. Crucially, the precise dynamics of frequencies and phases across a wide range of stimulus conditions was predicted from a physics theory that describes how weakly coupled oscillators influence each other's phase relations. Hence, the fundamental mathematical principle of synchronization through instantaneous frequency modulations applies to gamma in V1 and is likely generalizable to other brain regions and rhythms.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

Optimization of noise-induced synchronization of oscillator networks

We investigate common-noise-induced synchronization between two identical networks of coupled phase oscillators exhibiting fully locked collective oscillations. Using the collective phase description method for fully locked oscillators, we demonstrate that two noninteracting networks of coupled phase oscillators can exhibit in-phase synchronization between the networks when driven by weak common noise. We derive the Lyapunov exponent characterizing the relaxation time for synchronization and develop a method of obtaining the optimal input pattern of common noise to achieve fast synchronization. We illustrate the theory using three representative networks with heterogeneous, global, and local coupling. The theoretical results are validated by direct numerical simulations.

Criterion for noise-induced synchronization: Application to colloidal alignment

Colloidal bodies of irregular shape rotate as they descend under gravity in solution. This rotational response provides a means of bringing a dispersion of identical bodies into a synchronized rotation with the same orientation using programed forcing. We use the notion of statistical entropy to derive bounds on the rate of synchronization. These bounds apply generally to dynamical systems with stable periodic motion with a phase ϕ(t), when subjected to an impulsive perturbation. The impulse causes a change of phase expressible as a phase map ψ(ϕ). We derive an upper limit on the average change of entropy 〈ΔH〉 in terms of this phase map; when this limit is negative, alignment must occur. For systems that have achieved a low entropy, the 〈ΔH〉 approaches this upper limit.

The mechanics of state-dependent neural correlations

Simultaneous recordings from large neural populations are becoming increasingly common. An important feature of population activity is the trial-to-trial correlated fluctuation of spike train outputs from recorded neuron pairs. Similar to the firing rate of single neurons, correlated activity can be modulated by a number of factors, from changes in arousal and attentional state to learning and task engagement. However, the physiological mechanisms that underlie these changes are not fully understood. We review recent theoretical results that identify three separate mechanisms that modulate spike train correlations: changes in input correlations, internal fluctuations and the transfer function of single neurons. We first examine these mechanisms in feedforward pathways and then show how the same approach can explain the modulation of correlations in recurrent networks. Such mechanistic constraints on the modulation of population activity will be important in statistical analyses of high-dimensional neural data.

Stochastic Synchronization in Purkinje Cells with Feedforward Inhibition Could Be Studied with Equivalent Phase-Response Curves

Simple-spike synchrony between Purkinje cells projecting to a common neuron in the deep cerebellar nucleus is emerging as an important factor in the encoding of output information from cerebellar cortex. A phenomenon known as stochastic synchronization happens when uncoupled oscillators synchronize due to correlated inputs. Stochastic synchronization is a viable mechanism through which simple-spike synchrony could be generated, but it has received scarce attention, perhaps because the presence of feedforward inhibition in the input to Purkinje cells makes insights difficult. This paper presents a method to account for feedforward inhibition so the usual mathematical approaches to stochastic synchronization can be applied. The method consists in finding a single Phase Response Curve, called the equivalent PRC, that accounts for the effects of both excitatory inputs and delayed feedforward inhibition from molecular layer interneurons. The results suggest that a theory of stochastic synchronization for the case of feedforward inhibition may not be necessary, since this case can be approximately reduced to the case of inputs characterized by a single PRC. Moreover, feedforward inhibition could in many situations increase the level of synchrony experienced by Purkinje cells.

Phase-resetting as a tool of information transmission

Models of information transmission in the brain largely rely on firing rate codes. The abundance of oscillatory activity in the brain suggests that information may be also encoded using the phases of ongoing oscillations. Sensory perception, working memory and spatial navigation have been hypothesized to use phase codes, and cross-frequency coordination and phase synchronization between brain areas have been proposed to gate the flow of information. Phase codes generally require the phase of the oscillations to be reset at specific reference points for consistent coding, and coordination between oscillators requires favorable phase resetting characteristics. Recent evidence supports a role for neural oscillations in providing temporal reference windows that allow for correct parsing of phase-coded information.

Examination of rhythmicity of extracellularly recorded neurons in the entorhinal cortex

A number of studies have examined the theta-rhythmic modulation of neuronal firing in the hippocampal circuit. For extracellular recordings, this is often done by examining spectral properties of the spike-time autocorrelogram, most significantly, for validating the presence or absence of theta modulation across species. These techniques can show significant rhythmicity for high firing rate, highly rhythmic neurons; however, they are substantially biased by several factors including the peak firing rate of the neuron, the amount of time spent in the neuron's receptive field, and other temporal properties of the rhythmicity such as cycle-skipping. These limitations make it difficult to examine rhythmic modulation in neurons with low firing rates or when an animal has short dwell times within the firing field and difficult to compare rhythmicity under disparate experimental conditions when these factors frequently differ. Here, we describe in detail the challenges that researchers face when using these techniques and apply our findings to recent recordings from bat entorhinal grid cells, suggesting that they may have lacked enough data to examine theta rhythmicity robustly. We describe a more sensitive and statistically rigorous method using maximum likelihood estimation (MLE) of a parametric model of the lags within the autocorrelation window, which helps to alleviate some of the problems of traditional methods and was also unable to detect rhythmicity in bat grid cells. Using large batteries of simulated data, we explored the boundaries for which the MLE technique and the theta index can detect rhythmicity. The MLE technique is less sensitive to many features of the autocorrelogram and provides a framework for statistical testing to detect rhythmicity as well as changes in rhythmicity in individual sessions providing a substantial improvement over previous methods.