Intrinsic dendritic filtering gives low-pass power spectra of local field potentials

Large-scale interactions in predictive processing: oscillatory versus transient dynamics

How do the two main types of neural dynamics, aperiodic transients and oscillations, contribute to the interactions between feedforward (FF) and feedback (FB) pathways in sensory inference and predictive processing? We discuss three theoretical perspectives. First, we critically evaluate the theory that gamma and alpha/beta rhythms play a role in classic hierarchical predictive coding (HPC) by mediating FF and FB communication, respectively. Second, we outline an alternative functional model in which rapid sensory inference is mediated by aperiodic transients, whereas oscillations contribute to the stabilization of neural representations over time and plasticity processes. Third, we propose that the strong dependence of oscillations on predictability can be explained based on a biologically plausible alternative to classic HPC, namely dendritic HPC.

Strategies to decipher neuron identity from extracellular recordings in the cerebellum of behaving non-human primates

Identification of neuron type is critical to understand computation in neural circuits through extracellular recordings in awake, behaving animal subjects. Yet, modern recording probes have limited power to resolve neuron type. Here, we leverage the well-characterized architecture of the cerebellar circuit to perform expert identification of neuron type from extracellular recordings in behaving non-human primates. Using deep-learning classifiers we evaluate the information contained in readily accessible extracellular features for neuron identification. Waveform, discharge statistics, anatomical layer, and functional interactions each can inform neuron labels for a sizable fraction of cerebellar units. Together, as inputs to a deep-learning classifier, the features perform even better. Our tools and methodologies, validated during smooth pursuit eye movements in the cerebellar floccular complex of awake behaving monkeys, can guide expert identification of neuron type during cerebellar-dependent tasks in behaving animals across species. They lay the groundwork for characterization of information processing in the cerebellar cortex.

Single neuron contributions to the auditory brainstem EEG

The auditory brainstem response (ABR) is an acoustically evoked EEG potential that is an important diagnostic tool for hearing loss, especially in newborns. The ABR originates from the response sequence of auditory nerve and brainstem nuclei, and a click-evoked ABR typically shows three positive peaks ('waves') within the first six milliseconds. However, an assignment of the waves of the ABR to specific sources is difficult, and a quantification of contributions to the ABR waves is not available. Here, we exploit the large size and physical separation of the barn owl first-order cochlear nucleus magnocellularis (NM) to estimate single-cell contributions to the ABR. We simultaneously recorded NM neurons' spikes and the EEG, and found that ≳ 5, 000 spontaneous single-cell spikes are necessary to isolate a significant spike-triggered average response at the EEG electrode. An average single-neuron contribution to the ABR was predicted by convolving the spike-triggered average with the cell's peri-stimulus time histogram. Amplitudes of predicted contributions of single NM cells typically reached 32.9 ± 1.1 nV (mean ± SE, range: 2.5-162.7 nV), or 0.07±0.02% (median ± SE; range from 0.01% to 1%) of the ABR amplitude. The time of the predicted peak coincided best with the peak of the ABR wave II, independent of the click sound level. Our results suggest that individual neurons' contributions to an EEG can vary widely, and that wave II of the ABR is shaped by NM units.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. The contributions from different cortical layers within V1 could however not be robustly separated and identified with LPA. This is likely due to substantial synchrony in population firing rates across layers, which may be reduced with other stimulus protocols in the future. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Filtered Point Processes Tractably Capture Rhythmic And Broadband Power Spectral Structure in Neural Electrophysiological Recordings

Neural electrophysiological recordings arise from interacting rhythmic (oscillatory) and broadband (aperiodic) biological subprocesses. Both rhythmic and broadband processes contribute to the neural power spectrum, which decomposes the variance of a neural recording across frequencies. Although an extensive body of literature has successfully studied rhythms in various diseases and brain states, researchers only recently have systematically studied the characteristics of broadband effects in the power spectrum. Broadband effects can generally be categorized as 1) shifts in power across all frequencies, which correlate with changes in local firing rates and 2) changes in the overall shape of the power spectrum, such as the spectral slope or power law exponent. Shape changes are evident in various conditions and brain states, influenced by factors such as excitation to inhibition balance, age, and various diseases. It is increasingly recognized that broadband and rhythmic effects can interact on a sub-second timescale. For example, broadband power is time-locked to the phase of <1 Hz rhythms in propofol induced unconsciousness. Modeling tools that explicitly deal with both rhythmic and broadband contributors to the power spectrum and that capture their interactions are essential to help improve the interpretability of power spectral effects. Here, we introduce a tractable stochastic forward modeling framework designed to capture both narrowband and broadband spectral effects when prior knowledge or theory about the primary biophysical processes involved is available. Population-level neural recordings are modeled as the sum of filtered point processes (FPPs), each representing the contribution of a different biophysical process such as action potentials or postsynaptic potentials of different types. Our approach builds on prior neuroscience FPP work by allowing multiple interacting processes and time-varying firing rates and by deriving theoretical power spectra and cross-spectra. We demonstrate several properties of the models, including that they divide the power spectrum into frequency ranges dominated by rhythmic and broadband effects, and that they can capture spectral effects across multiple timescales, including sub-second cross-frequency coupling. The framework can be used to interpret empirically observed power spectra and cross-frequency coupling effects in biophysical terms, which bridges the gap between theoretical models and experimental results.

Convergent effects of different anesthetics on changes in phase alignment of cortical oscillations

Many anesthetics cause loss of responsiveness despite having diverse underlying molecular and circuit actions. To explore the convergent effects of these drugs, we examined how anesthetic doses of ketamine and dexmedetomidine affected oscillations in the prefrontal cortex of nonhuman primates. Both anesthetics caused increases in phase locking in the ventrolateral and dorsolateral prefrontal cortex, within and across hemispheres. However, the nature of the phase locking varied. Activity in different subregions within a hemisphere became more anti-phase with both drugs. Local analyses within a region suggested that this finding could be explained by broad cortical distance-based effects, such as large traveling waves. By contrast, homologous areas across hemispheres became more in-phase. Our results suggest that both anesthetics induce strong patterns of cortical phase alignment that are markedly different from those in the awake state, and that these patterns may be a common feature driving loss of responsiveness from different anesthetic drugs.

Reconciliation of weak pairwise spike-train correlations and highly coherent local field potentials across space

Multi-electrode arrays covering several square millimeters of neural tissue provide simultaneous access to population signals such as extracellular potentials and spiking activity of one hundred or more individual neurons. The interpretation of the recorded data calls for multiscale computational models with corresponding spatial dimensions and signal predictions. Multi-layer spiking neuron network models of local cortical circuits covering about $1\,{\text{mm}^{2}}$ have been developed, integrating experimentally obtained neuron-type-specific connectivity data and reproducing features of observed in-vivo spiking statistics. Local field potentials can be computed from the simulated spiking activity. We here extend a local network and local field potential model to an area of $4\times 4\,{\text{mm}^{2}}$, preserving the neuron density and introducing distance-dependent connection probabilities and conduction delays. We find that the upscaling procedure preserves the overall spiking statistics of the original model and reproduces asynchronous irregular spiking across populations and weak pairwise spike-train correlations in agreement with experimental recordings from sensory cortex. Also compatible with experimental observations, the correlation of local field potential signals is strong and decays over a distance of several hundred micrometers. Enhanced spatial coherence in the low-gamma band around $50\,\text{Hz}$ may explain the recent report of an apparent band-pass filter effect in the spatial reach of the local field potential.

Synaptic sensitization in the anterior cingulate cortex sustains the consciousness of pain via synchronized oscillating electromagnetic waves

A recent report showed that experiencing pain requires not only activities in the brain, but also the generation of electric fields in a defined area of the anterior cingulate cortex (ACC). The present manuscript presents evidence that electromagnetic (EM) waves are also necessary. Action potentials (APs) encoding information about an injury stimulate thousands synapses on pyramidal neurons within the ACC resulting in the generation of synchronized oscillating (EM) waves and the activation of NMDA receptors. The latter induces a long-term potentiation (LTP) in the pyramidal dendrites that is necessary to experience both neuropathic and visceral pain. The LTP sensitizes transmission across the synapses that sustains the duration of the waves and the pain, EM waves containing information about the injury travel throughout the brain and studies using transcranial stimulation indicate that they can induce NMDA-mediated LTP in distant neuronal circuits. What is ultimately experienced as pain depends on the almost instantaneous integration of information from numerous neuronal centers, such as the amygdala, that are widely separated in the brain. These centers also generate EM waves and I propose that the EM waves from these centers interact to rapidly adjust the intensity of the pain to accommodate past and present circumstances. Where the waves are transformed into a consciousness of pain is unknown. One possibility is the mind which, according to contemporary theories, is where conscious experiences arise. The hypothesis can be tested directly by blocking the waves from the ACC. If correct, the waves would open new avenues of research into the relationship between the brain, consciousness, and the mind.

Whole-brain simulation of interictal epileptic discharges for patient-specific interpretation of interictal SEEG data

In patients with refractory epilepsy, the clinical interpretation of stereoelectroencephalographic (SEEG) signals is crucial to delineate the epileptogenic network that should be targeted by surgery. We propose a pipeline of patient-specific computational modeling of interictal epileptic activity to improve the definition of regions of interest. Comparison between the computationally defined regions of interest and the resected region confirmed the efficiency of the pipeline. This result suggests that computational modeling can be used to reconstruct signals and aid clinical interpretation.

Multitask learning of a biophysically-detailed neuron model

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

Aperiodic components of local field potentials reflect inherent differences between cortical and subcortical activity

Information flow in brain networks is reflected in local field potentials that have both periodic and aperiodic components. The 1/fχ aperiodic component of the power spectra tracks arousal and correlates with other physiological and pathophysiological states. Here we explored the aperiodic activity in the human thalamus and basal ganglia in relation to simultaneously recorded cortical activity. We elaborated on the parameterization of the aperiodic component implemented by specparam (formerly known as FOOOF) to avoid parameter unidentifiability and to obtain independent and more easily interpretable parameters. This allowed us to seamlessly fit spectra with and without an aperiodic knee, a parameter that captures a change in the slope of the aperiodic component. We found that the cortical aperiodic exponent χ, which reflects the decay of the aperiodic component with frequency, is correlated with Parkinson's disease symptom severity. Interestingly, no aperiodic knee was detected from the thalamus, the pallidum, or the subthalamic nucleus, which exhibited an aperiodic exponent significantly lower than in cortex. These differences were replicated in epilepsy patients undergoing intracranial monitoring that included thalamic recordings. The consistently lower aperiodic exponent and lack of an aperiodic knee from all subcortical recordings may reflect cytoarchitectonic and/or functional differences.

SIGNIFICANCE STATEMENT

The aperiodic component of local field potentials can be modeled to produce useful and reproducible indices of neural activity. Here we refined a widely used phenomenological model for extracting aperiodic parameters (namely the exponent, offset and knee), with which we fit cortical, basal ganglia, and thalamic intracranial local field potentials, recorded from unique cohorts of movement disorders and epilepsy patients. We found that the aperiodic exponent in motor cortex is higher in Parkinson's disease patients with more severe motor symptoms, suggesting that aperiodic features may have potential as electrophysiological biomarkers for movement disorders symptoms. Remarkably, we found conspicuous differences in the aperiodic parameters of basal ganglia and thalamic signals compared to those from neocortex.

Aperiodic EEG and 7T MRSI evidence for maturation of E/I balance supporting the development of working memory through adolescence

Adolescence has been hypothesized to be a critical period for the development of human association cortex and higher-order cognition. A defining feature of critical period development is a shift in the excitation: inhibition (E/I) balance of neural circuitry, however how changes in E/I may enhance cortical circuit function to support maturational improvements in cognitive capacities is not known. Harnessing ultra-high field 7 T MR spectroscopy and EEG in a large, longitudinal cohort of youth (N = 164, ages 10-32 years old, 347 neuroimaging sessions), we delineate biologically specific associations between age-related changes in excitatory glutamate and inhibitory GABA neurotransmitters and EEG-derived measures of aperiodic neural activity reflective of E/I balance in prefrontal association cortex. Specifically, we find that developmental increases in E/I balance reflected in glutamate:GABA balance are linked to changes in E/I balance assessed by the suppression of prefrontal aperiodic activity, which in turn facilitates robust improvements in working memory. These findings indicate a role for E/I-engendered changes in prefrontal signaling mechanisms in the maturation of cognitive maintenance. More broadly, this multi-modal imaging study provides evidence that human association cortex undergoes physiological changes consistent with critical period plasticity during adolescence.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

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.

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.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

Macroscopic resting-state brain dynamics are best described by linear models

It is typically assumed that large networks of neurons exhibit a large repertoire of nonlinear behaviours. Here we challenge this assumption by leveraging mathematical models derived from measurements of local field potentials via intracranial electroencephalography and of whole-brain blood-oxygen-level-dependent brain activity via functional magnetic resonance imaging. We used state-of-the-art linear and nonlinear families of models to describe spontaneous resting-state activity of 700 participants in the Human Connectome Project and 122 participants in the Restoring Active Memory project. We found that linear autoregressive models provide the best fit across both data types and three performance metrics: predictive power, computational complexity and the extent of the residual dynamics unexplained by the model. To explain this observation, we show that microscopic nonlinear dynamics can be counteracted or masked by four factors associated with macroscopic dynamics: averaging over space and over time, which are inherent to aggregated macroscopic brain activity, and observation noise and limited data samples, which stem from technological limitations. We therefore argue that easier-to-interpret linear models can faithfully describe macroscopic brain dynamics during resting-state conditions.

Magnetic seizure therapy and electroconvulsive therapy increase aperiodic activity

Major depressive disorder (MDD) is a leading cause of disability worldwide. One of the most efficacious treatments for treatment-resistant MDD is electroconvulsive therapy (ECT). Recently, magnetic seizure therapy (MST) was developed as an alternative to ECT due to its more favorable side effect profile. While these approaches have been very successful clinically, the neural mechanisms underlying their therapeutic effects are unknown. For example, clinical "slowing" of the electroencephalogram beginning in the postictal state and extending days to weeks post-treatment has been observed in both treatment modalities. However, a recent longitudinal study of a small cohort of ECT patients revealed that, rather than delta oscillations, clinical slowing was better explained by increases in aperiodic activity, an emerging EEG signal linked to neural inhibition. Here we investigate the role of aperiodic activity in a cohort of patients who received ECT and a cohort of patients who received MST treatment. We find that aperiodic neural activity increases significantly in patients receiving either ECT or MST. Although not directly related to clinical efficacy in this dataset, increased aperiodic activity is linked to greater amounts of neural inhibition, which is suggestive of a potential shared neural mechanism of action across ECT and MST.

Metamodelling of a two-population spiking neural network

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

Magnetic seizure therapy and electroconvulsive therapy increase aperiodic activity

Major depressive disorder (MDD) is a leading cause of disability worldwide. One of the most efficacious treatments for treatment-resistant MDD is electroconvulsive therapy (ECT). Recently, magnetic seizure therapy (MST) was developed as an alternative to ECT due to its more favorable side effect profile. While these approaches have been very successful clinically, the neural mechanisms underlying their therapeutic effects are unknown. For example, clinical "slowing" of the electroencephalogram beginning in the postictal state and extending days to weeks post-treatment has been observed in both treatment modalities. However, a recent longitudinal study of a small cohort of ECT patients revealed that, rather than delta oscillations, clinical slowing was better explained by increases in aperiodic activity, an emerging EEG signal linked to neural inhibition. Here we investigate the role of aperiodic activity in a cohort of patients who received ECT and a cohort of patients who received MST treatment. We find that aperiodic neural activity increases significantly in patients receiving either ECT or MST. Although not directly related to clinical efficacy in this dataset, increased aperiodic activity is linked to greater amounts of neural inhibition, which is suggestive of a potential shared neural mechanism of action across ECT and MST.

In vivo ephaptic coupling allows memory network formation

It is increasingly clear that memories are distributed across multiple brain areas. Such "engram complexes" are important features of memory formation and consolidation. Here, we test the hypothesis that engram complexes are formed in part by bioelectric fields that sculpt and guide the neural activity and tie together the areas that participate in engram complexes. Like the conductor of an orchestra, the fields influence each musician or neuron and orchestrate the output, the symphony. Our results use the theory of synergetics, machine learning, and data from a spatial delayed saccade task and provide evidence for in vivo ephaptic coupling in memory representations.

Comparative analysis of multifaceted neural effects associated with varying endogenous cognitive load

Contemporary neuroscience has firmly established that mental state variation concurs with changes in neural dynamic activity in a complex way that a one-to-one mapping cannot describe. To explore the scenario of the multifaceted changes in neural dynamics associated with simple mental state variation, we took cognitive load - a common cognitive manipulation in psychology - as a venue to characterize how multiple neural dynamic features are simultaneously altered by the manipulation and how their sensitivity differs. Electroencephalogram was collected from 152 participants performing stimulus-free tasks with different demands. The results show that task demand alters wide-ranging neural dynamic features, including band-specific oscillations across broad frequency bands, scale-free dynamics, and cross-frequency phase-amplitude coupling. The scale-free dynamics outperformed others in indexing cognitive load variation. This study demonstrates a complex relationship between cognitive dynamics and neural dynamics, which points to a necessity to integrate multifaceted neural dynamic features when studying mind-brain relationship in the future.

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.

Subtractive adaptation is a more effective and general mechanism in binocular rivalry than divisive adaptation

The activity of neurons is influenced by random fluctuations and can be strongly modulated by firing rate adaptation, particularly in sensory systems. Still, there is ongoing debate about the characteristics of neuronal noise and the mechanisms of adaptation, and even less is known about how exactly they affect perception. Noise and adaptation are critical in binocular rivalry, a visual phenomenon where two images compete for perceptual dominance. Here, we investigated the effects of different noise processes and adaptation mechanisms on visual perception by simulating a model of binocular rivalry with Gaussian white noise, Ornstein-Uhlenbeck noise, and pink noise, in variants with divisive adaptation, subtractive adaptation, and without adaptation. By simulating the nine models in parameter space, we find that white noise only produces rivalry when paired with subtractive adaptation and that subtractive adaptation reduces the influence of noise intensity on rivalry strength and introduces convergence of the mean percept duration, an important metric of binocular rivalry, across all noise processes. In sum, our results show that white noise is an insufficient description of background activity in the brain and that subtractive adaptation is a stronger and more general switching mechanism in binocular rivalry than divisive adaptation, with important noise-filtering properties.

Visual cortical LFP in relation to the hippocampal theta rhythm in track running rats

Theta oscillations in the primary visual cortex (VC) have been observed during running tasks, but the mechanism behind their generation is not well understood. Some studies have suggested that theta in the VC is locally generated, while others have proposed that it is volume conducted from the hippocampus. The present study aimed to investigate the relationship between hippocampal and VC LFP dynamics. Analysis of power spectral density revealed that LFP in the VC was similar to that in the hippocampus, but with lower overall magnitude. As running velocity increased, both the power and frequency of theta and its harmonics increased in the VC, similarly to what is observed in the hippocampus. Current source density analysis triggered to theta did not identify distinct current sources and sinks in the VC, supporting the idea that theta in the VC is conducted from the adjacent hippocampus. Phase coupling between theta, its harmonics, and gamma is a notable feature in the hippocampus, particularly in the lacunosum moleculare. While some evidence of coupling between theta and its harmonics in the VC was found, bicoherence estimates did not reveal significant phase coupling between theta and gamma. Similar results were seen in the cross-region bicoherence analysis, where theta showed strong coupling with its harmonics with increasing velocity. Thus, theta oscillations observed in the VC during running tasks are likely due to volume conduction from the hippocampus.

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.

Decreased aperiodic neural activity in Parkinson's disease and dementia with Lewy bodies

Neural oscillations and signal complexity have been widely studied in neurodegenerative diseases, whereas aperiodic activity has not been explored yet in those disorders. Here, we assessed whether the study of aperiodic activity brings new insights relating to disease as compared to the conventional spectral and complexity analyses. Eyes-closed resting-state electroencephalography (EEG) was recorded in 21 patients with dementia with Lewy bodies (DLB), 28 patients with Parkinson's disease (PD), 27 patients with mild cognitive impairment (MCI) and 22 age-matched healthy controls. Spectral power was differentiated into its oscillatory and aperiodic components using the Irregularly Resampled Auto-Spectral Analysis. Signal complexity was explored using the Lempel-Ziv algorithm (LZC). We found that DLB patients showed steeper slopes of the aperiodic power component with large effect sizes compared to the controls and MCI and with a moderate effect size compared to PD. PD patients showed steeper slopes with a moderate effect size compared to controls and MCI. Oscillatory power and LZC differentiated only between DLB and other study groups and were not sensitive enough to detect differences between PD, MCI, and controls. In conclusion, both DLB and PD are characterized by alterations in aperiodic dynamics, which are more sensitive in detecting disease-related neural changes than the traditional spectral and complexity analyses. Our findings suggest that steeper aperiodic slopes may serve as a marker of network dysfunction in DLB and PD features.

Glutamatergic and GABAergic Receptor Modulation Present Unique Electrophysiological Fingerprints in a Concentration-Dependent and Region-Specific Manner

Brain function depends on complex circuit interactions between excitatory and inhibitory neurons embedded in local and long-range networks. Systemic GABAA-receptor (GABAAR) or NMDA-receptor (NMDAR) modulation alters the excitatory-inhibitory balance (EIB), measurable with electroencephalography (EEG). However, EEG signatures are complex in localization and spectral composition. We developed and applied analytical tools to investigate the effects of two EIB modulators, MK801 (NMDAR antagonist) and diazepam (GABAAR modulator), on periodic and aperiodic EEG features in freely-moving male Sprague Dawley rats. We investigated how, across three brain regions, EEG features are correlated with EIB modulation. We found that the periodic component was composed of seven frequency bands that presented region-dependent and compound-dependent changes. The aperiodic component was also different between compounds and brain regions. Importantly, the parametrization into periodic and aperiodic components unveiled correlations between quantitative EEG and plasma concentrations of pharmacological compounds. MK-801 exposures were positively correlated with the slope of the aperiodic component. Concerning the periodic component, MK-801 exposures correlated negatively with the peak frequency of low-γ oscillations but positively with those of high-γ and high-frequency oscillations (HFOs). As for the power, θ and low-γ oscillations correlated negatively with MK-801, whereas mid-γ correlated positively. Diazepam correlated negatively with the knee of the aperiodic component, positively to β and negatively to low-γ oscillatory power, and positively to the modal frequency of θ, low-γ, mid-γ, and high-γ. In conclusion, correlations between exposures and pharmacodynamic effects can be better-understood thanks to the parametrization of EEG into periodic and aperiodic components. Such parametrization could be key in functional biomarker discovery.

Adaptive Stimulations in a Biophysical Network Model of Parkinson’s Disease

Deep brain stimulation (DBS)—through a surgically implanted electrode to the subthalamic nucleus (STN)—has become a widely used therapeutic option for the treatment of Parkinson’s disease and other neurological disorders. The standard conventional high-frequency stimulation (HF) that is currently used has several drawbacks. To overcome the limitations of HF, researchers have been developing closed-loop and demand-controlled, adaptive stimulation protocols wherein the amount of current that is delivered is turned on and off in real-time in accordance with a biophysical signal. Computational modeling of DBS in neural network models is an increasingly important tool in the development of new protocols that aid researchers in animal and clinical studies. In this computational study, we seek to implement a novel technique of DBS where we stimulate the STN in an adaptive fashion using the interspike time of the neurons to control stimulation. Our results show that our protocol eliminates bursts in the synchronized bursting neuronal activity of the STN, which is hypothesized to cause the failure of thalamocortical neurons (TC) to respond properly to excitatory cortical inputs. Further, we are able to significantly decrease the TC relay errors, representing potential therapeutics for Parkinson’s disease.

Broadband aperiodic components of local field potentials reflect inherent differences between cortical and subcortical activity

Information flow in brain networks is reflected in intracerebral local field potential (LFP) measurements that have both periodic and aperiodic components. The 1/fχ broadband aperiodic component of the power spectra has been shown to track arousal level and to correlate with other physiological and pathophysiological states, with consistent patterns across cortical regions. Previous studies have focused almost exclusively on cortical neurophysiology. Here we explored the aperiodic activity of subcortical nuclei from the human thalamus and basal ganglia, in relation to simultaneously recorded cortical activity. We elaborated on the FOOOF (fitting of one over f) method by creating a new parameterization of the aperiodic component with independent and more easily interpretable parameters, which allows seamlessly fitting spectra with and without an aperiodic knee, a component of the signal that reflects the dominant timescale of aperiodic fluctuations. First, we found that the aperiodic exponent from sensorimotor cortex in Parkinson's disease (PD) patients correlated with disease severity. Second, although the aperiodic knee frequency changed across cortical regions as previously reported, no aperiodic knee was detected from subcortical regions across movement disorders patients, including the ventral thalamus (VIM), globus pallidus internus (GPi) and subthalamic nucleus (STN). All subcortical region studied exhibited a relatively low aperiodic exponent (χSTN=1.3±0.2, χVIM=1.4±0.1, χGPi =1.4±0.1) that differed markedly from cortical values (χCortex=3.2±0.4, fkCortex=17±5 Hz). These differences were replicated in a second dataset from epilepsy patients undergoing intracranial monitoring that included thalamic recordings. The consistently lower aperiodic exponent and lack of an aperiodic knee from all subcortical recordings may reflect cytoarchitectonic and/or functional differences between subcortical nuclei and the cortex.

Sampling effects and measurement overlap can bias the inference of neuronal avalanches

To date, it is still impossible to sample the entire mammalian brain with single-neuron precision. This forces one to either use spikes (focusing on few neurons) or to use coarse-sampled activity (averaging over many neurons, e.g. LFP). Naturally, the sampling technique impacts inference about collective properties. Here, we emulate both sampling techniques on a simple spiking model to quantify how they alter observed correlations and signatures of criticality. We describe a general effect: when the inter-electrode distance is small, electrodes sample overlapping regions in space, which increases the correlation between the signals. For coarse-sampled activity, this can produce power-law distributions even for non-critical systems. In contrast, spike recordings do not suffer this particular bias and underlying dynamics can be identified. This may resolve why coarse measures and spikes have produced contradicting results in the past.

Local networks from different parts of the human cerebral cortex generate and share the same population dynamic

A major goal of neuroscience is to reveal mechanisms supporting collaborative actions of neurons in local and larger-scale networks. However, no clear overall principle of operation has emerged despite decades-long experimental efforts. Here, we used an unbiased method to extract and identify the dynamics of local postsynaptic network states contained in the cortical field potential. Field potentials were recorded by depth electrodes targeting a wide selection of cortical regions during spontaneous activities, and sensory, motor, and cognitive experimental tasks. Despite different architectures and different activities, all local cortical networks generated the same type of dynamic confined to one region only of state space. Surprisingly, within this region, state trajectories expanded and contracted continuously during all brain activities and generated a single expansion followed by a contraction in a single trial. This behavior deviates from known attractors and attractor networks. The state-space contractions of particular subsets of brain regions cross-correlated during perceptive, motor, and cognitive tasks. Our results imply that the cortex does not need to change its dynamic to shift between different activities, making task-switching inherent in the dynamic of collective cortical operations. Our results provide a mathematically described general explanation of local and larger scale cortical dynamic.

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

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

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

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

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

Using noise for the better: The effects of transcranial random noise stimulation on the brain and behavior

Van der Groen, O., Potok, W., Wenderoth, N., Edwards, G., Mattingley, J.B. and Edwards, D. Using noise for the better: The effects of transcranial random noise stimulation on the brain and behavior. NEUROSCI BIOBEHAV REV X (X) XXX-XXX 2021.- Transcranial random noise stimulation (tRNS) is a non-invasive electrical brain stimulation method that is increasingly employed in studies of human brain function and behavior, in health and disease. tRNS is effective in modulating perception acutely and can improve learning. By contrast, its effectiveness for modulating higher cognitive processes is variable. Prolonged stimulation with tRNS, either as one longer application, or multiple shorter applications, may engage plasticity mechanisms that can result in long-term benefits. Here we provide an overview of the current understanding of the effects of tRNS on the brain and behavior and provide some specific recommendations for future research.

Dopaminergic Modulation of Local Non-oscillatory Activity and Global-Network Properties in Parkinson's Disease: An EEG Study

Dopaminergic medication for Parkinson's disease (PD) modulates neuronal oscillations and functional connectivity (FC) across the basal ganglia-thalamic-cortical circuit. However, the non-oscillatory component of the neuronal activity, potentially indicating a state of excitation/inhibition balance, has not yet been investigated and previous studies have shown inconsistent changes of cortico-cortical connectivity as a response to dopaminergic medication. To further elucidate changes of regional non-oscillatory component of the neuronal power spectra, FC, and to determine which aspects of network organization obtained with graph theory respond to dopaminergic medication, we analyzed a resting-state electroencephalography (EEG) dataset including 15 PD patients during OFF and ON medication conditions. We found that the spectral slope, typically used to quantify the broadband non-oscillatory component of power spectra, steepened particularly in the left central region in the ON compared to OFF condition. In addition, using lagged coherence as a FC measure, we found that the FC in the beta frequency range between centro-parietal and frontal regions was enhanced in the ON compared to the OFF condition. After applying graph theory analysis, we observed that at the lower level of topology the node degree was increased, particularly in the centro-parietal area. Yet, results showed no significant difference in global topological organization between the two conditions: either in global efficiency or clustering coefficient for measuring global and local integration, respectively. Interestingly, we found a close association between local/global spectral slope and functional network global efficiency in the OFF condition, suggesting a crucial role of local non-oscillatory dynamics in forming the functional global integration which characterizes PD. These results provide further evidence and a more complete picture for the engagement of multiple cortical regions at various levels in response to dopaminergic medication in PD.

Alpha blocking and 1/fβ spectral scaling in resting EEG can be accounted for by a sum of damped alpha band oscillatory processes

The dynamical and physiological basis of alpha band activity and 1/f^β noise in the EEG are the subject of continued speculation. Here we conjecture, on the basis of empirical data analysis, that both of these features may be economically accounted for through a single process if the resting EEG is conceived of being the sum of multiple stochastically perturbed alpha band damped linear oscillators with a distribution of dampings (relaxation rates). The modulation of alpha-band and 1/f^β noise activity by changes in damping is explored in eyes closed (EC) and eyes open (EO) resting state EEG. We aim to estimate the distribution of dampings by solving an inverse problem applied to EEG power spectra. The characteristics of the damping distribution are examined across subjects, sensors and recording condition (EC/EO). We find that there are robust changes in the damping distribution between EC and EO recording conditions across participants. The estimated damping distributions are found to be predominantly bimodal, with the number and position of the modes related to the sharpness of the alpha resonance and the scaling (β) of the power spectrum (1/f^β). The results suggest that there exists an intimate relationship between resting state alpha activity and 1/f^β noise with changes in both governed by changes to the damping of the underlying alpha oscillatory processes. In particular, alpha-blocking is observed to be the result of the most weakly damped distribution mode becoming more heavily damped. The results suggest a novel way of characterizing resting EEG power spectra and provides new insight into the central role that damped alpha-band activity may play in characterising the spatio-temporal features of resting state EEG.

The resting human electroencephalogram (EEG) exhibits two dominant spectral features: the alpha rhythm (8–13 Hz) and its associated attenuation between eyes-closed and eyes-open resting state (alpha blocking), and the 1/f^β scaling of the power spectrum. While these phenomena are well studied a thorough understanding of their respective generative processes remains elusive. By employing a theoretical approach that follows from neural population models of EEG we demonstrate that it is possible to economically account for both of these phenomena using a singular mechanistic framework: resting EEG is assumed to arise from the summed activity of multiple uncorrelated, stochastically driven, damped alpha band linear oscillatory processes having a distribution of relaxation rates or dampings. By numerically estimating these damping distributions from eyes-closed and eyes-open EEG data, in a total of 136 participants, it is found that such damping distributions are predominantly bimodal in shape. The most weakly damped mode is found to account for alpha band power, with alpha blocking being driven by an increase in the damping of this weakly damped mode, whereas the second, and more heavily damped mode, is able to explain 1/f^β scaling present in the resting state EEG spectra.

Discovering Microcircuit Secrets With Multi-Spot Imaging and Electrophysiological Recordings: The Example of Cerebellar Network Dynamics

The cerebellar cortex microcircuit is characterized by a highly ordered neuronal architecture having a relatively simple and stereotyped connectivity pattern. For a long time, this structural simplicity has incorrectly led to the idea that anatomical considerations would be sufficient to understand the dynamics of the underlying circuitry. However, recent experimental evidence indicates that cerebellar operations are much more complex than solely predicted by anatomy, due to the crucial role played by neuronal and synaptic properties. To be able to explore neuronal and microcircuit dynamics, advanced imaging, electrophysiological techniques and computational models have been combined, allowing us to investigate neuronal ensembles activity and to connect microscale to mesoscale phenomena. Here, we review what is known about cerebellar network organization, neural dynamics and synaptic plasticity and point out what is still missing and would require experimental assessments. We consider the available experimental techniques that allow a comprehensive assessment of circuit dynamics, including voltage and calcium imaging and extracellular electrophysiological recordings with multi-electrode arrays (MEAs). These techniques are proving essential to investigate the spatiotemporal pattern of activity and plasticity in the cerebellar network, providing new clues on how circuit dynamics contribute to motor control and higher cognitive functions.

A Simulation on Relation between Power Distribution of Low-Frequency Field Potentials and Conducting Direction of Rhythm Generator Flowing through 3D Asymmetrical Brain Tissue

Although the power of low-frequency oscillatory field potentials (FP) has been extensively applied previously, few studies have investigated the influence of conducting direction of deep-brain rhythm generator on the power distribution of low-frequency oscillatory FPs on the head surface. To address this issue, a simulation was designed based on the principle of electroencephalogram (EEG) generation of equivalent dipole current in deep brain, where a single oscillatory dipole current represented the rhythm generator, the dipole moment for the rhythm generator’s conducting direction (which was orthogonal and rotating every 30 degrees and at pointing to or parallel to the frontal lobe surface) and the (an)isotropic conduction medium for the 3D (a)symmetrical brain tissue. Both the power above average (significant power value, SP value) and its space (SP area) of low-frequency oscillatory FPs were employed to respectively evaluate the strength and the space of the influence. The computation was conducted using the finite element method (FEM) and Hilbert transform. The finding was that either the SP value or the SP area could be reduced or extended, depending on the conducting direction of deep-brain rhythm generator flowing in the (an)isotropic medium, suggesting that the 3D (a)symmetrical brain tissue could decay or strengthen the spatial spread of a rhythm generator conducting in a different direction.

Flexible Frequency Switching in Adult Mouse Visual Cortex Is Mediated by Competition Between Parvalbumin and Somatostatin Expressing Interneurons

Neuronal networks in rodent primary visual cortex (V1) can generate oscillations in different frequency bands depending on the network state and the level of visual stimulation. High-frequency gamma rhythms, for example, dominate the network's spontaneous activity in adult mice but are attenuated upon visual stimulation, during which the network switches to the beta band instead. The spontaneous local field potential (LFP) of juvenile mouse V1, however, mainly contains beta rhythms and presenting a stimulus does not elicit drastic changes in network oscillations. We study, in a spiking neuron network model, the mechanism in adult mice allowing for flexible switches between multiple frequency bands and contrast this to the network structure in juvenile mice that lack this flexibility. The model comprises excitatory pyramidal cells (PCs) and two types of interneurons: the parvalbumin-expressing (PV) and the somatostatinexpressing (SOM) interneuron. In accordance with experimental findings, the pyramidal-PV and pyramidal-SOM cell subnetworks are associated with gamma and beta oscillations, respectively. In our model, they are both generated via a pyramidal-interneuron gamma (PING) mechanism, wherein the PCs drive the oscillations. Furthermore, we demonstrate that large but not small visual stimulation activates SOM cells, which shift the frequency of resting-state gamma oscillations produced by the pyramidal-PV cell subnetwork so that beta rhythms emerge. Finally, we show that this behavior is obtained for only a subset of PV and SOM interneuron projection strengths, indicating that their influence on the PCs should be balanced so that they can compete for oscillatory control of the PCs. In sum, we propose a mechanism by which visual beta rhythms can emerge from spontaneous gamma oscillations in a network model of the mouse V1; for this mechanism to reproduce V1 dynamics in adult mice, balance between the effective strengths of PV and SOM cells is required.

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.

Significance and Translational Value of High-Frequency Cortico-Basal Ganglia Oscillations in Parkinson's Disease

The mechanisms and significance of basal ganglia oscillations is a fundamental research question engaging both clinical and basic investigators. In Parkinson’s disease (PD), neural activity in basal ganglia nuclei is characterized by oscillatory patterns that are believed to disrupt the dynamic processing of movement-related information and thus generate motor symptoms. Beta-band oscillations associated with hypokinetic states have been reviewed in several excellent previous articles. Here we focus on faster oscillatory phenomena that have been reported in association with a diverse range of motor states. We review the occurrence of different types of fast oscillations and the evidence supporting their pathophysiological role. We also provide a general discussion on the definition, possible mechanisms, and translational value of synchronized oscillations of different frequencies in cortico-basal ganglia structures. Revealing how oscillatory phenomena are caused and spread in cortico-basal ganglia-thalamocortical networks will offer a key to unlock the neural codes of both motor and non-motor symptoms in PD. In preclinical therapeutic research, recording of oscillatory neural activities holds the promise to unravel mechanisms of action of current and future treatments.

Standard multiscale entropy reflects neural dynamics at mismatched temporal scales: What's signal irregularity got to do with it?

Multiscale Entropy (MSE) is used to characterize the temporal irregularity of neural time series patterns. Due to its' presumed sensitivity to non-linear signal characteristics, MSE is typically considered a complementary measure of brain dynamics to signal variance and spectral power. However, the divergence between these measures is often unclear in application. Furthermore, it is commonly assumed (yet sparingly verified) that entropy estimated at specific time scales reflects signal irregularity at those precise time scales of brain function. We argue that such assumptions are not tenable. Using simulated and empirical electroencephalogram (EEG) data from 47 younger and 52 older adults, we indicate strong and previously underappreciated associations between MSE and spectral power, and highlight how these links preclude traditional interpretations of MSE time scales. Specifically, we show that the typical definition of temporal patterns via "similarity bounds" biases coarse MSE scales-that are thought to reflect slow dynamics-by high-frequency dynamics. Moreover, we demonstrate that entropy at fine time scales-presumed to indicate fast dynamics-is highly sensitive to broadband spectral power, a measure dominated by low-frequency contributions. Jointly, these issues produce counterintuitive reflections of frequency-specific content on MSE time scales. We emphasize the resulting inferential problems in a conceptual replication of cross-sectional age differences at rest, in which scale-specific entropy age effects could be explained by spectral power differences at mismatched temporal scales. Furthermore, we demonstrate how such problems may be alleviated, resulting in the indication of scale-specific age differences in rhythmic irregularity. By controlling for narrowband contributions, we indicate that spontaneous alpha rhythms during eyes open rest transiently reduce broadband signal irregularity. Finally, we recommend best practices that may better permit a valid estimation and interpretation of neural signal irregularity at time scales of interest.

An Indexing Theory for Working Memory Based on Fast Hebbian Plasticity

Working memory (WM) is a key component of human memory and cognition. Computational models have been used to study the underlying neural mechanisms, but neglected the important role of short-term memory (STM) and long-term memory (LTM) interactions for WM. Here, we investigate these using a novel multiarea spiking neural network model of prefrontal cortex (PFC) and two parietotemporal cortical areas based on macaque data. We propose a WM indexing theory that explains how PFC could associate, maintain, and update multimodal LTM representations. Our simulations demonstrate how simultaneous, brief multimodal memory cues could build a temporary joint memory representation as an “index” in PFC by means of fast Hebbian synaptic plasticity. This index can then reactivate spontaneously and thereby also the associated LTM representations. Cueing one LTM item rapidly pattern completes the associated uncued item via PFC. The PFC–STM network updates flexibly as new stimuli arrive, thereby gradually overwriting older representations.

Basal ganglia oscillations as biomarkers for targeting circuit dysfunction in Parkinson's disease

Oscillations are a naturally occurring phenomenon in highly interconnected dynamical systems. However, it is thought that excessive synchronized oscillations in brain circuits can be detrimental for many brain functions by disrupting neuronal information processing. Because synchronized basal ganglia oscillations are a hallmark of Parkinson's disease (PD), it has been suggested that aberrant rhythmic activity associated with symptoms of the disease could be used as a physiological biomarker to guide pharmacological and electrical neuromodulatory interventions. We here briefly review the various manifestations of basal ganglia oscillations observed in human subjects and in animal models of PD. In this context, we also review the evidence supporting a pathophysiological role of different oscillations for the suppression of voluntary movements as well as for the induction of excessive motor activity. In light of these findings, it is discussed how oscillations could be used to guide a more precise targeting of dysfunctional circuits to obtain improved symptomatic treatment of PD.

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

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

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

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

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