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 (≳ 106 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.

Genome-wide analyses reveal novel opioid use disorder loci and genetic overlap with schizophrenia, bipolar disorder, and major depression.

Opioid use disorder (OUD) and mental disorders are often comorbid, with increased morbidity and mortality. The causes underlying this relationship are poorly understood. Although these conditions are highly heritable, their shared genetic vulnerabilities remain unaccounted for. We applied the conditional/conjunctional false discovery rate (cond/conjFDR) approach to analyse summary statistics from independent genome wide association studies of OUD, schizophrenia (SCZ), bipolar disorder (BD) and major depression (MD) of European ancestry. Next, we characterized the identified shared loci using biological annotation resources. OUD data were obtained from the Million Veteran Program, Yale-Penn and Study of Addiction: Genetics and Environment (SAGE) (15 756 cases, 99 039 controls). SCZ (53 386 cases, 77 258 controls), BD (41 917 cases, 371 549 controls) and MD (170 756 cases, 329 443 controls) data were provided by the Psychiatric Genomics Consortium. We discovered genetic enrichment for OUD conditional on associations with SCZ, BD, MD and vice versa, indicating polygenic overlap with identification of 14 novel OUD loci at condFDR < 0.05 and 7 unique loci shared between OUD and SCZ (n = 2), BD (n = 2) and MD (n = 7) at conjFDR < 0.05 with concordant effect directions, in line with estimated positive genetic correlations. Two loci were novel for OUD, one for BD and one for MD. Three OUD risk loci were shared with more than one psychiatric disorder, at DRD2 on chromosome 11 (BD and MD), at FURIN on chromosome 15 (SCZ, BD and MD) and at the major histocompatibility complex region (SCZ and MD). Our findings provide new insights into the shared genetic architecture between OUD and SCZ, BD and MD, indicating a complex genetic relationship, suggesting overlapping neurobiological pathways.

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.

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.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations.

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Focal Local Field Potential Signature of the Single-Axon Monosynaptic Thalamocortical Connection.

A resurgence has taken place in recent years in the use of the extracellularly recorded local field potential (LFP) to investigate neural network activity. To probe monosynaptic thalamic activation of cortical postsynaptic target cells, so called spike-trigger-averaged LFP (stLFP) signatures have been measured. In these experiments, the cortical LFP is measured by multielectrodes covering several cortical lamina and averaged on spontaneous spikes of thalamocortical (TC) cells. Using a well established forward-modeling scheme, we investigated the biophysical origin of this stLFP signature with simultaneous synaptic activation of cortical layer-4 neurons, mimicking the effect of a single afferent spike from a single TC neuron. Constrained by previously measured intracellular responses of the main postsynaptic target cell types and with biologically plausible assumptions regarding the spatial distribution of thalamic synaptic inputs into layer 4, the model predicted characteristic contributions to monosynaptic stLFP signatures both for the regular-spiking (RS) excitatory neurons and the fast-spiking (FS) inhibitory interneurons. In particular, the FS cells generated stLFP signatures of shorter temporal duration than the RS cells. Added together, a sum of the stLFP signatures of these two principal synaptic targets of TC cells were observed to resemble experimentally measured stLFP signatures. Outside the volume targeted by TC afferents, the resulting postsynaptic LFP signals were found to be sharply attenuated. This implies that such stLFP signatures provide a very local measure of TC synaptic activation, and that newly developed inverse current-source density (CSD)-estimation methods are needed for precise assessment of the underlying spatiotemporal CSD profiles.SIGNIFICANCE STATEMENT Despite its long history and prevalent use, the proper interpretation of the extracellularly recorded local field potential (LFP) is still not fully established. Here we investigate by biophysical modeling the origin of the focal LFP signature of the single-axon monosynaptic thalamocortical connection as measured by spike-trigger-averaging of cortical LFPs on spontaneous spikes of thalamocortical neurons. We find that this LFP signature is well accounted for by a model assuming thalamic projections to two cortical layer-4 cell populations: one excitatory (putatively regular-spiking cells) and one inhibitory (putatively fast-spiking cells). The LFP signature is observed to decay sharply outside the cortical region receiving the thalamocortical projection, implying that it indeed provides a very local measure of thalamocortical synaptic activation.

Biophysical Network Modelling of the dLGN Circuit: Different Effects of Triadic and Axonal Inhibition on Visual Responses of Relay Cells.

Despite its prominent placement between the retina and primary visual cortex in the early visual pathway, the role of the dorsal lateral geniculate nucleus (dLGN) in molding and regulating the visual signals entering the brain is still poorly understood. A striking feature of the dLGN circuit is that relay cells (RCs) and interneurons (INs) form so-called triadic synapses, where an IN dendritic terminal can be simultaneously postsynaptic to a retinal ganglion cell (GC) input and presynaptic to an RC dendrite, allowing for so-called triadic inhibition. Taking advantage of a recently developed biophysically detailed multicompartmental model for an IN, we here investigate putative effects of these different inhibitory actions of INs, i.e., triadic inhibition and standard axonal inhibition, on the response properties of RCs. We compute and investigate so-called area-response curves, that is, trial-averaged visual spike responses vs. spot size, for circular flashing spots in a network of RCs and INs. The model parameters are grossly tuned to give results in qualitative accordance with previous in vivo data of responses to such stimuli for cat GCs and RCs. We particularly investigate how the model ingredients affect salient response properties such as the receptive-field center size of RCs and INs, maximal responses and center-surround antagonisms. For example, while triadic inhibition not involving firing of IN action potentials was found to provide only a non-linear gain control of the conversion of input spikes to output spikes by RCs, axonal inhibition was in contrast found to substantially affect the receptive-field center size: the larger the inhibition, the more the RC center size shrinks compared to the GC providing the feedforward excitation. Thus, a possible role of the different inhibitory actions from INs to RCs in the dLGN circuit is to provide separate mechanisms for overall gain control (direct triadic inhibition) and regulation of spatial resolution (axonal inhibition) of visual signals sent to cortex.

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

With rapidly advancing multi-electrode recording technology, the local field potential (LFP) has again become a popular measure of neuronal activity in both research and clinical applications. Proper understanding of the LFP requires detailed mathematical modeling incorporating the anatomical and electrophysiological features of neurons near the recording electrode, as well as synaptic inputs from the entire network. Here we propose a hybrid modeling scheme combining efficient point-neuron network models with biophysical principles underlying LFP generation by real neurons. The LFP predictions rely on populations of network-equivalent multicompartment neuron models with layer-specific synaptic connectivity, can be used with an arbitrary number of point-neuron network populations, and allows for a full separation of simulated network dynamics and LFPs. We apply the scheme to a full-scale cortical network model for a ∼1 mm2 patch of primary visual cortex, predict laminar LFPs for different network states, assess the relative LFP contribution from different laminar populations, and investigate effects of input correlations and neuron density on the LFP. The generic nature of the hybrid scheme and its public implementation in hybridLFPy form the basis for LFP predictions from other and larger point-neuron network models, as well as extensions of the current application with additional biological detail.

Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

COSGAP: COntainerized Statistical Genetics Analysis Pipelines.

Summary: The collection and analysis of sensitive data in large-scale consortia for statistical genetics is hampered by multiple challenges, due to their non-shareable nature. Time-consuming issues in installing software frequently arise due to different operating systems, software dependencies, and limited internet access. For federated analysis across sites, it can be challenging to resolve different problems, including format requirements, data wrangling, setting up analysis on high-performance computing (HPC) facilities, etc. Easier, more standardized, automated protocols and pipelines can be solutions to overcome these issues. We have developed one such solution for statistical genetic data analysis using software container technologies. This solution, named COSGAP: "COntainerized Statistical Genetics Analysis Pipelines," consists of already established software tools placed into Singularity containers, alongside corresponding code and instructions on how to perform statistical genetic analyses, such as genome-wide association studies, polygenic scoring, LD score regression, Gaussian Mixture Models, and gene-set analysis. Using provided helper scripts written in Python, users can obtain auto-generated scripts to conduct the desired analysis either on HPC facilities or on a personal computer. COSGAP is actively being applied by users from different countries and projects to conduct genetic data analyses without spending much effort on software installation, converting data formats, and other technical requirements. Availability and implementation: COSGAP is freely available on GitHub (https://github.com/comorment/containers) under the GPLv3 license.

Dissecting the Shared Genetic Architecture of Common Epilepsies With Cortical Brain Morphology.

Background and Objectives: Epilepsies are associated with differences in cortical thickness (TH) and surface area (SA). However, the mechanisms underlying these relationships remain elusive. We investigated the extent to which these phenotypes share genetic influences. Methods: We analyzed genome-wide association study data on common epilepsies (n = 69,995) and TH and SA (n = 32,877) using Gaussian mixture modeling MiXeR and conjunctional false discovery rate (conjFDR) analysis to quantify their shared genetic architecture and identify overlapping loci. We biologically interrogated the loci using a variety of resources and validated in independent samples. Results: The epilepsies (2.4 k-2.9 k variants) were more polygenic than both SA (1.8 k variants) and TH (1.3 k variants). Despite absent genome-wide genetic correlations, there was a substantial genetic overlap between SA and genetic generalized epilepsy (GGE) (1.1 k), all epilepsies (1.1 k), and juvenile myoclonic epilepsy (JME) (0.7 k), as well as between TH and GGE (0.8 k), all epilepsies (0.7 k), and JME (0.8 k), estimated with MiXeR. Furthermore, conjFDR analysis identified 15 GGE loci jointly associated with SA and 15 with TH, 3 loci shared between SA and childhood absence epilepsy, and 6 loci overlapping between SA and JME. 23 loci were novel for epilepsies and 11 for cortical morphology. We observed a high degree of sign concordance in the independent samples. Discussion: Our findings show extensive genetic overlap between generalized epilepsies and cortical morphology, indicating a complex genetic relationship with mixed-effect directions. The results suggest that shared genetic influences may contribute to cortical abnormalities in epilepsies.

Dissecting the genetic overlap between three complex phenotypes with trivariate MiXeR.

Comorbidities are an increasing global health challenge. Accumulating evidence suggests overlapping genetic architectures underlying comorbid complex human traits and disorders. The bivariate causal mixture model (MiXeR) can quantify the polygenic overlap between complex phenotypes beyond global genetic correlation. Still, the pattern of genetic overlap between three distinct phenotypes, which is important to better characterize multimorbidities, has previously not been possible to quantify. Here, we present and validate the trivariate MiXeR tool, which disentangles the pattern of genetic overlap between three phenotypes using summary statistics from genome-wide association studies (GWAS). Our simulations show that the trivariate MiXeR can reliably reconstruct different patterns of genetic overlap. We further demonstrate how the tool can be used to estimate the proportions of genetic overlap between three phenotypes using real GWAS data, providing examples of complex patterns of genetic overlap between diverse human traits and diseases that could not be deduced from bivariate analyses. This contributes to a better understanding of the etiology of complex phenotypes and the nature of their relationship, which may aid in dissecting comorbidity patterns and their biological underpinnings.

Unraveling the shared genetics of common epilepsies and general cognitive ability.

Objective: Cognitive impairment is prevalent among individuals with epilepsy, and it is possible that genetic factors can underlie this relationship. Here, we investigated the potential shared genetic basis of common epilepsies and general cognitive ability (COG). Methods: We applied linkage disequilibrium score (LDSC) regression, MiXeR and conjunctional false discovery rate (conjFDR) to analyze different aspects of genetic overlap between COG and epilepsies. We used the largest available genome-wide association study data on COG (n = 269,867) and common epilepsies (n = 27,559 cases, 42,436 controls), including the broad phenotypes 'all epilepsy', focal epilepsies and genetic generalized epilepsies (GGE), and as well as specific subtypes. We functionally annotated the identified loci using a variety of biological resources and validated the results in independent samples. Results: Using MiXeR, COG (11.2k variants) was estimated to be almost four times more polygenic than 'all epilepsy', GGE, juvenile myoclonic epilepsy (JME), and childhood absence epilepsy (CAE) (2.5k - 2.9k variants). The other epilepsy phenotypes were insufficiently powered for analysis. We show extensive genetic overlap between COG and epilepsies with significant negative genetic correlations (-0.23 to -0.04). COG was estimated to share 2.9k variants with both GGE and 'all epilepsy', and 2.3k variants with both JME and CAE. Using conjFDR, we identified 66 distinct loci shared between COG and epilepsies, including novel associations for GGE (27), 'all epilepsy' (5), JME (5) and CAE (5). The implicated genes were significantly expressed in multiple brain regions. The results were validated in independent samples (COG: p = 1.0 × 10-14; 'all epilepsy': p = 5.6 × 10-3). Significance: Our study demonstrates a substantial genetic basis shared between epilepsies and COG and identifies novel overlapping genomic loci. Enhancing our understanding of the relationship between epilepsies and COG may lead to the development of novel comorbidity-targeted epilepsy treatments.

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.

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.

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.

Estimation of population firing rates and current source densities from laminar electrode recordings.

This model study investigates the validity of methods used to interpret linear (laminar) multielectrode recordings. In computer experiments extracellular potentials from a synaptically activated population of about 1,000 pyramidal neurons are calculated using biologically realistic compartmental neuron models combined with electrostatic forward modeling. The somas of the pyramidal neurons are located in a 0.4 mm high and wide columnar cylinder, mimicking a stimulus-evoked layer-5 population in a neocortical column. Current-source density (CSD) analysis of the low-frequency part (<500 Hz) of the calculated potentials (local field potentials, LFP) based on the 'inverse' CSD method is, in contrast to the 'standard' CSD method, seen to give excellent estimates of the true underlying CSD. The high-frequency part (>750 Hz) of the potentials (multi-unit activity, MUA) is found to scale approximately as the population firing rate to the power 3/4 and to give excellent estimates of the underlying population firing rate for trial-averaged data. The MUA signal is found to decay much more sharply outside the columnar populations than the LFP.

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.

Independent Component Analysis for Fully Automated Multi-Electrode Array Spike Sorting.

In neural electrophysiology, spike sorting allows to separate different neurons from extracellularly measured recordings. It is an essential processing step in order to understand neural activity and it is an unsupervised problem in nature, since no ground truth information is available. There are several available spike sorting packages, but many of them require a manual intervention to curate the results, which makes the process time consuming and hard to reproduce. Here, we focus on high-density Multi-Electrode Array (MEA) recordings and we present a fully automated pipeline based on Independent Component Analysis (ICA). While ICA has been previously investigated for spike sorting, it has never been compared with fully automated state-of-the-art algorithms. We use realistic simulated datasets to compare the spike sorting performance in terms of complexity, signal-to-noise ratio, and recording duration. We show that an ICA-based fully automated spike sorting approach can be a viable alternative approach due to its precision and robustness, but it needs to be optimized for time constraints and requires sufficient density of electrodes to cover active neurons in the proximity of the MEA.

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.