Transition between Functional Regimes in an Integrate-And-Fire Network Model of the Thalamus

Mal de Debarquement Syndrome explained by a vestibulo-cerebellar oscillator

Mal de Debarquement Syndrome (MdDS) is a puzzling central vestibular disorder characterized by a long-lasting perception of oscillatory postural instability that may occur after sea travels or flights. We have postulated that MdDS originates from the post-disembarking persistence of an adaptive internal oscillator consisting of a loop system, involving the right and left vestibular nuclei, and the Purkinje cells of the right and left flocculonodular cerebellar cortex, connected by GABAergic and glutamatergic fibers. We have formulated here a mathematical model of the vestibulo-cerebellar loop system and carried out a computational analysis based on a set of differential equations describing the interactions among the loop elements and containing Hill functions that model input-output firing rates relationships among neurons. The analysis indicates that the system acquires a spontaneous and permanent oscillatory behavior for a decrease of threshold and an increase of sensitivity in neuronal input-output responses. These results suggest a role for synaptic plasticity in MdDS pathophysiology, thus reinforcing our previous hypothesis that MdDS may be the result of excessive synaptic plasticity acting on the vestibulo-cerebellar network during its entraining to an oscillatory environment. Hence, our study points to neuroendocrine pathways that lead to increased synaptic response as possible new therapeutic targets for the clinical treatment of the disorder.

Thalamocortical bistable switch as a theoretical model of fibromyalgia pathogenesis inferred from a literature survey

Fibromyalgia (FM) is an unsolved central pain processing disturbance. We aim to provide a unifying model for FM pathogenesis based on a loop network involving thalamocortical regions, i.e., the ventroposterior lateral thalamus (VPL), the somatosensory cortex (SC), and the thalamic reticular nucleus (TRN). The dynamics of the loop have been described by three differential equations having neuron mean firing rates as variables and containing Hill functions to model mutual interactions among the loop elements. A computational analysis conducted with MATLAB has shown a transition from monostability to bistability of the loop behavior for a weakening of GABAergic transmission between TRN and VPL. This involves the appearance of a high-firing-rate steady state, which becomes dominant and is assumed to represent pathogenic pain processing giving rise to chronic pain. Our model is consistent with a bulk of literature evidence, such as neuroimaging and pharmacological data collected on FM patients, and with correlations between FM and immunoendocrine conditions, such as stress, perimenopause, chronic inflammation, obesity, and chronic dizziness. The model suggests that critical targets for FM treatment are to be found among immunoendocrine pathways leading to GABA/glutamate imbalance having an impact on the thalamocortical system.

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.

Computational Models in Electroencephalography

Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses in silico and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by “computational model” is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses in silico. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.

Cortical Responses to Input From Distant Areas are Modulated by Local Spontaneous Alpha/Beta Oscillations

Any given area in human cortex may receive input from multiple, functionally heterogeneous areas, potentially representing different processing threads. Alpha (8-13 Hz) and beta oscillations (13-20 Hz) have been hypothesized by other investigators to gate local cortical processing, but their influence on cortical responses to input from other cortical areas is unknown. To study this, we measured the effect of local oscillatory power and phase on cortical responses elicited by single-pulse electrical stimulation (SPES) at distant cortical sites, in awake human subjects implanted with intracranial electrodes for epilepsy surgery. In 4 out of 5 subjects, the amplitudes of corticocortical evoked potentials (CCEPs) elicited by distant SPES were reproducibly modulated by the power, but not the phase, of local oscillations in alpha and beta frequencies. Specifically, CCEP amplitudes were higher when average oscillatory power just before distant SPES (-110 to -10 ms) was high. This effect was observed in only a subset (0-33%) of sites with CCEPs and, like the CCEPs themselves, varied with stimulation at different distant sites. Our results suggest that although alpha and beta oscillations may gate local processing, they may also enhance the responsiveness of cortex to input from distant cortical sites.

Integrate-and-fire network model of activity propagation from thalamus to cortex

The thalamus plays a crucial role in modulating the cortical activity underlying sensory and cognitive processes. In particular, recent experimental findings highlighted that the thalamus does not merely act as a binary gate for sensory stimuli, but rather participates to the processing of sensory information. Clarifying such thalamic influence on cortical dynamics is also important as the thalamus is the target of therapies such as DBS for Tourette patients. In this perspective, various computational models have been proposed in the last decades. However, a detailed description of the propagation of thalamic activity to the cortex is missing. Here we present a simple computational model of thalamocortical connectivity accounting for the propagation of activity from the thalamus to the cortex. The model includes both the single-neuron scale and the mesoscopic level of Local Field Potential (LFP) signals. Numerical simulations at both levels reproduce typical thalamocortical dynamics which are consistent with experimental measurements and robust to parameters changes. In particular, our model correctly reproduces locally generated rhythms as spindle oscillations in the thalamus and gamma oscillations in the cortex. Our model paves the way to deeper investigations of the thalamic influence on cortical dynamics, with and without sensory inputs or therapeutic electrical stimulation.

A Proposed Mechanism for Spontaneous Transitions between Interictal and Ictal Activity

Epileptic networks are characterized by two outputs: brief interictal spikes and rarer, more prolonged seizures. Although either output state is readily modeled in silico and induced experimentally, the transition mechanisms are unknown, in part because no models exhibit both output states spontaneously. In silico small-world neural networks were built using single-compartment neurons whose physiological parameters were derived from dual whole-cell recordings of pyramidal cells in organotypic hippocampal slice cultures that were generating spontaneous seizure-like activity. In silico, neurons were connected by abundant local synapses and rare long-distance synapses. Activity-dependent synaptic depression and gradual recovery delimited synchronous activity. Full synaptic recovery engendered interictal population spikes that spread via long-distance synapses. When synaptic recovery was incomplete, postsynaptic neurons required coincident activation of multiple presynaptic terminals to reach firing threshold. Only local connections were sufficiently dense to spread activity under these conditions. This coalesced network activity into traveling waves whose velocity varied with synaptic recovery. Seizures were comprised of sustained traveling waves that were similar to those recorded during experimental and human neocortical seizures. Sustained traveling waves occurred only when wave velocity, network dimensions, and the rate of synaptic recovery enabled wave reentry into previously depressed areas at precisely ictogenic levels of synaptic recovery. Wide-field, cellular-resolution GCamP7b calcium imaging demonstrated similar initial patterns of activation in the hippocampus, although the anatomical distribution of traveling waves of synaptic activation was altered by the pattern of synaptic connectivity in the organotypic hippocampal cultures.SIGNIFICANCE STATEMENT When computerized distributed neural network models are required to generate both features of epileptic networks (i.e., spontaneous interictal population spikes and seizures), the network structure is substantially constrained. These constraints provide important new hypotheses regarding the nature of epileptic networks and mechanisms of seizure onset.

Firing-rate based network modeling of the dLGN circuit: Effects of cortical feedback on spatiotemporal response properties of relay cells

Visually evoked signals in the retina pass through the dorsal geniculate nucleus (dLGN) on the way to the visual cortex. This is however not a simple feedforward flow of information: there is a significant feedback from cortical cells back to both relay cells and interneurons in the dLGN. Despite four decades of experimental and theoretical studies, the functional role of this feedback is still debated. Here we use a firing-rate model, the extended difference-of-Gaussians (eDOG) model, to explore cortical feedback effects on visual responses of dLGN relay cells. For this model the responses are found by direct evaluation of two- or three-dimensional integrals allowing for fast and comprehensive studies of putative effects of different candidate organizations of the cortical feedback. Our analysis identifies a special mixed configuration of excitatory and inhibitory cortical feedback which seems to best account for available experimental data. This configuration consists of (i) a slow (long-delay) and spatially widespread inhibitory feedback, combined with (ii) a fast (short-delayed) and spatially narrow excitatory feedback, where (iii) the excitatory/inhibitory ON-ON connections are accompanied respectively by inhibitory/excitatory OFF-ON connections, i.e. following a phase-reversed arrangement. The recent development of optogenetic and pharmacogenetic methods has provided new tools for more precise manipulation and investigation of the thalamocortical circuit, in particular for mice. Such data will expectedly allow the eDOG model to be better constrained by data from specific animal model systems than has been possible until now for cat. We have therefore made the Python tool pyLGN which allows for easy adaptation of the eDOG model to new situations.

On route from the retina to primary visual cortex, visually evoked signals have to pass through the dorsal lateral geniculate nucleus (dLGN). However, this is not an exclusive feedforward flow of information as feedback exists from neurons in the cortex back to both relay cells and interneurons in the dLGN. The functional role of this feedback remains mostly unresolved. Here, we use a firing-rate model, the extended difference-of-Gaussians (eDOG) model, to explore cortical feedback effects on visual responses of dLGN relay cells. Our analysis indicates that a particular mix of excitatory and inhibitory cortical feedback agrees best with available experimental observations. In this configuration ON-center relay cells receive both excitatory and (indirect) inhibitory feedback from ON-center cortical cells (ON-ON feedback) where the excitatory feedback is fast and spatially narrow while the inhibitory feedback is slow and spatially widespread. In addition to the ON-ON feedback, the connections are accompanied by OFF-ON connections following a so-called phase-reversed (push-pull) arrangement. To facilitate further applications of the model, we have made the Python tool pyLGN which allows for easy modification and evaluation of the a priori quite general eDOG model to new situations.