Origin of slow spontaneous resting-state neuronal fluctuations in brain networks

Chronic full-band recordings with graphene microtransistors as neural interfaces for discrimination of brain states

Brain states such as sleep, anesthesia, wakefulness, or coma are characterized by specific patterns of cortical activity dynamics, from local circuits to full-brain emergent properties. We previously demonstrated that full-spectrum signals, including the infraslow component (DC, direct current-coupled), can be recorded acutely in multiple sites using flexible arrays of graphene solution-gated field-effect transistors (gSGFETs). Here, we performed chronic implantation of 16-channel gSGFET arrays over the rat cerebral cortex and recorded full-band neuronal activity with two objectives: (1) to test the long-term stability of implanted devices; and (2) to investigate full-band activity during the transition across different levels of anesthesia. First, we demonstrate it is possible to record full-band signals with stability, fidelity, and spatiotemporal resolution for up to 5.5 months using chronic epicortical gSGFET implants. Second, brain states generated by progressive variation of levels of anesthesia could be identified as traditionally using the high-pass filtered (AC, alternating current-coupled) spectrogram: from synchronous slow oscillations in deep anesthesia through to asynchronous activity in the awake state. However, the DC signal introduced a highly significant improvement for brain-state discrimination: the DC band provided an almost linear information prediction of the depth of anesthesia, with about 85% precision, using a trained algorithm. This prediction rose to about 95% precision when the full-band (AC + DC) spectrogram was taken into account. We conclude that recording infraslow activity using gSGFET interfaces is superior for the identification of brain states, and further supports the preclinical and clinical use of graphene neural interfaces for long-term recordings of cortical activity.

Spatiotemporal dedifferentiation of the global brain signal topography along the adult lifespan

Age-related variations in many regions and/or networks of the human brain have been uncovered using resting-state functional magnetic resonance imaging. However, these findings did not account for the dynamical effect the brain's global activity (global signal [GS]) causes on local characteristics, which is measured by GS topography. To address this gap, we tested GS topography including its correlation with age using a large-scale cross-sectional adult lifespan dataset (n = 492). Both GS topography and its variation with age showed frequency-specific patterns, reflecting the spatiotemporal characteristics of the dynamic change of GS topography with age. A general trend toward dedifferentiation of GS topography with age was observed in both spatial (i.e., less differences of GS between different regions) and temporal (i.e., less differences of GS between different frequencies) dimensions. Further, methodological control analyses suggested that although most age-related dedifferentiation effects remained across different preprocessing strategies, some were triggered by neuro-vascular coupling and physiological noises. Together, these results provide the first evidence for age-related effects on global brain activity and its topographic-dynamic representation in terms of spatiotemporal dedifferentiation.

Spontaneous Brain Activity Alterations in First-Episode Psychosis: A Meta-analysis of Functional Magnetic Resonance Imaging Studies

BACKGROUND AND HYPOTHESIS

Several studies have shown that spontaneous brain activity, including the total and fractional amplitude of low-frequency fluctuations (LFF) and regional homogeneity (ReHo), is altered in psychosis. Nonetheless, neuroimaging results show a high heterogeneity. For this reason, we gathered the extant literature on spontaneous brain activity in first-episode psychosis (FEP), where the effects of long-term treatment and chronic disease are minimal.

STUDY DESIGN

A systematic research was conducted on PubMed, Scopus, and Web of Science to identify studies exploring spontaneous brain activity and local connectivity in FEP estimated using functional magnetic resonance imaging. 20 LFF and 15 ReHo studies were included. Coordinate-Based Activation Likelihood Estimation Meta-Analyses stratified by brain measures, age (adolescent vs adult), and drug-naïve status were performed to identify spatially-convergent alterations in spontaneous brain activity in FEP.

STUDY RESULTS

We found a significant increase in LFF in FEP compared to healthy controls (HC) in the right striatum and in ReHo in the left striatum. When pooling together all studies on LFF and ReHo, spontaneous brain activity was increased in the bilateral striatum and superior and middle frontal gyri and decreased in the right precentral gyrus and the right inferior frontal gyrus compared to HC. These results were also replicated in the adult and drug-naïve samples.

CONCLUSIONS

Abnormalities in the frontostriatal circuit are present in early psychosis independently of treatment status. Our findings support the view that altered frontostriatal can represent a core neural alteration of the disorder and could be a target of treatment.

The Changes of Histone Methylation Induced by Adolescent Social Stress Regulate the Resting-State Activity in mPFC

Early-life stress can lead to sustained alterations in regional resting-state brain functions, but the underlying molecular mechanism remains unclear. Stress can also induce sustained changes in epigenetic modifications across brain regions, which are not limited to a few genes; rather, they often tend to produce global levels of change. The functional implication of these changes also remains to be elucidated. We hypothesize that global epigenetic changes may partly modulate the resting-state functions of brain regions to influence behavior. To test this hypothesis, we used an adolescent social stress (ASS) model in mice and examined the relationship between epigenetic modifications and regional resting-state brain activity using resting-state functional magnetic resonance imaging (rs-fMRI). The results showed that, compared to the control mice, the stressed mice showed increased anxiety and social avoidance behaviors and greater levels of dimethylation of histone H3 at lysine 9 (H3K9me2) in the medial prefrontal cortex (mPFC). In addition, the resting-state activity represented by the amplitude of low-frequency fluctuation (ALFF) was significantly lower in the mPFC of stressed mice. To verify the relationship of H3K9me2 and ALFF, the specific inhibition of H3Kme2 was performed by using the drug UNC0642, which reversed the anxiety behavior induced by ASS and significantly increase the ALFF value of mPFC in both normal and ASS animals. Our study is the first to report an association between histone modifications and rs-fMRI findings, providing a new perspective for understanding of the significance of regional brain epigenetic changes and a possible molecular explanation for rs-fMRI findings.

The decoupling between hemodynamic parameters and neural activity implies a complex origin of spontaneous brain oscillations

Introduction

Spontaneous low-frequency oscillations play a key role in brain activity. However, the underlying mechanism and origin of low-frequency oscillations remain under debate.

Methods

Optical imaging and an electrophysiological recording system were combined to investigate spontaneous oscillations in the hemodynamic parameters and neuronal activity of awake and anesthetized mice after Nω-nitro-L-arginine methyl ester (L-NAME) administration.

Results

The spectrum of local field potential (LFP) signals was significantly changed by L-NAME, which was further corroborated by the increase in energy and spatial synchronization. The important finding was that L-NAME triggered regular oscillations in both LFP signals and hemodynamic signals. Notably, the frequency peak of hemodynamic signals can be different from that of LFP oscillations in awake mice.

Discussion

A model of the neurovascular system was proposed to interpret this mismatch of peak frequencies, supporting the view that spontaneous low-frequency oscillations arise from multiple sources.

Astrocytic chloride is brain state dependent and modulates inhibitory neurotransmission in mice

Information transfer within neuronal circuits depends on the balance and recurrent activity of excitatory and inhibitory neurotransmission. Chloride (Cl-) is the major central nervous system (CNS) anion mediating inhibitory neurotransmission. Astrocytes are key homoeostatic glial cells populating the CNS, although the role of these cells in regulating excitatory-inhibitory balance remains unexplored. Here we show that astrocytes act as a dynamic Cl- reservoir regulating Cl- homoeostasis in the CNS. We found that intracellular chloride concentration ([Cl-]i) in astrocytes is high and stable during sleep. In awake mice astrocytic [Cl-]i is lower and exhibits large fluctuation in response to both sensory input and motor activity. Optogenetic manipulation of astrocytic [Cl-]i directly modulates neuronal activity during locomotion or whisker stimulation. Astrocytes thus serve as a dynamic source of extracellular Cl- available for GABAergic transmission in awake mice, which represents a mechanism for modulation of the inhibitory tone during sustained neuronal activity.

Brain-wide network analysis of resting-state neuromagnetic data

The present study performed a brain-wide network analysis of resting-state magnetoencephalograms recorded from 53 healthy participants to visualize elaborate brain maps of phase- and amplitude-derived graph-theory metrics at different frequencies. To achieve this, we conducted a vertex-wise computation of threshold-independent graph metrics by combining proportional thresholding and a conjunction analysis and applied them to a correlation analysis of age and brain networks. Source power showed a frequency-dependent cortical distribution. Threshold-independent graph metrics derived from phase- and amplitude-based connectivity showed similar or different distributions depending on frequency. Vertex-wise age-brain correlation maps revealed that source power at the beta band and the amplitude-based degree at the alpha band changed with age in local regions. The present results indicate that a brain-wide analysis of neuromagnetic data has the potential to reveal neurophysiological network features in the human brain in a resting state.

Aberrant protein S-nitrosylation contributes to hyperexcitability-induced synaptic damage in Alzheimer's disease: Mechanistic insights and potential therapies

Alzheimer's disease (AD) is arguably the most common cause of dementia in the elderly and is marked by progressive synaptic degeneration, which in turn leads to cognitive decline. Studies in patients and in various AD models have shown that one of the early signatures of AD is neuronal hyperactivity. This excessive electrical activity contributes to dysregulated neural network function and synaptic damage. Mechanistically, evidence suggests that hyperexcitability accelerates production of reactive oxygen species (ROS) and reactive nitrogen species (RNS) that contribute to neural network impairment and synapse loss. This review focuses on the pathways and molecular changes that cause hyperexcitability and how RNS-dependent posttranslational modifications, represented predominantly by protein S-nitrosylation, mediate, at least in part, the deleterious effects of hyperexcitability on single neurons and the neural network, resulting in synaptic loss in AD.

Computational models of Idling brain activity for memory processing

Studying the underlying neural mechanisms of cognitive functions of the brain is one of the central questions in modern biology. Moreover, it has significantly impacted the development of novel technologies in artificial intelligence. Spontaneous activity is a unique feature of the brain and is currently lacking in many artificially constructed intelligent machines. Spontaneous activity may represent the brain's idling states, which are internally driven by neuronal networks and possibly participate in offline processing during awake, sleep, and resting states. Evidence is accumulating that the brain's spontaneous activity is not mere noise but part of the mechanisms to process information about previous experiences. A bunch of literature has shown how previous sensory and behavioral experiences influence the subsequent patterns of brain activity with various methods in various animals. It seems, however, that the patterns of neural activity and their computational roles differ significantly from area to area and from function to function. In this article, I review the various forms of the brain's spontaneous activity, especially those observed during memory processing, and some attempts to model the generation mechanisms and computational roles of such activities.

Infraslow closed-loop brain training for anxiety and depression (ISAD): a protocol for a randomized, double-blind, sham-controlled pilot trial in adult females with internalizing disorders

Background

The core intrinsic connectivity networks (core-ICNs), encompassing the default-mode network (DMN), salience network (SN) and central executive network (CEN), have been shown to be dysfunctional in individuals with internalizing disorders (IDs, e.g. major depressive disorder, MDD; generalized anxiety disorder, GAD; social anxiety disorder, SOC). As such, source-localized, closed-loop brain training of electrophysiological signals, also known as standardized low-resolution electromagnetic tomography (sLORETA) neurofeedback (NFB), targeting key cortical nodes within these networks has the potential to reduce symptoms associated with IDs and restore normal core ICN function. We intend to conduct a randomized, double-blind (participant and assessor), sham-controlled, parallel-group (3-arm) trial of sLORETA infraslow (<0.1 Hz) fluctuation neurofeedback (sLORETA ISF-NFB) 3 times per week over 4 weeks in participants (n=60) with IDs. Our primary objectives will be to examine patient-reported outcomes (PROs) and neurophysiological measures to (1) compare the potential effects of sham ISF-NFB to either genuine 1-region ISF-NFB or genuine 2-region ISF-NFB, and (2) assess for potential associations between changes in PRO scores and modifications of electroencephalographic (EEG) activity/connectivity within/between the trained regions of interest (ROIs). As part of an exploratory analysis, we will investigate the effects of additional training sessions and the potential for the potentiation of the effects over time.

Methods

We will randomly assign participants who meet the criteria for MDD, GAD, and/or SOC per the MINI (Mini International Neuropsychiatric Interview for DSM-5) to one of three groups: (1) 12 sessions of posterior cingulate cortex (PCC) ISF-NFB up-training (n=15), (2) 12 sessions of concurrent PCC ISF up-training and dorsal anterior cingulate cortex (dACC) ISF-NFB down-training (n=15), or (3) 6 sessions of yoked-sham training followed by 6 sessions genuine ISF-NFB (n=30). Transdiagnostic PROs (Hospital Anxiety and Depression Scale, HADS; Inventory of Depression and Anxiety Symptoms – Second Version, IDAS-II; Multidimensional Emotional Disorder Inventory, MEDI; Intolerance of Uncertainty Scale – Short Form, IUS-12; Repetitive Thinking Questionnaire, RTQ-10) as well as resting-state neurophysiological measures (full-band EEG and ECG) will be collected from all subjects during two baseline sessions (approximately 1 week apart) then at post 6 sessions, post 12 sessions, and follow-up (1 month later). We will employ Bayesian methods in R and advanced source-localisation software (i.e. exact low-resolution brain electromagnetic tomography; eLORETA) in our analysis.

Discussion

This protocol will outline the rationale and research methodology for a clinical pilot trial of sLORETA ISF-NFB targeting key nodes within the core-ICNs in a female ID population with the primary aims being to assess its potential efficacy via transdiagnostic PROs and relevant neurophysiological measures.

Trial registration

Our study was prospectively registered with the Australia New Zealand Clinical Trials Registry (ANZCTR; Trial ID: ACTRN12619001428156). Registered on October 15, 2019.

Transition behavior of the seizure dynamics modulated by the astrocyte inositol triphosphate noise

Epilepsy is a neurological disorder with recurrent seizures, which convey complex dynamical characteristics including chaos and randomness. Until now, the underlying mechanism has not been fully elucidated, especially the bistable property beneath the epileptic random induction phenomena in certain conditions. Inspired by the recent finding that astrocyte GTPase-activating protein (G-protein)-coupled receptors could be involved in stochastic epileptic seizures, we proposed a neuron-astrocyte network model, incorporating the noise of the astrocytic second messenger, inositol triphosphate (IP3) that is modulated by G-protein-coupled receptor activation. Based on this model, we have statistically analyzed the transitions of epileptic seizures by performing repeatable simulation trials. Our simulation results show that the increase in the IP3 noise intensity induces depolarization-block epileptic seizures together with an increase in neuronal firing frequency, consistent with corresponding experiments. Meanwhile, the bistable states of the seizure dynamics were present under certain noise intensities, during which the neuronal firing pattern switches between regular sparse spiking and epileptic seizure states. This random presence of epileptic seizures is absent when the noise intensity continues to increase, accompanying with an increase in the epileptic depolarization block duration. The simulation results also shed light on the fact that calcium signals in astrocytes play significant roles in the pattern formations of the epileptic seizure. Our results provide a potential pathway for understanding the epileptic randomness in certain conditions.

Anomalous and heterogeneous characteristics of the BOLD hemodynamic response function in white matter

Detailed knowledge of the BOLD hemodynamic response function (HRF) is crucial for accurate analyses and interpretation of functional MRI data. Considerable efforts have been made to characterize the HRF in gray matter (GM), but much less attention has been paid to BOLD effects in white matter (WM). However, several recent reports have demonstrated reliable detection and analyses of WM BOLD signals both after stimulation and in a resting state. WM and GM differ in composition, energy requirements, and blood flow, so their neurovascular couplings also may well be different. We aimed to derive a comprehensive characterization of the HRF in WM across a population, including accurate measurements of its shape and its variation along and between WM pathways, using resting-state fMRI acquisitions. Our results show that the HRF is significantly different between WM and GM. Features of the HRF, such as a prominent initial dip, show strong relationships with features of the tissue microstructure derived from diffusion imaging, and these relationships differ between WM and GM, consistent with BOLD signal fluctuations reflecting different energy demands and neurovascular couplings in tissues of different composition and function. We also show that the HRF varies in shape significantly along WM pathways and is different between different WM pathways, suggesting the temporal evolution of BOLD signals after an event vary in different parts of the WM. These features of the HRF in WM are especially relevant for interpretation of the biophysical basis of BOLD effects in WM.

Inferring excitation-inhibition dynamics using a maximum entropy model unifying brain structure and function

Neural activity coordinated across different scales from neuronal circuits to large-scale brain networks gives rise to complex cognitive functions. Bridging the gap between micro- and macroscale processes, we present a novel framework based on the maximum entropy model to infer a hybrid resting-state structural connectome, representing functional interactions constrained by structural connectivity. We demonstrate that the structurally informed network outperforms the unconstrained model in simulating brain dynamics, wherein by constraining the inference model with the network structure we may improve the estimation of pairwise BOLD signal interactions. Further, we simulate brain network dynamics using Monte Carlo simulations with the new hybrid connectome to probe connectome-level differences in excitation-inhibition balance between apolipoprotein E (APOE)-ε4 carriers and noncarriers. Our results reveal sex differences among APOE-ε4 carriers in functional dynamics at criticality; specifically, female carriers appear to exhibit a lower tolerance to network disruptions resulting from increased excitatory interactions. In sum, the new multimodal network explored here enables analysis of brain dynamics through the integration of structure and function, providing insight into the complex interactions underlying neural activity such as the balance of excitation and inhibition.

Associations between polygenic risk scores and amplitude of low-frequency fluctuation of inferior frontal gyrus in schizophrenia

Schizophrenia (SCZ) is a serious and complex mental disorder with high heritability. Polygenic risk score (PRS) is a useful tool calculating the accumulating effects of multiple common genetic variants of schizophrenia. The amplitude of low-frequency fluctuation (ALFF) is an efficient index to reflect spontaneous, intrinsic neuronal activity. Aberrant ALFF of brain regions were reported in schizophrenia frequently, but the relationship between PRS and ALFF has not been studied. In the present study, we compared PRS and ALFF in 101 schizophrenia patients and 106 age-matched healthy controls to test their associations with schizophrenia. Then, the correlation of PRS with ALFF was measured to reveal the effect of polygenic risk on brain activity in schizophrenia. We found that schizophrenia patients showed significant differences in PRS and ALFF compared with controls. Twenty-six brain regions showed significant difference of ALFF between schizophrenia cases and controls, of which left inferior frontal gyrus, triangular part (IFGtriang.L) showed increased activity in schizophrenia. PRS-SCZ was positively correlated with ALFF in IFGtriang.L in 57 non-chronic patients. Genes involved in synaptic organization and transmission, especially in glutamatergic synapse, were highly enriched in PRS-SCZ genes, suggesting the dysfunction of synapses in schizophrenia. These results help to understand the molecular mechanism underlying schizophrenia and related brain dysfunction.

Impact of DC-Coupled Electrophysiological Recordings for Translational Neuroscience: Case Study of Tracking Neural Dynamics in Rodent Models of Seizures

We propose that to fully understand biological mechanisms underlying pathological brain activity with transitions (e.g., into and out of seizures), wide-bandwidth electrophysiological recordings are important. We demonstrate the importance of ultraslow potential shifts and infraslow oscillations for reliable tracking of synaptic physiology, within a neural mass model, from brain recordings that undergo pathological phase transitions. We use wide-bandwidth data (direct current (DC) to high-frequency activity), recorded using epidural and penetrating graphene micro-transistor arrays in a rodent model of acute seizures. Using this technological approach, we capture the dynamics of infraslow changes that contribute to seizure initiation (active pre-seizure DC shifts) and progression (passive DC shifts). By employing a continuous-discrete unscented Kalman filter, we track biological mechanisms from full-bandwidth data with and without active pre-seizure DC shifts during paroxysmal transitions. We then apply the same methodological approach for tracking the same parameters after application of high-pass-filtering >0.3Hz to both data sets. This approach reveals that ultraslow potential shifts play a fundamental role in the transition to seizure, and the use of high-pass-filtered data results in the loss of key information in regard to seizure onset and termination dynamics.

Grey-box modeling and hypothesis testing of functional near-infrared spectroscopy-based cerebrovascular reactivity to anodal high-definition tDCS in healthy humans

Transcranial direct current stimulation (tDCS) has been shown to evoke hemodynamics response; however, the mechanisms have not been investigated systematically using systems biology approaches. Our study presents a grey-box linear model that was developed from a physiologically detailed multi-compartmental neurovascular unit model consisting of the vascular smooth muscle, perivascular space, synaptic space, and astrocyte glial cell. Then, model linearization was performed on the physiologically detailed nonlinear model to find appropriate complexity (Akaike information criterion) to fit functional near-infrared spectroscopy (fNIRS) based measure of blood volume changes, called cerebrovascular reactivity (CVR), to high-definition (HD) tDCS. The grey-box linear model was applied on the fNIRS-based CVR during the first 150 seconds of anodal HD-tDCS in eleven healthy humans. The grey-box linear models for each of the four nested pathways starting from tDCS scalp current density that perturbed synaptic potassium released from active neurons for Pathway 1, astrocytic transmembrane current for Pathway 2, perivascular potassium concentration for Pathway 3, and voltage-gated ion channel current on the smooth muscle cell for Pathway 4 were fitted to the total hemoglobin concentration (tHb) changes from optodes in the vicinity of 4x1 HD-tDCS electrodes as well as on the contralateral sensorimotor cortex. We found that the tDCS perturbation Pathway 3 presented the least mean square error (MSE, median <2.5%) and the lowest Akaike information criterion (AIC, median -1.726) from the individual grey-box linear model fitting at the targeted-region. Then, minimal realization transfer function with reduced-order approximations of the grey-box model pathways was fitted to the ensemble average tHb time series. Again, Pathway 3 with nine poles and two zeros (all free parameters), provided the best Goodness of Fit of 0.0078 for Chi-Square difference test of nested pathways. Therefore, our study provided a systems biology approach to investigate the initial transient hemodynamic response to tDCS based on fNIRS tHb data. Future studies need to investigate the steady-state responses, including steady-state oscillations found to be driven by calcium dynamics, where transcranial alternating current stimulation may provide frequency-dependent physiological entrainment for system identification. We postulate that such a mechanistic understanding from system identification of the hemodynamics response to transcranial electrical stimulation can facilitate adequate delivery of the current density to the neurovascular tissue under simultaneous portable imaging in various cerebrovascular diseases.

NitroSynapsin ameliorates hypersynchronous neural network activity in Alzheimer hiPSC models

Beginning at early stages, human Alzheimer's disease (AD) brains manifest hyperexcitability, contributing to subsequent extensive synapse loss, which has been linked to cognitive dysfunction. No current therapy for AD is disease-modifying. Part of the problem with AD drug discovery is that transgenic mouse models have been poor predictors of potential human treatment. While it is undoubtedly important to test drugs in these animal models, additional evidence for drug efficacy in a human context might improve our chances of success. Accordingly, in order to test drugs in a human context, we have developed a platform of physiological assays using patch-clamp electrophysiology, calcium imaging, and multielectrode array (MEA) experiments on human (h)iPSC-derived 2D cortical neuronal cultures and 3D cerebral organoids. We compare hiPSCs bearing familial AD mutations vs. their wild-type (WT) isogenic controls in order to characterize the aberrant electrical activity in such a human context. Here, we show that these AD neuronal cultures and organoids manifest increased spontaneous action potentials, slow oscillatory events (~1 Hz), and hypersynchronous network activity. Importantly, the dual-allosteric NMDAR antagonist NitroSynapsin, but not the FDA-approved drug memantine, abrogated this hyperactivity. We propose a novel model of synaptic plasticity in which aberrant neural networks are rebalanced by NitroSynapsin. We propose that hiPSC models may be useful for screening drugs to treat hyperexcitability and related synaptic damage in AD.

Progressive Stabilization of Brain Network Dynamics during Childhood and Adolescence

Functional brain networks require dynamic reconfiguration to support flexible cognitive function. However, the developmental principles shaping brain network dynamics remain poorly understood. Here, we report the longitudinal development of large-scale brain network dynamics during childhood and adolescence, and its connection with gene expression profiles. Using a multilayer network model, we show the temporally varying modular architecture of child brain networks, with higher network switching primarily in the association cortex and lower switching in the primary regions. This topographical profile exhibits progressive maturation, which manifests as reduced modular dynamics, particularly in the transmodal (e.g., default-mode and frontoparietal) and sensorimotor regions. These developmental refinements mediate age-related enhancements of global network segregation and are linked with the expression profiles of genes associated with the enrichment of ion transport and nucleobase-containing compound transport. These results highlight a progressive stabilization of brain dynamics, which expand our understanding of the neural mechanisms that underlie cognitive development.

Global waves synchronize the brain's functional systems with fluctuating arousal

We propose and empirically support a parsimonious account of intrinsic, brain-wide spatiotemporal organization arising from traveling waves linked to arousal. We hypothesize that these waves are the predominant physiological process reflected in spontaneous functional magnetic resonance imaging (fMRI) signal fluctuations. The correlation structure ("functional connectivity") of these fluctuations recapitulates the large-scale functional organization of the brain. However, a unifying physiological account of this structure has so far been lacking. Here, using fMRI in humans, we show that ongoing arousal fluctuations are associated with global waves of activity that slowly propagate in parallel throughout the neocortex, thalamus, striatum, and cerebellum. We show that these waves can parsimoniously account for many features of spontaneous fMRI signal fluctuations, including topographically organized functional connectivity. Last, we demonstrate similar, cortex-wide propagation of neural activity measured with electrocorticography in macaques. These findings suggest that traveling waves spatiotemporally pattern brain-wide excitability in relation to arousal.

Functional Connectivity of the Cognitive Cerebellum

Anatomical tracing, human clinical data, and stimulation functional imaging have firmly established the major role of the (neo-)cerebellum in cognition and emotion. Telencephalization characterized by the great expansion of associative cortices, especially the prefrontal one, has been associated with parallel expansion of the neocerebellar cortex, especially the lobule VII, and by an increased number of interconnections between these two cortical structures. These anatomical modifications underlie the implication of the neocerebellum in cognitive control of complex motor and non-motor tasks. In humans, resting state functional connectivity has been used to determine a thorough anatomo-functional parcellation of the neocerebellum. This technique has identified central networks involving the neocerebellum and subserving its cognitive function. Neocerebellum participates in all intrinsic connected networks such as central executive, default mode, salience, dorsal and ventral attentional, and language-dedicated networks. The central executive network constitutes the main circuit represented within the neocerebellar cortex. Cerebellar zones devoted to these intrinsic networks appear multiple, interdigitated, and spatially ordered in three gradients. Such complex neocerebellar organization enables the neocerebellum to monitor and synchronize the main networks involved in cognition and emotion, likely by computing internal models.

Spontaneous Afferent Activity Carves Olfactory Circuits

Electrical activity has a key role in shaping neuronal circuits during development. In most sensory modalities, early in development, internally generated spontaneous activity sculpts the initial layout of neuronal wiring. With the maturation of the sense organs, the system relies more on sensory-evoked electrical activity. Stimuli-driven neuronal discharge is required for the transformation of immature circuits in the specific patterns of neuronal connectivity that subserve normal brain function. The olfactory system (OS) differs from this organizational plan. Despite the important role of odorant receptors (ORs) in shaping olfactory topography, odor-evoked activity does not have a prominent role in refining neuronal wiring. On the contrary, afferent spontaneous discharge is required to achieve and maintain the specific diagram of connectivity that defines the topography of the olfactory bulb (OB). Here, we provide an overview of the development of olfactory topography, with a focus on the role of afferent spontaneous discharge in the formation and maintenance of the specific synaptic contacts that result in the topographic organization of the OB.

Resting cerebral oxygen metabolism exhibits archetypal network features

Standard magnetic resonance imaging approaches offer high‐resolution but indirect measures of neural activity, limiting understanding of the physiological processes associated with imaging findings. Here, we used calibrated functional magnetic resonance imaging during the resting state to recover low‐frequency fluctuations of the cerebral metabolic rate of oxygen (CMRO₂). We tested whether functional connections derived from these fluctuations exhibited organization properties similar to those established by previous standard functional and anatomical connectivity studies. Seventeen participants underwent 20 min of resting imaging during dual‐echo, pseudocontinuous arterial spin labeling, and blood‐oxygen‐level dependent (BOLD) signal acquisition. Participants also underwent a 10 min normocapnic and hypercapnic procedure. Brain‐wide, CMRO₂ low‐frequency fluctuations were subjected to graph‐based and voxel‐wise functional connectivity analyses. Results demonstrated that connections derived from resting CMRO₂ fluctuations exhibited complex, small‐world topological properties (i.e., high integration and segregation, cost efficiency) consistent with those observed in previous studies using functional and anatomical connectivity approaches. Voxel‐wise CMRO₂ connectivity also exhibited spatial patterns consistent with four targeted resting‐state subnetworks: two association (i.e., frontoparietal and default mode) and two perceptual (i.e., auditory and occipital‐visual). These are the first findings to support the use of calibration‐derived CMRO₂ low‐frequency fluctuations for detecting brain‐wide organizational properties typical of healthy participants. We discuss interpretations, advantages, and challenges in using calibration‐derived oxygen metabolism signals for examining the intrinsic organization of the human brain.

Graphene active sensor arrays for long-term and wireless mapping of wide frequency band epicortical brain activity

Graphene active sensors have demonstrated promising capabilities for the detection of electrophysiological signals in the brain. Their functional properties, together with their flexibility as well as their expected stability and biocompatibility have raised them as a promising building block for large-scale sensing neural interfaces. However, in order to provide reliable tools for neuroscience and biomedical engineering applications, the maturity of this technology must be thoroughly studied. Here, we evaluate the performance of 64-channel graphene sensor arrays in terms of homogeneity, sensitivity and stability using a wireless, quasi-commercial headstage and demonstrate the biocompatibility of epicortical graphene chronic implants. Furthermore, to illustrate the potential of the technology to detect cortical signals from infra-slow to high-gamma frequency bands, we perform proof-of-concept long-term wireless recording in a freely behaving rodent. Our work demonstrates the maturity of the graphene-based technology, which represents a promising candidate for chronic, wide frequency band neural sensing interfaces.

Identifying and validating subtypes within major psychiatric disorders based on frontal-posterior functional imbalance via deep learning

Converging evidence increasingly implicates shared etiologic and pathophysiological characteristics among major psychiatric disorders (MPDs), such as schizophrenia (SZ), bipolar disorder (BD), and major depressive disorder (MDD). Examining the neurobiology of the psychotic-affective spectrum may greatly advance biological determination of psychiatric diagnosis, which is critical for the development of more effective treatments. In this study, ensemble clustering was developed to identify subtypes within a trans-diagnostic sample of MPDs. Whole brain amplitude of low-frequency fluctuations (ALFF) was used to extract the low-dimensional features for clustering in a total of 944 participants: 581 psychiatric patients (193 with SZ, 171 with BD, and 217 with MDD) and 363 healthy controls (HC). We identified two subtypes with differentiating patterns of functional imbalance between frontal and posterior brain regions, as compared to HC: (1) Archetypal MPDs (60% of MPDs) had increased frontal and decreased posterior ALFF, and decreased cortical thickness and white matter integrity in multiple brain regions that were associated with increased polygenic risk scores and enriched risk gene expression in brain tissues; (2) Atypical MPDs (40% of MPDs) had decreased frontal and increased posterior ALFF with no associated alterations in validity measures. Medicated Archetypal MPDs had lower symptom severity than their unmedicated counterparts; whereas medicated and unmedicated Atypical MPDs had no differences in symptom scores. Our findings suggest that frontal versus posterior functional imbalance as measured by ALFF is a novel putative trans-diagnostic biomarker differentiating subtypes of MPDs that could have implications for precision medicine.

On the origin of ultraslow spontaneous Na+ fluctuations in neurons of the neonatal forebrain

Spontaneous neuronal and astrocytic activity in the neonate forebrain is believed to drive the maturation of individual cells and their integration into complex brain-region-specific networks. The previously reported forms include bursts of electrical activity and oscillations in intracellular Ca²⁺ concentration. Here, we use ratiometric Na⁺ imaging to demonstrate spontaneous fluctuations in the intracellular Na⁺ concentration of CA1 pyramidal neurons and astrocytes in tissue slices obtained from the hippocampus of mice at postnatal days 2-4 (P2-4). These occur at very low frequency (∼2/h), can last minutes with amplitudes up to several millimolar, and mostly disappear after the first postnatal week. To further investigate their mechanisms, we model a network consisting of pyramidal neurons and interneurons. Experimentally observed Na⁺ fluctuations are mimicked when GABAergic inhibition in the simulated network is made depolarizing. Both our experiments and computational model show that blocking voltage-gated Na⁺ channels or GABAergic signaling significantly diminish the neuronal Na⁺ fluctuations. On the other hand, blocking a variety of other ion channels, receptors, or transporters including glutamatergic pathways does not have significant effects. Our model also shows that the amplitude and duration of Na⁺ fluctuations decrease as we increase the strength of glial K⁺ uptake. Furthermore, neurons with smaller somatic volumes exhibit fluctuations with higher frequency and amplitude. As opposed to this, larger extracellular to intracellular volume ratio observed in neonatal brain exerts a dampening effect. Finally, our model predicts that these periods of spontaneous Na⁺ influx leave neonatal neuronal networks more vulnerable to seizure-like states when compared with mature brain.NEW & NOTEWORTHY Spontaneous activity in the neonate forebrain plays a key role in cell maturation and brain development. We report spontaneous, ultraslow, asynchronous fluctuations in the intracellular Na⁺ concentration of neurons and astrocytes. We show that this activity is not correlated with the previously reported synchronous neuronal population bursting or Ca²⁺ oscillations, both of which occur at much faster timescales. Furthermore, extracellular K⁺ concentration remains nearly constant. The spontaneous Na⁺ fluctuations disappear after the first postnatal week.

Frequency-dependent circuits anchored in the dorsal and ventral left anterior insula

The hub role of the right anterior insula (AI) has been emphasized in cognitive neurosciences and been demonstrated to be frequency-dependently organized. However, the functional organization of left AI (LAI) has not been systematically investigated. Here we used 100 unrelated datasets from the Human Connectome Project to study the frequency-dependent organization of LAI along slow 6 to slow 1 bands. The broadband functional connectivity of LAI was similar to previous findings. In slow 6-slow 3 bands, both dorsal and ventral seeds in LAI were correlated to the salience network (SN) and language network (LN) and anti-correlated to the default mode network (DMN). However, these seeds were only correlated to the LAI in slow 2-slow 1 bands. These findings indicate that broadband and narrow band functional connections reflect different functional organizations of the LAI. Furthermore, the dorsal seed had a stronger connection with the LN and anti-correlation with DMN while the ventral seed had a stronger connection within the SN in slow 6-slow 3 bands. In slow 2-slow 1 bands, both seeds had stronger connections with themselves. These observations indicate distinctive functional organizations for the two parts of LAI. Significant frequency effect and frequency by seed interaction were also found, suggesting different frequency characteristics of these two seeds. The functional integration and functional segregation of LDAI and LVAI were further supported by their cognitive associations. The frequency- and seed-dependent functional organizations of LAI may enlighten future clinical and cognitive investigations.

Interstitial ions: A key regulator of state-dependent neural activity?

Throughout the nervous system, ion gradients drive fundamental processes. Yet, the roles of interstitial ions in brain functioning is largely forgotten. Emerging literature is now revitalizing this area of neuroscience by showing that interstitial cations (K+, Ca2+ and Mg2+) are not static quantities but change dynamically across states such as sleep and locomotion. In turn, these state-dependent changes are capable of sculpting neuronal activity; for example, changing the local interstitial ion composition in the cortex is sufficient for modulating the prevalence of slow-frequency neuronal oscillations, or potentiating the gain of visually evoked responses. Disturbances in interstitial ionic homeostasis may also play a central role in the pathogenesis of central nervous system diseases. For example, impairments in K+ buffering occur in a number of neurodegenerative diseases, and abnormalities in neuronal activity in disease models disappear when interstitial K+ is normalized. Here we provide an overview of the roles of interstitial ions in physiology and pathology. We propose the brain uses interstitial ion signaling as a global mechanism to coordinate its complex activity patterns, and ion homeostasis failure contributes to central nervous system diseases affecting cognitive functions and behavior.

Distinct Structure of Cortical Population Activity on Fast and Infraslow Timescales

Cortical activity is organized across multiple spatial and temporal scales. Most research on the dynamics of neuronal spiking is concerned with timescales of 1 ms–1 s, and little is known about spiking dynamics on timescales of tens of seconds and minutes. Here, we used frequency domain analyses to study the structure of individual neurons’ spiking activity and its coupling to local population rate and to arousal level across 0.01–100 Hz frequency range. In mouse medial prefrontal cortex, the spiking dynamics of individual neurons could be quantitatively captured by a combination of interspike interval and firing rate power spectrum distributions. The relative strength of coherence with local population often differed across timescales: a neuron strongly coupled to population rate on fast timescales could be weakly coupled on slow timescales, and vice versa. On slow but not fast timescales, a substantial proportion of neurons showed firing anticorrelated with the population. Infraslow firing rate changes were largely determined by arousal rather than by local factors, which could explain the timescale dependence of individual neurons’ population coupling strength. These observations demonstrate how neurons simultaneously partake in fast local dynamics, and slow brain-wide dynamics, extending our understanding of infraslow cortical activity beyond the mesoscale resolution of fMRI.

Brain states govern the spatio-temporal dynamics of resting-state functional connectivity

Previously, using simultaneous resting-state functional magnetic resonance imaging (fMRI) and photometry-based neuronal calcium recordings in the anesthetized rat, we identified blood oxygenation level-dependent (BOLD) responses directly related to slow calcium waves, revealing a cortex-wide and spatially organized correlate of locally recorded neuronal activity (Schwalm et al., 2017). Here, using the same techniques, we investigate two distinct cortical activity states: persistent activity, in which compartmentalized network dynamics were observed; and slow wave activity, dominated by a cortex-wide BOLD component, suggesting a strong functional coupling of inter-cortical activity. During slow wave activity, we find a correlation between the occurring slow wave events and the strength of functional connectivity between different cortical areas. These findings suggest that down-up transitions of neuronal excitability can drive cortex-wide functional connectivity. This study provides further evidence that changes in functional connectivity are dependent on the brain's current state, directly linked to the generation of slow waves.

Decoding Task-Specific Cognitive States with Slow, Directed Functional Networks in the Human Brain

Flexible functional interactions among brain regions mediate critical cognitive functions. Such interactions can be measured using functional magnetic resonance imaging (fMRI) data either with instantaneous (zero-lag) or lag-based (time-lagged) functional connectivity. Because the fMRI hemodynamic response is slow, and is sampled at a timescale (seconds) several orders of magnitude slower than the underlying neural dynamics (milliseconds), simulation studies have shown that lag-based fMRI functional connectivity, measured with approaches like Granger–Geweke causality (GC), provides spurious and unreliable estimates of underlying neural interactions. Experimental verification of this claim is challenging because neural ground truth connectivity is often unavailable concurrently with fMRI recordings. Here we demonstrate that, despite these widely held caveats, GC networks estimated from fMRI recordings contain useful information for classifying task-specific cognitive states. We estimated instantaneous and lag-based GC functional connectivity networks using fMRI data from 1000 participants (Human Connectome Project database). A linear classifier, trained on either instantaneous or lag-based GC, reliably discriminated among seven different task and resting brain states, with >80% cross-validation accuracy. With network simulations, we demonstrate that instantaneous and lag-based GC exploited interactions at fast and slow timescales, respectively, to achieve robust classification. With human fMRI data, instantaneous and lag-based GC identified complementary, task–core networks. Finally, variations in GC connectivity explained inter-individual variations in a variety of cognitive scores. Our findings show that instantaneous and lag-based methods reveal complementary aspects of functional connectivity in the brain, and suggest that slow, directed functional interactions, estimated with fMRI, may provide useful markers of behaviorally relevant cognitive states.

Dissociated neuronal phase- and amplitude-coupling patterns in the human brain

Coupling of neuronal oscillations may reflect and facilitate the communication between neuronal populations. Two primary neuronal coupling modes have been described: phase-coupling and amplitude-coupling. Theoretically, both coupling modes are independent, but so far, their neuronal relationship remains unclear. Here, we combined MEG, source-reconstruction and simulations to systematically compare cortical amplitude-coupling and phase-coupling patterns in the human brain. Importantly, we took into account a critical bias of amplitude-coupling measures due to phase-coupling. We found differences between both coupling modes across a broad frequency range and most of the cortex. Furthermore, by combining empirical measurements and simulations we ruled out that these results were caused by methodological biases, but instead reflected genuine neuronal amplitude coupling. Our results show that cortical phase- and amplitude-coupling patterns are non-redundant, which may reflect at least partly distinct neuronal mechanisms. Furthermore, our findings highlight and clarify the compound nature of amplitude coupling measures.

Diet modulates brain network stability, a biomarker for brain aging, in young adults

To better understand how diet influences brain aging, we focus here on the presymptomatic period during which prevention may be most effective. Large-scale life span neuroimaging datasets show functional communication between brain regions destabilizes with age, typically starting in the late 40s, and that destabilization correlates with poorer cognition and accelerates with insulin resistance. Targeted experiments show that this biomarker for brain aging is reliably modulated with consumption of different fuel sources: Glucose decreases, and ketones increase the stability of brain networks. This effect replicated across both changes to total diet as well as fuel-specific calorie-matched bolus, producing changes in overall brain activity that suggest that network “switching” may reflect the brain’s adaptive response to conserve energy under resource constraint.

Epidemiological studies suggest that insulin resistance accelerates progression of age-based cognitive impairment, which neuroimaging has linked to brain glucose hypometabolism. As cellular inputs, ketones increase Gibbs free energy change for ATP by 27% compared to glucose. Here we test whether dietary changes are capable of modulating sustained functional communication between brain regions (network stability) by changing their predominant dietary fuel from glucose to ketones. We first established network stability as a biomarker for brain aging using two large-scale (n = 292, ages 20 to 85 y; n = 636, ages 18 to 88 y) 3 T functional MRI (fMRI) datasets. To determine whether diet can influence brain network stability, we additionally scanned 42 adults, age < 50 y, using ultrahigh-field (7 T) ultrafast (802 ms) fMRI optimized for single-participant-level detection sensitivity. One cohort was scanned under standard diet, overnight fasting, and ketogenic diet conditions. To isolate the impact of fuel type, an independent overnight fasted cohort was scanned before and after administration of a calorie-matched glucose and exogenous ketone ester (d-β-hydroxybutyrate) bolus. Across the life span, brain network destabilization correlated with decreased brain activity and cognitive acuity. Effects emerged at 47 y, with the most rapid degeneration occurring at 60 y. Networks were destabilized by glucose and stabilized by ketones, irrespective of whether ketosis was achieved with a ketogenic diet or exogenous ketone ester. Together, our results suggest that brain network destabilization may reflect early signs of hypometabolism, associated with dementia. Dietary interventions resulting in ketone utilization increase available energy and thus may show potential in protecting the aging brain.

Cortex-wide Changes in Extracellular Potassium Ions Parallel Brain State Transitions in Awake Behaving Mice

Brain state fluctuations modulate sensory processing, but the factors governing state-dependent neural activity remain unclear. Here, we tracked the dynamics of cortical extracellular K+ concentrations ([K+]o) during awake state transitions and manipulated [K+]o in slices, during visual processing, and during skilled motor execution. When mice transitioned from quiescence to locomotion, [K+]o increased by 0.6-1.0 mM in all cortical areas analyzed, and this preceded locomotion by 1 s. Emulating the state-dependent [K+]o increase in cortical slices caused neuronal depolarization and enhanced input-output transformation. In vivo, locomotion increased the gain of visually evoked responses in layer 2/3 of visual cortex; this effect was recreated by imposing a [K+]o increase. Elevating [K+]o in the motor cortex increased movement-induced neuronal spiking in layer 5 and improved motor performance. Thus, [K+]o increases in a cortex-wide state-dependent manner, and this [K+]o increase affects both sensory and motor processing through the dynamic modulation of neural activity.

Slow cortical potential signal classification using concave-convex feature

BACKGROUND

The classification of the slow cortical potential (SCP) signals plays a key role in a variety of research areas, including disease diagnostics, human-machine interaction, and education. The widely used classification methods, which combine multiple kinds of features, can be unsuitable in practical applications due to their low robustness to scenario changes.

NEW METHOD

A flexible concave-convex (C-C) feature is reported. The C-C feature is extracted by two steps: (1) the low-frequency node coefficients of the SCP signals are first extracted using wavelet packet decomposition; (2) then the underlying trend of the low-frequency node coefficients is estimated using third-order polynomial fitting, and the feature is constructed using the minimum and maximum second derivative values of the trend curve as |y_min| - y_max where y is the second derivative value.

RESULTS

Experimental results on real datasets reveal that our method with the single C-C feature exhibits high average classification accuracies (92.5% and 84.9% on the BCI competition II dataset Ia and the TJU dataset). The accuracy can be further improved (94.5% and 85.9%) by adding the commonly used mean voltage feature and using the naive Bayesian classifier, indicating the flexibility and scalability of the proposed method.

COMPARISON WITH EXISTING METHODS

The proposed C-C feature based method outperforms state-of-the-art (SOTA) multi-feature classification method from the perspective of classification accuracy.

CONCLUSIONS

The effectiveness of the C-C feature for SCP classification is validated. The proposed feature will represent a useful contribution to the SCP classification, balancing the strengths of traditional features and the proposed one.

Ionic and synaptic mechanisms of seizure generation and epileptogenesis

The biophysical mechanisms underlying epileptogenesis and the generation of seizures remain to be better understood. Among many factors triggering epileptogenesis are traumatic brain injury breaking normal synaptic homeostasis and genetic mutations disrupting ionic concentration homeostasis. Impairments in these mechanisms, as seen in various brain diseases, may push the brain network to a pathological state characterized by increased susceptibility to unprovoked seizures. Here, we review recent computational studies exploring the roles of ionic concentration dynamics in the generation, maintenance, and termination of seizures. We further discuss how ionic and synaptic homeostatic mechanisms may give rise to conditions which prime brain networks to exhibit recurrent spontaneous seizures and epilepsy.

Biophysical Modeling Suggests Optimal Drug Combinations for Improving the Efficacy of GABA Agonists after Traumatic Brain Injuries

Traumatic brain injuries (TBI) lead to dramatic changes in the surviving brain tissue. Altered ion concentrations, coupled with changes in the expression of membrane-spanning proteins, create a post-TBI brain state that can lead to further neuronal loss caused by secondary excitotoxicity. Several GABA receptor agonists have been tested in the search for neuroprotection immediately after an injury, with paradoxical results. These drugs not only fail to offer neuroprotection, but can also slow down functional recovery after TBI. Here, using computational modeling, we provide a biophysical hypothesis to explain these observations. We show that the accumulation of intracellular chloride ions caused by a transient upregulation of Na+-K+-2Cl- (NKCC1) co-transporters as observed following TBI, causes GABA receptor agonists to lead to excitation and depolarization block, rather than the expected hyperpolarization. The likelihood of prolonged, excitotoxic depolarization block is further exacerbated by the extremely high levels of extracellular potassium seen after TBI. Our modeling results predict that the neuroprotective efficacy of GABA receptor agonists can be substantially enhanced when they are combined with NKCC1 co-transporter inhibitors. This suggests a rational, biophysically principled method for identifying drug combinations for neuroprotection after TBI.

Simulating human sleep spindle MEG and EEG from ion channel and circuit level dynamics

BACKGROUND

Although they form a unitary phenomenon, the relationship between extracranial M/EEG and transmembrane ion flows is understood only as a general principle rather than as a well-articulated and quantified causal chain.

METHOD

We present an integrated multiscale model, consisting of a neural simulation of thalamus and cortex during stage N2 sleep and a biophysical model projecting cortical current densities to M/EEG fields. Sleep spindles were generated through the interactions of local and distant network connections and intrinsic currents within thalamocortical circuits. 32,652 cortical neurons were mapped onto the cortical surface reconstructed from subjects' MRI, interconnected based on geodesic distances, and scaled-up to current dipole densities based on laminar recordings in humans. MRIs were used to generate a quasi-static electromagnetic model enabling simulated cortical activity to be projected to the M/EEG sensors.

RESULTS

The simulated M/EEG spindles were similar in amplitude and topography to empirical examples in the same subjects. Simulated spindles with more core-dominant activity were more MEG weighted.

COMPARISON WITH EXISTING METHODS

Previous models lacked either spindle-generating thalamic neural dynamics or whole head biophysical modeling; the framework presented here is the first to simultaneously capture these disparate scales.

CONCLUSIONS

This multiscale model provides a platform for the principled quantitative integration of existing information relevant to the generation of sleep spindles, and allows the implications of future findings to be explored. It provides a proof of principle for a methodological framework allowing large-scale integrative brain oscillations to be understood in terms of their underlying channels and synapses.