Hybrid Scheme for Modeling Local Field Potentials from Point-Neuron Networks

Multitask learning of a biophysically-detailed neuron model

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials.

Neurobiological Causal Models of Language Processing

The language faculty is physically realized in the neurobiological infrastructure of the human brain. Despite significant efforts, an integrated understanding of this system remains a formidable challenge. What is missing from most theoretical accounts is a specification of the neural mechanisms that implement language function. Computational models that have been put forward generally lack an explicit neurobiological foundation. We propose a neurobiologically informed causal modeling approach which offers a framework for how to bridge this gap. A neurobiological causal model is a mechanistic description of language processing that is grounded in, and constrained by, the characteristics of the neurobiological substrate. It intends to model the generators of language behavior at the level of implementational causality. We describe key features and neurobiological component parts from which causal models can be built and provide guidelines on how to implement them in model simulations. Then we outline how this approach can shed new light on the core computational machinery for language, the long-term storage of words in the mental lexicon and combinatorial processing in sentence comprehension. In contrast to cognitive theories of behavior, causal models are formulated in the "machine language" of neurobiology which is universal to human cognition. We argue that neurobiological causal modeling should be pursued in addition to existing approaches. Eventually, this approach will allow us to develop an explicit computational neurobiology of language.

A neurophysiological basis for aperiodic EEG and the background spectral trend

Electroencephalograms (EEGs) display a mixture of rhythmic and broadband fluctuations, the latter manifesting as an apparent 1/f spectral trend. While network oscillations are known to generate rhythmic EEG, the neural basis of broadband EEG remains unexplained. Here, we use biophysical modelling to show that aperiodic neural activity can generate detectable scalp potentials and shape broadband EEG features, but that these aperiodic signals do not significantly perturb brain rhythm quantification. Further model analysis demonstrated that rhythmic EEG signals are profoundly corrupted by shifts in synapse properties. To examine this scenario, we recorded EEGs of human subjects being administered propofol, a general anesthetic and GABA receptor agonist. Drug administration caused broadband EEG changes that quantitatively matched propofol's known effects on GABA receptors. We used our model to correct for these confounding broadband changes, which revealed that delta power, uniquely, increased within seconds of individuals losing consciousness. Altogether, this work details how EEG signals are shaped by neurophysiological factors other than brain rhythms and elucidates how these signals can undermine traditional EEG interpretation.

Multiscale co-simulation design pattern for neuroscience applications

Integration of information across heterogeneous sources creates added scientific value. Interoperability of data, tools and models is, however, difficult to accomplish across spatial and temporal scales. Here we introduce the toolbox Parallel Co-Simulation, which enables the interoperation of simulators operating at different scales. We provide a software science co-design pattern and illustrate its functioning along a neuroscience example, in which individual regions of interest are simulated on the cellular level allowing us to study detailed mechanisms, while the remaining network is efficiently simulated on the population level. A workflow is illustrated for the use case of The Virtual Brain and NEST, in which the CA1 region of the cellular-level hippocampus of the mouse is embedded into a full brain network involving micro and macro electrode recordings. This new tool allows integrating knowledge across scales in the same simulation framework and validating them against multiscale experiments, thereby largely widening the explanatory power of computational models.

Laminar neural dynamics of auditory evoked responses: Computational modeling of local field potentials in auditory cortex of non-human primates

Evoked neural responses to sensory stimuli have been extensively investigated in humans and animal models both to enhance our understanding of brain function and to aid in clinical diagnosis of neurological and neuropsychiatric conditions. Recording and imaging techniques such as electroencephalography (EEG), magnetoencephalography (MEG), local field potentials (LFPs), and calcium imaging provide complementary information about different aspects of brain activity at different spatial and temporal scales. Modeling and simulations provide a way to integrate these different types of information to clarify underlying neural mechanisms. In this study, we aimed to shed light on the neural dynamics underlying auditory evoked responses by fitting a rate-based model to LFPs recorded via multi-contact electrodes which simultaneously sampled neural activity across cortical laminae. Recordings included neural population responses to best-frequency (BF) and non-BF tones at four representative sites in primary auditory cortex (A1) of awake monkeys. The model considered major neural populations of excitatory, parvalbumin-expressing (PV), and somatostatin-expressing (SOM) neurons across layers 2/3, 4, and 5/6. Unknown parameters, including the connection strength between the populations, were fitted to the data. Our results revealed similar population dynamics, fitted model parameters, predicted equivalent current dipoles (ECD), tuning curves, and lateral inhibition profiles across recording sites and animals, in spite of quite different extracellular current distributions. We found that PV firing rates were higher in BF than in non-BF responses, mainly due to different strengths of tonotopic thalamic input, whereas SOM firing rates were higher in non-BF than in BF responses due to lateral inhibition. In conclusion, we demonstrate the feasibility of the model-fitting approach in identifying the contributions of cell-type specific population activity to stimulus-evoked LFPs across cortical laminae, providing a foundation for further investigations into the dynamics of neural circuits underlying cortical sensory processing.

Metamodelling of a two-population spiking neural network

In computational neuroscience, hypotheses are often formulated as bottom-up mechanistic models of the systems in question, consisting of differential equations that can be numerically integrated forward in time. Candidate models can then be validated by comparison against experimental data. The model outputs of neural network models depend on both neuron parameters, connectivity parameters and other model inputs. Successful model fitting requires sufficient exploration of the model parameter space, which can be computationally demanding. Additionally, identifying degeneracy in the parameters, i.e. different combinations of parameter values that produce similar outputs, is of interest, as they define the subset of parameter values consistent with the data. In this computational study, we apply metamodels to a two-population recurrent spiking network of point-neurons, the so-called Brunel network. Metamodels are data-driven approximations to more complex models with more desirable computational properties, which can be run considerably faster than the original model. Specifically, we apply and compare two different metamodelling techniques, masked autoregressive flows (MAF) and deep Gaussian process regression (DGPR), to estimate the power spectra of two different signals; the population spiking activities and the local field potential. We find that the metamodels are able to accurately model the power spectra in the asynchronous irregular regime, and that the DGPR metamodel provides a more accurate representation of the simulator compared to the MAF metamodel. Using the metamodels, we estimate the posterior probability distributions over parameters given observed simulator outputs separately for both LFP and population spiking activities. We find that these distributions correctly identify parameter combinations that give similar model outputs, and that some parameters are significantly more constrained by observing the LFP than by observing the population spiking activities.

Neuronal connected burst cascades bridge macroscale adaptive signatures across arousal states

The human brain displays a rich repertoire of states that emerge from the microscopic interactions of cortical and subcortical neurons. Difficulties inherent within large-scale simultaneous neuronal recording limit our ability to link biophysical processes at the microscale to emergent macroscopic brain states. Here we introduce a microscale biophysical network model of layer-5 pyramidal neurons that display graded coarse-sampled dynamics matching those observed in macroscale electrophysiological recordings from macaques and humans. We invert our model to identify the neuronal spike and burst dynamics that differentiate unconscious, dreaming, and awake arousal states and provide insights into their functional signatures. We further show that neuromodulatory arousal can mediate different modes of neuronal dynamics around a low-dimensional energy landscape, which in turn changes the response of the model to external stimuli. Our results highlight the promise of multiscale modelling to bridge theories of consciousness across spatiotemporal scales.

Uncovering circuit mechanisms of current sinks and sources with biophysical simulations of primary visual cortex

Local field potential (LFP) recordings reflect the dynamics of the current source density (CSD) in brain tissue. The synaptic, cellular, and circuit contributions to current sinks and sources are ill-understood. We investigated these in mouse primary visual cortex using public Neuropixels recordings and a detailed circuit model based on simulating the Hodgkin-Huxley dynamics of >50,000 neurons belonging to 17 cell types. The model simultaneously captured spiking and CSD responses and demonstrated a two-way dissociation: firing rates are altered with minor effects on the CSD pattern by adjusting synaptic weights, and CSD is altered with minor effects on firing rates by adjusting synaptic placement on the dendrites. We describe how thalamocortical inputs and recurrent connections sculpt specific sinks and sources early in the visual response, whereas cortical feedback crucially alters them in later stages. These results establish quantitative links between macroscopic brain measurements (LFP/CSD) and microscopic biophysics-based understanding of neuron dynamics and show that CSD analysis provides powerful constraints for modeling beyond those from considering spikes.

Theta oscillations represent collective dynamics of multineuronal membrane potentials of murine hippocampal pyramidal cells

Theta (θ) oscillations are one of the characteristic local field potentials (LFPs) in the hippocampus that emerge during spatial navigation, exploratory sniffing, and rapid eye movement sleep. LFPs are thought to summarize multineuronal events, including synaptic currents and action potentials. However, no in vivo study to date has directly interrelated θ oscillations with the membrane potentials (Vm) of multiple neurons, and it remains unclear whether LFPs can be predicted from multineuronal Vms. Here, we simultaneously patch-clamp up to three CA1 pyramidal neurons in awake or anesthetized mice and find that the temporal evolution of the power and frequency of θ oscillations in Vms (θ_Vms) are weakly but significantly correlate with LFP θ oscillations (θ_LFP) such that a deep neural network could predict the θ_LFP waveforms based on the θ_Vm traces of three neurons. Therefore, individual neurons are loosely interdependent to ensure freedom of activity, but they partially share information to collectively produce θ_LFP.

Uncovering the organization of neural circuits with Generalized Phase Locking Analysis

Despite the considerable progress of in vivo neural recording techniques, inferring the biophysical mechanisms underlying large scale coordination of brain activity from neural data remains challenging. One obstacle is the difficulty to link high dimensional functional connectivity measures to mechanistic models of network activity. We address this issue by investigating spike-field coupling (SFC) measurements, which quantify the synchronization between, on the one hand, the action potentials produced by neurons, and on the other hand mesoscopic "field" signals, reflecting subthreshold activities at possibly multiple recording sites. As the number of recording sites gets large, the amount of pairwise SFC measurements becomes overwhelmingly challenging to interpret. We develop Generalized Phase Locking Analysis (GPLA) as an interpretable dimensionality reduction of this multivariate SFC. GPLA describes the dominant coupling between field activity and neural ensembles across space and frequencies. We show that GPLA features are biophysically interpretable when used in conjunction with appropriate network models, such that we can identify the influence of underlying circuit properties on these features. We demonstrate the statistical benefits and interpretability of this approach in various computational models and Utah array recordings. The results suggest that GPLA, used jointly with biophysical modeling, can help uncover the contribution of recurrent microcircuits to the spatio-temporal dynamics observed in multi-channel experimental recordings.

Towards a mesoscale physical modeling framework for stereotactic-EEG recordings

Objective.Stereotactic-electroencephalography (SEEG) and scalp EEG recordings can be modeled using mesoscale neural mass population models (NMMs). However, the relationship between those mathematical models and the physics of the measurements is unclear. In addition, it is challenging to represent SEEG data by combining NMMs and volume conductor models due to the intermediate spatial scale represented by these measurements.Approach.We provide a framework combining the multi-compartmental modeling formalism and a detailed geometrical model to simulate the transmembrane currents that appear in layer 3, 5 and 6 pyramidal cells due to a synaptic input. With this approach, it is possible to realistically simulate the current source density (CSD) depth profile inside a cortical patch due to inputs localized into a single cortical layer and the induced voltage measured by two SEEG contacts using a volume conductor model. Based on this approach, we built a framework to connect the activity of a NMM with a volume conductor model and we simulated an example of SEEG signal as a proof of concept.Main results.CSD depends strongly on the distribution of the synaptic inputs onto the different cortical layers and the equivalent current dipole strengths display substantial differences (of up to a factor of four in magnitude in our example). Thus, the inputs coming from different neural populations do not contribute equally to the electrophysiological recordings. A direct consequence of this is that the raw output of NMMs is not a good proxy for electrical recordings. We also show that the simplest CSD model that can accurately reproduce SEEG measurements can be constructed from discrete monopolar sources (one per cortical layer).Significance.Our results highlight the importance of including a physical model in NMMs to represent measurements. We provide a framework connecting microscale neuron models with the neural mass formalism and with physical models of the measurement process that can improve the accuracy of predicted electrophysiological recordings.

Cortical Pyramidal and Parvalbumin Cells Exhibit Distinct Spatiotemporal Extracellular Electric Potentials

Brain circuits are composed of diverse cell types with distinct morphologies, connections, and distributions of ion channels. Modeling suggests that the spatial distribution of the extracellular voltage during a spike depends on cellular morphology, connectivity, and identity. However, experimental evidence from the intact brain is lacking. Here, we combined high-density recordings from hippocampal region CA1 and neocortex of freely moving mice with optogenetic tagging of parvalbumin-immunoreactive (PV) cells. We used ground truth tagging of the recorded pyramidal cells (PYR) and PV cells to construct binary classification models. Features derived from single-channel waveforms or from spike timing alone allowed near-perfect classification of PYR and PV cells. To determine whether there is unique information in the spatial distribution of the extracellular potentials, we removed all single-channel waveform information from the multichannel waveforms using an event-based delta-transformation. We found that spatiotemporal features derived from the transformed waveforms yield accurate classification. The extracellular analog of the spatial distribution of the initial depolarization phase provided the highest contribution to the spatially based prediction. Compared with PV cell spikes, PYR spikes exhibited higher spatial synchrony at the beginning of the extracellular spike and lower synchrony at the trough. The successful classification of PYR and PV cells based on purely spatial features provides direct experimental evidence that spikes of distinct cell types are associated with distinct spatial distributions of extracellular potentials.

Mean-field based framework for forward modeling of LFP and MEG signals

The use of mean-field models to describe the activity of large neuronal populations has become a very powerful tool for large-scale or whole brain simulations. However, the calculation of brain signals from mean-field models, such as the electric and magnetic fields, is still under development. Thus, the emergence of new methods for an accurate and efficient calculation of such brain signals is currently of great relevance. In this paper we propose a novel method to calculate the local field potentials (LFP) and magnetic fields from mean-field models. The calculation of LFP is done via a kernel method based on unitary LFP's (the LFP generated by a single axon) that was recently introduced for spiking-networks simulations and that we adapt here for mean-field models. The calculation of the magnetic field is based on current-dipole and volume-conductor models, where the secondary currents (due to the conducting extracellular medium) are estimated using the LFP calculated via the kernel method and the effects of medium-inhomogeneities are incorporated. We provide an example of the application of our method for the calculation of LFP and MEG under slow-waves of neuronal activity generated by a mean-field model of a network of Adaptive-Exponential Integrate-and-Fire (AdEx) neurons. We validate our method via comparison with results obtained from the corresponding spiking neuronal networks. Finally we provide an example of our method for whole brain simulations performed with The Virtual Brain (TVB), a recently developed tool for large scale simulations of the brain. Our method provides an efficient way of calculating electric and magnetic fields from mean-field models. This method exhibits a great potential for its application in large-scale or whole-brain simulations, where calculations via detailed biological models are not feasible.

Brain signal predictions from multi-scale networks using a linearized framework

Simulations of neural activity at different levels of detail are ubiquitous in modern neurosciences, aiding the interpretation of experimental data and underlying neural mechanisms at the level of cells and circuits. Extracellular measurements of brain signals reflecting transmembrane currents throughout the neural tissue remain commonplace. The lower frequencies (≲ 300Hz) of measured signals generally stem from synaptic activity driven by recurrent interactions among neural populations and computational models should also incorporate accurate predictions of such signals. Due to limited computational resources, large-scale neuronal network models (≳ 10⁶ neurons or so) often require reducing the level of biophysical detail and account mainly for times of action potentials (‘spikes’) or spike rates. Corresponding extracellular signal predictions have thus poorly accounted for their biophysical origin.

Here we propose a computational framework for predicting spatiotemporal filter kernels for such extracellular signals stemming from synaptic activity, accounting for the biophysics of neurons, populations, and recurrent connections. Signals are obtained by convolving population spike rates by appropriate kernels for each connection pathway and summing the contributions. Our main results are that kernels derived via linearized synapse and membrane dynamics, distributions of cells, conduction delay, and volume conductor model allow for accurately capturing the spatiotemporal dynamics of ground truth extracellular signals from conductance-based multicompartment neuron networks. One particular observation is that changes in the effective membrane time constants caused by persistent synapse activation must be accounted for.

The work also constitutes a major advance in computational efficiency of accurate, biophysics-based signal predictions from large-scale spike and rate-based neuron network models drastically reducing signal prediction times compared to biophysically detailed network models. This work also provides insight into how experimentally recorded low-frequency extracellular signals of neuronal activity may be approximately linearly dependent on spiking activity. A new software tool LFPykernels serves as a reference implementation of the framework.

Understanding the brain’s function and activity in healthy and pathological states across spatial scales and times spanning entire lives is one of humanity’s great undertakings. In experimental and clinical work probing the brain’s activity, a variety of electric and magnetic measurement techniques are routinely applied. However interpreting the extracellularly measured signals remains arduous due to multiple factors, mainly the large number of neurons contributing to the signals and complex interactions occurring in recurrently connected neuronal circuits. To understand how neurons give rise to such signals, mechanistic modeling combined with forward models derived using volume conductor theory has proven to be successful, but this approach currently does not scale to the systems level (encompassing millions of neurons or more) where simplified or abstract neuron representations typically are used. Motivated by experimental findings implying approximately linear relationships between times of neuronal action potentials and extracellular population signals, we provide a biophysics-based method for computing causal filters relating spikes and extracellular signals that can be applied with spike times or rates of large-scale neuronal network models for predictions of population signals without relying on ad hoc approximations.

NNMT: Mean-Field Based Analysis Tools for Neuronal Network Models

Mean-field theory of neuronal networks has led to numerous advances in our analytical and intuitive understanding of their dynamics during the past decades. In order to make mean-field based analysis tools more accessible, we implemented an extensible, easy-to-use open-source Python toolbox that collects a variety of mean-field methods for the leaky integrate-and-fire neuron model. The Neuronal Network Mean-field Toolbox (NNMT) in its current state allows for estimating properties of large neuronal networks, such as firing rates, power spectra, and dynamical stability in mean-field and linear response approximation, without running simulations. In this article, we describe how the toolbox is implemented, show how it is used to reproduce results of previous studies, and discuss different use-cases, such as parameter space explorations, or mapping different network models. Although the initial version of the toolbox focuses on methods for leaky integrate-and-fire neurons, its structure is designed to be open and extensible. It aims to provide a platform for collecting analytical methods for neuronal network model analysis, such that the neuroscientific community can take maximal advantage of them.

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.

Active Dendrites and Local Field Potentials: Biophysical Mechanisms and Computational Explorations

Neurons and glial cells are endowed with membranes that express a rich repertoire of ion channels, transporters, and receptors. The constant flux of ions across the neuronal and glial membranes results in voltage fluctuations that can be recorded from the extracellular matrix. The high frequency components of this voltage signal contain information about the spiking activity, reflecting the output from the neurons surrounding the recording location. The low frequency components of the signal, referred to as the local field potential (LFP), have been traditionally thought to provide information about the synaptic inputs that impinge on the large dendritic trees of various neurons. In this review, we discuss recent computational and experimental studies pointing to a critical role of several active dendritic mechanisms that can influence the genesis and the location-dependent spectro-temporal dynamics of LFPs, spanning different brain regions. We strongly emphasize the need to account for the several fast and slow dendritic events and associated active mechanisms - including gradients in their expression profiles, inter- and intra-cellular spatio-temporal interactions spanning neurons and glia, heterogeneities and degeneracy across scales, neuromodulatory influences, and activitydependent plasticity - towards gaining important insights about the origins of LFP under different behavioral states in health and disease. We provide simple but essential guidelines on how to model LFPs taking into account these dendritic mechanisms, with detailed methodology on how to account for various heterogeneities and electrophysiological properties of neurons and synapses while studying LFPs.

NetPyNE Implementation and Scaling of the Potjans-Diesmann Cortical Microcircuit Model

The Potjans-Diesmann cortical microcircuit model is a widely used model originally implemented in NEST. Here, we reimplemented the model using NetPyNE, a high-level Python interface to the NEURON simulator, and reproduced the findings of the original publication. We also implemented a method for scaling the network size that preserves first- and second-order statistics, building on existing work on network theory. Our new implementation enabled the use of more detailed neuron models with multicompartmental morphologies and multiple biophysically realistic ion channels. This opens the model to new research, including the study of dendritic processing, the influence of individual channel parameters, the relation to local field potentials, and other multiscale interactions. The scaling method we used provides flexibility to increase or decrease the network size as needed when running these CPU-intensive detailed simulations. Finally, NetPyNE facilitates modifying or extending the model using its declarative language; optimizing model parameters; running efficient, large-scale parallelized simulations; and analyzing the model through built-in methods, including local field potential calculation and information flow measures.

Hybrid modelling for stroke care: Review and suggestions of new approaches for risk assessment and simulation of scenarios

Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.

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.

Brain-to-brain communication: the possible role of brain electromagnetic fields (As a Potential Hypothesis)

Up now, the communication between brains of different humans or animals has been confirmed and confined by the sensory medium and motor facilities of body. Recently, direct brain-to-brain communication (DBBC) outside the conventional five senses has been verified between animals and humans. Nevertheless, no empirical studies or serious discussion have been performed to elucidate the mechanism behind this process. The validation of DBBC has been documented via recording similar pattern of action potentials occurring in the brain cortex of two animals. With regard to action potentials in brain neurons, the magnetic field resulting from the action potentials created in neurons is one of the tools where the brain of one animal can affect the brain of another. It has been shown that different animals, even humans, have the power to understand the magnetic field. Cryptochrome, which exists in the retina and in different regions of the brain, has been confirmed to be able to perceive magnetic fields and convert magnetic fields to action potentials. Recently, iron particles (Fe₃O₄) believed to be functioning as magnets have been found in various parts of the brain, and are postulated as magnetic field receptors. Newly developed supersensitive magnetic sensors made of iron magnets that can sense the brain's magnetic field have suggested the idea that these Fe₃O₄ particles or magnets may be capable of perceiving the brain's extremely weak magnetic field. The present study suggests that it is possible the extremely week magnetic field in one animal's brain to transmit vital and accurate information to another animal's brain.

Empirically constrained network models for contrast-dependent modulation of gamma rhythm in V1

Gamma oscillations are thought to play a key role in neuronal network function and neuronal communication, yet the underlying generating mechanisms have not been fully elucidated to date. At least partly, this may be due to the fact that even in simple network models of interconnected inhibitory (I) and excitatory (E) neurons, many parameters remain unknown and are set based on practical considerations or by convention. Here, we mitigate this problem by requiring PING (Pyramidal Interneuron Network Gamma) models to simultaneously satisfy a broad set of criteria for realistic behaviour based on empirical data spanning both the single unit (spikes) and local population (LFP) levels while unknown parameters are varied. By doing so, we were able to constrain the parameter ranges and select empirically valid models. The derived model constraints implied weak rather than strong PING as the generating mechanism for gamma, connectivity between E and I neurons within specific bounds, and variations of the external input to E but not I neurons. Constrained models showed valid behaviours, including gamma frequency increases with contrast and power saturation or decay at high contrasts. Using an empirically-validated model we studied the route to gamma instability at high contrasts. This involved increased heterogeneity of E neurons with increasing input triggering a breakdown of I neuron pacemaker function. Further, we illustrate the model's capacity to resolve disputes in the literature concerning gamma oscillation properties and GABA conductance proxies. We propose that the models derived in our study will be useful for other modelling studies, and that our approach to the empirical constraining of PING models can be expanded when richer empirical datasets become available. As local gamma networks are the building blocks of larger networks that aim to understand complex cognition through their interactions, there is considerable value in improving our models of these building blocks.

RippleNet: a Recurrent Neural Network for Sharp Wave Ripple (SPW-R) Detection

Hippocampal sharp wave ripples (SPW-R) have been identified as key bio-markers of important brain functions such as memory consolidation and decision making. Understanding their underlying mechanisms in healthy and pathological brain function and behaviour rely on accurate SPW-R detection. In this multidisciplinary study, we propose a novel, self-improving artificial intelligence (AI) detection method in the form of deep Recurrent Neural Networks (RNN) with Long Short-Term memory (LSTM) layers that can learn features of SPW-R events from raw, labeled input data. The approach contrasts conventional routines that typically relies on hand-crafted, heuristic feature extraction and often laborious manual curation. The algorithm is trained using supervised learning on hand-curated data sets with SPW-R events obtained under controlled conditions. The input to the algorithm is the local field potential (LFP), the low-frequency part of extracellularly recorded electric potentials from the CA1 region of the hippocampus. Its output predictions can be interpreted as time-varying probabilities of SPW-R events for the duration of the inputs. A simple thresholding applied to the output probabilities is found to identify times of SPW-R events with high precision. The non-causal, or bidirectional variant of the proposed algorithm demonstrates consistently better accuracy compared to the causal, or unidirectional counterpart. Reference implementations of the algorithm, named ‘RippleNet’, are open source, freely available, and implemented using a common open-source framework for neural networks (tensorflow.keras) and can be easily incorporated into existing data analysis workflows for processing experimental data.

Modeling nucleus accumbens : A Computational Model from Single Cell to Circuit Level

Nucleus accumbens is part of the neural structures required for reward based learning and cognitive processing of motivation. Understanding its cellular dynamics and its role in basal ganglia circuits is important not only in diagnosing behavioral disorders and psychiatric problems as addiction and depression but also for developing therapeutic treatments for them. Building a computational model would expand our comprehension of nucleus accumbens. In this work, we are focusing on establishing a model of nucleus accumbens which has not been considered as much as dorsal striatum in computational neuroscience. We will begin by modeling the behavior of single cells and then build a holistic model of nucleus accumbens considering the effect of synaptic currents. We will verify the validity of the model by showing the consistency of simulation results with the empirical data. Furthermore, the simulation results reveal the joint effect of cortical stimulation and dopaminergic modulation on the activity of medium spiny neurons. This effect differentiates with the type of dopamine receptors.

Biophysically detailed forward modeling of the neural origin of EEG and MEG signals

Electroencephalography (EEG) and magnetoencephalography (MEG) are among the most important techniques for non-invasively studying cognition and disease in the human brain. These signals are known to originate from cortical neural activity, typically described in terms of current dipoles. While the link between cortical current dipoles and EEG/MEG signals is relatively well understood, surprisingly little is known about the link between different kinds of neural activity and the current dipoles themselves. Detailed biophysical modeling has played an important role in exploring the neural origin of intracranial electric signals, like extracellular spikes and local field potentials. However, this approach has not yet been taken full advantage of in the context of exploring the neural origin of the cortical current dipoles that are causing EEG/MEG signals. Here, we present a method for reducing arbitrary simulated neural activity to single current dipoles. We find that the method is applicable for calculating extracranial signals, but less suited for calculating intracranial electrocorticography (ECoG) signals. We demonstrate that this approach can serve as a powerful tool for investigating the neural origin of EEG/MEG signals. This is done through example studies of the single-neuron EEG contribution, the putative EEG contribution from calcium spikes, and from calculating EEG signals from large-scale neural network simulations. We also demonstrate how the simulated current dipoles can be used directly in combination with detailed head models, allowing for simulated EEG signals with an unprecedented level of biophysical details. In conclusion, this paper presents a framework for biophysically detailed modeling of EEG and MEG signals, which can be used to better our understanding of non-inasively measured neural activity in humans.

A kernel-based method to calculate local field potentials from networks of spiking neurons

BACKGROUND

The local field potential (LFP) is usually calculated from current sources arising from transmembrane currents, in particular in asymmetric cellular morphologies such as pyramidal neurons.

NEW METHOD

Here, we adopt a different point of view and relate the spiking of neurons to the LFP through efferent synaptic connections and provide a method to calculate LFPs.

RESULTS

We show that the so-called unitary LFPs (uLFP) provide the key to such a calculation. We show experimental measurements and simulations of uLFPs in neocortex and hippocampus, for both excitatory and inhibitory neurons. We fit a "kernel" function to measurements of uLFPs, and we estimate its spatial and temporal spread by using simulations of morphologically detailed reconstructions of hippocampal pyramidal neurons. Assuming that LFPs are the sum of uLFPs generated by every neuron in the network, the LFP generated by excitatory and inhibitory neurons can be calculated by convolving the trains of action potentials with the kernels estimated from uLFPs. This provides a method to calculate the LFP from networks of spiking neurons, even for point neurons for which the LFP is not easily defined. We show examples of LFPs calculated from networks of point neurons and compare to the LFP calculated from synaptic currents.

CONCLUSIONS

The kernel-based method provides a practical way to calculate LFPs from networks of point neurons.

Modelling unitary fields and the single-neuron contribution to local field potentials in the hippocampus

Key points

We simulate the unitary local field potential (uLFP) generated in the hippocampus CA3, using morphologically detailed models.

The model suggests that cancelling effects between apical and basal dendritic synapses explain the low amplitude of excitatory uLFPs.

Inhibitory synapses around the soma do not cancel and could explain the high‐amplitude inhibitory uLFPs.

These results suggest that somatic inhibition constitutes a strong component of LFPs, which may explain a number of experimental observations.

Abstract

Synaptic currents represent a major contribution to the local field potential (LFP) in brain tissue, but the respective contribution of excitatory and inhibitory synapses is not known. Here, we provide estimates of this contribution by using computational models of hippocampal pyramidal neurons, constrained by in vitro recordings. We focus on the unitary LFP (uLFP) generated by single neurons in the CA3 region of the hippocampus. We first reproduce experimental results for hippocampal basket cells, and in particular how inhibitory uLFP are distributed within hippocampal layers. Next, we calculate the uLFP generated by pyramidal neurons, using morphologically reconstructed CA3 pyramidal cells. The model shows that the excitatory uLFP is of small amplitude, smaller than inhibitory uLFPs. Indeed, when the two are simulated together, inhibitory uLFPs mask excitatory uLFPs, which might create the illusion that the inhibitory field is generated by pyramidal cells. These results provide an explanation for the observation that excitatory and inhibitory uLFPs are of the same polarity, in vivo and in vitro. These results suggest that somatic inhibitory currents are large contributors to the LFP, which is important information for interpreting this signal. Finally, the results of our model might form the basis of a simple method to compute the LFP, which could be applied to point neurons for each cell type, thus providing a simple biologically grounded method for calculating LFPs from neural networks.

In conclusion, computational models constrained by in vitro recordings suggest that: (1) Excitatory uLFPs are of smaller amplitude than inhibitory uLFPs. (2) Inhibitory uLFPs form the major contribution to LFPs. (3) uLFPs can be used as a simple model to generate LFPs from spiking networks.

We simulate the unitary local field potential (uLFP) generated in the hippocampus CA3, using morphologically detailed models.

The model suggests that cancelling effects between apical and basal dendritic synapses explain the low amplitude of excitatory uLFPs.

Inhibitory synapses around the soma do not cancel and could explain the high‐amplitude inhibitory uLFPs.

These results suggest that somatic inhibition constitutes a strong component of LFPs, which may explain a number of experimental observations.

Estimation of neural network model parameters from local field potentials (LFPs)

Most modeling in systems neuroscience has been descriptive where neural representations such as ‘receptive fields’, have been found by statistically correlating neural activity to sensory input. In the traditional physics approach to modelling, hypotheses are represented by mechanistic models based on the underlying building blocks of the system, and candidate models are validated by comparing with experiments. Until now validation of mechanistic cortical network models has been based on comparison with neuronal spikes, found from the high-frequency part of extracellular electrical potentials. In this computational study we investigated to what extent the low-frequency part of the signal, the local field potential (LFP), can be used to validate and infer properties of mechanistic cortical network models. In particular, we asked the question whether the LFP can be used to accurately estimate synaptic connection weights in the underlying network. We considered the thoroughly analysed Brunel network comprising an excitatory and an inhibitory population of recurrently connected integrate-and-fire (LIF) neurons. This model exhibits a high diversity of spiking network dynamics depending on the values of only three network parameters. The LFP generated by the network was computed using a hybrid scheme where spikes computed from the point-neuron network were replayed on biophysically detailed multicompartmental neurons. We assessed how accurately the three model parameters could be estimated from power spectra of stationary ‘background’ LFP signals by application of convolutional neural nets (CNNs). All network parameters could be very accurately estimated, suggesting that LFPs indeed can be used for network model validation.

Most of what we have learned about brain networks in vivo has come from the measurement of spikes (action potentials) recorded by extracellular electrodes. The low-frequency part of these signals, the local field potential (LFP), contains unique information about how dendrites in neuronal populations integrate synaptic inputs, but has so far played a lesser role. To investigate whether the LFP can be used to validate network models, we computed LFP signals for a recurrent network model (the Brunel network) for which the ground-truth parameters are known. By application of convolutional neural nets (CNNs) we found that the synaptic weights indeed could be accurately estimated from ‘background’ LFP signals, suggesting a future key role for LFP in development of network models.

Fast simulation of extracellular action potential signatures based on a morphological filtering approximation

Simulating extracellular recordings of neuronal populations is an important and challenging task both for understanding the nature and relationships between extracellular field potentials at different scales, and for the validation of methodological tools for signal analysis such as spike detection and sorting algorithms. Detailed neuronal multicompartmental models with active or passive compartments are commonly used in this objective. Although using such realistic NEURON models could lead to realistic extracellular potentials, it may require a high computational burden making the simulation of large populations difficult without a workstation. We propose in this paper a novel method to simulate extracellular potentials of firing neurons, taking into account the NEURON geometry and the relative positions of the electrodes. The simulator takes the form of a linear geometry based filter that models the shape of an action potential by taking into account its generation in the cell body / axon hillock and its propagation along the axon. The validity of the approach for different NEURON morphologies is assessed. We demonstrate that our method is able to reproduce realistic extracellular action potentials in a given range of axon/dendrites surface ratio, with a time-efficient computational burden.

The asynchronous state's relation to large-scale potentials in cortex

Understanding the relation between large-scale potentials (M/EEG) and their underlying neural activity can improve the precision of research and clinical diagnosis. Recent insights into cortical dynamics highlighted a state of strongly reduced spike count correlations, termed the asynchronous state (AS). The AS has received considerable attention from experimenters and theorists alike, regarding its implications for cortical dynamics and coding of information. However, how reconcilable are these vanishing correlations in the AS with large-scale potentials such as M/EEG observed in most experiments? Typically the latter are assumed to be based on underlying correlations in activity, in particular between subthreshold potentials. We survey the occurrence of the AS across brain states, regions, and layers and argue for a reconciliation of this seeming disparity: large-scale potentials are either observed, first, at transitions between cortical activity states, which entail transient changes in population firing rate, as well as during the AS, and, second, on the basis of sufficiently large, asynchronous populations that only need to exhibit weak correlations in activity. Cells with no or little spiking activity can contribute to large-scale potentials via their subthreshold currents, while they do not contribute to the estimation of spiking correlations, defining the AS. Furthermore, third, the AS occurs only within particular cortical regions and layers associated with the currently selected modality, allowing for correlations at other times and between other areas and layers.

Mind the last spike - firing rate models for mesoscopic populations of spiking neurons

The dominant modeling framework for understanding cortical computations are heuristic firing rate models. Despite their success, these models fall short to capture spike synchronization effects, to link to biophysical parameters and to describe finite-size fluctuations. In this opinion article, we propose that the refractory density method (RDM), also known as age-structured population dynamics or quasi-renewal theory, yields a powerful theoretical framework to build rate-based models for mesoscopic neural populations from realistic neuron dynamics at the microscopic level. We review recent advances achieved by the RDM to obtain efficient population density equations for networks of generalized integrate-and-fire (GIF) neurons - a class of neuron models that has been successfully fitted to various cell types. The theory not only predicts the nonstationary dynamics of large populations of neurons but also permits an extension to finite-size populations and a systematic reduction to low-dimensional rate dynamics. The new types of rate models will allow a re-examination of models of cortical computations under biological constraints.

CoreNEURON : An Optimized Compute Engine for the NEURON Simulator

The NEURON simulator has been developed over the past three decades and is widely used by neuroscientists to model the electrical activity of neuronal networks. Large network simulation projects using NEURON have supercomputer allocations that individually measure in the millions of core hours. Supercomputer centers are transitioning to next generation architectures and the work accomplished per core hour for these simulations could be improved by an order of magnitude if NEURON was able to better utilize those new hardware capabilities. In order to adapt NEURON to evolving computer architectures, the compute engine of the NEURON simulator has been extracted and has been optimized as a library called CoreNEURON. This paper presents the design, implementation, and optimizations of CoreNEURON. We describe how CoreNEURON can be used as a library with NEURON and then compare performance of different network models on multiple architectures including IBM BlueGene/Q, Intel Skylake, Intel MIC and NVIDIA GPU. We show how CoreNEURON can simulate existing NEURON network models with 4-7x less memory usage and 2-7x less execution time while maintaining binary result compatibility with NEURON.

An architectonic type principle integrates macroscopic cortico-cortical connections with intrinsic cortical circuits of the primate brain

The connections linking neurons within and between cerebral cortical areas form a multiscale network for communication. We review recent work relating essential features of cortico-cortical connections, such as their existence and laminar origins and terminations, to fundamental structural parameters of cortical areas, such as their distance, similarity in cytoarchitecture, defined by lamination or neuronal density, and other macroscopic and microscopic structural features. These analyses demonstrate the presence of an architectonic type principle. Across species and cortices, the essential features of cortico-cortical connections vary consistently and strongly with the cytoarchitectonic similarity of cortical areas. By contrast, in multivariate analyses such relations were not found consistently for distance, similarity of cortical thickness, or cellular morphology. Gradients of laminar cortical differentiation, as reflected in overall neuronal density, also correspond to regional variations of cellular features, forming a spatially ordered natural axis of concerted architectonic and connectional changes across the cortical sheet. The robustness of findings across mammalian brains allows cross-species predictions of the existence and laminar patterns of projections, including estimates for the human brain that are not yet available experimentally. The architectonic type principle integrates cortical connectivity and architecture across scales, with implications for computational explorations of cortical physiology and developmental mechanisms.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences

The basolateral nucleus of the amygdala (BL) is thought to support numerous emotional behaviors through specific microcircuits. These are often thought to be comprised of feedforward networks of principal cells (PNs) and interneurons. Neither well-understood nor often considered are recurrent and feedback connections, which likely engender oscillatory dynamics within BL. Indeed, oscillations in the gamma frequency range (40 − 100 Hz) are known to occur in the BL, and yet their origin and effect on local circuits remains unknown. To address this, we constructed a biophysically and anatomically detailed model of the rat BL and its local field potential (LFP) based on the physiological and anatomical literature, along with in vivo and in vitro data we collected on the activities of neurons within the rat BL. Remarkably, the model produced intermittent gamma oscillations (∼50 − 70 Hz) whose properties matched those recorded in vivo, including their entrainment of spiking. BL gamma-band oscillations were generated by the intrinsic circuitry, depending upon reciprocal interactions between PNs and fast-spiking interneurons (FSIs), while connections within these cell types affected the rhythm’s frequency. The model allowed us to conduct experimentally impossible tests to characterize the synaptic and spatial properties of gamma. The entrainment of individual neurons to gamma depended on the number of afferent connections they received, and gamma bursts were spatially restricted in the BL. Importantly, the gamma rhythm synchronized PNs and mediated competition between ensembles. Together, these results indicate that the recurrent connectivity of BL expands its computational and communication repertoire.

The Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX) Version 2.0: Modelling in vitro electrical stimulation of brain tissue

Neuronal circuits can be modelled in detail allowing us to predict the effects of stimulation on individual neurons. Electrical stimulation of neuronal circuits in vitro and in vivo excites a range of neurons within the tissue and measurements of neural activity, e.g the local field potential (LFP), are again an aggregate of a large pool of cells. The previous version of our Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX) allowed for the simulation of the LFP generated by a patch of brain tissue. Here, we extend VERTEX to simulate the effect of electrical stimulation through a focal electric field. We observe both direct changes in neural activity and changes in synaptic plasticity. Testing our software in a model of a rat neocortical slice, we determine the currents contributing to the LFP, the effects of paired pulse stimulation to induce short term plasticity (STP), and the effect of theta burst stimulation (TBS) to induce long term potentiation (LTP).

Multimodal Modeling of Neural Network Activity: Computing LFP, ECoG, EEG, and MEG Signals With LFPy 2.0

Recordings of extracellular electrical, and later also magnetic, brain signals have been the dominant technique for measuring brain activity for decades. The interpretation of such signals is however nontrivial, as the measured signals result from both local and distant neuronal activity. In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons. Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed. This allows for the development of computational tools implementing forward models grounded in the biophysics underlying electrical and magnetic measurement modalities. LFPy (LFPy.readthedocs.io) incorporated a well-established scheme for predicting extracellular potentials of individual neurons with arbitrary levels of biological detail. It relies on NEURON (neuron.yale.edu) to compute transmembrane currents of multicompartment neurons which is then used in combination with an electrostatic forward model. Its functionality is now extended to allow for modeling of networks of multicompartment neurons with concurrent calculations of extracellular potentials and current dipole moments. The current dipole moments are then, in combination with suitable volume-conductor head models, used to compute non-invasive measures of neuronal activity, like scalp potentials (electroencephalographic recordings; EEG) and magnetic fields outside the head (magnetoencephalographic recordings; MEG). One such built-in head model is the four-sphere head model incorporating the different electric conductivities of brain, cerebrospinal fluid, skull and scalp. We demonstrate the new functionality of the software by constructing a network of biophysically detailed multicompartment neuron models from the Neocortical Microcircuit Collaboration (NMC) Portal (bbp.epfl.ch/nmc-portal) with corresponding statistics of connections and synapses, and compute in vivo-like extracellular potentials (local field potentials, LFP; electrocorticographical signals, ECoG) and corresponding current dipole moments. From the current dipole moments we estimate corresponding EEG and MEG signals using the four-sphere head model. We also show strong scaling performance of LFPy with different numbers of message-passing interface (MPI) processes, and for different network sizes with different density of connections. The open-source software LFPy is equally suitable for execution on laptops and in parallel on high-performance computing (HPC) facilities and is publicly available on GitHub.com.

VIOLA-A Multi-Purpose and Web-Based Visualization Tool for Neuronal-Network Simulation Output

Neuronal network models and corresponding computer simulations are invaluable tools to aid the interpretation of the relationship between neuron properties, connectivity, and measured activity in cortical tissue. Spatiotemporal patterns of activity propagating across the cortical surface as observed experimentally can for example be described by neuronal network models with layered geometry and distance-dependent connectivity. In order to cover the surface area captured by today's experimental techniques and to achieve sufficient self-consistency, such models contain millions of nerve cells. The interpretation of the resulting stream of multi-modal and multi-dimensional simulation data calls for integrating interactive visualization steps into existing simulation-analysis workflows. Here, we present a set of interactive visualization concepts called views for the visual analysis of activity data in topological network models, and a corresponding reference implementation VIOLA (VIsualization Of Layer Activity). The software is a lightweight, open-source, web-based, and platform-independent application combining and adapting modern interactive visualization paradigms, such as coordinated multiple views, for massively parallel neurophysiological data. For a use-case demonstration we consider spiking activity data of a two-population, layered point-neuron network model incorporating distance-dependent connectivity subject to a spatially confined excitation originating from an external population. With the multiple coordinated views, an explorative and qualitative assessment of the spatiotemporal features of neuronal activity can be performed upfront of a detailed quantitative data analysis of specific aspects of the data. Interactive multi-view analysis therefore assists existing data analysis workflows. Furthermore, ongoing efforts including the European Human Brain Project aim at providing online user portals for integrated model development, simulation, analysis, and provenance tracking, wherein interactive visual analysis tools are one component. Browser-compatible, web-technology based solutions are therefore required. Within this scope, with VIOLA we provide a first prototype.

Deciphering the Contribution of Oriens-Lacunosum/Moleculare (OLM) Cells to Intrinsic θ Rhythms Using Biophysical Local Field Potential (LFP) Models

Oscillations in local field potentials (LFPs) are prevalent and contribute to brain function. An understanding of the cellular correlates and pathways affecting LFPs is needed, but many overlapping pathways in vivo make this difficult to achieve. A prevalent LFP rhythm in the hippocampus associated with memory processing and spatial navigation is the θ (3–12 Hz) oscillation. θ rhythms emerge intrinsically in an in vitro whole hippocampus preparation and this reduced preparation makes it possible to assess the contribution of different cell types to LFP generation. We focus on oriens-lacunosum/moleculare (OLM) cells as a major class of interneurons in the hippocampus. OLM cells can influence pyramidal (PYR) cells through two distinct pathways: by direct inhibition of PYR cell distal dendrites, and by indirect disinhibition of PYR cell proximal dendrites. We use previous inhibitory network models and build biophysical LFP models using volume conductor theory. We examine the effect of OLM cells to ongoing intrinsic LFP θ rhythms by directly comparing our model LFP features with experiment. We find that OLM cell inputs regulate the robustness of LFP responses without affecting their average power and that this robust response depends on coactivation of distal inhibition and basal excitation. We use our models to estimate the spatial extent of the region generating LFP θ rhythms, leading us to predict that about 22,000 PYR cells participate in intrinsic θ generation. Besides obtaining an understanding of OLM cell contributions to intrinsic LFP θ rhythms, our work can help decipher cellular correlates of in vivo LFPs.

BioNet: A Python interface to NEURON for modeling large-scale networks

There is a significant interest in the neuroscience community in the development of large-scale network models that would integrate diverse sets of experimental data to help elucidate mechanisms underlying neuronal activity and computations. Although powerful numerical simulators (e.g., NEURON, NEST) exist, data-driven large-scale modeling remains challenging due to difficulties involved in setting up and running network simulations. We developed a high-level application programming interface (API) in Python that facilitates building large-scale biophysically detailed networks and simulating them with NEURON on parallel computer architecture. This tool, termed "BioNet", is designed to support a modular workflow whereby the description of a constructed model is saved as files that could be subsequently loaded for further refinement and/or simulation. The API supports both NEURON's built-in as well as user-defined models of cells and synapses. It is capable of simulating a variety of observables directly supported by NEURON (e.g., spikes, membrane voltage, intracellular [Ca++]), as well as plugging in modules for computing additional observables (e.g. extracellular potential). The high-level API platform obviates the time-consuming development of custom code for implementing individual models, and enables easy model sharing via standardized files. This tool will help refocus neuroscientists on addressing outstanding scientific questions rather than developing narrow-purpose modeling code.

h-Type Membrane Current Shapes the Local Field Potential from Populations of Pyramidal Neurons

In cortex, the local field potential (LFP) is thought to mainly stem from correlated synaptic input to populations of geometrically aligned neurons. Computer models of single cortical pyramidal neurons showed that subthreshold voltage-dependent membrane conductances can also shape the LFP signal, in particular the hyperpolarization-activated cation current (I_h; h-type). This ion channel is prominent in various types of pyramidal neurons, typically showing an increasing density gradient along the apical dendrites. Here, we investigate how I_h affects the LFP generated by a model of a population of cortical pyramidal neurons. We find that the LFP from populations of neurons that receive uncorrelated synaptic input can be well predicted by the LFP from single neurons. In this case, when input impinges on the distal dendrites, where most h-type channels are located, a strong resonance in the LFP was measured near the soma, whereas the opposite configuration does not reveal an I_h contribution to the LFP. Introducing correlations in the synaptic inputs to the pyramidal cells strongly amplifies the LFP, while maintaining the differential effects of I_h for distal dendritic versus perisomatic input. Previous theoretical work showed that input correlations do not amplify LFP power when neurons receive synaptic input uniformly across the cell. We find that this crucially depends on the membrane conductance distribution: the asymmetric distribution of I_h results in a strong amplification of the LFP when synaptic inputs to the cell population are correlated. In conclusion, we find that the h-type current is particularly suited to shape the LFP signal in cortical populations.SIGNIFICANCE STATEMENT The local field potential (LFP), the low-frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. While the cortical LFP is thought to mainly reflect synaptic inputs onto pyramidal neurons, little is known about the role of subthreshold active conductances in shaping the LFP. By means of biophysical modeling we obtain a comprehensive, qualitative understanding of how LFPs generated by populations of cortical pyramidal neurons depend on active subthreshold currents, and identify the key importance of the h-type channel. Our results show that LFPs can give information about the active properties of neurons and that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Performance Comparison of the Digital Neuromorphic Hardware SpiNNaker and the Neural Network Simulation Software NEST for a Full-Scale Cortical Microcircuit Model

The digital neuromorphic hardware SpiNNaker has been developed with the aim of enabling large-scale neural network simulations in real time and with low power consumption. Real-time performance is achieved with 1 ms integration time steps, and thus applies to neural networks for which faster time scales of the dynamics can be neglected. By slowing down the simulation, shorter integration time steps and hence faster time scales, which are often biologically relevant, can be incorporated. We here describe the first full-scale simulations of a cortical microcircuit with biological time scales on SpiNNaker. Since about half the synapses onto the neurons arise within the microcircuit, larger cortical circuits have only moderately more synapses per neuron. Therefore, the full-scale microcircuit paves the way for simulating cortical circuits of arbitrary size. With approximately 80, 000 neurons and 0.3 billion synapses, this model is the largest simulated on SpiNNaker to date. The scale-up is enabled by recent developments in the SpiNNaker software stack that allow simulations to be spread across multiple boards. Comparison with simulations using the NEST software on a high-performance cluster shows that both simulators can reach a similar accuracy, despite the fixed-point arithmetic of SpiNNaker, demonstrating the usability of SpiNNaker for computational neuroscience applications with biological time scales and large network size. The runtime and power consumption are also assessed for both simulators on the example of the cortical microcircuit model. To obtain an accuracy similar to that of NEST with 0.1 ms time steps, SpiNNaker requires a slowdown factor of around 20 compared to real time. The runtime for NEST saturates around 3 times real time using hybrid parallelization with MPI and multi-threading. However, achieving this runtime comes at the cost of increased power and energy consumption. The lowest total energy consumption for NEST is reached at around 144 parallel threads and 4.6 times slowdown. At this setting, NEST and SpiNNaker have a comparable energy consumption per synaptic event. Our results widen the application domain of SpiNNaker and help guide its development, showing that further optimizations such as synapse-centric network representation are necessary to enable real-time simulation of large biological neural networks.

The brain as a "hyper-network": the key role of neural networks as main producers of the integrated brain actions especially via the "broadcasted" neuroconnectomics

Investigations of brain complex integrative actions should consider beside neural networks, glial, extracellular molecular, and fluid channels networks. The present paper proposes that all these networks are assembled into the brain hyper-network that has as fundamental components, the tetra-partite synapses, formed by neural, glial, and extracellular molecular networks. Furthermore, peri-synaptic astrocytic processes by modulating the perviousness of extracellular fluid channels control the signals impinging on the tetra-partite synapses. It has also been surmised that global signalling via astrocytes networks and highly pervasive signals, such as electromagnetic fields (EMFs), allow the appropriate integration of the various networks especially at crucial nodes level, the tetra-partite synapses. As a matter of fact, it has been shown that astrocytes can form gap-junction-coupled syncytia allowing intercellular communication characterised by a rapid and possibly long-distance transfer of signals. As far as the EMFs are concerned, the concept of broadcasted neuroconnectomics (BNC) has been introduced to describe highly pervasive signals involved in resetting the information handling of brain networks at various miniaturisation levels. In other words, BNC creates, thanks to the EMFs, generated especially by neurons, different assemblages among the various networks forming the brain hyper-network. Thus, it is surmised that neuronal networks are the "core components" of the brain hyper-network that has as special "nodes" the multi-facet tetra-partite synapses. Furthermore, it is suggested that investigations on the functional plasticity of multi-partite synapses in response to BNC can be the background for a new understanding and perhaps a new modelling of brain morpho-functional organisation and integrative actions.

The roadmap for estimation of cell-type-specific neuronal activity from non-invasive measurements

The computational properties of the human brain arise from an intricate interplay between billions of neurons connected in complex networks. However, our ability to study these networks in healthy human brain is limited by the necessity to use non-invasive technologies. This is in contrast to animal models where a rich, detailed view of cellular-level brain function with cell-type-specific molecular identity has become available due to recent advances in microscopic optical imaging and genetics. Thus, a central challenge facing neuroscience today is leveraging these mechanistic insights from animal studies to accurately draw physiological inferences from non-invasive signals in humans. On the essential path towards this goal is the development of a detailed 'bottom-up' forward model bridging neuronal activity at the level of cell-type-specific populations to non-invasive imaging signals. The general idea is that specific neuronal cell types have identifiable signatures in the way they drive changes in cerebral blood flow, cerebral metabolic rate of O2 (measurable with quantitative functional Magnetic Resonance Imaging), and electrical currents/potentials (measurable with magneto/electroencephalography). This forward model would then provide the 'ground truth' for the development of new tools for tackling the inverse problem-estimation of neuronal activity from multimodal non-invasive imaging data.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Functional roles of three cues that provide nonsynaptic modes of communication in the brain: electromagnetic field, oxygen, and carbon dioxide

The integrative actions of the brain depend on the exchange of information among its computational elements. Hence, this phenomenon plays the key role in driving the complex dynamics of the central nervous system, in which true computations interact with noncomputational dynamical processes to generate brain representations of the body and of the body in the external world, and hence the finalistic behavior of the organism. In this context, it should be pointed out that, besides the intercellular interactions mediated by classical electrochemical signals, other types of interactions, namely, "cues" and "coercions," also appear to be exploited by the system to achieve its function. The present review focuses mainly on cues present in the environment and on those produced by cells of the body, which "pervade" the brain and contribute to its dynamics. These cues can also be metabolic substrates, and, in most cases, they are of fundamental importance to brain function and the survival of the entire organism. Three of these highly pervasive cues will be analyzed in greater detail, namely, oxygen, carbon dioxide, and electromagnetic fields (EMF). Special emphasis will be placed on EMF, since several authors have suggested that these highly pervasive energy fluctuations may play an important role in the global integrative actions of the brain; hence, EMF signaling may transcend classical connectionist models of brain function. Thus the new concept of "broadcasted neuroconnectomics" has been introduced, which transcends the current connectomics view of the brain.