The role of cortical feedback in the generation of the temporal receptive field responses of lateral geniculate nucleus neurons: a computational modelling study

A systematic analysis of the joint effects of ganglion cells, lagged LGN cells, and intercortical inhibition on spatiotemporal processing and direction selectivity

Simple cells in the visual cortex process spatial as well as temporal information of the visual stream and enable the perception of motion information. Previous work suggests different mechanisms associated with direction selectivity, such as a temporal offset in thalamocortical input stream through lagged and non-lagged cells of the lateral geniculate nucleus (LGN), or solely from intercortical inhibition, or through a baseline selectivity provided by the thalamocortical connection tuned by intercortical inhibition. While there exists a large corpus of models for spatiotemporal receptive fields, the majority of them built-in the spatiotemporal dynamics by utilizing a combination of spatial and temporal functions and thus, do not explain the emergence of spatiotemporal dynamics on basis of network dynamics emerging in the retina and the LGN. In order to better comprehend the emergence of spatiotemporal processing and direction selectivity, we used a spiking neural network to implement the visual pathway from the retina to the primary visual cortex. By varying different functional parts in our network, we demonstrate how the direction selectivity of simple cells emerges through the interplay between two components: tuned intercortical inhibition and a temporal offset in the feedforward path through lagged LGN cells. In contrast to previous findings, our model simulations suggest an alternative dynamic between these two mechanisms: While intercortical inhibition alone leads to bidirectional selectivity, a temporal shift in the thalamocortical pathway breaks this symmetry in favor of one direction, leading to unidirectional selectivity.

Resonating neurons stabilize heterogeneous grid-cell networks

A central theme that governs the functional design of biological networks is their ability to sustain stable function despite widespread parametric variability. Here, we investigated the impact of distinct forms of biological heterogeneities on the stability of a two-dimensional continuous attractor network (CAN) implicated in grid-patterned activity generation. We show that increasing degrees of biological heterogeneities progressively disrupted the emergence of grid-patterned activity and resulted in progressively large perturbations in low-frequency neural activity. We postulated that targeted suppression of low-frequency perturbations could ameliorate heterogeneity-induced disruptions of grid-patterned activity. To test this, we introduced intrinsic resonance, a physiological mechanism to suppress low-frequency activity, either by adding an additional high-pass filter (phenomenological) or by incorporating a slow negative feedback loop (mechanistic) into our model neurons. Strikingly, CAN models with resonating neurons were resilient to the incorporation of heterogeneities and exhibited stable grid-patterned firing. We found CAN models with mechanistic resonators to be more effective in targeted suppression of low-frequency activity, with the slow kinetics of the negative feedback loop essential in stabilizing these networks. As low-frequency perturbations (1/f noise) are pervasive across biological systems, our analyses suggest a universal role for mechanisms that suppress low-frequency activity in stabilizing heterogeneous biological networks.

Rebound or Entrainment? The Influence of Alternating Current Stimulation on Individual Alpha

Alternating current stimulation (ACS) is an established means to manipulate intrinsic cortical oscillations. While working towards clinical impact, ACS mechanisms of action remain unclear. For ACS’s well-documented influence on occipital alpha, hypotheses include neuronal entrainment as well as rebound phenomena. As a retinal origin is also discussed, we employed a novel form of ACS with the advantage that it specifically targets occipital alpha-oscillations via retinofugal pathways retinofugal ACS (rACS). We aimed to confirm alpha-enhancement outlasting the duration of stimulation with 10 Hz rACS. To distinguish entrainment from rebound effects, we investigated the correlation between alpha peak frequency change and alpha-enhancement strength. We quantified the alpha band power before and after 10 Hz rACS in 15 healthy subjects. Alpha power enhancement and alpha peak frequency change were assessed over the occipital electrodes and compared to sham stimulation. RACS significantly enhanced occipital alpha power in comparison to sham stimulation (p < 0.05). Alpha peak frequency changed by a mean 0.02 Hz (± 0.04). A greater change in alpha peak frequency did not correlate with greater effects on alpha power. Our findings show an alpha-enhancement consistent with studies conducted for transcranial ACS (tACS) and contribute evidence for a retinal involvement in tACS effects on occipital alpha. Furthermore, the lack of correlation between alpha peak frequency change and alpha-enhancement strength provides an argument against entrainment effects and in favor of a rebound phenomenon.

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.

Biophysical network modeling of the dLGN circuit: Effects of cortical feedback on spatial response properties of relay cells

Despite half-a-century of research since the seminal work of Hubel and Wiesel, the role of the dorsal lateral geniculate nucleus (dLGN) in shaping the visual signals is not properly understood. Placed on route from retina to primary visual cortex in the early visual pathway, a striking feature of the dLGN circuit is that both the relay cells (RCs) and interneurons (INs) not only receive feedforward input from retinal ganglion cells, but also a prominent feedback from cells in layer 6 of visual cortex. This feedback has been proposed to affect synchronicity and other temporal properties of the RC firing. It has also been seen to affect spatial properties such as the center-surround antagonism of thalamic receptive fields, i.e., the suppression of the response to very large stimuli compared to smaller, more optimal stimuli. Here we explore the spatial effects of cortical feedback on the RC response by means of a a comprehensive network model with biophysically detailed, single-compartment and multicompartment neuron models of RCs, INs and a population of orientation-selective layer 6 simple cells, consisting of pyramidal cells (PY). We have considered two different arrangements of synaptic feedback from the ON and OFF zones in the visual cortex to the dLGN: phase-reversed (‘push-pull’) and phase-matched (‘push-push’), as well as different spatial extents of the corticothalamic projection pattern. Our simulation results support that a phase-reversed arrangement provides a more effective way for cortical feedback to provide the increased center-surround antagonism seen in experiments both for flashing spots and, even more prominently, for patch gratings. This implies that ON-center RCs receive direct excitation from OFF-dominated cortical cells and indirect inhibitory feedback from ON-dominated cortical cells. The increased center-surround antagonism in the model is accompanied by spatial focusing, i.e., the maximum RC response occurs for smaller stimuli when feedback is present.

The functional role of the dorsal lateral geniculate nucleus (dLGN), placed on route from retina to primary visual cortex in the early visual pathway, is still poorly understood. A striking feature of the dLGN circuit is that dLGN cells not only receive feedforward input from the retina, but also a prominent feedback from cells in the visual cortex. It has been seen in experiments that cortical feedback modifies the spatial properties of dLGN cells in response to visual stimuli. In particular, it has been shown to increase the center-surround antagonism for flashing-spot and patch-grating visual stimuli, i.e., the suppression of responses to very large stimuli compared to smaller stimuli. Here we investigate the putative mechanisms behind this feature by means of a comprehensive network model of biophysically detailed neuron models for RCs and INs in the dLGN and orientation-selective cortical cells providing the feedback. Our results support that the experimentally observed feedback effects may be due to a phase-reversed (‘push-pull’) arrangement of the cortical feedback where ON-symmetry RCs receive (indirect) inhibitory feedback from ON-dominated cortical cell and excitation from OFF-dominated cortical cells, and vice versa for OFF-symmetry RCs.

A cross-species comparison of corticogeniculate structure and function

The corticogeniculate circuit is an evolutionarily conserved pathway linking the primary visual cortex with the visual thalamus in the feedback direction. While the corticogeniculate circuit is anatomically robust, the impact of corticogeniculate feedback on the visual response properties of visual thalamic neurons is subtle. Accordingly, discovering the function of corticogeniculate feedback in vision has been a particularly challenging task. In this review, the morphology, organization, physiology, and function of corticogeniculate feedback is compared across mammals commonly studied in visual neuroscience: primates, carnivores, rabbits, and rodents. Common structural and organizational motifs are present across species, including the organization of corticogeniculate feedback into parallel processing streams in highly visual mammals.

Corticogeniculate feedback sharpens the temporal precision and spatial resolution of visual signals in the ferret

The corticogeniculate (CG) pathway connects the visual cortex with the visual thalamus (LGN) in the feedback direction and enables the cortex to directly influence its own input. Despite numerous investigations, the role of this feedback circuit in visual perception remained elusive. To probe the function of CG feedback in a causal manner, we selectively and reversibly manipulated the activity of CG neurons in anesthetized ferrets in vivo using a combined viral-infection and optogenetics approach to drive expression of channelrhodopsin2 (ChR2) in CG neurons. We observed significant increases in temporal precision and spatial resolution of LGN neuronal responses to drifting grating and white noise stimuli when CG neurons expressing ChR2 were light activated. Enhancing CG feedback reduced visually evoked response latencies, increased spike-timing precision, and reduced classical receptive field size. Increased precision among LGN neurons led to increased spike-timing precision among granular layer V1 neurons as well. Together, our findings suggest that the function of CG feedback is to control the timing and precision of thalamic responses to incoming visual signals.

A Neurophysiological Perspective on a Preventive Treatment against Schizophrenia Using Transcranial Electric Stimulation of the Corticothalamic Pathway

Schizophrenia patients are waiting for a treatment free of detrimental effects. Psychotic disorders are devastating mental illnesses associated with dysfunctional brain networks. Ongoing brain network gamma frequency (30–80 Hz) oscillations, naturally implicated in integrative function, are excessively amplified during hallucinations, in at-risk mental states for psychosis and first-episode psychosis. So, gamma oscillations represent a bioelectrical marker for cerebral network disorders with prognostic and therapeutic potential. They accompany sensorimotor and cognitive deficits already present in prodromal schizophrenia. Abnormally amplified gamma oscillations are reproduced in the corticothalamic systems of healthy humans and rodents after a single systemic administration, at a psychotomimetic dose, of the glutamate N-methyl-d-aspartate receptor antagonist ketamine. These translational ketamine models of prodromal schizophrenia are thus promising to work out a preventive noninvasive treatment against first-episode psychosis and chronic schizophrenia. In the present essay, transcranial electric stimulation (TES) is considered an appropriate preventive therapeutic modality because it can influence cognitive performance and neural oscillations. Here, I highlight clinical and experimental findings showing that, together, the corticothalamic pathway, the thalamus, and the glutamatergic synaptic transmission form an etiopathophysiological backbone for schizophrenia and represent a potential therapeutic target for preventive TES of dysfunctional brain networks in at-risk mental state patients against psychotic disorders.

A Network Model of Local Field Potential Activity in Essential Tremor and the Impact of Deep Brain Stimulation

Essential tremor (ET), a movement disorder characterised by an uncontrollable shaking of the affected body part, is often professed to be the most common movement disorder, affecting up to one percent of adults over 40 years of age. The precise cause of ET is unknown, however pathological oscillations of a network of a number of brain regions are implicated in leading to the disorder. Deep brain stimulation (DBS) is a clinical therapy used to alleviate the symptoms of a number of movement disorders. DBS involves the surgical implantation of electrodes into specific nuclei in the brain. For ET the targeted region is the ventralis intermedius (Vim) nucleus of the thalamus. Though DBS is effective for treating ET, the mechanism through which the therapeutic effect is obtained is not understood. To elucidate the mechanism underlying the pathological network activity and the effect of DBS on such activity, we take a computational modelling approach combined with electrophysiological data. The pathological brain activity was recorded intra-operatively via implanted DBS electrodes, whilst simultaneously recording muscle activity of the affected limbs. We modelled the network hypothesised to underlie ET using the Wilson-Cowan approach. The modelled network exhibited oscillatory behaviour within the tremor frequency range, as did our electrophysiological data. By applying a DBS-like input we suppressed these oscillations. This study shows that the dynamics of the ET network support oscillations at the tremor frequency and the application of a DBS-like input disrupts this activity, which could be one mechanism underlying the therapeutic benefit.

Essential tremor (ET) is acknowledged to be the most common movement disorder affecting 1% of the population. Although the underlying mechanisms remain elusive, the thalamus, cortex and cerebellum are implicated in the underlying pathology. More recently, it has been shown that ET can be successfully treated by deep brain stimulation (DBS). This clinical treatment involves the surgical implantation of electrodes into the brain, through which current is applied. However, the mechanisms of how DBS achieves clinical benefit continue to be debated. A key question is whether ET can be modeled as a pathological network behavior as has been suggested previously. If so, we can then ask how DBS would modulate this brain activity. Our study combines: (i) simultaneous electrophysiological recordings from the brain and muscle; (ii) computational modelling; (iii) mathematical analysis. We found that the network supports oscillations in the tremor range, and the application of high frequency DBS switches this to low amplitude, high-frequency activity. We propose that our model can be used to predict DBS parameter settings that suppress pathological network activity and consequently tremor. In summary, we provide the first population level model of essential tremor including the effect of DBS on network behaviour.

Extended difference-of-Gaussians model incorporating cortical feedback for relay cells in the lateral geniculate nucleus of cat

A striking feature of the organization of the early visual pathway is the significant feedback from primary visual cortex to cells in the dorsal lateral geniculate nucleus (LGN). Despite numerous experimental and modeling studies, the functional role for this feedback remains elusive. We present a new firing-rate-based model for LGN relay cells in cat, explicitly accounting for thalamocortical loop effects. The established DOG model, here assumed to account for the spatial aspects of the feedforward processing of visual stimuli, is extended to incorporate the influence of thalamocortical loops including a full set of orientation-selective cortical cell populations. Assuming a phase-reversed push-pull arrangement of ON and OFF cortical feedback as seen experimentally, this extended DOG (eDOG) model exhibits linear firing properties despite non-linear firing characteristics of the corticothalamic cells. The spatiotemporal receptive field of the eDOG model has a simple algebraic structure in Fourier space, while the real-space receptive field, as well as responses to visual stimuli, are found by evaluation of an integral. As an example application we use the eDOG model to study effects of cortical feedback on responses to flashing circular spots and patch-grating stimuli and find that the eDOG model can qualitatively account for experimental findings.

A minimal mechanistic model for temporal signal processing in the lateral geniculate nucleus

The receptive fields of cells in the lateral geniculate nucleus (LGN) are shaped by their diverse set of impinging inputs: feedforward synaptic inputs stemming from retina, and feedback inputs stemming from the visual cortex and the thalamic reticular nucleus. To probe the possible roles of these feedforward and feedback inputs in shaping the temporal receptive-field structure of LGN relay cells, we here present and investigate a minimal mechanistic firing-rate model tailored to elucidate their disparate features. The model for LGN relay ON cells includes feedforward excitation and inhibition (via interneurons) from retinal ON cells and excitatory and inhibitory (via thalamic reticular nucleus cells and interneurons) feedback from cortical ON and OFF cells. From a general firing-rate model formulated in terms of Volterra integral equations, we derive a single delay differential equation with absolute delay governing the dynamics of the system. A freely available and easy-to-use GUI-based MATLAB version of this minimal mechanistic LGN circuit model is provided. We particularly investigate the LGN relay-cell impulse response and find through thorough explorations of the model's parameter space that both purely feedforward models and feedback models with feedforward excitation only, can account quantitatively for previously reported experimental results. We find, however, that the purely feedforward model predicts two impulse response measures, the time to first peak and the biphasic index (measuring the relative weight of the rebound phase) to be anticorrelated. In contrast, the models with feedback predict different correlations between these two measures. This suggests an experimental test assessing the relative importance of feedforward and feedback connections in shaping the impulse response of LGN relay cells.

Predictive feedback can account for biphasic responses in the lateral geniculate nucleus

Biphasic neural response properties, where the optimal stimulus for driving a neural response changes from one stimulus pattern to the opposite stimulus pattern over short periods of time, have been described in several visual areas, including lateral geniculate nucleus (LGN), primary visual cortex (V1), and middle temporal area (MT). We describe a hierarchical model of predictive coding and simulations that capture these temporal variations in neuronal response properties. We focus on the LGN-V1 circuit and find that after training on natural images the model exhibits the brain's LGN-V1 connectivity structure, in which the structure of V1 receptive fields is linked to the spatial alignment and properties of center-surround cells in the LGN. In addition, the spatio-temporal response profile of LGN model neurons is biphasic in structure, resembling the biphasic response structure of neurons in cat LGN. Moreover, the model displays a specific pattern of influence of feedback, where LGN receptive fields that are aligned over a simple cell receptive field zone of the same polarity decrease their responses while neurons of opposite polarity increase their responses with feedback. This phase-reversed pattern of influence was recently observed in neurophysiology. These results corroborate the idea that predictive feedback is a general coding strategy in the brain.

For many neurons in the early visual brain the optimal stimulation for driving a response changes from one stimulus pattern to the opposite stimulus pattern over short periods of time. For example, many neurons in the lateral geniculate nucleus (LGN) respond to a bright stimulus initially but prefer a dark stimulus only 20 milliseconds later in time, and similar changes in response preference have been found for neurons in other areas. What would be the computational reason for these biphasic response dynamics? We describe a hierarchical model of predictive coding that explains these response properties. In the model, higher-level neurons attempt to predict their lower-level input, while lower-level neurons signal the difference between actual input and the higher-level predictions. In our simulations we focus on the LGN and area V1 and find that after training on natural images the layout of model connections resembles the brain's LGN-V1 connectivity structure. In addition, the responses of model LGN neurons are biphasic in time, resembling biphasic responses as found in neurophysiology. Moreover, the model displays a specific pattern of influence of feedback from higher-level areas that was recently observed in neurophysiology. These results corroborate the idea that predictive feedback is a general coding strategy in the brain.