Human Neocortical Neurosolver (HNN), a new software tool for interpreting the cellular and network origin of human MEG/EEG data

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Towards the definition of a standard in TMS-EEG data preprocessing

Combining Non-Invasive Brain Stimulation (NIBS) techniques with the recording of brain electrophysiological activity is an increasingly widespread approach in neuroscience. Particularly successful has been the simultaneous combination of Transcranial Magnetic Stimulation (TMS) and Electroencephalography (EEG). Unfortunately, the strong magnetic pulse required to effectively interact with brain activity inevitably induces artifacts in the concurrent EEG acquisition. Therefore, a careful but aggressive pre-processing is required to efficiently remove artifacts. Unfortunately, as already reported in the literature, different preprocessing approaches can introduce variability in the results. Here we aim at characterizing the three main TMS-EEG preprocessing pipelines currently available, namely ARTIST (Wu et al., 2018), TESA (Rogasch et al., 2017) and SOUND/SSP-SIR (Mutanen et al., 2018, 2016), providing an insight to researchers who need to choose between different approaches. Differently from previous works, we tested the pipelines using a synthetic TMS-EEG signal with a known ground-truth (the artifacts-free to-be-reconstructed signal). In this way, it was possible to assess the reliability of each pipeline precisely and quantitatively, providing a more robust reference for future research. In summary, we found that all pipelines performed well, but with differences in terms of the spatio-temporal precision of the ground-truth reconstruction. Crucially, the three pipelines impacted differently on the inter-trial variability, with ARTIST introducing inter-trial variability not already intrinsic to the ground-truth signal.

Simulation-based Inference of Developmental EEG Maturation with the Spectral Graph Model

The spectral content of macroscopic neural activity evolves throughout development, yet how this maturation relates to underlying brain network formation and dynamics remains unknown. Here, we assess the developmental maturation of electroencephalogram spectra via Bayesian model inversion of the spectral graph model, a parsimonious whole-brain model of spatiospectral neural activity derived from linearized neural field models coupled by the structural connectome. Simulation-based inference was used to estimate age-varying spectral graph model parameter posterior distributions from electroencephalogram spectra spanning the developmental period. This model-fitting approach accurately captures observed developmental electroencephalogram spectral maturation via a neurobiologically consistent progression of key neural parameters: long-range coupling, axonal conduction speed, and excitatory:inhibitory balance. These results suggest that the spectral maturation of macroscopic neural activity observed during typical development is supported by age-dependent functional adaptations in localized neural dynamics and their long-range coupling across the macroscopic structural network.

Bridging model and experiment in systems neuroscience with Cleo: the Closed-Loop, Electrophysiology, and Optophysiology simulation testbed

Systems neuroscience has experienced an explosion of new tools for reading and writing neural activity, enabling exciting new experiments such as all-optical or closed-loop control that effect powerful causal interventions. At the same time, improved computational models are capable of reproducing behavior and neural activity with increasing fidelity. Unfortunately, these advances have drastically increased the complexity of integrating different lines of research, resulting in the missed opportunities and untapped potential of suboptimal experiments. Experiment simulation can help bridge this gap, allowing model and experiment to better inform each other by providing a low-cost testbed for experiment design, model validation, and methods engineering. Specifically, this can be achieved by incorporating the simulation of the experimental interface into our models, but no existing tool integrates optogenetics, two-photon calcium imaging, electrode recording, and flexible closed-loop processing with neural population simulations. To address this need, we have developed Cleo: the Closed-Loop, Electrophysiology, and Optophysiology experiment simulation testbed. Cleo is a Python package enabling injection of recording and stimulation devices as well as closed-loop control with realistic latency into a Brian spiking neural network model. It is the only publicly available tool currently supporting two-photon and multi-opsin/wavelength optogenetics. To facilitate adoption and extension by the community, Cleo is open-source, modular, tested, and documented, and can export results to various data formats. Here we describe the design and features of Cleo, validate output of individual components and integrated experiments, and demonstrate its utility for advancing optogenetic techniques in prospective experiments using previously published systems neuroscience models.

Neuronal Modeling of Cross-Sensory Visual Evoked Magnetoencephalography Responses in the Auditory Cortex

Previous studies have demonstrated that auditory cortex activity can be influenced by cross-sensory visual inputs. Intracortical laminar recordings in nonhuman primates have suggested a feedforward (FF) type profile for auditory evoked but feedback (FB) type for visual evoked activity in the auditory cortex. To test whether cross-sensory visual evoked activity in the auditory cortex is associated with FB inputs also in humans, we analyzed magnetoencephalography (MEG) responses from eight human subjects (six females) evoked by simple auditory or visual stimuli. In the estimated MEG source waveforms for auditory cortex regions of interest, auditory evoked response showed peaks at 37 and 90 ms and visual evoked response at 125 ms. The inputs to the auditory cortex were modeled through FF- and FB-type connections targeting different cortical layers using the Human Neocortical Neurosolver (HNN), which links cellular- and circuit-level mechanisms to MEG signals. HNN modeling suggested that the experimentally observed auditory response could be explained by an FF input followed by an FB input, whereas the cross-sensory visual response could be adequately explained by just an FB input. Thus, the combined MEG and HNN results support the hypothesis that cross-sensory visual input in the auditory cortex is of FB type. The results also illustrate how the dynamic patterns of the estimated MEG source activity can provide information about the characteristics of the input into a cortical area in terms of the hierarchical organization among areas.

intensity- and frequency-specific effects of transcranial alternating current stimulation are explained by network dynamics

Objective. Transcranial alternating current stimulation (tACS) can be used to non-invasively entrain neural activity and thereby cause changes in local neural oscillatory power. Despite its increased use in cognitive and clinical neuroscience, the fundamental mechanisms of tACS are still not fully understood.Approach. We developed a computational neuronal network model of two-compartment pyramidal neurons (PY) and inhibitory interneurons, which mimic the local cortical circuits. We modeled tACS with electric field strengths that are achievable in human applications. We then simulated intrinsic network activity and measured neural entrainment to investigate how tACS modulates ongoing endogenous oscillations.Main results. The intensity-specific effects of tACS are non-linear. At low intensities (<0.3 mV mm-1), tACS desynchronizes neural firing relative to the endogenous oscillations. At higher intensities (>0.3 mV mm-1), neurons are entrained to the exogenous electric field. We then further explore the stimulation parameter space and find that the entrainment of ongoing cortical oscillations also depends on stimulation frequency by following an Arnold tongue. Moreover, neuronal networks can amplify the tACS-induced entrainment via synaptic coupling and network effects. Our model shows that PY are directly entrained by the exogenous electric field and drive the inhibitory neurons.Significance. The results presented in this study provide a mechanistic framework for understanding the intensity- and frequency-specific effects of oscillating electric fields on neuronal networks. This is crucial for rational parameter selection for tACS in cognitive studies and clinical applications.

In-silico testing of new pharmacology for restoring inhibition and human cortical function in depression

Reduced inhibition by somatostatin-expressing interneurons is associated with depression. Administration of positive allosteric modulators of α5 subunit-containing GABAA receptor (α5-PAM) that selectively target this lost inhibition exhibit antidepressant and pro-cognitive effects in rodent models of chronic stress. However, the functional effects of α5-PAM on the human brain in vivo are unknown, and currently cannot be assessed experimentally. We modeled the effects of α5-PAM on tonic inhibition as measured in human neurons, and tested in silico α5-PAM effects on detailed models of human cortical microcircuits in health and depression. We found that α5-PAM effectively recovered impaired cortical processing as quantified by stimulus detection metrics, and also recovered the power spectral density profile of the microcircuit EEG signals. We performed an α5-PAM dose-response and identified simulated EEG biomarker candidates. Our results serve to de-risk and facilitate α5-PAM translation and provide biomarkers in non-invasive brain signals for monitoring target engagement and drug efficacy.

Neuronal modeling of magnetoencephalography responses in auditory cortex to auditory and visual stimuli

Previous studies have demonstrated that auditory cortex activity can be influenced by crosssensory visual inputs. Intracortical recordings in non-human primates (NHP) have suggested a bottom-up feedforward (FF) type laminar profile for auditory evoked but top-down feedback (FB) type for cross-sensory visual evoked activity in the auditory cortex. To test whether this principle applies also to humans, we analyzed magnetoencephalography (MEG) responses from eight human subjects (six females) evoked by simple auditory or visual stimuli. In the estimated MEG source waveforms for auditory cortex region of interest, auditory evoked responses showed peaks at 37 and 90 ms and cross-sensory visual responses at 125 ms. The inputs to the auditory cortex were then modeled through FF and FB type connections targeting different cortical layers using the Human Neocortical Neurosolver (HNN), which consists of a neocortical circuit model linking the cellular- and circuit-level mechanisms to MEG. The HNN models suggested that the measured auditory response could be explained by an FF input followed by an FB input, and the crosssensory visual response by an FB input. Thus, the combined MEG and HNN results support the hypothesis that cross-sensory visual input in the auditory cortex is of FB type. The results also illustrate how the dynamic patterns of the estimated MEG/EEG source activity can provide information about the characteristics of the input into a cortical area in terms of the hierarchical organization among areas.

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.

Human Sensorimotor Beta Event Characteristics and Aperiodic Signal Are Highly Heritable

Individuals' phenotypes, including the brain's structure and function, are largely determined by genes and their interplay. The resting brain generates salient rhythmic patterns that can be characterized noninvasively using functional neuroimaging such as magnetoencephalography (MEG). One of these rhythms, the somatomotor (rolandic) beta rhythm, shows intermittent high amplitude "events" that predict behavior across tasks and species. Beta rhythm is altered in neurological disease. The aperiodic (1/f) signal present in electrophysiological recordings is also modulated by some neurological conditions and aging. Both sensorimotor beta and aperiodic signal could thus serve as biomarkers of sensorimotor function. Knowledge about the extent to which these brain functional measures are heritable could shed light on the mechanisms underlying their generation. We investigated the heritability and variability of human spontaneous sensorimotor beta rhythm events and aperiodic activity in 210 healthy male and female adult siblings' spontaneous MEG activity. The most heritable trait was the aperiodic 1/f signal, with a heritability of 0.87 in the right hemisphere. Time-resolved beta event amplitude parameters were also highly heritable, whereas the heritabilities for overall beta power, peak frequency, and measures of event duration remained nonsignificant. Human sensorimotor neural activity can thus be dissected into different components with variable heritability. We postulate that these differences partially reflect different underlying signal-generating mechanisms. The 1/f signal and beta event amplitude measures may depend more on fixed, anatomical parameters, whereas beta event duration and its modulation reflect dynamic characteristics, guiding their use as potential disease biomarkers.

One out of ten independent components shows flipped polarity with poorer data quality: EEG database study

Independent component analysis (ICA) is widely used today for scalp-recorded EEG analysis. One of the limitations of ICA-based analysis is polarity indeterminacy. It is not easy to find detailed documentations that explains engineering solutions of how the polarity indeterminacy is addressed in a given implementation. We investigated how it is implemented in the case of EEGLAB and also the relation between the outcome of the polarity determination and classification of independent components (ICs) in terms of the estimated nature of the sources (brain, muscle, eye, etc.) using an open database of n = 212 EEG dataset of resting state recordings. We found that (1) about 91% of ICs showed positive-dominant IC scalp topographies; (2) positive-dominant ICs were more associated with brain-originated signals; (3) positive-dominant ICs showed more radial (peaked at 10-30 degrees deviations from the radial axis) dipolar projection pattern with less residual variance from fitting the equivalent current dipole. In conclusion, using the EEGLAB's default ICA algorithm, one out of 10 ICs results in flipping its polarity to negative, which is associated with non-radial dipole orientation with higher residual variance. Thus, we determined EEGLAB biases toward positive polarity in decomposing high-quality brain ICs.

HNN-core: A Python software for cellular and circuit-level interpretation of human MEG/EEG

HNN-core is a library for circuit and cellular level interpretation of non-invasive human magneto-/electro-encephalography (MEG/EEG) data. It is based on the Human Neocortical Neurosolver (HNN) software (Neymotin et al., 2020), a modeling tool designed to simulate multiscale neural mechanisms generating current dipoles in a localized patch of neocortex. HNN's foundation is a biophysically detailed neural network representing a canonical neocortical column containing populations of pyramidal and inhibitory neurons together with layer-specific exogenous synaptic drive (Figure 1 left). In addition to simulating network-level interactions, HNN produces the intracellular currents in the long apical dendrites of pyramidal cells across the cortical layers known to be responsible for macroscopic current dipole generation.

Data-driven multiscale model of macaque auditory thalamocortical circuits reproduces in vivo dynamics

We developed a detailed model of macaque auditory thalamocortical circuits, including primary auditory cortex (A1), medial geniculate body (MGB), and thalamic reticular nucleus, utilizing the NEURON simulator and NetPyNE tool. The A1 model simulates a cortical column with over 12,000 neurons and 25 million synapses, incorporating data on cell-type-specific neuron densities, morphology, and connectivity across six cortical layers. It is reciprocally connected to the MGB thalamus, which includes interneurons and core and matrix-layer-specific projections to A1. The model simulates multiscale measures, including physiological firing rates, local field potentials (LFPs), current source densities (CSDs), and electroencephalography (EEG) signals. Laminar CSD patterns, during spontaneous activity and in response to broadband noise stimulus trains, mirror experimental findings. Physiological oscillations emerge spontaneously across frequency bands comparable to those recorded in vivo. We elucidate population-specific contributions to observed oscillation events and relate them to firing and presynaptic input patterns. The model offers a quantitative theoretical framework to integrate and interpret experimental data and predict its underlying cellular and circuit mechanisms.

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.

Biophysical modeling of frontocentral ERP generation links circuit-level mechanisms of action-stopping to a behavioral race model

Human frontocentral event-related potentials (FC-ERPs) are ubiquitous neural correlates of cognition and control, but their generating multiscale mechanisms remain mostly unknown. We used the Human Neocortical Neurosolver(HNN)'s biophysical model of a canonical neocortical circuit under exogenous thalamic and cortical drive to simulate the cell and circuit mechanisms underpinning the P2, N2, and P3 features of the FC-ERP observed after Stop-Signals in the Stop-Signal task (SST). We demonstrate that a sequence of simulated external thalamocortical and cortico-cortical drives can produce the FC-ERP, similar to what has been shown for primary sensory cortices. We used this model of the FC-ERP to examine likely circuit-mechanisms underlying FC-ERP features that distinguish between successful and failed action-stopping. We also tested their adherence to the predictions of the horse-race model of the SST, with specific hypotheses motivated by theoretical links between the P3 and Stop process. These simulations revealed that a difference in P3 onset between successful and failed Stops is most likely due to a later arrival of thalamocortical drive in failed Stops, rather than, for example, a difference in effective strength of the input. In contrast, the same model predicted that early thalamocortical drives underpinning the P2 and N2 differed in both strength and timing across stopping accuracy conditions. Overall, this model generates novel testable predictions of the thalamocortical dynamics underlying FC-ERP generation during action-stopping. Moreover, it provides a detailed cellular and circuit-level interpretation that supports links between these macroscale signatures and predictions of the behavioral race model.

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.

Diverse beta burst waveform motifs characterize movement-related cortical dynamics

Classical analyses of induced, frequency-specific neural activity typically average band-limited power over trials. More recently, it has become widely appreciated that in individual trials, beta band activity occurs as transient bursts rather than amplitude-modulated oscillations. Most studies of beta bursts treat them as unitary, and having a stereotyped waveform. However, we show there is a wide diversity of burst shapes. Using a biophysical model of burst generation, we demonstrate that waveform variability is predicted by variability in the synaptic drives that generate beta bursts. We then use a novel, adaptive burst detection algorithm to identify bursts from human MEG sensor data recorded during a joystick-based reaching task, and apply principal component analysis to burst waveforms to define a set of dimensions, or motifs, that best explain waveform variance. Finally, we show that bursts with a particular range of waveform motifs, ones not fully accounted for by the biophysical model, differentially contribute to movement-related beta dynamics. Sensorimotor beta bursts are therefore not homogeneous events and likely reflect distinct computational processes.

A guide towards optimal detection of transient oscillatory bursts with unknown parameters

Objectives. Recent event-based analyses of transient neural activities have characterized the oscillatory bursts as a neural signature that bridges dynamic neural states to cognition and behaviors. Following this insight, our study aimed to (1) compare the efficacy of common burst detection algorithms under varying signal-to-noise ratios and event durations using synthetic signals and (2) establish a strategic guideline for selecting the optimal algorithm for real datasets with undefined properties.Approach.We tested the robustness of burst detection algorithms using a simulation dataset comprising bursts of multiple frequencies. To systematically assess their performance, we used a metric called 'detection confidence', quantifying classification accuracy and temporal precision in a balanced manner. Given that burst properties in empirical data are often unknown in advance, we then proposed a selection rule to identify an optimal algorithm for a given dataset and validated its application on local field potentials of basolateral amygdala recorded from male mice (n=8) exposed to a natural threat.Main Results.Our simulation-based evaluation demonstrated that burst detection is contingent upon event duration, whereas accurately pinpointing burst onsets is more susceptible to noise level. For real data, the algorithm chosen based on the selection rule exhibited superior detection and temporal accuracy, although its statistical significance differed across frequency bands. Notably, the algorithm chosen by human visual screening differed from the one recommended by the rule, implying a potential misalignment between human priors and mathematical assumptions of the algorithms.Significance.Therefore, our findings underscore that the precise detection of transient bursts is fundamentally influenced by the chosen algorithm. The proposed algorithm-selection rule suggests a potentially viable solution, while also emphasizing the inherent limitations originating from algorithmic design and volatile performances across datasets. Consequently, this study cautions against relying solely on heuristic-based approaches, advocating for a careful algorithm selection in burst detection studies.

Age-related differences in transient gamma band activity during working memory maintenance through adolescence

Adolescence is a stage of development characterized by neurodevelopmental specialization of cognitive processes. In particular, working memory continues to improve through adolescence, with increases in response accuracy and decreases in response latency continuing well into the twenties. Human electroencephalogram (EEG) studies indicate that gamma oscillations (35-65 Hz) during the working memory delay period support the maintenance of mnemonic information guiding subsequent goal-driven behavior, which decrease in power with development. Importantly, recent electrophysiological studies have shown that gamma events, more so than sustained activity, may underlie working memory maintenance during the delay period. However, developmental differences in gamma events during working memory have not been studied. Here, we used EEG in conjunction with a novel spectral event processing approach to investigate age-related differences in transient gamma band activity during a memory guided saccade (MGS) task in 164 10- to 30-year-olds. Total gamma power was found to significantly decrease through adolescence, replicating prior findings. Results from the spectral event pipeline showed age-related decreases in the mean power of gamma events and trial-by-trial power variability across both the delay period and fixation epochs of the MGS task. In addition, we found that while event number decreased with age during the fixation period, the developmental decrease during the delay period was more dramatic, resulting in an increase in event spiking from fixation to delay in adolescence but not adulthood. While average power of the transient gamma events was found to mediate age-related differences in total gamma power in the fixation and delay periods, the number of gamma events was related to total power in only the delay period, suggesting that the power of gamma events may underlie the sustained gamma activity seen in EEG literature while the number of events may directly support age-related improvements in working memory maintenance. Our findings provide compelling new evidence for mechanistic changes in neural processing characterized by refinements in neural function as behavior becomes optimized in adulthood.

Laminar specificity of the auditory perceptual awareness negativity: A biophysical modeling study

How perception of sensory stimuli emerges from brain activity is a fundamental question of neuroscience. To date, two disparate lines of research have examined this question. On one hand, human neuroimaging studies have helped us understand the large-scale brain dynamics of perception. On the other hand, work in animal models (mice, typically) has led to fundamental insight into the micro-scale neural circuits underlying perception. However, translating such fundamental insight from animal models to humans has been challenging. Here, using biophysical modeling, we show that the auditory awareness negativity (AAN), an evoked response associated with perception of target sounds in noise, can be accounted for by synaptic input to the supragranular layers of auditory cortex (AC) that is present when target sounds are heard but absent when they are missed. This additional input likely arises from cortico-cortical feedback and/or non-lemniscal thalamic projections and targets the apical dendrites of layer-5 (L5) pyramidal neurons. In turn, this leads to increased local field potential activity, increased spiking activity in L5 pyramidal neurons, and the AAN. The results are consistent with current cellular models of conscious processing and help bridge the gap between the macro and micro levels of perception-related brain activity.

Relative Contributions of Different Neural Sources to the EEG

Dogma dictates that the EEG signal is generated by postsynaptic currents (PSCs) because there are an enormous number of synapses in the brain, and PSCs have relatively long durations. However, PSCs are not the only potential source of electric fields in the brain. Action potentials, afterpolarizations, and presynaptic activity can also generate electric fields. Experimentally it is exceedingly difficult to delineate the contributions of different sources because they are casually linked. However, using computational modeling, we can interrogate the relative contributions of different neural elements to the EEG. We used a library of neuron models with morphologically realistic axonal arbors to quantify the relative contributions of PSCs, action potentials, and presynaptic activity to the EEG signal. Consistent with prior assertions, PSCs were the largest contributor to the EEG, but action potentials and afterpolarizations can also make appreciable contributions. For a population of neurons generating simultaneous PSCs and action potentials, we found that the action potentials accounted for up to 20% of the source strength while PSCs accounted for the other 80% and presynaptic activity negligibly contributed. Additionally, L5 PCs generated the largest PSC and action potential signals indicating that they the dominant EEG signal generator. Further, action potentials and afterpolarizations were sufficient to generate physiological oscillations, indicating that they are valid source contributors to the EEG. The EEG emerges from a combination of multiple different source, and, while PSCs are the largest contributor, other sources are non-negligible and should be included in modeling, analysis and interpretation of the EEG.

Intensity- and frequency-specific effects of transcranial alternating current stimulation are explained by network dynamics

Transcranial alternating current stimulation (tACS) can be used to non-invasively entrain neural activity, and thereby cause changes in local neural oscillatory power. Despite an increased use in cognitive and clinical neuroscience, the fundamental mechanisms of tACS are still not fully understood. Here, we develop a computational neuronal network model of two-compartment pyramidal neurons and inhibitory interneurons which mimic the local cortical circuits. We model tACS with electric field strengths that are achievable in human applications. We then simulate intrinsic network activity and measure neural entrainment to investigate how tACS modulates ongoing endogenous oscillations. First, we show that intensity-specific effects of tACS are non-linear. At low intensities (<0.3 mV/mm), tACS desynchronizes neural firing relative to the endogenous oscillations. At higher intensities (>0.3 mV/mm), neurons are entrained to the exogenous electric field. We then further explore the stimulation parameter space and find that entrainment of ongoing cortical oscillations also depends on frequency by following an Arnold tongue. Moreover, neuronal networks can amplify the tACS induced entrainment via excitation-inhibition balance. Our model shows that pyramidal neurons are directly entrained by the exogenous electric field and drive the inhibitory neurons. Our findings can thus provide a mechanistic framework for understanding the intensity- and frequency- specific effects of oscillating electric fields on neuronal networks. This is crucial for rational parameters selection for tACS in cognitive studies and clinical applications.

Fronto-central resting-state 15-29 Hz transient beta events change with therapeutic transcranial magnetic stimulation for posttraumatic stress disorder and major depressive disorder

Repetitive transcranial magnetic stimulation (rTMS) is an established treatment for major depressive disorder (MDD) and shows promise for posttraumatic stress disorder (PTSD), yet effectiveness varies. Electroencephalography (EEG) can identify rTMS-associated brain changes. EEG oscillations are often examined using averaging approaches that mask finer time-scale dynamics. Recent advances show some brain oscillations emerge as transient increases in power, a phenomenon termed "Spectral Events," and that event characteristics correspond with cognitive functions. We applied Spectral Event analyses to identify potential EEG biomarkers of effective rTMS treatment. Resting 8-electrode EEG was collected from 23 patients with MDD and PTSD before and after 5 Hz rTMS targeting the left dorsolateral prefrontal cortex. Using an open-source toolbox ( https://github.com/jonescompneurolab/SpectralEvents ), we quantified event features and tested for treatment associated changes. Spectral Events in delta/theta (1-6 Hz), alpha (7-14 Hz), and beta (15-29 Hz) bands occurred in all patients. rTMS-induced improvement in comorbid MDD PTSD were associated with pre- to post-treatment changes in fronto-central electrode beta event features, including frontal beta event frequency spans and durations, and central beta event maxima power. Furthermore, frontal pre-treatment beta event duration correlated negatively with MDD symptom improvement. Beta events may provide new biomarkers of clinical response and advance the understanding of rTMS.

Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference

Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.

Spatiotemporal organisation of human sensorimotor beta burst activity

Beta oscillations in human sensorimotor cortex are hallmark signatures of healthy and pathological movement. In single trials, beta oscillations include bursts of intermittent, transient periods of high-power activity. These burst events have been linked to a range of sensory and motor processes, but their precise spatial, spectral, and temporal structure remains unclear. Specifically, a role for beta burst activity in information coding and communication suggests spatiotemporal patterns, or travelling wave activity, along specific anatomical gradients. We here show in human magnetoencephalography recordings that burst activity in sensorimotor cortex occurs in planar spatiotemporal wave-like patterns that dominate along two axes either parallel or perpendicular to the central sulcus. Moreover, we find that the two propagation directions are characterised by distinct anatomical and physiological features. Finally, our results suggest that sensorimotor beta bursts occurring before and after a movement can be distinguished by their anatomical, spectral, and spatiotemporal characteristics, indicating distinct functional roles.

Developmental trajectory of transmission speed in the human brain

The structure of the human connectome develops from childhood throughout adolescence to middle age, but how these structural changes affect the speed of neuronal signaling is not well described. In 74 subjects, we measured the latency of cortico-cortical evoked responses across association and U-fibers and calculated their corresponding transmission speeds. Decreases in conduction delays until at least 30 years show that the speed of neuronal communication develops well into adulthood.

Laminar Specificity of the Auditory Perceptual Awareness Negativity: A Biophysical Modeling Study

How perception of sensory stimuli emerges from brain activity is a fundamental question of neuroscience. To date, two disparate lines of research have examined this question. On one hand, human neuroimaging studies have helped us understand the large-scale brain dynamics of perception. On the other hand, work in animal models (mice, typically) has led to fundamental insight into the micro-scale neural circuits underlying perception. However, translating such fundamental insight from animal models to humans has been challenging. Here, using biophysical modeling, we show that the auditory awareness negativity (AAN), an evoked response associated with perception of target sounds in noise, can be accounted for by synaptic input to the supragranular layers of auditory cortex (AC) that is present when target sounds are heard but absent when they are missed. This additional input likely arises from cortico-cortical feedback and/or non-lemniscal thalamic projections and targets the apical dendrites of layer-V pyramidal neurons (PNs). In turn, this leads to increased local field potential activity, increased spiking activity in layer-V PNs, and the AAN. The results are consistent with current cellular models of conscious processing and help bridge the gap between the macro and micro levels of perception-related brain activity.

Author Summary

To date, our understanding of the brain basis of conscious perception has mostly been restricted to large-scale, network-level activity that can be measured non-invasively in human subjects. However, we lack understanding of how such network-level activity is supported by individual neurons and neural circuits. This is at least partially because conscious perception is difficult to study in experimental animals, where such detailed characterization of neural activity is possible. To address this gap, we used biophysical modeling to gain circuit-level insight into an auditory brain response known as the auditory awareness negativity (AAN). This response can be recorded non-invasively in humans and is associated with perceptual awareness of sounds of interest. Our model shows that the AAN likely arises from specific cortical layers and cell types. These data help bridge the gap between circuit- and network-level theories of consciousness, and could lead to new, targeted treatments for perceptual dysfunction and disorders of consciousness.

Full-scale scaffold model of the human hippocampus CA1 area

The increasing availability of quantitative data on the human brain is opening new avenues to study neural function and dysfunction, thus bringing us closer and closer to the implementation of digital twin applications for personalized medicine. Here we provide a resource to the neuroscience community: a computational method to generate full-scale scaffold model of human brain regions starting from microscopy images. We have benchmarked the method to reconstruct the CA1 region of a right human hippocampus, which accounts for about half of the entire right hippocampal formation. Together with 3D soma positioning we provide a connectivity matrix generated using a morpho-anatomical connection strategy based on axonal and dendritic probability density functions accounting for morphological properties of hippocampal neurons. The data and algorithms are supplied in a ready-to-use format, suited to implement computational models at different scales and detail.

NRN-EZ: an application to streamline biophysical modeling of synaptic integration using NEURON

One of the fundamental goals in neuroscience is to determine how the brain processes information and ultimately controls the execution of complex behaviors. Over the past four decades, there has been a steady growth in our knowledge of the morphological and functional diversity of neurons, the building blocks of the brain. These cells clearly differ not only for their anatomy and ion channel distribution, but also for the type, strength, location, and temporal pattern of activity of the many synaptic inputs they receive. Compartmental modeling programs like NEURON have become widely used in the neuroscience community to address a broad range of research questions, including how neurons integrate synaptic inputs and propagate information through complex neural networks. One of the main strengths of NEURON is its ability to incorporate user-defined information about the realistic morphology and biophysical properties of different cell types. Although the graphical user interface of the program can be used to run initial exploratory simulations, introducing a stochastic representation of synaptic weights, locations and activation times typically requires users to develop their own codes, a task that can be overwhelming for some beginner users. Here we describe NRN-EZ, an interactive application that allows users to specify complex patterns of synaptic input activity that can be integrated as part of NEURON simulations. Through its graphical user interface, NRN-EZ aims to ease the learning curve to run computational models in NEURON, for users that do not necessarily have a computer science background.

FiNN: A toolbox for neurophysiological network analysis

Recently, neuroscience has seen a shift from localist approaches to network-wide investigations of brain function. Neurophysiological signals across different spatial and temporal scales provide insight into neural communication. However, additional methodological considerations arise when investigating network-wide brain dynamics rather than local effects. Specifically, larger amounts of data, investigated across a higher dimensional space, are necessary. Here, we present FiNN (Find Neurophysiological Networks), a novel toolbox for the analysis of neurophysiological data with a focus on functional and effective connectivity. FiNN provides a wide range of data processing methods and statistical and visualization tools to facilitate inspection of connectivity estimates and the resulting metrics of brain dynamics. The Python toolbox and its documentation are freely available as Supporting Information. We evaluated FiNN against a number of established frameworks on both a conceptual and an implementation level. We found FiNN to require much less processing time and memory than other toolboxes. In addition, FiNN adheres to a design philosophy of easy access and modifiability, while providing efficient data processing implementations. Since the investigation of network-level neural dynamics is experiencing increasing interest, we place FiNN at the disposal of the neuroscientific community as open-source software.

Non-zero mean alpha oscillations revealed with computational model and empirical data

Ongoing oscillations and evoked responses are two main types of neuronal activity obtained with diverse electrophysiological recordings (EEG/MEG/iEEG/LFP). Although typically studied separately, they might in fact be closely related. One possibility to unite them is to demonstrate that neuronal oscillations have non-zero mean which predicts that stimulus- or task-triggered amplitude modulation of oscillations can contribute to the generation of evoked responses. We validated this mechanism using computational modelling and analysis of a large EEG data set. With a biophysical model, we indeed demonstrated that intracellular currents in the neuron are asymmetric and, consequently, the mean of alpha oscillations is non-zero. To understand the effect that neuronal currents exert on oscillatory mean, we varied several biophysical and morphological properties of neurons in the network, such as voltage-gated channel densities, length of dendrites, and intensity of incoming stimuli. For a very large range of model parameters, we observed evidence for non-zero mean of oscillations. Complimentary, we analysed empirical rest EEG recordings of 90 participants (50 young, 40 elderly) and, with spatio-spectral decomposition, detected at least one spatially-filtred oscillatory component of non-zero mean alpha oscillations in 93% of participants. In order to explain a complex relationship between the dynamics of amplitude-envelope and corresponding baseline shifts, we performed additional simulations with simple oscillators coupled with different time delays. We demonstrated that the extent of spatial synchronisation may obscure macroscopic estimation of alpha rhythm modulation while leaving baseline shifts unchanged. Overall, our results predict that amplitude modulation of neural oscillations should at least partially explain the generation of evoked responses. Therefore, inference about changes in evoked responses with respect to cognitive conditions, age or neuropathologies should be constructed while taking into account oscillatory neuronal dynamics.

Detecting Spontaneous Neural Oscillation Events in Primate Auditory Cortex

Electrophysiological oscillations in the brain have been shown to occur as multicycle events, with onset and offset dependent on behavioral and cognitive state. To provide a baseline for state-related and task-related events, we quantified oscillation features in resting-state recordings. We developed an open-source wavelet-based tool to detect and characterize such oscillation events (OEvents) and exemplify the use of this tool in both simulations and two invasively-recorded electrophysiology datasets: one from human, and one from nonhuman primate (NHP) auditory system. After removing incidentally occurring event-related potentials (ERPs), we used OEvents to quantify oscillation features. We identified ∼2 million oscillation events, classified within traditional frequency bands: δ, θ, α, β, low γ, γ, and high γ. Oscillation events of 1-44 cycles could be identified in at least one frequency band 90% of the time in human and NHP recordings. Individual oscillation events were characterized by nonconstant frequency and amplitude. This result necessarily contrasts with prior studies which assumed frequency constancy, but is consistent with evidence from event-associated oscillations. We measured oscillation event duration, frequency span, and waveform shape. Oscillations tended to exhibit multiple cycles per event, verifiable by comparing filtered to unfiltered waveforms. In addition to the clear intraevent rhythmicity, there was also evidence of interevent rhythmicity within bands, demonstrated by finding that coefficient of variation of interval distributions and Fano factor (FF) measures differed significantly from a Poisson distribution assumption. Overall, our study provides an easy-to-use tool to study oscillation events at the single-trial level or in ongoing recordings, and demonstrates that rhythmic, multicycle oscillation events dominate auditory cortical dynamics.

Modernizing the NEURON Simulator for Sustainability, Portability, and Performance

The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.

Methodological considerations for studying neural oscillations

Neural oscillations are ubiquitous across recording methodologies and species, broadly associated with cognitive tasks, and amenable to computational modelling that investigates neural circuit generating mechanisms and neural population dynamics. Because of this, neural oscillations offer an exciting potential opportunity for linking theory, physiology and mechanisms of cognition. However, despite their prevalence, there are many concerns-new and old-about how our analysis assumptions are violated by known properties of field potential data. For investigations of neural oscillations to be properly interpreted, and ultimately developed into mechanistic theories, it is necessary to carefully consider the underlying assumptions of the methods we employ. Here, we discuss seven methodological considerations for analysing neural oscillations. The considerations are to (1) verify the presence of oscillations, as they may be absent; (2) validate oscillation band definitions, to address variable peak frequencies; (3) account for concurrent non-oscillatory aperiodic activity, which might otherwise confound measures; measure and account for (4) temporal variability and (5) waveform shape of neural oscillations, which are often bursty and/or nonsinusoidal, potentially leading to spurious results; (6) separate spatially overlapping rhythms, which may interfere with each other; and (7) consider the required signal-to-noise ratio for obtaining reliable estimates. For each topic, we provide relevant examples, demonstrate potential errors of interpretation, and offer suggestions to address these issues. We primarily focus on univariate measures, such as power and phase estimates, though we discuss how these issues can propagate to multivariate measures. These considerations and recommendations offer a helpful guide for measuring and interpreting neural oscillations.

Training a spiking neuronal network model of visual-motor cortex to play a virtual racket-ball game using reinforcement learning

Recent models of spiking neuronal networks have been trained to perform behaviors in static environments using a variety of learning rules, with varying degrees of biological realism. Most of these models have not been tested in dynamic visual environments where models must make predictions on future states and adjust their behavior accordingly. The models using these learning rules are often treated as black boxes, with little analysis on circuit architectures and learning mechanisms supporting optimal performance. Here we developed visual/motor spiking neuronal network models and trained them to play a virtual racket-ball game using several reinforcement learning algorithms inspired by the dopaminergic reward system. We systematically investigated how different architectures and circuit-motifs (feed-forward, recurrent, feedback) contributed to learning and performance. We also developed a new biologically-inspired learning rule that significantly enhanced performance, while reducing training time. Our models included visual areas encoding game inputs and relaying the information to motor areas, which used this information to learn to move the racket to hit the ball. Neurons in the early visual area relayed information encoding object location and motion direction across the network. Neuronal association areas encoded spatial relationships between objects in the visual scene. Motor populations received inputs from visual and association areas representing the dorsal pathway. Two populations of motor neurons generated commands to move the racket up or down. Model-generated actions updated the environment and triggered reward or punishment signals that adjusted synaptic weights so that the models could learn which actions led to reward. Here we demonstrate that our biologically-plausible learning rules were effective in training spiking neuronal network models to solve problems in dynamic environments. We used our models to dissect the circuit architectures and learning rules most effective for learning. Our model shows that learning mechanisms involving different neural circuits produce similar performance in sensory-motor tasks. In biological networks, all learning mechanisms may complement one another, accelerating the learning capabilities of animals. Furthermore, this also highlights the resilience and redundancy in biological systems.

A visual encoding model links magnetoencephalography signals to neural synchrony in human cortex

Synchronization of neuronal responses over large distances is hypothesized to be important for many cortical functions. However, no straightforward methods exist to estimate synchrony non-invasively in the living human brain. MEG and EEG measure the whole brain, but the sensors pool over large, overlapping cortical regions, obscuring the underlying neural synchrony. Here, we developed a model from stimulus to cortex to MEG sensors to disentangle neural synchrony from spatial pooling of the instrument. We find that synchrony across cortex has a surprisingly large and systematic effect on predicted MEG spatial topography. We then conducted visual MEG experiments and separated responses into stimulus-locked and broadband components. The stimulus-locked topography was similar to model predictions assuming synchronous neural sources, whereas the broadband topography was similar to model predictions assuming asynchronous sources. We infer that visual stimulation elicits two distinct types of neural responses, one highly synchronous and one largely asynchronous across cortex.

Technical note: A fast and robust integrator of delay differential equations in DCM for electrophysiological data

Dynamic causal models (DCMs) of electrophysiological data allow, in principle, for inference on hidden, bulk synaptic function in neural circuits. The directed influences between the neuronal elements of modeled circuits are subject to delays due to the finite transmission speed of axonal connections. Ordinary differential equations are therefore not adequate to capture the ensuing circuit dynamics, and delay differential equations (DDEs) are required instead. Previous work has illustrated that the integration of DDEs in DCMs benefits from sophisticated integration schemes in order to ensure rigorous parameter estimation and correct model identification. However, integration schemes that have been proposed for DCMs either emphasize speed (at the possible expense of accuracy) or robustness (but with computational costs that are problematic in practice). In this technical note, we propose an alternative integration scheme that overcomes these shortcomings and offers high computational efficiency while correctly preserving the nature of delayed effects. This integration scheme is available as open-source code in the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) toolbox and can be easily integrated into existing software (SPM) for the analysis of DCMs for electrophysiological data. While this paper focuses on its application to the convolution-based formalism of DCMs, the new integration scheme can be equally applied to more advanced formulations of DCMs (e.g. conductance based models). Our method provides a new option for electrophysiological DCMs that offers the speed required for scientific projects, but also the accuracy required for rigorous translational applications, e.g. in computational psychiatry.

Laminar dynamics of high amplitude beta bursts in human motor cortex

Motor cortical activity in the beta frequency range is one of the strongest and most studied movement-related neural signals. At the single trial level, beta band activity is often characterized by transient, high amplitude, bursting events rather than slowly modulating oscillations. The timing of these bursting events is tightly linked to behavior, suggesting a more dynamic functional role for beta activity than previously believed. However, the neural mechanisms underlying beta bursts in sensorimotor circuits are poorly understood. To address this, we here leverage and extend recent developments in high precision MEG for temporally resolved laminar analysis of burst activity, combined with a neocortical circuit model that simulates the biophysical generators of the electrical currents which drive beta bursts. This approach pinpoints the generation of beta bursts in human motor cortex to distinct excitatory synaptic inputs to deep and superficial cortical layers, which drive current flow in opposite directions. These laminar dynamics of beta bursts in motor cortex align with prior invasive animal recordings within the somatosensory cortex, and suggest a conserved mechanism for somatosensory and motor cortical beta bursts. More generally, we demonstrate the ability for uncovering the laminar dynamics of event-related neural signals in human non-invasive recordings. This provides important constraints to theories about the functional role of burst activity for movement control in health and disease, and crucial links between macro-scale phenomena measured in humans and micro-circuit activity recorded from animal models.

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.

Targeting neural oscillations with transcranial alternating current stimulation

Neural oscillations at the network level synchronize activity between regions and temporal scales. Transcranial alternating current stimulation (tACS), the delivery of low-amplitude electric current to the scalp, provides a tool for investigating the causal role of neural oscillations in cognition. The parameter space for tACS is vast and optimization is required in terms of temporal and spatial targeting. We review emerging techniques and suggest novel approaches that capitalize on the non-sinusoidal and transient nature of neural oscillations and leverage the flexibility provided by a customizable electrode montage and electrical waveform. The customizability and safety profile of tACS make it a promising tool for precision intervention in psychiatric illnesses.

The impact of artifact removal approaches on TMS-EEG signal

Transcranial magnetic stimulation (TMS)-evoked potentials (TEPs) allow one to assess cortical excitability and effective connectivity in clinical and basic research. However, obtaining clean TEPs is challenging due to the various TMS-related artifacts that contaminate the electroencephalographic (EEG) signal when the TMS pulse is delivered. Different preprocessing approaches have been employed to remove the artifacts, but the degree of artifact reduction or signal distortion introduced in this phase of analysis is still unknown. Knowing and controlling this potential source of uncertainty will increase the inter-rater reliability of TEPs and improve the comparability between TMS-EEG studies. The goal of this study was to assess the variability in TEP waveforms due to of the use of different preprocessing pipelines. To accomplish this aim, we preprocessed the same TMS-EEG data with four different pipelines and compared the results. The dataset was obtained from 16 subjects in two identical recording sessions, each session consisting of both left dorsolateral prefrontal cortex and left inferior parietal lobule stimulation at 100% of the resting motor threshold. Considerable differences in TEP amplitudes and global mean field power (GMFP) were found between the preprocessing pipelines. Topographies of TEPs from the different pipelines were all highly correlated (ρ>0.8) at latencies over 100 ms. By contrast, waveforms at latencies under 100 ms showed a variable level of correlation, with ρ ranging between 0.2 and 0.9. Moreover, the test-retest reliability of TEPs depended on the preprocessing pipeline. Taken together, these results take us to suggest that the choice of the preprocessing approach has a marked impact on the final TEP, and that further studies are needed to understand advantages and disadvantages of the different approaches.

Whole-head OPM-MEG enables noninvasive assessment of functional connectivity

A recent paper by Boto et al. established reliability of noninvasive functional connectivity measurements with a new whole-head optically pumped magnetometer magnetoencephalography (OPM-MEG) system. This rapidly developing technology enables a conformal sensor array to accommodate different head sizes and opens up new avenues for experiments in more naturalistic settings.

Neural Mechanisms Underlying Human Auditory Evoked Responses Revealed By Human Neocortical Neurosolver

Auditory evoked fields (AEFs) are commonly studied, yet their underlying neural mechanisms remain poorly understood. Here, we used the biophysical modelling software Human Neocortical Neurosolver (HNN) whose foundation is a canonical neocortical circuit model to interpret the cell and network mechanisms contributing to macroscale AEFs elicited by a simple tone, measured with magnetoencephalography. We found that AEFs can be reproduced by activating the neocortical circuit through a layer specific sequence of feedforward and feedback excitatory synaptic drives, similar to prior simulation of somatosensory evoked responses, supporting the notion that basic structures and activation patterns are preserved across sensory regions. We also applied the modeling framework to develop and test predictions on neural mechanisms underlying AEF differences in the left and right hemispheres, as well as in hemispheres contralateral and ipsilateral to the presentation of the auditory stimulus. We found that increasing the strength of the excitatory synaptic cortical feedback inputs to supragranular layers simulates the commonly observed right hemisphere dominance, while decreasing the input latencies and simultaneously increasing the number of cells contributing to the signal accounted for the contralateral dominance. These results provide a direct link between human data and prior animal studies and lay the foundation for future translational research examining the mechanisms underlying alteration in this fundamental biomarker of auditory processing in healthy cognition and neuropathology.

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.

Decoding and interpreting cortical signals with a compact convolutional neural network

OBJECTIVE

Brain-computer interfaces (BCIs) decode information from neural activity and send it to external devices. The use of Deep Learning approaches for decoding allows for automatic feature engineering within the specific decoding task. Physiologically plausible interpretation of the network parameters ensures the robustness of the learned decision rules and opens the exciting opportunity for automatic knowledge discovery.

APPROACH

We describe a compact convolutional network-based architecture for adaptive decoding of electrocorticographic (ECoG) data into finger kinematics. We also propose a novel theoretically justified approach to interpreting the spatial and temporal weights in the architectures that combine adaptation in both space and time. The obtained spatial and frequency patterns characterizing the neuronal populations pivotal to the specific decoding task can then be interpreted by fitting appropriate spatial and dynamical models.

MAIN RESULTS

We first tested our solution using realistic Monte-Carlo simulations. Then, when applied to the ECoG data from Berlin BCI competition IV dataset, our architecture performed comparably to the competition winners without requiring explicit feature engineering. Using the proposed approach to the network weights interpretation we could unravel the spatial and the spectral patterns of the neuronal processes underlying the successful decoding of finger kinematics from an ECoG dataset. Finally we have also applied the entire pipeline to the analysis of a 32-channel EEG motor-imagery dataset and observed physiologically plausible patterns specific to the task.

SIGNIFICANCE

We described a compact and interpretable CNN architecture derived from the basic principles and encompassing the knowledge in the field of neural electrophysiology. For the first time in the context of such multibranch architectures with factorized spatial and temporal processing we presented theoretically justified weights interpretation rules. We verified our recipes using simulations and real data and demonstrated that the proposed solution offers a good decoder and a tool for investigating motor control neural mechanisms.

Brain Modeling ToolKit: An open source software suite for multiscale modeling of brain circuits

Experimental studies in neuroscience are producing data at a rapidly increasing rate, providing exciting opportunities and formidable challenges to existing theoretical and modeling approaches. To turn massive datasets into predictive quantitative frameworks, the field needs software solutions for systematic integration of data into realistic, multiscale models. Here we describe the Brain Modeling ToolKit (BMTK), a software suite for building models and performing simulations at multiple levels of resolution, from biophysically detailed multi-compartmental, to point-neuron, to population-statistical approaches. Leveraging the SONATA file format and existing software such as NEURON, NEST, and others, BMTK offers a consistent user experience across multiple levels of resolution. It permits highly sophisticated simulations to be set up with little coding required, thus lowering entry barriers to new users. We illustrate successful applications of BMTK to large-scale simulations of a cortical area. BMTK is an open-source package provided as a resource supporting modeling-based discovery in the community.

Models of communication and control for brain networks: distinctions, convergence, and future outlook

Recent advances in computational models of signal propagation and routing in the human brain have underscored the critical role of white-matter structure. A complementary approach has utilized the framework of network control theory to better understand how white matter constrains the manner in which a region or set of regions can direct or control the activity of other regions. Despite the potential for both of these approaches to enhance our understanding of the role of network structure in brain function, little work has sought to understand the relations between them. Here, we seek to explicitly bridge computational models of communication and principles of network control in a conceptual review of the current literature. By drawing comparisons between communication and control models in terms of the level of abstraction, the dynamical complexity, the dependence on network attributes, and the interplay of multiple spatiotemporal scales, we highlight the convergence of and distinctions between the two frameworks. Based on the understanding of the intertwined nature of communication and control in human brain networks, this work provides an integrative perspective for the field and outlines exciting directions for future work.

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.