Computing the Local Field Potential (LFP) from Integrate-and-Fire Network Models

Revisiting the role of computational neuroimaging in the era of integrative neuroscience

Computational models have become integral to human neuroimaging research, providing both mechanistic insights and predictive tools for human cognition and behavior. However, concerns persist regarding the ecological validity of lab-based neuroimaging studies and whether their spatiotemporal resolution is not sufficient for capturing neural dynamics. This review aims to re-examine the utility of computational neuroimaging, particularly in light of the growing prominence of alternative neuroscientific methods and the growing emphasis on more naturalistic behaviors and paradigms. Specifically, we will explore how computational modeling can both enhance the analysis of high-dimensional imaging datasets and, conversely, how neuroimaging, in conjunction with other data modalities, can inform computational models through the lens of neurobiological plausibility. Collectively, this evidence suggests that neuroimaging remains critical for human neuroscience research, and when enhanced by computational models, imaging can serve an important role in bridging levels of analysis and understanding. We conclude by proposing key directions for future research, emphasizing the development of standardized paradigms and the integrative use of computational modeling across neuroimaging techniques.

Horizontal cortical connections shape intrinsic traveling waves into feature-selective motifs that regulate perceptual sensitivity

Intrinsic cortical activity forms traveling waves that modulate sensory-evoked responses and perceptual sensitivity. These intrinsic traveling waves (iTWs) may arise from the coordination of synaptic activity through long-range feature-dependent horizontal connectivity within cortical areas. In a spiking network model that incorporates feature-selective patchy connections, we observe iTW motifs that result from shifts in excitatory/inhibitory balance as action potentials traverse these patchy connections. To test whether feature-selective motifs occur in vivo, we examined data recorded in the middle temporal visual area (Area MT) of marmosets performing a visual detection task. We find that some iTWs form motifs that are feature selective, exhibiting direction-selective modulations in spiking activity. Further, motifs modulate the gain of target-evoked responses and perceptual sensitivity if the target matches the preference of the motif. These results suggest that iTWs are shaped by the patchy horizontal fiber projections in the cortex and can regulate neural and perceptual sensitivity in a feature-selective manner.

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.

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling

The major anxiety disorders (ANX; including generalized anxiety disorder, panic disorder, and phobias) are highly prevalent, often onset early, persist throughout life, and cause substantial global disability. Although distinct in their clinical presentations, they likely represent differential expressions of a dysregulated threat-response system. Here we present a genome-wide association meta-analysis comprising 122,341 European ancestry ANX cases and 729,881 controls. We identified 58 independent genome-wide significant ANX risk variants and 66 genes with robust biological support. In an independent sample of 1,175,012 self-report ANX cases and 1,956,379 controls, 51 of the 58 associated variants were replicated. As predicted by twin studies, we found substantial genetic correlation between ANX and depression, neuroticism, and other internalizing phenotypes. Follow-up analyses demonstrated enrichment in all major brain regions and highlighted GABAergic signaling as one potential mechanism underlying ANX genetic risk. These results advance our understanding of the genetic architecture of ANX and prioritize genes for functional follow-up studies.

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.

Map plasticity following noise exposure in auditory cortex of rats: implications for disentangling neural correlates of tinnitus and hyperacusis

Introduction

Both tinnitus and hyperacusis, likely triggered by hearing loss, can be attributed to maladaptive plasticity in auditory perception. However, owing to their co-occurrence, disentangling their neural mechanisms proves difficult. We hypothesized that the neural correlates of tinnitus are associated with neural activities triggered by low-intensity tones, while hyperacusis is linked to responses to moderate- and high-intensity tones.

Methods

To test these hypotheses, we conducted behavioral and electrophysiological experiments in rats 2 to 8 days after traumatic tone exposure.

Results

In the behavioral experiments, prepulse and gap inhibition tended to exhibit different frequency characteristics (although not reaching sufficient statistical levels), suggesting that exposure to traumatic tones led to acute symptoms of hyperacusis and tinnitus at different frequency ranges. When examining the auditory cortex at the thalamocortical recipient layer, we observed that tinnitus symptoms correlated with a disorganized tonotopic map, typically characterized by responses to low-intensity tones. Neural correlates of hyperacusis were found in the cortical recruitment function at the multi-unit activity (MUA) level, but not at the local field potential (LFP) level, in response to moderate- and high-intensity tones. This shift from LFP to MUA was associated with a loss of monotonicity, suggesting a crucial role for inhibitory synapses.

Discussion

Thus, in acute symptoms of traumatic tone exposure, our experiments successfully disentangled the neural correlates of tinnitus and hyperacusis at the thalamocortical recipient layer of the auditory cortex. They also suggested that tinnitus is linked to central noise, whereas hyperacusis is associated with aberrant gain control. Further interactions between animal experiments and clinical studies will offer insights into neural mechanisms, diagnosis and treatments of tinnitus and hyperacusis, specifically in terms of long-term plasticity of chronic symptoms.

A Translaminar Spacetime Code Supports Touch-Evoked Traveling Waves

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked cortical traveling waves and their underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a variable late wave that lasted hundreds of milliseconds post-stimulus. Strikingly, late-wave dynamics were modulated by stimulus value and correlated with task performance. Mechanistically, the late wave component was i) modulated by motor feedback, ii) complemented by a sparse ensemble pattern across layer 2/3, which a balanced-state network model reconciled via inhibitory stabilization, and iii) aligned to regenerative Layer-5 apical dendritic Ca 2+ events. Our results reveal a translaminar spacetime pattern organized by cortical feedback in the sensory cortex that supports touch-evoked traveling waves.

GRAPHICAL ABSTRACT AND HIGHLIGHTS

Whisker touch evokes both early- and late-traveling waves in the barrel cortex over 100's of millisecondsReward reinforcement modulates wave dynamics Late wave emergence coincides with network sparsity in L23 and time-locked L5 dendritic Ca 2+ spikes Experimental and computational results link motor feedback to distinct translaminar spacetime patterns.

Distributed and dynamical communication: a mechanism for flexible cortico-cortical interactions and its functional roles in visual attention

Perceptual and cognitive processing relies on flexible communication among cortical areas; however, the underlying neural mechanism remains unclear. Here we report a mechanism based on the realistic spatiotemporal dynamics of propagating wave patterns in neural population activity. Using a biophysically plausible, multiarea spiking neural circuit model, we demonstrate that these wave patterns, characterized by their rich and complex dynamics, can account for a wide variety of empirically observed neural processes. The coordinated interactions of these wave patterns give rise to distributed and dynamic communication (DDC) that enables flexible and rapid routing of neural activity across cortical areas. We elucidate how DDC unifies the previously proposed oscillation synchronization-based and subspace-based views of interareal communication, offering experimentally testable predictions that we validate through the analysis of Allen Institute Neuropixels data. Furthermore, we demonstrate that DDC can be effectively modulated during attention tasks through the interplay of neuromodulators and cortical feedback loops. This modulation process explains many neural effects of attention, underscoring the fundamental functional role of DDC in cognition.

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.

Neurobiological Causal Models of Language Processing

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

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Rodents' visual gamma as a biomarker of pathological neural conditions

Neural gamma oscillations (indicatively 30-100 Hz) are ubiquitous: they are associated with a broad range of functions in multiple cortical areas and across many animal species. Experimental and computational works established gamma rhythms as a global emergent property of neuronal networks generated by the balanced and coordinated interaction of excitation and inhibition. Coherently, gamma activity is strongly influenced by the alterations of synaptic dynamics which are often associated with pathological neural dysfunctions. We argue therefore that these oscillations are an optimal biomarker for probing the mechanism of cortical dysfunctions. Gamma oscillations are also highly sensitive to external stimuli in sensory cortices, especially the primary visual cortex (V1), where the stimulus dependence of gamma oscillations has been thoroughly investigated. Gamma manipulation by visual stimuli tuning is particularly easy in rodents, which have become a standard animal model for investigating the effects of network alterations on gamma oscillations. Overall, gamma in the rodents' visual cortex offers an accessible probe on dysfunctional information processing in pathological conditions. Beyond vision-related dysfunctions, alterations of gamma oscillations in rodents were indeed also reported in neural deficits such as migraine, epilepsy and neurodegenerative or neuropsychiatric conditions such as Alzheimer's, schizophrenia and autism spectrum disorders. Altogether, the connections between visual cortical gamma activity and physio-pathological conditions in rodent models underscore the potential of gamma oscillations as markers of neuronal (dys)functioning.

Sleep and wake in a model of the thalamocortical system with Martinotti cells

The mechanisms leading to the alternation between active (UP) and silent (DOWN) states during sleep slow waves (SWs) remain poorly understood. Previous models have explained the transition to the DOWN state by a progressive failure of excitation because of the build-up of adaptation currents or synaptic depression. However, these models are at odds with recent studies suggesting a role for presynaptic inhibition by Martinotti cells (MaCs) in generating SWs. Here, we update a classical large-scale model of sleep SWs to include MaCs and propose a different mechanism for the generation of SWs. In the wake mode, the network exhibits irregular and selective activity with low firing rates (FRs). Following an increase in the strength of background inputs and a modulation of synaptic strength and potassium leak potential mimicking the reduced effect of acetylcholine during sleep, the network enters a sleep-like regime in which local increases of network activity trigger bursts of MaC activity, resulting in strong disfacilitation of the local network via presynaptic GABAB1a -type inhibition. This model replicates findings on slow wave activity (SWA) during sleep that challenge previous models, including low and skewed FRs that are comparable between the wake and sleep modes, higher synchrony of transitions to DOWN states than to UP states, the possibility of triggering SWs by optogenetic stimulation of MaCs, and the local dependence of SWA on synaptic strength. Overall, this work points to a role for presynaptic inhibition by MaCs in the generation of DOWN states during sleep.

A developmental increase of inhibition promotes the emergence of hippocampal ripples

Sharp wave-ripples (SPW-Rs) are a hippocampal network phenomenon critical for memory consolidation and planning. SPW-Rs have been extensively studied in the adult brain, yet their developmental trajectory is poorly understood. While SPWs have been recorded in rodents shortly after birth, the time point and mechanisms of ripple emergence are still unclear. Here, we combine in vivo electrophysiology with optogenetics and chemogenetics in 4 to 12-day-old mice to address this knowledge gap. We show that ripples are robustly detected and induced by light stimulation of channelrhodopsin-2-transfected CA1 pyramidal neurons only from postnatal day 10 onwards. Leveraging a spiking neural network model, we mechanistically link the maturation of inhibition and ripple emergence. We corroborate these findings by reducing ripple rate upon chemogenetic silencing of CA1 interneurons. Finally, we show that early SPW-Rs elicit a more robust prefrontal cortex response than SPWs lacking ripples. Thus, development of inhibition promotes ripples emergence.

Information representation in an oscillating neural field model modulated by working memory signals

We study how stimulus information can be represented in the dynamical signatures of an oscillatory model of neural activity-a model whose activity can be modulated by input akin to signals involved in working memory (WM). We developed a neural field model, tuned near an oscillatory instability, in which the WM-like input can modulate the frequency and amplitude of the oscillation. Our neural field model has a spatial-like domain in which an input that preferentially targets a point-a stimulus feature-on the domain will induce feature-specific activity changes. These feature-specific activity changes affect both the mean rate of spikes and the relative timing of spiking activity to the global field oscillation-the phase of the spiking activity. From these two dynamical signatures, we define both a spike rate code and an oscillatory phase code. We assess the performance of these two codes to discriminate stimulus features using an information-theoretic analysis. We show that global WM input modulations can enhance phase code discrimination while simultaneously reducing rate code discrimination. Moreover, we find that the phase code performance is roughly two orders of magnitude larger than that of the rate code defined for the same model solutions. The results of our model have applications to sensory areas of the brain, to which prefrontal areas send inputs reflecting the content of WM. These WM inputs to sensory areas have been established to induce oscillatory changes similar to our model. Our model results suggest a mechanism by which WM signals may enhance sensory information represented in oscillatory activity beyond the comparatively weak representations based on the mean rate activity.

Horizontal cortical connections shape intrinsic traveling waves into feature-selective motifs that regulate perceptual sensitivity

Intrinsic, ongoing fluctuations of cortical activity form traveling waves that modulate the gain of sensory-evoked responses and perceptual sensitivity. Several lines of evidence suggest that intrinsic traveling waves (iTWs) may arise, in part, from the coordination of synaptic activity through the recurrent horizontal connectivity within cortical areas, which include long range patchy connections that link neurons with shared feature preferences. In a spiking network model with anatomical topology that incorporates feature-selective patchy connections, which we call the Balanced Patchy Network (BPN), we observe repeated iTWs, which we refer to as motifs. In the model, motifs stem from fluctuations in the excitability of like-tuned neurons that result from shifts in E/I balance as action potentials traverse these patchy connections. To test if feature-selective motifs occur in vivo, we examined data previously recorded using multielectrode arrays in Area MT of marmosets trained to perform a threshold visual detection task. Using a newly developed method for comparing the similarity of wave patterns we found that some iTWs can be grouped into motifs. As predicted by the BPN, many of these motifs are feature selective, exhibiting direction-selective modulations in ongoing spiking activity. Further, motifs modulate the gain of the response evoked by a target and perceptual sensitivity to the target if the target matches the preference of the motif. These results provide evidence that iTWs are shaped by the patterns of horizontal fiber projections in the cortex and that patchy connections enable iTWs to regulate neural and perceptual sensitivity in a feature selective manner.

Criticality supports cross-frequency cortical-thalamic information transfer during conscious states

Consciousness is thought to be regulated by bidirectional information transfer between the cortex and thalamus, but the nature of this bidirectional communication - and its possible disruption in unconsciousness - remains poorly understood. Here, we present two main findings elucidating mechanisms of corticothalamic information transfer during conscious states. First, we identify a highly preserved spectral channel of cortical-thalamic communication that is present during conscious states, but which is diminished during the loss of consciousness and enhanced during psychedelic states. Specifically, we show that in humans, mice, and rats, information sent from either the cortex or thalamus via δ/θ/α waves (∼1-13 Hz) is consistently encoded by the other brain region by high γ waves (52-104 Hz); moreover, unconsciousness induced by propofol anesthesia or generalized spike-and-wave seizures diminishes this cross-frequency communication, whereas the psychedelic 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) enhances this low-to-high frequency interregional communication. Second, we leverage numerical simulations and neural electrophysiology recordings from the thalamus and cortex of human patients, rats, and mice to show that these changes in cross-frequency cortical-thalamic information transfer may be mediated by excursions of low-frequency thalamocortical electrodynamics toward/away from edge-of-chaos criticality, or the phase transition from stability to chaos. Overall, our findings link thalamic-cortical communication to consciousness, and further offer a novel, mathematically well-defined framework to explain the disruption to thalamic-cortical information transfer during unconscious states.

Dynamic evolution of the anterior cingulate-insula network during seizures

OBJECTIVES

In physiological situations, the anterior cingulate cortex (ACC) and anterior insular cortex (AIC) are prone to coactivation. The functional connectivity and interaction between ACC and AIC in the context of epilepsy remain unclear. This study aimed to investigate the dynamic coupling between these two brain regions during seizures.

METHODS

Patients who underwent stereoelectroencephalography (SEEG) recording were included in this study. The SEEG data were visually inspected and quantitatively analyzed. The narrowband oscillations and aperiodic components at seizure onset were parameterized. The frequency-specific non-linear correlation analysis was applied to the functional connectivity. The excitation/inhibition ratio (E:I ratio) reflected by the aperiodic slope was performed to evaluate the excitability.

RESULTS

Twenty patients were included in the study, with 10 diagnosed with anterior cingulate epilepsy and 10 with anterior insular epilepsy. In both types of epilepsy, the correlation coefficient (h2 ) between the ACC and AIC at seizure onset exhibited a significantly higher value than that during interictal and preictal periods (p < 0.05). The direction index (D) showed a significant increase at seizure onset, serving as an indicator for the direction of information flow between these two brain regions with up to 90% accuracy. The E:I ratio increased significantly at seizure onset, with the seizure-onset zone (SOZ) demonstrating a more pronounced increase compared to non-SOZ (p < 0.05). For seizures originating from AIC, the E:I ratio was significantly higher in the AIC than in the ACC (p = 0.0364).

CONCLUSIONS

In the context of epilepsy, the ACC and AIC are dynamically coupled during seizures. The functional connectivity and excitability exhibit a significant increase at seizure onset. By analyzing connectivity and excitability, the SOZ in ACC and AIC can be identified. The direction index (D) serves as an indicator for the direction of information flow from SOZ to non-SOZ. Notably, the excitability of SOZ changes more significantly than that of non-SOZ.

γ And β Band Oscillation in Working Memory Given Sequential or Concurrent Multiple Items: A Spiking Network Model

Working memory (WM) can maintain sequential and concurrent information, and the load enhances the γ band oscillation during the delay period. To provide a unified account for these phenomena in working memory, we investigated a continuous network model consisting of pyramidal cells, high-threshold fast-spiking interneurons (FS), and low-threshold nonfast-spiking interneurons (nFS) for working memory of sequential and concurrent directional cues. Our model exhibits the γ (30-100 Hz) and β (10-30 Hz) band oscillation during the retention of both concurrent cues and sequential cues. We found that the β oscillation results from the interaction between pyramidal cells and nFS, whereas the γ oscillation emerges from the interaction between pyramidal cells and FS because of the strong excitation elicited by cue presentation, shedding light on the mechanism underlying the enhancement of γ power in many cognitive executions.

A supervised data-driven spatial filter denoising method for speech artifacts in intracranial electrophysiological recordings

Neurosurgical procedures that enable direct brain recordings in awake patients offer unique opportunities to explore the neurophysiology of human speech. The scarcity of these opportunities and the altruism of participating patients compel us to apply the highest rigor to signal analysis. Intracranial electroencephalography (iEEG) signals recorded during overt speech can contain a speech artifact that tracks the fundamental frequency (F0) of the participant's voice, involving the same high-gamma frequencies that are modulated during speech production and perception. To address this artifact, we developed a spatial-filtering approach to identify and remove acoustic-induced contaminations of the recorded signal. We found that traditional reference schemes jeopardized signal quality, whereas our data-driven method denoised the recordings while preserving underlying neural activity.

Origins and Evolution of Novel Bacteroides in Captive Apes

Bacterial strains evolve in response to the gut environment of their hosts, with genomic changes that influence their interactions with hosts as well as with other members of the gut community. Great apes in captivity have acquired strains of Bacteroides xylanisolvens, which are common within gut microbiome of humans but not typically found other apes, thereby enabling characterization of strain evolution following colonization. Here, we isolate, sequence and reconstruct the history of gene gain and loss events in numerous captive-ape-associated strains since their divergence from their closest human-associated strains. We show that multiple captive-ape-associated B. xylanisolvens lineages have independently acquired gene complexes that encode functions related to host mucin metabolism. Our results support the finding of high genome fluidity in Bacteroides, in that several strains, in moving from humans to captive apes, have rapidly gained large genomic regions that augment metabolic properties not previously present in their relatives.

Temporal characteristics of gamma rhythm constrain properties of noise in an inhibition-stabilized network model

Gamma rhythm refers to oscillatory neural activity between 30 and 80 Hz, induced in visual cortex by stimuli such as iso-luminant hues or gratings. The power and peak frequency of gamma depend on the properties of the stimulus such as size and contrast. Gamma waveform is typically arch-shaped, with narrow troughs and broad peaks, and can be replicated in a self-oscillating Wilson-Cowan (WC) model operating in an appropriate regime. However, oscillations in this model are infinitely long, unlike physiological gamma that occurs in short bursts. Further, unlike the model, gamma is faster after stimulus onset and slows down over time. Here, we first characterized gamma burst duration in local field potential data recorded from two monkeys as they viewed full screen iso-luminant hues. We then added different types of noise in the inputs to the WC model and tested how that affected duration and temporal dynamics of gamma. While the model failed with the often-used Poisson noise, Ornstein-Uhlenbeck noise applied to both the excitatory and the inhibitory populations replicated the duration and slowing of gamma and replicated the shape and stimulus dependencies. Thus, the temporal dynamics of gamma oscillations put constraints on the type and properties of underlying neural noise.

Mechanisms underlying pathological cortical bursts during metabolic depletion

Cortical activity depends upon a continuous supply of oxygen and other metabolic resources. Perinatal disruption of oxygen availability is a common clinical scenario in neonatal intensive care units, and a leading cause of lifelong disability. Pathological patterns of brain activity including burst suppression and seizures are a hallmark of the recovery period, yet the mechanisms by which these patterns arise remain poorly understood. Here, we use computational modeling of coupled metabolic-neuronal activity to explore the mechanisms by which oxygen depletion generates pathological brain activity. We find that restricting oxygen supply drives transitions from normal activity to several pathological activity patterns (isoelectric, burst suppression, and seizures), depending on the potassium supply. Trajectories through parameter space track key features of clinical electrophysiology recordings and reveal how infants with good recovery outcomes track toward normal parameter values, whereas the parameter values for infants with poor outcomes dwell around the pathological values. These findings open avenues for studying and monitoring the metabolically challenged infant brain, and deepen our understanding of the link between neuronal and metabolic activity.

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

INTRODUCTION

Accumulation and interaction of amyloid-beta (Aβ) and tau proteins during progression of Alzheimer's disease (AD) are shown to tilt neuronal circuits away from balanced excitation/inhibition (E/I). Current available techniques for noninvasive interrogation of E/I in the intact human brain, for example, magnetic resonance spectroscopy (MRS), are highly restrictive (i.e., limited spatial extent), have low temporal and spatial resolution and suffer from the limited ability to distinguish accurately between different neurotransmitters complicating its interpretation. As such, these methods alone offer an incomplete explanation of E/I. Recently, the aperiodic component of neural power spectrum, often referred to in the literature as the '1/f slope', has been described as a promising and scalable biomarker that can track disruptions in E/I potentially underlying a spectrum of clinical conditions, such as autism, schizophrenia, or epilepsy, as well as developmental E/I changes as seen in aging.

METHODS

Using 1/f slopes from resting-state spectral data and computational modeling, we developed a new method for inferring E/I alterations in AD.

RESULTS

We tested our method on recent freely and publicly available electroencephalography (EEG) and magnetoencephalography (MEG) datasets of patients with AD or prodromal disease and demonstrated the method's potential for uncovering regional patterns of abnormal excitatory and inhibitory parameters.

DISCUSSION

Our results provide a general framework for investigating circuit-level disorders in AD and developing therapeutic interventions that aim to restore the balance between excitation and inhibition.

Spiking Recurrent Neural Networks Represent Task-Relevant Neural Sequences in Rule-Dependent Computation

Background

Prefrontal cortical neurons play essential roles in performing rule-dependent tasks and working memory-based decision making.

Methods

Motivated by PFG recordings of task-performing mice, we developed an excitatory-inhibitory spiking recurrent neural network (SRNN) to perform a rule-dependent two-alternative forced choice (2AFC) task. We imposed several important biological constraints onto the SRNN, and adapted spike frequency adaptation (SFA) and SuperSpike gradient methods to train the SRNN efficiently.

Results

The trained SRNN produced emergent rule-specific tunings in single-unit representations, showing rule-dependent population dynamics that resembled experimentally observed data. Under varying test conditions, we manipulated the SRNN parameters or configuration in computer simulations, and we investigated the impacts of rule-coding error, delay duration, recurrent weight connectivity and sparsity, and excitation/inhibition (E/I) balance on both task performance and neural representations.

Conclusions

Overall, our modeling study provides a computational framework to understand neuronal representations at a fine timescale during working memory and cognitive control, and provides new experimentally testable hypotheses in future experiments.

Neurogenetic and multi-omic sources of overlap among sensation seeking, alcohol consumption, and alcohol use disorder

Sensation seeking is bidirectionally associated with levels of alcohol consumption in both adult and adolescent samples and shared neurobiological and genetic influences may in part explain this association. Links between sensation seeking and alcohol use disorder (AUD) may primarily manifest via increased alcohol consumption rather than through direct effects on increasing problems and consequences. Here the overlap between sensation seeking, alcohol consumption, and AUD was examined using multivariate modeling approaches for genome-wide association study (GWAS) summary statistics in conjunction with neurobiologically-informed analyses at multiple levels of investigation. Meta-analytic and genomic structural equation modeling (GenomicSEM) approaches were used to conduct GWAS of sensation seeking, alcohol consumption, and AUD. Resulting summary statistics were used in downstream analyses to examine shared brain tissue enrichment of heritability and genome-wide evidence of overlap (e.g., stratified GenomicSEM, RRHO, genetic correlations with neuroimaging phenotypes) and to identify genomic regions likely contributing to observed genetic overlap across traits (e.g., HMAGMA, LAVA). Across approaches, results supported shared neurogenetic architecture between sensation seeking and alcohol consumption characterized by overlapping enrichment of genes expressed in midbrain and striatal tissues and variants associated with increased cortical surface area. Alcohol consumption and AUD evidenced overlap in relation to variants associated with decreased frontocortical thickness. Finally, genetic mediation models provided evidence of alcohol consumption mediating associations between sensation seeking and AUD. This study extends previous research by examining critical sources of neurogenetic and multi-omic overlap among sensation seeking, alcohol consumption, and AUD which may underlie observed phenotypic associations.

Attractor-Like Dynamics Extracted from Human Electrocorticographic Recordings Underlie Computational Principles of Auditory Bistable Perception

In bistable perception, observers experience alternations between two interpretations of an unchanging stimulus. Neurophysiological studies of bistable perception typically partition neural measurements into stimulus-based epochs and assess neuronal differences between epochs based on subjects' perceptual reports. Computational studies replicate statistical properties of percept durations with modeling principles like competitive attractors or Bayesian inference. However, bridging neuro-behavioral findings with modeling theory requires the analysis of single-trial dynamic data. Here, we propose an algorithm for extracting nonstationary timeseries features from single-trial electrocorticography (ECoG) data. We applied the proposed algorithm to 5-min ECoG recordings from human primary auditory cortex obtained during perceptual alternations in an auditory triplet streaming task (six subjects: four male, two female). We report two ensembles of emergent neuronal features in all trial blocks. One ensemble consists of periodic functions that encode a stereotypical response to the stimulus. The other comprises more transient features and encodes dynamics associated with bistable perception at multiple time scales: minutes (within-trial alternations), seconds (duration of individual percepts), and milliseconds (switches between percepts). Within the second ensemble, we identified a slowly drifting rhythm that correlates with the perceptual states and several oscillators with phase shifts near perceptual switches. Projections of single-trial ECoG data onto these features establish low-dimensional attractor-like geometric structures invariant across subjects and stimulus types. These findings provide supporting neural evidence for computational models with oscillatory-driven attractor-based principles. The feature extraction techniques described here generalize across recording modality and are appropriate when hypothesized low-dimensional dynamics characterize an underlying neural system.SIGNIFICANCE STATEMENT Irrespective of the sensory modality, neurophysiological studies of multistable perception have typically investigated events time-locked to the perceptual switching rather than the time course of the perceptual states per se. Here, we propose an algorithm that extracts neuronal features of bistable auditory perception from largescale single-trial data while remaining agnostic to the subject's perceptual reports. The algorithm captures the dynamics of perception at multiple timescales, minutes (within-trial alternations), seconds (durations of individual percepts), and milliseconds (timing of switches), and distinguishes attributes of neural encoding of the stimulus from those encoding the perceptual states. Finally, our analysis identifies a set of latent variables that exhibit alternating dynamics along a low-dimensional manifold, similar to trajectories in attractor-based models for perceptual bistability.

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

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

Uncovering the organization of neural circuits with Generalized Phase Locking Analysis

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

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.

Efficient suppression of parkinsonian beta oscillations in a closed-loop model of deep brain stimulation with amplitude modulation

Introduction

Parkinson's disease (PD) is a movement disorder characterized by the pathological beta band (15-30 Hz) neural oscillations within the basal ganglia (BG). It is shown that the suppression of abnormal beta oscillations is correlated with the improvement of PD motor symptoms, which is a goal of standard therapies including deep brain stimulation (DBS). To overcome the stimulation-induced side effects and inefficiencies of conventional DBS (cDBS) and to reduce the administered stimulation current, closed-loop adaptive DBS (aDBS) techniques were developed. In this method, the frequency and/or amplitude of stimulation are modulated based on various disease biomarkers.

Methods

Here, by computational modeling of a cortico-BG-thalamic network in normal and PD conditions, we show that closed-loop aDBS of the subthalamic nucleus (STN) with amplitude modulation leads to a more effective suppression of pathological beta oscillations within the parkinsonian BG.

Results

Our results show that beta band neural oscillations are restored to their normal range and the reliability of the response of the thalamic neurons to motor cortex commands is retained due to aDBS with amplitude modulation. Furthermore, notably less stimulation current is administered during aDBS compared with cDBS due to a closed-loop control of stimulation amplitude based on the STN local field potential (LFP) beta activity.

Discussion

Efficient models of closed-loop stimulation may contribute to the clinical development of optimized aDBS techniques designed to reduce potential stimulation-induced side effects of cDBS in PD patients while leading to a better therapeutic outcome.

Improving the validity of neuroimaging decoding tests of invariant and configural neural representation

Many research questions in sensory neuroscience involve determining whether the neural representation of a stimulus property is invariant or specific to a particular stimulus context (e.g., Is object representation invariant to translation? Is the representation of a face feature specific to the context of other face features?). Between these two extremes, representations may also be context-tolerant or context-sensitive. Most neuroimaging studies have used operational tests in which a target property is inferred from a significant test against the null hypothesis of the opposite property. For example, the popular cross-classification test concludes that representations are invariant or tolerant when the null hypothesis of specificity is rejected. A recently developed neurocomputational theory suggests two insights regarding such tests. First, tests against the null of context-specificity, and for the alternative of context-invariance, are prone to false positives due to the way in which the underlying neural representations are transformed into indirect measurements in neuroimaging studies. Second, jointly performing tests against the nulls of invariance and specificity allows one to reach more precise and valid conclusions about the underlying representations, particularly when the null of invariance is tested using the fine-grained information from classifier decision variables rather than only accuracies (i.e., using the decoding separability test). Here, we provide empirical and computational evidence supporting both of these theoretical insights. In our empirical study, we use encoding of orientation and spatial position in primary visual cortex as a case study, as previous research has established that these properties are encoded in a context-sensitive way. Using fMRI decoding, we show that the cross-classification test produces false-positive conclusions of invariance, but that more valid conclusions can be reached by jointly performing tests against the null of invariance. The results of two simulations further support both of these conclusions. We conclude that more valid inferences about invariance or specificity of neural representations can be reached by jointly testing against both hypotheses, and using neurocomputational theory to guide the interpretation of results.

Towards a mesoscale physical modeling framework for stereotactic-EEG recordings

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

Rhythms of human attention and memory: An embedded process perspective

It remains a dogma in cognitive neuroscience to separate human attention and memory into distinct modules and processes. Here we propose that brain rhythms reflect the embedded nature of these processes in the human brain, as evident from their shared neural signatures: gamma oscillations (30-90 Hz) reflect sensory information processing and activated neural representations (memory items). The theta rhythm (3-8 Hz) is a pacemaker of explicit control processes (central executive), structuring neural information processing, bit by bit, as reflected in the theta-gamma code. By representing memory items in a sequential and time-compressed manner the theta-gamma code is hypothesized to solve key problems of neural computation: (1) attentional sampling (integrating and segregating information processing), (2) mnemonic updating (implementing Hebbian learning), and (3) predictive coding (advancing information processing ahead of the real time to guide behavior). In this framework, reduced alpha oscillations (8-14 Hz) reflect activated semantic networks, involved in both explicit and implicit mnemonic processes. Linking recent theoretical accounts and empirical insights on neural rhythms to the embedded-process model advances our understanding of the integrated nature of attention and memory - as the bedrock of human cognition.

Synaptic alterations in visual cortex reshape contrast-dependent gamma oscillations and inhibition-excitation ratio in a genetic mouse model of migraine

Background

Migraine affects a significant fraction of the world population, yet its etiology is not completely understood. In vitro results highlighted thalamocortical and intra-cortical glutamatergic synaptic gain-of-function associated with a monogenic form of migraine (familial-hemiplegic-migraine-type-1: FHM1). However, how these alterations reverberate on cortical activity remains unclear. As altered responsivity to visual stimuli and abnormal processing of visual sensory information are common hallmarks of migraine, herein we investigated the effects of FHM1-driven synaptic alterations in the visual cortex of awake mice.

Methods

We recorded extracellular field potentials from the primary visual cortex (V1) of head-fixed awake FHM1 knock-in (n = 12) and wild type (n = 12) mice in response to square-wave gratings with different visual contrasts. Additionally, we reproduced in silico the obtained experimental results with a novel spiking neurons network model of mouse V1, by implementing in the model both the synaptic alterations characterizing the FHM1 genetic mouse model adopted.

Results

FHM1 mice displayed similar amplitude but slower temporal evolution of visual evoked potentials. Visual contrast stimuli induced a lower increase of multi-unit activity in FHM1 mice, while the amount of information content about contrast level remained, however, similar to WT.

Spectral analysis of the local field potentials revealed an increase in the β/low γ range of WT mice following the abrupt reversal of contrast gratings. Such frequency range transitioned to the high γ range in FHM1 mice. Despite this change in the encoding channel, these oscillations preserved the amount of information conveyed about visual contrast. The computational model showed how these network effects may arise from a combination of changes in thalamocortical and intra-cortical synaptic transmission, with the former inducing a lower cortical activity and the latter inducing the higher frequencies ɣ oscillations.

Conclusions

Contrast-driven ɣ modulation in V1 activity occurs at a much higher frequency in FHM1. This is likely to play a role in the altered processing of visual information. Computational studies suggest that this shift is specifically due to enhanced cortical excitatory transmission. Our network model can help to shed light on the relationship between cellular and network levels of migraine neural alterations.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s10194-022-01495-9.

A microcircuit model involving parvalbumin, somatostatin, and vasoactive intestinal polypeptide inhibitory interneurons for the modulation of neuronal oscillation during visual processing

Various subtypes of inhibitory interneurons contact one another to organize cortical networks. Most cortical inhibitory interneurons express 1 of 3 genes: parvalbumin (PV), somatostatin (SOM), or vasoactive intestinal polypeptide (VIP). This diversity of inhibition allows the flexible regulation of neuronal responses within and between cortical areas. However, the exact roles of these interneuron subtypes and of excitatory pyramidal (Pyr) neurons in regulating neuronal network activity and establishing perception (via interactions between feedforward sensory and feedback attentional signals) remain largely unknown. To explore the regulatory roles of distinct neuronal types in cortical computation, we developed a computational microcircuit model with biologically plausible visual cortex layers 2/3 that combined Pyr neurons and the 3 inhibitory interneuron subtypes to generate network activity. In simulations with our model, inhibitory signals from PV and SOM neurons preferentially induced neuronal firing at gamma (30-80 Hz) and beta (20-30 Hz) frequencies, respectively, in agreement with observed physiological results. Furthermore, our model indicated that rapid inhibition from VIP to SOM subtypes underlies marked attentional modulation for low-gamma frequency (30-50 Hz) in Pyr neuron responses. Our results suggest the distinct but cooperative roles of inhibitory interneuron subtypes in the establishment of visual perception.

A personalizable autonomous neural mass model of epileptic seizures

Work in the last two decades has shown that neural mass models (NMM) can realistically reproduce and explain epileptic seizure transitions as recorded by electrophysiological methods (EEG, SEEG). In previous work, advances were achieved by increasing excitation and heuristically varying network inhibitory coupling parameters in the models. Based on these early studies, we provide a laminar NMM capable of realistically reproducing the electrical activity recorded by SEEG in the epileptogenic zone during interictal to ictal states. With the exception of the external noise input into the pyramidal cell population, the model dynamics are autonomous. By setting the system at a point close to bifurcation, seizure-like transitions are generated, including pre-ictal spikes, low voltage fast activity, and ictal rhythmic activity. A novel element in the model is a physiologically motivated algorithm for chloride dynamics: the gain of GABAergic post-synaptic potentials is modulated by the pathological accumulation of chloride in pyramidal cells due to high inhibitory input and/or dysfunctional chloride transport. In addition, in order to simulate SEEG signals for comparison with real seizure recordings, the NMM is embedded first in a layered model of the neocortex and then in a realistic physical model. We compare modeling results with data from four epilepsy patient cases. By including key pathophysiological mechanisms, the proposed framework captures succinctly the electrophysiological phenomenology observed in ictal states, paving the way for robust personalization methods based on NMMs.

Cell to network computational model of the epileptic human hippocampus suggests specific roles of network and channel dysfunctions in the ictal and interictal oscillations

The mechanisms underlying the generation of hippocampal epileptic seizures and interictal events and their interactions with the sleep-wake cycle are not yet fully understood. Indeed, medial temporal lobe epilepsy is associated with hippocampal abnormalities both at the neuronal (channelopathies, impaired potassium and chloride dynamics) and network level (neuronal and axonal loss, mossy fiber sprouting), with more frequent seizures during wakefulness compared with slow-wave sleep. In this article, starting from our previous computational modeling work of the hippocampal formation based on realistic topology and synaptic connectivity, we study the role of micro- and mesoscale pathological conditions of the epileptic hippocampus in the generation and maintenance of seizure-like theta and interictal oscillations. We show, through the simulations of hippocampal activity during slow-wave sleep and wakefulness that: (i) both mossy fiber sprouting and sclerosis account for seizure-like theta activity, (ii) but they have antagonist effects (seizure-like activity occurrence increases with sprouting but decreases with sclerosis), (iii) though impaired potassium and chloride dynamics have little influence on the generation of seizure-like activity, they do play a role on the generation of interictal patterns, and (iv) seizure-like activity and fast ripples are more likely to occur during wakefulness and interictal spikes during sleep.

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.

Extracting synchronized neuronal activity from local field potentials based on a marked point process framework

OBJECTIVE

Brain-machine interfaces (BMIs) translate neural activity into motor commands to restore motor functions for people with paralysis. Local field potentials (LFPs) are promising for long-term BMIs, since the quality of the recording lasts longer than single neuronal spikes. Inferring neuronal spike activity from population activities such as LFPs is challenging, because LFPs stem from synaptic currents flowing in the neural tissue produced by various neuronal ensembles and reflect neural synchronization. Existing studies that combine LFPs with spikes leverage the spectrogram of LFPs, which can neither detect the transient characteristics of LFP features (here, neuromodulation in a specific frequency band) with high accuracy, nor correlate them with relevant neuronal activity with a sufficient time resolution.

APPROACH

We propose a feature extraction and validation framework to directly extract LFP neuromodulations related to synchronized spike activity using recordings from the primary motor cortex of six Sprague Dawley (SD) rats during a lever-press task. We first select important LFP frequency bands relevant to behavior, and then implement a marked point process (MPP) methodology to extract transient LFP neuromodulations. We validate the LFP feature extraction by examining the correlation with the pairwise synchronized firing probability of important neurons, which are selected according to their contribution to behavioral decoding. The highly correlated synchronized firings identified by the LFP neuromodulations are fed into a decoder to check whether they can serve as a reliable neural data source for movement decoding.

MAIN RESULTS

We find that the gamma band (30-80Hz) LFP neuromodulations demonstrate significant correlation with synchronized firings. Compared with traditional spectrogram-based method, the higher-temporal resolution MPP method captures the synchronized firing patterns with fewer false alarms, and demonstrates significantly higher correlation than single neuron spikes. The decoding performance using the synchronized neuronal firings identified by the LFP neuromodulations can reach 90% compared to the full recorded neuronal ensembles.

SIGNIFICANCE

Our proposed framework successfully extracts the sparse LFP neuromodulations that can identify temporal synchronized neuronal spikes with high correlation. The identified neuronal spike pattern demonstrates high decoding performance, which reveals the possibility of using LFP as an effective modality for long-term BMI decoding.

Waveform detection by deep learning reveals multi-area spindles that are selectively modulated by memory load

Sleep is generally considered to be a state of large-scale synchrony across thalamus and neocortex; however, recent work has challenged this idea by reporting isolated sleep rhythms such as slow oscillations and spindles. What is the spatial scale of sleep rhythms? To answer this question, we adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves in high-noise settings for analysis of neural recordings in sleep. We then studied sleep spindles in non-human primate electrocorticography (ECoG), human electroencephalogram (EEG), and clinical intracranial electroencephalogram (iEEG) recordings in the human. Within each recording type, we find widespread spindles occur much more frequently than previously reported. We then analyzed the spatiotemporal patterns of these large-scale, multi-area spindles and, in the EEG recordings, how spindle patterns change following a visual memory task. Our results reveal a potential role for widespread, multi-area spindles in consolidation of memories in networks widely distributed across primate cortex.

The brain processes memories as we sleep, generating rhythms of electrical activity called ‘sleep spindles’. Sleep spindles were long thought to be a state where the entire brain was fully synchronized by this rhythm. This was based on EEG recordings, short for electroencephalogram, a technique that uses electrodes on the scalp to measure electrical activity in the outermost layer of the brain, the cortex. But more recent intracranial recordings of people undergoing brain surgery have challenged this idea and suggested that sleep spindles may not be a state of global brain synchronization, but rather localised to specific areas.

Mofrad et al. sought to clarify the extent to which spindles co-occur at multiple sites in the brain, which could shed light on how networks of neurons coordinate memory storage during sleep. To analyse highly variable brain wave recordings, Mofrad et al. adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves. The resulting algorithm, designed to more sensitively detect spindles amongst other brain activity, was then applied to a range of sleep recordings from humans and macaque monkeys.

The analyses revealed that widespread and complex patterns of spindle rhythms, spanning multiple areas in the cortex of the brain, actually appear much more frequently than previously thought. This finding was consistent across all the recordings analysed, even recordings under the skull, which provide the clearest window into brain circuits. Further analyses found that these multi-area spindles occurred more often in sleep after people had completed tasks that required holding many visual scenes in memory, as opposed to control conditions with fewer visual scenes.

In summary, Mofrad et al. show that neuroscientists had previously not appreciated the complex and dynamic patterns in this sleep rhythm. These patterns in sleep spindles may be able to adapt based on the demands needed for memory storage, and this will be the subject of future work. Moreover, the findings support the idea that sleep spindles help coordinate the consolidation of memories in brain circuits that stretch across the cortex. Understanding this mechanism may provide insights into how memory falters in aging and sleep-related diseases, such as Alzheimer’s disease. Lastly, the algorithm developed by Mofrad et al. stands to be a useful tool for analysing other rhythmic waveforms in noisy recordings.

Spike-Gamma Phase Relationship in the Visual Cortex

Gamma oscillations (30-70 Hz) have been hypothesized to play a role in cortical function. Most of the proposed mechanisms involve rhythmic modulation of neuronal excitability at gamma frequencies, leading to modulation of spike timing relative to the rhythm. I first show that the gamma band could be more privileged than other frequencies in observing spike-field interactions even in the absence of genuine gamma rhythmicity and discuss several biases in spike-gamma phase estimation. I then discuss the expected spike-gamma phase according to several hypotheses. Inconsistent with the phase-coding hypothesis (but not with others), the spike-gamma phase does not change with changes in stimulus intensity or attentional state, with spikes preferentially occurring 2-4 ms before the trough, but with substantial variability. However, this phase relationship is expected even when gamma is a byproduct of excitatory-inhibitory interactions. Given that gamma occurs in short bursts, I argue that the debate over the role of gamma is a matter of semantics. Expected final online publication date for the Annual Review of Vision Science, Volume 8 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Columnar Localization and Laminar Origin of Cortical Surface Electrical Potentials

Electrocorticography (ECoG) methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. To address this gap, we recorded from rat auditory cortex using customized μECoG, and simulated cortical surface electrical potentials with a full-scale, biophysically detailed cortical column model. Experimentally, μECoG-derived auditory representations were tonotopically organized and signals were anisotropically localized to less than or equal to ±200 μm, that is, a single cortical column. Biophysical simulations reproduce experimental findings and indicate that neurons in cortical layers V and VI contribute ∼85% of evoked high-gamma signal recorded at the surface. Cell number and synchrony were the primary biophysical properties determining laminar contributions to evoked μECoG signals, whereas distance was only a minimal factor. Thus, evoked μECoG signals primarily originate from neurons in the infragranular layers of a single cortical column.SIGNIFICANCE STATEMENT ECoG methodologically bridges basic neuroscience and understanding of human brains in health and disease. However, the localization of ECoG signals across the surface of the brain and the spatial distribution of their generating neuronal sources are poorly understood. We investigated the localization and origins of sensory-evoked ECoG responses. We experimentally found that ECoG responses were anisotropically localized to a cortical column. Biophysically detailed simulations revealed that neurons in layers V and VI were the primary sources of evoked ECoG responses. These results indicate that evoked ECoG high-gamma responses are primarily generated by the population spike rate of pyramidal neurons in layers V and VI of single cortical columns and highlight the possibility of understanding how microscopic sources produce mesoscale signals.

Modeling multiscale causal interactions between spiking and field potential signals during behavior

Objective.Brain recordings exhibit dynamics at multiple spatiotemporal scales, which are measured with spike trains and larger-scale field potential signals. To study neural processes, it is important to identify and model causal interactions not only at a single scale of activity, but also across multiple scales, i.e. between spike trains and field potential signals. Standard causality measures are not directly applicable here because spike trains are binary-valued but field potentials are continuous-valued. It is thus important to develop computational tools to recover multiscale neural causality during behavior, assess their performance on neural datasets, and study whether modeling multiscale causalities can improve the prediction of neural signals beyond what is possible with single-scale causality.Approach.We design a multiscale model-based Granger-like causality method based on directed information and evaluate its success both in realistic biophysical spike-field simulations and in motor cortical datasets from two non-human primates (NHP) performing a motor behavior. To compute multiscale causality, we learn point-process generalized linear models that predict the spike events at a given time based on the history of both spike trains and field potential signals. We also learn linear Gaussian models that predict the field potential signals at a given time based on their own history as well as either the history of binary spike events or that of latent firing rates.Main results.We find that our method reveals the true multiscale causality network structure in biophysical simulations despite the presence of model mismatch. Further, models with the identified multiscale causalities in the NHP neural datasets lead to better prediction of both spike trains and field potential signals compared to just modeling single-scale causalities. Finally, we find that latent firing rates are better predictors of field potential signals compared with the binary spike events in the NHP datasets.Significance.This multiscale causality method can reveal the directed functional interactions across spatiotemporal scales of brain activity to inform basic science investigations and neurotechnologies.

Shape analysis of gamma rhythm supports a superlinear inhibitory regime in an inhibition-stabilized network

Visual inspection of stimulus-induced gamma oscillations (30–70 Hz) often reveals a non-sinusoidal shape. Such distortions are a hallmark of non-linear systems and are also observed in mean-field models of gamma oscillations. A thorough characterization of the shape of the gamma cycle can therefore provide additional constraints on the operating regime of such models. However, the gamma waveform has not been quantitatively characterized, partially because the first harmonic of gamma, which arises because of the non-sinusoidal nature of the signal, is typically weak and gets masked due to a broadband increase in power related to spiking. To address this, we recorded local field potential (LFP) from the primary visual cortex (V1) of two awake female macaques while presenting full-field gratings or iso-luminant chromatic hues that produced huge gamma oscillations with prominent peaks at harmonic frequencies in the power spectra. We found that gamma and its first harmonic always maintained a specific phase relationship, resulting in a distinctive shape with a sharp trough and a shallow peak. Interestingly, a Wilson-Cowan (WC) model operating in an inhibition stabilized mode could replicate this shape, but only when the inhibitory population operated in the super-linear regime, as predicted recently. However, another recently developed model of gamma that operates in a linear regime driven by stochastic noise failed to produce salient harmonics or the observed shape. Our results impose additional constraints on models that generate gamma oscillations and their operating regimes.

Gamma rhythm is not sinusoidal. Understanding these distortions could provide clues about the cortical network that generates the rhythm. Here, we use harmonic phase analysis to describe these waveforms quantitatively and show that gamma rhythm recorded from the primary visual cortex of macaques has a signature arch shaped waveform, with a sharp trough and a shallow peak, when visual stimuli such as full-screen plain hues and achromatic gratings are presented. This arch shaped waveform is observed over a wide range of stimuli, despite the variation in power and frequency of the rhythm. We then compare two population rate models that have been used to accurately describe the stimulus dependencies of gamma rhythm and show that this arch shaped waveform is obtained only in one of those models. Further, the waveform shape is dependent on the operating domain of the system. Therefore, shape analysis provides additional constraints on cortical models and their operating regimes.

Coordinated hippocampal-thalamic-cortical communication crucial for engram dynamics underneath systems consolidation

Systems consolidation refers to the time-dependent reorganization of memory representations or engrams across brain regions. Despite recent advancements in unravelling this process, the exact mechanisms behind engram dynamics and the role of associated pathways remain largely unknown. Here we propose a biologically-plausible computational model to address this knowledge gap. By coordinating synaptic plasticity timescales and incorporating a hippocampus-thalamus-cortex circuit, our model is able to couple engram reactivations across these regions and thereby reproduce key dynamics of cortical and hippocampal engram cells along with their interdependencies. Decoupling hippocampal-thalamic-cortical activity disrupts systems consolidation. Critically, our model yields testable predictions regarding hippocampal and thalamic engram cells, inhibitory engrams, thalamic inhibitory input, and the effect of thalamocortical synaptic coupling on retrograde amnesia induced by hippocampal lesions. Overall, our results suggest that systems consolidation emerges from coupled reactivations of engram cells in distributed brain regions enabled by coordinated synaptic plasticity timescales in multisynaptic subcortical-cortical circuits.

Synchronization Through Uncorrelated Noise in Excitatory-Inhibitory Networks

Gamma rhythms play a major role in many different processes in the brain, such as attention, working memory, and sensory processing. While typically considered detrimental, counterintuitively noise can sometimes have beneficial effects on communication and information transfer. Recently, Meng and Riecke showed that synchronization of interacting networks of inhibitory neurons in the gamma band (i.e., gamma generated through an ING mechanism) increases while synchronization within these networks decreases when neurons are subject to uncorrelated noise. However, experimental and modeling studies point towardz an important role of the pyramidal-interneuronal network gamma (PING) mechanism in the cortex. Therefore, we investigated the effect of uncorrelated noise on the communication between excitatory-inhibitory networks producing gamma oscillations via a PING mechanism. Our results suggest that, at least in a certain range of noise strengths and natural frequency differences between the regions, synaptic noise can have a supporting role in facilitating inter-regional communication, similar to the ING case for a slightly larger parameter range. Furthermore, the noise-induced synchronization between networks is generated via a different mechanism than when synchronization is mediated by strong synaptic coupling. Noise-induced synchronization is achieved by lowering synchronization within networks which allows the respective other network to impose its own gamma rhythm resulting in synchronization between networks.

Criticality enhances the multilevel reliability of stimulus responses in cortical neural networks

Cortical neural networks exhibit high internal variability in spontaneous dynamic activities and they can robustly and reliably respond to external stimuli with multilevel features–from microscopic irregular spiking of neurons to macroscopic oscillatory local field potential. A comprehensive study integrating these multilevel features in spontaneous and stimulus–evoked dynamics with seemingly distinct mechanisms is still lacking. Here, we study the stimulus–response dynamics of biologically plausible excitation–inhibition (E–I) balanced networks. We confirm that networks around critical synchronous transition states can maintain strong internal variability but are sensitive to external stimuli. In this dynamical region, applying a stimulus to the network can reduce the trial-to-trial variability and shift the network oscillatory frequency while preserving the dynamical criticality. These multilevel features widely observed in different experiments cannot simultaneously occur in non-critical dynamical states. Furthermore, the dynamical mechanisms underlying these multilevel features are revealed using a semi-analytical mean-field theory that derives the macroscopic network field equations from the microscopic neuronal networks, enabling the analysis by nonlinear dynamics theory and linear noise approximation. The generic dynamical principle revealed here contributes to a more integrative understanding of neural systems and brain functions and incorporates multimodal and multilevel experimental observations. The E–I balanced neural network in combination with the effective mean-field theory can serve as a mechanistic modeling framework to study the multilevel neural dynamics underlying neural information and cognitive processes.

The complexity and variability of brain dynamical activity range from neuronal spiking and neural avalanches to oscillatory local field potentials of local neural circuits in both spontaneous and stimulus-evoked states. Such multilevel variable brain dynamics are functionally and behaviorally relevant and are principal components of the underlying circuit organization. To more comprehensively clarify their neural mechanisms, we use a bottom-up approach to study the stimulus–response dynamics of neural circuits. Our model assumes the following key biologically plausible components: excitation–inhibition (E–I) neuronal interaction and chemical synaptic coupling. We show that the circuits with E–I balance have a special dynamic sub-region, the critical region. Circuits around this region could account for the emergence of multilevel brain response patterns, both ongoing and stimulus-induced, observed in different experiments, including the reduction of trial-to-trial variability, effective modulation of gamma frequency, and preservation of criticality in the presence of a stimulus. We further analyze the corresponding nonlinear dynamical principles using a novel and highly generalizable semi-analytical mean-field theory. Our computational and theoretical studies explain the cross-level brain dynamical organization of spontaneous and evoked states in a more integrative manner.

Directed functional and structural connectivity in a large-scale model for the mouse cortex

Inferring the structural connectivity from electrophysiological measurements is a fundamental challenge in systems neuroscience. Directed functional connectivity measures, such as the generalized partial directed coherence (GPDC), provide estimates of the causal influence between areas. However, the relation between causality estimates and structural connectivity is still not clear. We analyzed this problem by evaluating the effectiveness of GPDC to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. The model contains 19 cortical areas composed of spiking neurons, with areas connected by long-range projections with weights obtained from a tract-tracing cortical connectome. We show that GPDC values provide a reasonable estimate of structural connectivity, with an average Pearson correlation over simulations of 0.74. Moreover, even in a typical electrophysiological recording scenario containing five areas, the mean correlation was above 0.6. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

We analyzed the relationship between structural and directed functional connectivity by evaluating the effectiveness of generalized partial directed coherence (GPDC) to estimate the connectivity of a ground-truth, data-constrained computational model of a large-scale network model of the mouse cortex. We show that GPDC values provide a reasonable estimate of structural connectivity even in a typical electrophysiological recording scenario containing few areas. These results suggest that it may be possible to empirically estimate structural connectivity from functional connectivity even when detailed whole-brain recordings are not achievable.

Plasticity impairment exposes CA3 vulnerability in a hippocampal network model of mild traumatic brain injury

Proper function of the hippocampus is critical for executing cognitive tasks such as learning and memory. Traumatic brain injury (TBI) and other neurological disorders are commonly associated with cognitive deficits and hippocampal dysfunction. Although there are many existing models of individual subregions of the hippocampus, few models attempt to integrate the primary areas into one system. In this work, we developed a computational model of the hippocampus, including the dentate gyrus, CA3, and CA1. The subregions are represented as an interconnected neuronal network, incorporating well-characterized ex vivo slice electrophysiology into the functional neuron models and well-documented anatomical connections into the network structure. In addition, since plasticity is foundational to the role of the hippocampus in learning and memory as well as necessary for studying adaptation to injury, we implemented spike-timing-dependent plasticity among the synaptic connections. Our model mimics key features of hippocampal activity, including signal frequencies in the theta and gamma bands and phase-amplitude coupling in area CA1. We also studied the effects of spike-timing-dependent plasticity impairment, a potential consequence of TBI, in our model and found that impairment decreases broadband power in CA3 and CA1 and reduces phase coherence between these two subregions, yet phase-amplitude coupling in CA1 remains intact. Altogether, our work demonstrates characteristic hippocampal activity with a scaled network model of spiking neurons and reveals the sensitive balance of plasticity mechanisms in the circuit through one manifestation of mild traumatic injury.