Combining aperiodic 1/f slopes and brain simulation: An EEG/MEG proxy marker of excitation/inhibition imbalance in Alzheimer's disease.

INTRODUCTION: Accumulation and interaction of amyloid-beta (Aß) and tau proteins during progression of Alzheimer's disease (AD) are shown to tilt neuronal circuits away from balanced excitation/inhibition (E/I). Current available techniques for noninvasive interrogation of E/I in the intact human brain, for example, magnetic resonance spectroscopy (MRS), are highly restrictive (i.e., limited spatial extent), have low temporal and spatial resolution and suffer from the limited ability to distinguish accurately between different neurotransmitters complicating its interpretation. As such, these methods alone offer an incomplete explanation of E/I. Recently, the aperiodic component of neural power spectrum, often referred to in the literature as the '1/f slope', has been described as a promising and scalable biomarker that can track disruptions in E/I potentially underlying a spectrum of clinical conditions, such as autism, schizophrenia, or epilepsy, as well as developmental E/I changes as seen in aging. METHODS: Using 1/f slopes from resting-state spectral data and computational modeling, we developed a new method for inferring E/I alterations in AD. RESULTS: We tested our method on recent freely and publicly available electroencephalography (EEG) and magnetoencephalography (MEG) datasets of patients with AD or prodromal disease and demonstrated the method's potential for uncovering regional patterns of abnormal excitatory and inhibitory parameters. DISCUSSION: Our results provide a general framework for investigating circuit-level disorders in AD and developing therapeutic interventions that aim to restore the balance between excitation and inhibition.

Methods for inferring neural circuit interactions and neuromodulation from local field potential and electroencephalogram measures.

Electrical recordings of neural mass activity, such as local field potentials (LFPs) and electroencephalograms (EEGs), have been instrumental in studying brain function. However, these aggregate signals lack cellular resolution and thus are not easy to be interpreted directly in terms of parameters of neural microcircuits. Developing tools for a reliable estimation of key neural parameters from these signals, such as the interaction between excitation and inhibition or the level of neuromodulation, is important for both neuroscientific and clinical applications. Over the years, we have developed tools based on neural network modeling and computational analysis of empirical data to estimate neural parameters from aggregate neural signals. This review article gives an overview of the main computational tools that we have developed and employed to invert LFPs and EEGs in terms of circuit-level neural phenomena, and outlines future challenges and directions for future research.

Computation of the electroencephalogram (EEG) from network models of point neurons.

The electroencephalogram (EEG) is a major tool for non-invasively studying brain function and dysfunction. Comparing experimentally recorded EEGs with neural network models is important to better interpret EEGs in terms of neural mechanisms. Most current neural network models use networks of simple point neurons. They capture important properties of cortical dynamics, and are numerically or analytically tractable. However, point neurons cannot generate an EEG, as EEG generation requires spatially separated transmembrane currents. Here, we explored how to compute an accurate approximation of a rodent's EEG with quantities defined in point-neuron network models. We constructed different approximations (or proxies) of the EEG signal that can be computed from networks of leaky integrate-and-fire (LIF) point neurons, such as firing rates, membrane potentials, and combinations of synaptic currents. We then evaluated how well each proxy reconstructed a ground-truth EEG obtained when the synaptic currents of the LIF model network were fed into a three-dimensional network model of multicompartmental neurons with realistic morphologies. Proxies based on linear combinations of AMPA and GABA currents performed better than proxies based on firing rates or membrane potentials. A new class of proxies, based on an optimized linear combination of time-shifted AMPA and GABA currents, provided the most accurate estimate of the EEG over a wide range of network states. The new linear proxies explained 85-95% of the variance of the ground-truth EEG for a wide range of network configurations including different cell morphologies, distributions of presynaptic inputs, positions of the recording electrode, and spatial extensions of the network. Non-linear EEG proxies using a convolutional neural network (CNN) on synaptic currents increased proxy performance by a further 2-8%. Our proxies can be used to easily calculate a biologically realistic EEG signal directly from point-neuron simulations thus facilitating a quantitative comparison between computational models and experimental EEG recordings.

Intrinsic excitation-inhibition imbalance affects medial prefrontal cortex differently in autistic men versus women.

Excitation-inhibition (E:I) imbalance is theorized as an important pathophysiological mechanism in autism. Autism affects males more frequently than females and sex-related mechanisms (e.g., X-linked genes, androgen hormones) can influence E:I balance. This suggests that E:I imbalance may affect autism differently in males versus females. With a combination of in-silico modeling and in-vivo chemogenetic manipulations in mice, we first show that a time-series metric estimated from fMRI BOLD signal, the Hurst exponent (H), can be an index for underlying change in the synaptic E:I ratio. In autism we find that H is reduced, indicating increased excitation, in the medial prefrontal cortex (MPFC) of autistic males but not females. Increasingly intact MPFC H is also associated with heightened ability to behaviorally camouflage social-communicative difficulties, but only in autistic females. This work suggests that H in BOLD can index synaptic E:I ratio and that E:I imbalance affects autistic males and females differently.

Autism is a condition that is usually diagnosed early in life that affects how a person communicates and socializes, and is often characterized by repetitive behaviors. One key theory of autism is that it reflects an imbalance in levels of excitation and inhibition in the brain. Excitatory signals are those that make other brain cells more likely to become active; inhibitory signals have the opposite effect. In non-autistic individuals, inhibitory activity outweighs excitatory activity. In people with autism, by contrast, an increase in excitatory activity is believed to produce an imbalance in excitation and inhibition. Most of the evidence to support this excitation-inhibition imbalance theory has come from studies of rare mutations that cause autism. Many of these mutations occur on the sex chromosomes or are influenced by androgen hormones (hormones that usually play a role on typically male traits). However, most people with autism do not possess these particular mutations. It was thus unclear whether the theory could apply to everyone with autism or, for example, whether it may better apply to specific groups of individuals based on their sex or gender. This is especially important given that about four times as many men and boys compared to women and girls are diagnosed with autism. Trakoshis, Martínez-Cañada et al. have now found a way to ask whether any imbalance in excitation and inhibition in the brain occurs differently in men and women. Using computer modeling, they identified a signal in brain scans that corresponds to an imbalance of excitation and inhibition. After showing that the technique works to identify real increases in excitation in the brain scans of mice, Trakoshis, Martínez-Cañada et al. looked for this signal, or biomarker, in brain scans of people with and without autism. All the people in the study identified with the gender that matched the sex they were assigned at birth. The results revealed differences between the men and women with autism. Men with autism showed an imbalance in excitation and inhibition in specific 'social brain' regions including the medial prefrontal cortex, but women with autism did not. Notably, many of these brain regions are strongly affected by androgen hormones. Previous studies have found that women with autism are sometimes better at hiding or 'camouflaging' their difficulties when socializing or communicating than men with autism. Trakoshis, Martínez-Cañada et al. showed that the better a woman was at camouflaging her autism, the more her brain activity in this region resembled that of non-autistic women. Excitation-inhibition imbalance may thus affect specific brain regions involved in socializing and communication more in men who have autism than in women with the condition. Balanced excitation and inhibition in these brain areas may enable some women with autism to camouflage their difficulties socializing or communicating. Being able to detect imbalances in activity using standard brain imaging could be useful for clinical trials. Future studies could use this biomarker to monitor responses to drug treatments that aim to adjust the balance between excitation and inhibition.

A Neuronal Network Model of the Primate Visual System: Color Mechanisms in the Retina, LGN and V1.

Color plays a key role in human vision but the neural machinery that underlies the transformation from stimulus to perception is not well understood. Here, we implemented a two-dimensional network model of the first stages in the primate parvocellular pathway (retina, lateral geniculate nucleus and layer 4C ß in V1) consisting of conductance-based point neurons. Model parameters were tuned based on physiological and anatomical data from the primate foveal and parafoveal vision, the most relevant visual field areas for color vision. We exhaustively benchmarked the model against well-established chromatic and achromatic visual stimuli, showing spatial and temporal responses of the model to disk- and ring-shaped light flashes, spatially uniform squares and sine-wave gratings of varying spatial frequency. The spatiotemporal patterns of parvocellular cells and cortical cells are consistent with their classification into chromatically single-opponent and double-opponent groups, and nonopponent cells selective for luminance stimuli. The model was implemented in the widely used neural simulation tool NEST and released as open source software. The aim of our modeling is to provide a biologically realistic framework within which a broad range of neuronal interactions can be examined at several different levels, with a focus on understanding how color information is processed.

A Minimal-Model Approach to Analyze Neuronal Circuit Dynamics from multifocal ERG (mERG).

The multifocal electroretinogram (mERG) is a valuable non-invasive tool used by clinicians for studying both the local physiology of the normal human retina, as well as retinal function affected by disease. However, the neural basis of the major components of the mERG is not well understood. Computational modeling of neural circuits has lead to important insights into the structure and workings of nervous systems. A better understanding of the different neural factors that contribute to the retina function and progression of the pathology can help in developing an effective clinical intervention. In this paper, we propose a minimal computational model as an analysis tool for inferring the relationship between mERG data of the human retina and its underlying neural activity and interpreting the data in terms of the specific neural features that contribute to it. In this preliminary study, the rigor of this analysis technique was assessed by fitting the model with a genetic algorithm (GA) to collected mERG data from ten healthy patients.

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.

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.

Connecting Artificial Brains to Robots in a Comprehensive Simulation Framework: The Neurorobotics Platform.

Combined efforts in the fields of neuroscience, computer science, and biology allowed to design biologically realistic models of the brain based on spiking neural networks. For a proper validation of these models, an embodiment in a dynamic and rich sensory environment, where the model is exposed to a realistic sensory-motor task, is needed. Due to the complexity of these brain models that, at the current stage, cannot deal with real-time constraints, it is not possible to embed them into a real-world task. Rather, the embodiment has to be simulated as well. While adequate tools exist to simulate either complex neural networks or robots and their environments, there is so far no tool that allows to easily establish a communication between brain and body models. The Neurorobotics Platform is a new web-based environment that aims to fill this gap by offering scientists and technology developers a software infrastructure allowing them to connect brain models to detailed simulations of robot bodies and environments and to use the resulting neurorobotic systems for in silico experimentation. In order to simplify the workflow and reduce the level of the required programming skills, the platform provides editors for the specification of experimental sequences and conditions, environments, robots, and brain-body connectors. In addition to that, a variety of existing robots and environments are provided. This work presents the architecture of the first release of the Neurorobotics Platform developed in subproject 10 "Neurorobotics" of the Human Brain Project (HBP). At the current state, the Neurorobotics Platform allows researchers to design and run basic experiments in neurorobotics using simulated robots and simulated environments linked to simplified versions of brain models. We illustrate the capabilities of the platform with three example experiments: a Braitenberg task implemented on a mobile robot, a sensory-motor learning task based on a robotic controller, and a visual tracking embedding a retina model on the iCub humanoid robot. These use-cases allow to assess the applicability of the Neurorobotics Platform for robotic tasks as well as in neuroscientific experiments.

A Computational Framework for Realistic Retina Modeling.

Computational simulations of the retina have led to valuable insights about the biophysics of its neuronal activity and processing principles. A great number of retina models have been proposed to reproduce the behavioral diversity of the different visual processing pathways. While many of these models share common computational stages, previous efforts have been more focused on fitting specific retina functions rather than generalizing them beyond a particular model. Here, we define a set of computational retinal microcircuits that can be used as basic building blocks for the modeling of different retina mechanisms. To validate the hypothesis that similar processing structures may be repeatedly found in different retina functions, we implemented a series of retina models simply by combining these computational retinal microcircuits. Accuracy of the retina models for capturing neural behavior was assessed by fitting published electrophysiological recordings that characterize some of the best-known phenomena observed in the retina: adaptation to the mean light intensity and temporal contrast, and differential motion sensitivity. The retinal microcircuits are part of a new software platform for efficient computational retina modeling from single-cell to large-scale levels. It includes an interface with spiking neural networks that allows simulation of the spiking response of ganglion cells and integration with models of higher visual areas.