Tracking brain arousal fluctuations with fMRI

Integration of Simultaneous Resting-State Electroencephalography, Functional Magnetic Resonance Imaging, and Eye-Tracker Methods to Determine and Verify Electroencephalography Vigilance Measure

Background/Introduction: Concurrent electroencephalography and resting-state functional magnetic resonance imaging (rsfMRI) have been widely used for studying the (presumably) awake and alert human brain with high temporal/spatial resolution. Although rsfMRI scans are typically collected while individuals are instructed to focus their eyes on a fixated cross, objective and verified experimental measures to quantify degree of vigilance are not readily available. Electroencephalography (EEG) is the modality extensively used for estimating vigilance, especially during eyes-closed resting state. However, pupil size measured using an eye-tracker device could provide an indirect index of vigilance. Methods: Three 12-min resting scans (eyes open, fixating on the cross) were collected from 10 healthy control participants. We simultaneously collected EEG, fMRI, physiological, and eye-tracker data and investigated the correlation between EEG features, pupil size, and heart rate. Furthermore, we used pupil size and EEG features as regressors to find their correlations with blood-oxygen-level-dependent fMRI measures. Results: EEG frontal and occipital beta power (FOBP) correlates with pupil size changes, an indirect index for locus coeruleus activity implicated in vigilance regulation (r = 0.306, p < 0.001). Moreover, FOBP also correlated with heart rate (r = 0.255, p < 0.001), as well as several brain regions in the anticorrelated network, including the bilateral insula and inferior parietal lobule. Discussion: In this study, we investigated whether simultaneous EEG-fMRI combined with eye-tracker measurements can be used to determine EEG signal feature associated with vigilance measures during eyes-open rsfMRI. Our results support the conclusion that FOBP is an objective measure of vigilance in healthy human subjects. Impact statement We revealed an association between electroencephalography frontal and occipital beta power (FOBP) and pupil size changes during an eyes-open resting state, which supports the conclusion that FOBP could serve as an objective measure of vigilance in healthy human subjects. The results were validated by using simultaneously recorded heart rate and functional magnetic resonance imaging (fMRI). Interestingly, independently verified heart rate changes can also provide an easy-to-determine measure of vigilance during resting-state fMRI. These findings have important implications for an analysis and interpretation of dynamic resting-state fMRI connectivity studies in health and disease.

Neurovascular coupling and bilateral connectivity during NREM and REM sleep

To understand how arousal state impacts cerebral hemodynamics and neurovascular coupling, we monitored neural activity, behavior, and hemodynamic signals in un-anesthetized, head-fixed mice. Mice frequently fell asleep during imaging, and these sleep events were interspersed with periods of wake. During both NREM and REM sleep, mice showed large increases in cerebral blood volume ([HbT]) and arteriole diameter relative to the awake state, two to five times larger than those evoked by sensory stimulation. During NREM, the amplitude of bilateral low-frequency oscillations in [HbT] increased markedly, and coherency between neural activity and hemodynamic signals was higher than the awake resting and REM states. Bilateral correlations in neural activity and [HbT] were highest during NREM, and lowest in the awake state. Hemodynamic signals in the cortex are strongly modulated by arousal state, and changes during sleep are substantially larger than sensory-evoked responses.

Predicting the fMRI Signal Fluctuation with Recurrent Neural Networks Trained on Vascular Network Dynamics

Resting-state functional MRI (rs-fMRI) studies have revealed specific low-frequency hemodynamic signal fluctuations (<0.1 Hz) in the brain, which could be related to neuronal oscillations through the neurovascular coupling mechanism. Given the vascular origin of the fMRI signal, it remains challenging to separate the neural correlates of global rs-fMRI signal fluctuations from other confounding sources. However, the slow-oscillation detected from individual vessels by single-vessel fMRI presents strong correlation to neural oscillations. Here, we use recurrent neural networks (RNNs) to predict the future temporal evolution of the rs-fMRI slow oscillation from both rodent and human brains. The RNNs trained with vessel-specific rs-fMRI signals encode the unique brain oscillatory dynamic feature, presenting more effective prediction than the conventional autoregressive model. This RNN-based predictive modeling of rs-fMRI datasets from the Human Connectome Project (HCP) reveals brain state-specific characteristics, demonstrating an inverse relationship between the global rs-fMRI signal fluctuation with the internal default-mode network (DMN) correlation. The RNN prediction method presents a unique data-driven encoding scheme to specify potential brain state differences based on the global fMRI signal fluctuation, but not solely dependent on the global variance.

Altered Global Brain Signal during Physiologic, Pharmacologic, and Pathologic States of Unconsciousness in Humans and Rats

BACKGROUND

Consciousness is supported by integrated brain activity across widespread functionally segregated networks. The functional magnetic resonance imaging-derived global brain signal is a candidate marker for a conscious state, and thus the authors hypothesized that unconsciousness would be accompanied by a loss of global temporal coordination, with specific patterns of decoupling between local regions and global activity differentiating among various unconscious states.

METHODS

Functional magnetic resonance imaging global signals were studied in physiologic, pharmacologic, and pathologic states of unconsciousness in human natural sleep (n = 9), propofol anesthesia (humans, n = 14; male rats, n = 12), and neuropathological patients (n = 21). The global signal amplitude as well as the correlation between global signal and signals of local voxels were quantified. The former reflects the net strength of global temporal coordination, and the latter yields global signal topography.

RESULTS

A profound reduction of global signal amplitude was seen consistently across the various unconscious states: wakefulness (median [1st, 3rd quartile], 0.46 [0.21, 0.50]) versus non-rapid eye movement stage 3 of sleep (0.30 [0.24, 0.32]; P = 0.035), wakefulness (0.36 [0.31, 0.42]) versus general anesthesia (0.25 [0.21, 0.28]; P = 0.001), healthy controls (0.30 [0.27, 0.37]) versus unresponsive wakefulness syndrome (0.22 [0.15, 0.24]; P < 0.001), and low dose (0.07 [0.06, 0.08]) versus high dose of propofol (0.04 [0.03, 0.05]; P = 0.028) in rats. Furthermore, non-rapid eye movement stage 3 of sleep was characterized by a decoupling of sensory and attention networks from the global network. General anesthesia and unresponsive wakefulness syndrome were characterized by a dissociation of the majority of functional networks from the global network. This decoupling, however, was dominated by distinct neuroanatomic foci (e.g., precuneus and anterior cingulate cortices).

CONCLUSIONS

The global temporal coordination of various modules across the brain may distinguish the coarse-grained state of consciousness versus unconsciousness, while the relationship between the global and local signals may define the particular qualities of a particular unconscious state.

Metastable neural dynamics underlies cognitive performance across multiple behavioural paradigms

Despite resting state networks being associated with a variety of cognitive abilities, it remains unclear how these local areas act in concert to express particular cognitive operations. Theoretical and empirical accounts indicate that large-scale resting state networks reconcile dual tendencies towards integration and segregation by operating in a metastable regime of their coordination dynamics. Metastability may confer important behavioural qualities by binding distributed local areas into large-scale neurocognitive networks. We tested this hypothesis by analysing fMRI data in a large cohort of healthy individuals (N = 566) and comparing the metastability of the brain's large-scale resting network architecture at rest and during the performance of several tasks. Metastability was estimated using a well-defined collective variable capturing the level of 'phase-locking' between large-scale networks over time. Task-based reasoning was principally characterised by high metastability in cognitive control networks and low metastability in sensory processing areas. Although metastability between resting state networks increased during task performance, cognitive ability was more closely linked to spontaneous activity. High metastability in the intrinsic connectivity of cognitive control networks was linked to novel problem solving or fluid intelligence, but was less important in tasks relying on previous experience or crystallised intelligence. Crucially, subjects with resting architectures similar or 'pre-configured' to a task-general arrangement demonstrated superior cognitive performance. Taken together, our findings support a key linkage between the spontaneous metastability of large-scale networks in the cerebral cortex and cognition.

Analysis of stimulus-induced brain dynamics during naturalistic paradigms

Naturalistic stimuli offer promising avenues for investigating brain function across the rich, realistic spectrum of human experiences. Functional magnetic resonance imaging (fMRI) studies of brain activity during naturalistic paradigms have provided new information about dynamic neural processing in ecologically valid contexts. Yet, the complex, uncontrolled nature of such stimuli -- and the resulting mixture of neuronal and physiological responses embedded within the fMRI signals -- present challenges with respect to data analysis and interpretation. In this brief commentary, we discuss methods and open challenges in naturalistic fMRI investigations, with a focus on extracting and interpreting stimulus-induced fMRI signals.

Rest-task modulation of fMRI-derived global signal topography is mediated by transient coactivation patterns

Recent resting-state functional MRI (fMRI) studies have revealed that the global signal (GS) exhibits a nonuniform spatial distribution across the gray matter. Whether this topography is informative remains largely unknown. We therefore tested rest-task modulation of GS topography by analyzing static GS correlation and dynamic coactivation patterns in a large sample of an fMRI dataset (n = 837) from the Human Connectome Project. The GS topography in the resting state and in seven different tasks was first measured by correlating the GS with the local time series (GSCORR). In the resting state, high GSCORR was observed mainly in the primary sensory and motor regions, whereas low GSCORR was seen in the association brain areas. This pattern changed during the seven tasks, with mainly decreased GSCORR in sensorimotor cortex. Importantly, this rest-task modulation of GSCORR could be traced to transient coactivation patterns at the peak period of GS (GS-peak). By comparing the topography of GSCORR and respiration effects, we observed that the topography of respiration mimicked the topography of GS in the resting state, whereas both differed during the task states; because of such partial dissociation, we assume that GSCORR could not be equated with a respiration effect. Finally, rest-task modulation of GS topography could not be exclusively explained by other sources of physiological noise. Together, we here demonstrate the informative nature of GS topography by showing its rest-task modulation, the underlying dynamic coactivation patterns, and its partial dissociation from respiration effects during task states.

Recent resting-state fMRI studies have shown that the global signal exhibits a nonuniform spatial distribution across gray matter, but is this informative? This neuroimaging study reveals novel insights into the informative nature of global signal by rest-task modulation of the global signal topography.

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.

Connectivity dynamics from wakefulness to sleep

Interest in time-resolved connectivity in fMRI has grown rapidly in recent years. The most widely used technique for studying connectivity changes over time utilizes a sliding windows approach. There has been some debate about the utility of shorter versus longer windows, the use of fixed versus adaptive windows, as well as whether observed resting state dynamics during wakefulness may be predominantly due to changes in sleep state and subject head motion. In this work we use an independent component analysis (ICA)-based pipeline applied to concurrent EEG/fMRI data collected during wakefulness and various sleep stages and show: 1) connectivity states obtained from clustering sliding windowed correlations of resting state functional network time courses well classify the sleep states obtained from EEG data, 2) using shorter sliding windows instead of longer non-overlapping windows improves the ability to capture transition dynamics even at windows as short as 30 ​s, 3) motion appears to be mostly associated with one of the states rather than spread across all of them 4) a fixed tapered sliding window approach outperforms an adaptive dynamic conditional correlation approach, and 5) consistent with prior EEG/fMRI work, we identify evidence of multiple states within the wakeful condition which are able to be classified with high accuracy. Classification of wakeful only states suggest the presence of time-varying changes in connectivity in fMRI data beyond sleep state or motion. Results also inform about advantageous technical choices, and the identification of different clusters within wakefulness that are separable suggest further studies in this direction.

Dynamic representations in networked neural systems

A group of neurons can generate patterns of activity that represent information about stimuli; subsequently, the group can transform and transmit activity patterns across synapses to spatially distributed areas. Recent studies in neuroscience have begun to independently address the two components of information processing: the representation of stimuli in neural activity and the transmission of information in networks that model neural interactions. Yet only recently are studies seeking to link these two types of approaches. Here we briefly review the two separate bodies of literature; we then review the recent strides made to address this gap. We continue with a discussion of how patterns of activity evolve from one representation to another, forming dynamic representations that unfold on the underlying network. Our goal is to offer a holistic framework for understanding and describing neural information representation and transmission while revealing exciting frontiers for future research.

Temporal fluctuations in the brain's modular architecture during movie-watching

Brain networks are flexible and reconfigure over time to support ongoing cognitive processes. However, tracking statistically meaningful reconfigurations across time has proven difficult. This has to do largely with issues related to sampling variability, making instantaneous estimation of network organization difficult, along with increased reliance on task-free (cognitively unconstrained) experimental paradigms, limiting the ability to interpret the origin of changes in network structure over time. Here, we address these challenges using time-varying network analysis in conjunction with a naturalistic viewing paradigm. Specifically, we developed a measure of inter-subject network similarity and used this measure as a coincidence filter to identify synchronous fluctuations in network organization across individuals. Applied to movie-watching data, we found that periods of high inter-subject similarity coincided with reductions in network modularity and increased connectivity between cognitive systems. In contrast, low inter-subject similarity was associated with increased system segregation and more rest-like architectures. We then used a data-driven approach to uncover clusters of functional connections that follow similar trajectories over time and are more strongly correlated during movie-watching than at rest. Finally, we show that synchronous fluctuations in network architecture over time can be linked to a subset of features in the movie. Our findings link dynamic fluctuations in network integration and segregation to patterns of inter-subject similarity, and suggest that moment-to-moment fluctuations in functional connectivity reflect shared cognitive processing across individuals.

Transient Arousal Modulations Contribute to Resting-State Functional Connectivity Changes Associated with Head Motion Parameters

Correlations of resting-state functional magnetic resonance imaging (rsfMRI) signals are being widely used for assessing the functional brain connectivity in health and disease. However, an association was recently observed between rsfMRI connectivity modulations and the head motion parameters and regarded as a causal relationship, which has raised serious concerns about the validity of many rsfMRI findings. Here, we studied the origin of this rsfMRI-motion association and its relationship to arousal modulations. By using a template-matching method to locate arousal-related fMRI changes, we showed that the effects of high motion time points on rsfMRI connectivity are largely due to their significant overlap with arousal-affected time points. The finding suggests that the association between rsfMRI connectivity and the head motion parameters arises from their comodulations at transient arousal modulations, and this information is critical not only for proper interpretation of motion-associated rsfMRI connectivity changes, but also for controlling the potential confounding effects of arousal modulation on rsfMRI metrics.

Vigilance Effects in Resting-State fMRI

Measures of resting-state functional magnetic resonance imaging (rsfMRI) activity have been shown to be sensitive to cognitive function and disease state. However, there is growing evidence that variations in vigilance can lead to pronounced and spatially widespread differences in resting-state brain activity. Unless properly accounted for, differences in vigilance can give rise to changes in resting-state activity that can be misinterpreted as primary cognitive or disease-related effects. In this paper, we examine in detail the link between vigilance and rsfMRI measures, such as signal variance and functional connectivity. We consider how state changes due to factors such as caffeine and sleep deprivation affect both vigilance and rsfMRI measures and review emerging approaches and methodological challenges for the estimation and interpretation of vigilance effects.

Concurrent Imaging of Markers of Current Flow and Neurophysiological Changes During tDCS

Despite being a popular neuromodulation technique, clinical translation of transcranial direct current stimulation (tDCS) is hampered by variable responses observed within treatment cohorts. Addressing this challenge has been difficult due to the lack of an effective means of mapping the neuromodulatory electromagnetic fields together with the brain's response. In this study, we present a novel imaging technique that provides the capability of concurrently mapping markers of tDCS currents, as well as the brain's response to tDCS. A dual-echo echo-planar imaging (DE-EPI) sequence is used, wherein the phase of the acquired MRI-signal encodes the tDCS current induced magnetic field, while the magnitude encodes the blood oxygenation level dependent (BOLD) contrast. The proposed technique was first validated in a custom designed phantom. Subsequent test-retest experiments in human participants showed that tDCS-induced magnetic fields can be detected reliably in vivo. The concurrently acquired BOLD data revealed large-scale networks characteristic of a brain in resting-state as well as a 'cathodal' and an 'anodal' resting-state component under each electrode. Moreover, 'cathodal's BOLD-signal was observed to significantly decrease with the applied current at the group level in all datasets. With its ability to image markers of electromagnetic cause as well as neurophysiological changes, the proposed technique may provide an effective means to visualize neural engagement in tDCS at the group level. Our technique also contributes to addressing confounding factors in applying BOLD fMRI concurrently with tDCS.

Indexing brain state-dependent pupil dynamics with simultaneous fMRI and optical fiber calcium recording

Pupillometry, a noninvasive measure of arousal, complements human functional MRI (fMRI) to detect periods of variable cognitive processing and identify networks that relate to particular attentional states. Even under anesthesia, pupil dynamics correlate with brain-state fluctuations, and extended dilations mark the transition to more arousable states. However, cross-scale neuronal activation patterns are seldom linked to brain state-dependent pupil dynamics. Here, we complemented resting-state fMRI in rats with cortical calcium recording (GCaMP-mediated) and pupillometry to tackle the linkage between brain-state changes and neural dynamics across different scales. This multimodal platform allowed us to identify a global brain network that covaried with pupil size, which served to generate an index indicative of the brain-state fluctuation during anesthesia. Besides, a specific correlation pattern was detected in the brainstem, at a location consistent with noradrenergic cell group 5 (A5), which appeared to be dependent on the coupling between different frequencies of cortical activity, possibly further indicating particular brain-state dynamics. The multimodal fMRI combining concurrent calcium recordings and pupillometry enables tracking brain state-dependent pupil dynamics and identifying unique cross-scale neuronal dynamic patterns under anesthesia.

Temporal circuit of macroscale dynamic brain activity supports human consciousness

The ongoing stream of human consciousness relies on two distinct cortical systems, the default mode network and the dorsal attention network, which alternate their activity in an anticorrelated manner. We examined how the two systems are regulated in the conscious brain and how they are disrupted when consciousness is diminished. We provide evidence for a "temporal circuit" characterized by a set of trajectories along which dynamic brain activity occurs. We demonstrate that the transitions between default mode and dorsal attention networks are embedded in this temporal circuit, in which a balanced reciprocal accessibility of brain states is characteristic of consciousness. Conversely, isolation of the default mode and dorsal attention networks from the temporal circuit is associated with unresponsiveness of diverse etiologies. These findings advance the foundational understanding of the functional role of anticorrelated systems in consciousness.

Identifying and removing widespread signal deflections from fMRI data: Rethinking the global signal regression problem

One of the most controversial procedures in the analysis of resting-state functional magnetic resonance imaging (rsfMRI) data is global signal regression (GSR): the removal, via linear regression, of the mean signal averaged over the entire brain. On one hand, the global mean signal contains variance associated with respiratory, scanner-, and motion-related artifacts, and its removal via GSR improves various quality-control metrics, enhances the anatomical specificity of functional-connectivity patterns, and can increase the behavioral variance explained by such patterns. On the other hand, GSR alters the distribution of regional signal correlations in the brain, can induce artifactual anticorrelations, may remove real neural signal, and can distort case-control comparisons of functional-connectivity measures. Global signal fluctuations can be identified visually from a matrix of colour-coded signal intensities, called a carpet plot, in which rows represent voxels and columns represent time. Prior to GSR, large, periodic bands of coherent signal changes that affect most of the brain are often apparent; after GSR, these apparently global changes are greatly diminished. Here, using three independent datasets, we show that reordering carpet plots to emphasize cluster structure in the data reveals a greater diversity of spatially widespread signal deflections (WSDs) than previously thought. Their precise form varies across time and participants, and GSR is only effective in removing specific kinds of WSDs. We present an alternative, iterative correction method called Diffuse Cluster Estimation and Regression (DiCER), that identifies representative signals associated with large clusters of coherent voxels. DiCER is more effective than GSR at removing diverse WSDs as visualized in carpet plots, reduces correlations between functional connectivity and head-motion estimates, reduces inter-individual variability in global correlation structure, and results in comparable or improved identification of canonical functional-connectivity networks. Using task fMRI data across 47 contrasts from 7 tasks in the Human Connectome Project, we also present evidence that DiCER is more successful than GSR in preserving the spatial structure of expected task-related activation patterns. Our findings indicate that care must be exercised when examining WSDs (and their possible removal) in rsfMRI data, and that DiCER is a viable alternative to GSR for removing anatomically widespread and temporally coherent signals. All code for implementing DiCER and replicating our results is available at https://github.com/BMHLab/DiCER.

Time of day is associated with paradoxical reductions in global signal fluctuation and functional connectivity

The brain exhibits substantial diurnal variation in physiology and function, but neuroscience studies rarely report or consider the effects of time of day. Here, we examined variation in resting-state functional MRI (fMRI) in around 900 individuals scanned between 8 AM and 10 PM on two different days. Multiple studies across animals and humans have demonstrated that the brain’s global signal (GS) amplitude (henceforth referred to as “fluctuation”) increases with decreased arousal. Thus, in accord with known circadian variation in arousal, we hypothesised that GS fluctuation would be lowest in the morning, increase in the midafternoon, and dip in the early evening. Instead, we observed a cumulative decrease in GS fluctuation as the day progressed. Although respiratory variation also decreased with time of day, control analyses suggested that this did not account for the reduction in GS fluctuation. Finally, time of day was associated with marked decreases in resting-state functional connectivity across the whole brain. The magnitude of decrease was significantly stronger than associations between functional connectivity and behaviour (e.g., fluid intelligence). These findings reveal time of day effects on global brain activity that are not easily explained by expected arousal state or physiological artefacts. We conclude by discussing potential mechanisms for the observed diurnal variation in resting brain activity and the importance of accounting for time of day in future studies.

The brain exhibits substantial diurnal variation in physiology and function. A large-scale fMRI study reveals that the brain’s global signal amplitude, typically elevated during drowsy states, unexpectedly reduces steadily as the day progresses.

Questions and controversies in the study of time-varying functional connectivity in resting fMRI

The brain is a complex, multiscale dynamical system composed of many interacting regions. Knowledge of the spatiotemporal organization of these interactions is critical for establishing a solid understanding of the brain's functional architecture and the relationship between neural dynamics and cognition in health and disease. The possibility of studying these dynamics through careful analysis of neuroimaging data has catalyzed substantial interest in methods that estimate time-resolved fluctuations in functional connectivity (often referred to as "dynamic" or time-varying functional connectivity; TVFC). At the same time, debates have emerged regarding the application of TVFC analyses to resting fMRI data, and about the statistical validity, physiological origins, and cognitive and behavioral relevance of resting TVFC. These and other unresolved issues complicate interpretation of resting TVFC findings and limit the insights that can be gained from this promising new research area. This article brings together scientists with a variety of perspectives on resting TVFC to review the current literature in light of these issues. We introduce core concepts, define key terms, summarize controversies and open questions, and present a forward-looking perspective on how resting TVFC analyses can be rigorously and productively applied to investigate a wide range of questions in cognitive and systems neuroscience.

Brain networks, dimensionality, and global signal averaging in resting-state fMRI: Hierarchical network structure results in low-dimensional spatiotemporal dynamics

One of the most controversial practices in resting-state fMRI functional connectivity studies is whether or not to regress out the global average brain signal (GS) during artifact removal. Some groups have argued that it is absolutely essential to regress out the GS in order to fully remove head motion, respiration, and other global imaging artifacts. Others have argued that removing the GS distorts the resulting correlation matrices and inappropriately alters the results of group comparisons and relationships to behavior. At the core of this argument is the assessment of dimensionality in terms of the number of brain networks with uncorrelated time series. If the dimensionality is high, then the distortions due to GS removal could be effectively negligible. In the current paper, we examine the dimensionality of resting-state fMRI data using principal component analyses (PCA) and network clustering analyses. In two independent datasets (Set 1: N = 62, Set 2: N = 32), scree plots of the eigenvalues level off at or prior to 10 principal components, with prominent elbows at 3 and 7 components. While network clustering analyses have previously demonstrated that numerous networks can be distinguished with high thresholding of the voxel-wise correlation matrices, lower thresholding reveals a lower-dimensional hierarchical structure, with the first prominent branch at 2 networks (corresponding to the previously described "task-positive"/"task-negative" distinction) and further stable subdivisions at 4, 7 and 17. Since inter-correlated time series within these larger branches do not cancel to zero when averaged, the hierarchical nature of the correlation structure results in low effective dimensionality. Consistent with this, partial correlation analyses revealed that network-specific variance remains present in the GS at each level of the hierarchy, accounting for at least 14-18% of the overall GS variance in each dataset. These results demonstrate that GS regression is expected to remove substantial portions of network-specific brain signals along with artifacts, not simply whole-brain signals corresponding to arousal levels. We highlight alternative means of controlling for residual global artifacts when not removing the GS.

Impaired vigilance networks in temporal lobe epilepsy: Mechanisms and clinical implications

Mesial temporal lobe epilepsy (mTLE) is a neurological disorder in which patients suffer from frequent consciousness-impairing seizures, broad neurocognitive deficits, and diminished quality of life. Although seizures in mTLE originate focally in the hippocampus or amygdala, mTLE patients demonstrate cognitive deficits that extend beyond temporal lobe function-such as decline in executive function, cognitive processing speed, and attention-as well as diffuse decreases in neocortical metabolism and functional connectivity. Given prior observations that mTLE patients exhibit impairments in vigilance, and that seizures may disrupt the activity and long-range connectivity of subcortical brain structures involved in vigilance regulation, we propose that subcortical activating networks underlying vigilance play a critical role in mediating the widespread neural and cognitive effects of focal mTLE. Here, we review evidence for impaired vigilance in mTLE, examine clinical implications and potential network underpinnings, and suggest neuroimaging strategies for determining the relationship between vigilance, brain connectivity, and neurocognition in patients and healthy controls.

Physiological changes in sleep that affect fMRI inference

fMRI relies on a localized cerebral blood flow (CBF) response to changes in cortical neuronal activity. An underappreciated aspect however is its sensitivity to contributions from autonomic physiology that may affect CBF through changes in vascular resistance and blood pressure. As is reviewed here, this is crucial to consider in fMRI studies of sleep, given the close linkage between the regulation of arousal state and autonomic physiology. Typical methods for separating these effects are based on the use of reference signals that may include physiological parameters such as heart rate and respiration; however, the use of time-invariant models may not be adequate due to the possibly changing relationship between reference and fMRI signals with arousal state. In addition, recent research indicates that additional physiological reference signals may be needed to accurately describe changes in systemic physiology, including sympathetic indicators such as finger skin vascular tone and blood pressure.

Modulation of simultaneously collected hemodynamic and electrophysiological functional connectivity by ketamine and midazolam

The pharmacological modulation of functional connectivity in the brain may underlie therapeutic efficacy for several neurological and psychiatric disorders. Functional magnetic resonance imaging (fMRI) provides a noninvasive method of assessing this modulation, however, the indirect nature of the blood‐oxygen level dependent signal restricts the discrimination of neural from physiological contributions. Here we followed two approaches to assess the validity of fMRI functional connectivity in developing drug biomarkers, using simultaneous electroencephalography (EEG)/fMRI in a placebo‐controlled, three‐way crossover design with ketamine and midazolam. First, we compared seven different preprocessing pipelines to determine their impact on the connectivity of common resting‐state networks. Independent components analysis (ICA)‐denoising resulted in stronger reductions in connectivity after ketamine, and weaker increases after midazolam, than pipelines employing physiological noise modelling or averaged signals from cerebrospinal fluid or white matter. This suggests that pipeline decisions should reflect a drug's unique noise structure, and if this is unknown then accepting possible signal loss when choosing extensive ICA denoising pipelines could engender more confidence in the remaining results. We then compared the temporal correlation structure of fMRI to that derived from two connectivity metrics of EEG, which provides a direct measure of neural activity. While electrophysiological estimates based on the power envelope were more closely aligned to BOLD signal connectivity than those based on phase consistency, no significant relationship between the change in electrophysiological and hemodynamic correlation structures was found, implying caution should be used when making cross‐modal comparisons of pharmacologically‐modulated functional connectivity.

Global brain signal in awake rats

Although often used as a nuisance in resting-state functional magnetic resonance imaging (rsfMRI), the global brain signal in humans and anesthetized animals has important neural basis. However, our knowledge of the global signal in awake rodents is sparse. To bridge this gap, we systematically analyzed rsfMRI data acquired with a conventional single-echo (SE) echo planar imaging (EPI) sequence in awake rats. The spatial pattern of rsfMRI frames during peaks of the global signal exhibited prominent co-activations in the thalamo-cortical and hippocampo-cortical networks, as well as in the basal forebrain, hinting that these neural networks might contribute to the global brain signal in awake rodents. To validate this concept, we acquired rsfMRI data using a multi-echo (ME) EPI sequence and removed non-neural components in the rsfMRI signal. Consistent co-activation patterns were obtained in extensively de-noised ME-rsfMRI data, corroborating the finding from SE-rsfMRI data. Furthermore, during rsfMRI experiments, we simultaneously recorded neural spiking activities in the hippocampus using GCaMP-based fiber photometry. The hippocampal calcium activity exhibited significant correspondence with the global rsfMRI signal. These data collectively suggest that the global rsfMRI signal contains significant neural components that involve coordinated activities in the thalamo-cortical and hippocampo-cortical networks. These results provide important insight into the neural substrate of the global brain signal in awake rodents.

Identifying reproducible individual differences in childhood functional brain networks: An ABCD study

The 21-site Adolescent Brain Cognitive Development (ABCD) study provides an unparalleled opportunity to characterize functional brain development via resting-state functional connectivity (RSFC) and to quantify relationships between RSFC and behavior. This multi-site data set includes potentially confounding sources of variance, such as differences between data collection sites and/or scanner manufacturers, in addition to those inherent to RSFC (e.g., head motion). The ABCD project provides a framework for characterizing and reproducing RSFC and RSFC-behavior associations, while quantifying the extent to which sources of variability bias RSFC estimates. We quantified RSFC and functional network architecture in 2,188 9-10-year old children from the ABCD study, segregated into demographically-matched discovery (N = 1,166) and replication datasets (N = 1,022). We found RSFC and network architecture to be highly reproducible across children. We did not observe strong effects of site; however, scanner manufacturer effects were large, reproducible, and followed a "short-to-long" association with distance between regions. Accounting for potential confounding variables, we replicated that RSFC between several higher-order networks was related to general cognition. In sum, we provide a framework for how to characterize RSFC-behavior relationships in a rigorous and reproducible manner using the ABCD dataset and other large multi-site projects.

Mapping optogenetically-driven single-vessel fMRI with concurrent neuronal calcium recordings in the rat hippocampus

Extensive in vivo imaging studies investigate the hippocampal neural network function, mainly focusing on the dorsal CA1 region given its optical accessibility. Multi-modality fMRI with simultaneous hippocampal electrophysiological recording reveal broad cortical correlation patterns, but the detailed spatial hippocampal functional map remains lacking given the limited fMRI resolution. In particular, hemodynamic responses linked to specific neural activity are unclear at the single-vessel level across hippocampal vasculature, which hinders the deciphering of the hippocampal malfunction in animal models and the translation to critical neurovascular coupling (NVC) patterns for human fMRI. We simultaneously acquired optogenetically-driven neuronal Ca2+ signals with single-vessel blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-fMRI from individual venules and arterioles. Distinct spatiotemporal patterns of hippocampal hemodynamic responses were correlated to optogenetically evoked and spreading depression-like calcium events. The calcium event-related single-vessel hemodynamic modeling revealed significantly reduced NVC efficiency upon spreading depression-like (SDL) events, providing a direct measure of the NVC function at various hippocampal states.

BOLD and EEG signal variability at rest differently relate to aging in the human brain

Variability of neural activity is regarded as a crucial feature of healthy brain function, and several neuroimaging approaches have been employed to assess it noninvasively. Studies on the variability of both evoked brain response and spontaneous brain signals have shown remarkable changes with aging but it is unclear if the different measures of brain signal variability - identified with either hemodynamic or electrophysiological methods - reflect the same underlying physiology. In this study, we aimed to explore age differences of spontaneous brain signal variability with two different imaging modalities (EEG, fMRI) in healthy younger (25 ​± ​3 years, N ​= ​135) and older (67 ​± ​4 years, N ​= ​54) adults. Consistent with the previous studies, we found lower blood oxygenation level dependent (BOLD) variability in the older subjects as well as less signal variability in the amplitude of low-frequency oscillations (1-12 ​Hz), measured in source space. These age-related reductions were mostly observed in the areas that overlap with the default mode network. Moreover, age-related increases of variability in the amplitude of beta-band frequency EEG oscillations (15-25 ​Hz) were seen predominantly in temporal brain regions. There were significant sex differences in EEG signal variability in various brain regions while no significant sex differences were observed in BOLD signal variability. Bivariate and multivariate correlation analyses revealed no significant associations between EEG- and fMRI-based variability measures. In summary, we show that both BOLD and EEG signal variability reflect aging-related processes but are likely to be dominated by different physiological origins, which relate differentially to age and sex.

Sympathetic activity contributes to the fMRI signal

The interpretation of functional magnetic resonance imaging (fMRI) studies of brain activity is often hampered by the presence of brain-wide signal variations that may arise from a variety of neuronal and non-neuronal sources. Recent work suggests a contribution from the sympathetic vascular innervation, which may affect the fMRI signal through its putative and poorly understood role in cerebral blood flow (CBF) regulation. By analyzing fMRI and (electro-) physiological signals concurrently acquired during sleep, we found that widespread fMRI signal changes often co-occur with electroencephalography (EEG) K-complexes, signatures of sub-cortical arousal, and episodic drops in finger skin vascular tone; phenomena that have been associated with intermittent sympathetic activity. These findings support the notion that the extrinsic sympathetic innervation of the cerebral vasculature contributes to CBF regulation and the fMRI signal. Accounting for this mechanism could help separate systemic from local signal contributions and improve interpretation of fMRI studies.

Arousal Contributions to Resting-State fMRI Connectivity and Dynamics

Resting-state functional magnetic resonance imaging (rsfMRI) is being widely used for charting brain connectivity and dynamics in healthy and diseased brains. However, the resting state paradigm allows an unconstrained fluctuation of brain arousal, which may have profound effects on resting-state fMRI signals and associated connectivity/dynamic metrics. Here, we review current understandings of the relationship between resting-state fMRI and brain arousal, in particular the effect of a recently discovered event of arousal modulation on resting-state fMRI. We further discuss potential implications of arousal-related fMRI modulation with a focus on its potential role in mediating spurious correlations between resting-state connectivity/dynamics with physiology and behavior. Multiple hypotheses are formulated based on existing evidence and remain to be tested by future studies.

A dual role of prestimulus spontaneous neural activity in visual object recognition

Vision relies on both specific knowledge of visual attributes, such as object categories, and general brain states, such as those reflecting arousal. We hypothesized that these phenomena independently influence recognition of forthcoming stimuli through distinct processes reflected in spontaneous neural activity. Here, we recorded magnetoencephalographic (MEG) activity in participants (N = 24) who viewed images of objects presented at recognition threshold. Using multivariate analysis applied to sensor-level activity patterns recorded before stimulus presentation, we identified two neural processes influencing subsequent subjective recognition: a general process, which disregards stimulus category and correlates with pupil size, and a specific process, which facilitates category-specific recognition. The two processes are doubly-dissociable: the general process correlates with changes in criterion but not in sensitivity, whereas the specific process correlates with changes in sensitivity but not in criterion. Our findings reveal distinct mechanisms of how spontaneous neural activity influences perception and provide a framework to integrate previous findings.

The effect of spontaneous variations in prestimulus neural activity on subsequent perception is incompletely understood. Here, using MEG, the authors identify two distinct neural processes that can influence object recognition in different ways.

Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration

Functional magnetic resonance imaging (fMRI) is widely viewed as the gold standard for studying brain function due to its high spatial resolution and non-invasive nature. However, it is well established that changes in breathing patterns and heart rate strongly influence the blood oxygen-level dependent (BOLD) fMRI signal and this, in turn, can have considerable effects on fMRI studies, particularly resting-state studies. The dynamic effects of physiological processes are often quantified by using convolution models along with simultaneously recorded physiological data. In this context, physiological response function (PRF) curves (cardiac and respiratory response functions), which are convolved with the corresponding physiological fluctuations, are commonly employed. While it has often been suggested that the PRF curves may be region- or subject-specific, it is still an open question whether this is the case. In the present study, we propose a novel framework for the robust estimation of PRF curves and use this framework to rigorously examine the implications of using population-, subject-, session- and scan-specific PRF curves. The proposed framework was tested on resting-state fMRI and physiological data from the Human Connectome Project. Our results suggest that PRF curves vary significantly across subjects and, to a lesser extent, across sessions from the same subject. These differences can be partly attributed to physiological variables such as the mean and variance of the heart rate during the scan. The proposed methodological framework can be used to obtain robust scan-specific PRF curves from data records with duration longer than 5 min, exhibiting significantly improved performance compared to previously defined canonical cardiac and respiration response functions. Besides removing physiological confounds from the BOLD signal, accurate modeling of subject- (or session-/scan-) specific PRF curves is of importance in studies that involve populations with altered vascular responses, such as aging subjects.

Imaging the spontaneous flow of thought: Distinct periods of cognition contribute to dynamic functional connectivity during rest

Brain functional connectivity (FC) changes have been measured across seconds using fMRI. This is true for both rest and task scenarios. Moreover, it is well accepted that task engagement alters FC, and that dynamic estimates of FC during and before task events can help predict their nature and performance. Yet, when it comes to dynamic FC (dFC) during rest, there is no consensus about its origin or significance. Some argue that rest dFC reflects fluctuations in on-going cognition, or is a manifestation of intrinsic brain maintenance mechanisms, which could have predictive clinical value. Conversely, others have concluded that rest dFC is mostly the result of sampling variability, head motion or fluctuating sleep states. Here, we present novel analyses suggesting that rest dFC is influenced by short periods of spontaneous cognitive-task-like processes, and that the cognitive nature of such mental processes can be inferred blindly from the data. As such, several different behaviorally relevant whole-brain FC configurations may occur during a single rest scan even when subjects were continuously awake and displayed minimal motion. In addition, using low dimensional embeddings as visualization aids, we show how FC states-commonly used to summarize and interpret resting dFC-can accurately and robustly reveal periods of externally imposed tasks; however, they may be less effective in capturing periods of distinct cognition during rest.

Spatiotemporal Empirical Mode Decomposition of Resting-State fMRI Signals: Application to Global Signal Regression

Resting-state functional connectivity MRI (rs-fcMRI) is a common method for mapping functional brain networks. However, estimation of these networks is affected by the presence of a common global systemic noise, or global signal (GS). Previous studies have shown that the common preprocessing steps of removing the GS may create spurious correlations between brain regions. In this paper, we decompose fMRI signals into 5 spatial and 3 temporal intrinsic mode functions (SIMF and TIMF, respectively) by means of the empirical mode decomposition (EMD), which is an adaptive data-driven method widely used to analyze non-linear and non-stationary phenomena. For each SIMF, functional connectivity matrices were computed by means of Pearson correlation between TIMFs of different brain areas. Thus, instead of a single connectivity matrix, we obtained 5 × 3 = 15 functional connectivity matrices. Given the high correlation and global efficiency values of the connectivity matrices related to the low spatial maps (SIMF3, SIMF4, and SIMF5), our results suggest that these maps can be considered as spatial global signal masks. Thus, by summing up the first two SIMFs extracted from the fMRI signals, we have automatically excluded the GS which is now voxel-specific. We compared the performance of our method with the conventional GS regression and to the results when the GS was not removed. While the correlation pattern identified by the other methods suffers from a low level of precision in identifying the correct brain network connectivity, our approach demonstrated expected connectivity patterns for the default mode network and task-positive network.

Twitches, Blinks, and Fidgets: Important Generators of Ongoing Neural Activity

Animals and humans continuously engage in small, spontaneous motor actions, such as blinking, whisking, and postural adjustments ("fidgeting"). These movements are accompanied by changes in neural activity in sensory and motor regions of the brain. The frequency of these motions varies in time, is affected by sensory stimuli, arousal levels, and pathology. These fidgeting behaviors can be entrained by sensory stimuli. Fidgeting behaviors will cause distributed, bilateral functional activation in the 0.01 to 0.1 Hz frequency range that will show up in functional magnetic resonance imaging and wide-field calcium neuroimaging studies, and will contribute to the observed functional connectivity among brain regions. However, despite the large potential of these behaviors to drive brain-wide activity, these fidget-like behaviors are rarely monitored. We argue that studies of spontaneous and evoked brain dynamics in awake animals and humans should closely monitor these fidgeting behaviors. Differences in these fidgeting behaviors due to arousal or pathology will "contaminate" ongoing neural activity, and lead to apparent differences in functional connectivity. Monitoring and accounting for the brain-wide activations by these behaviors is essential during experiments to differentiate fidget-driven activity from internally driven neural dynamics.

Brain Dynamics: Global Pulse and Brain State Switching

A major challenge in systems-level neuroscience is to understand the dynamic formation and succession of brain states. A new study has extracted reproducible brain states from mouse resting-state fMRI data, revealing interactions between occurrences of these states and the phase of global signal fluctuations and alterations of the states in a mouse model of autism.

Vascular and neural basis of the BOLD signal

Neural activity in the brain is usually coupled to increases in local cerebral blood flow, leading to the increase in oxygenation that generates the BOLD fMRI signal. Recent work has begun to elucidate the vascular and neural mechanisms underlying the BOLD signal. The dilatory response is distributed throughout the vascular network. Arteries actively dilate within a second following neural activity increases, while venous distensions are passive and have a time course that last tens of seconds. Vasodilation, and thus local blood flow, is controlled by the activity of both neurons and astrocytes via multiple different pathways. The relationship between sensory-driven neural activity and the vascular dynamics in sensory areas are well-captured with a linear convolution model. However, depending on the behavioral state or brain region, the coupling between neural activity and hemodynamic signals can be weak or even inverted.

Investigating the Intersession Reliability of Dynamic Brain-State Properties

Dynamic functional connectivity metrics have much to offer to the neuroscience of individual differences of cognition. Yet, despite the recent expansion in dynamic connectivity research, limited resources have been devoted to the study of the reliability of these connectivity measures. To address this, resting-state functional magnetic resonance imaging data from 100 Human Connectome Project subjects were compared across 2 scan days. Brain states (i.e., patterns of coactivity across regions) were identified by classifying each time frame using k means clustering. This was done with and without global signal regression (GSR). Multiple gauges of reliability indicated consistency in the brain-state properties across days and GSR attenuated the reliability of the brain states. Changes in the brain-state properties across the course of the scan were investigated as well. The results demonstrate that summary metrics describing the clustering of individual time frames have adequate test/retest reliability, and thus, these patterns of brain activation may hold promise for individual-difference research.

The resting-state fMRI arterial signal predicts differential blood transit time through the brain

Previous studies have found that aperiodic, systemic low-frequency oscillations (sLFOs) are present in blood-oxygen-level-dependent (BOLD) data. These signals are in the same low frequency band as the "resting state" signal; however, they are distinct signals which represent non-neuronal, physiological oscillations. The same sLFOs are found in the periphery (i.e. finger tips) as changes in oxy/deoxy-hemoglobin concentration using concurrent near-infrared spectroscopy. Together, this evidence points toward an extra-cerebral origin of these sLFOs. If this is the case, it is expected that these sLFO signals would be found in the carotid arteries with time delays that precede the signals found in the brain. To test this hypothesis, we employed the publicly available MyConnectome dataset (a two-year longitudinal study of a single subject) to extract the sLFOs in the internal carotid arteries (ICAs) with the help of the T1/T2-weighted images. Significant, but negative, correlations were found between the LFO BOLD signals from the ICAs and (1) the global signal (GS), (2) the superior sagittal sinus, and (3) the jugulars. We found the consistent time delays between the sLFO signals from ICAs, GS and veins which coincide with the blood transit time through the cerebral vascular tree.

Neuromodulatory Influences on Integration and Segregation in the Brain

Cognitive function relies on the dynamic cooperation of specialized regions of the brain; however, the elements of the system responsible for coordinating this interaction remain poorly understood. In this Opinion article I argue that this capacity is mediated in part by competitive and cooperative dynamic interactions between two prominent metabotropic neuromodulatory systems - the cholinergic basal forebrain and the noradrenergic locus coeruleus (LC). I assert that activity in these projection nuclei regulates the amount of segregation and integration within the whole brain network by modulating the activity of a diverse set of specialized regions of the brain on a timescale relevant for cognition and attention.

Task-related hemodynamic responses are modulated by reward and task engagement

Hemodynamic recordings from visual cortex contain powerful endogenous task-related responses that may reflect task-related arousal, or "task engagement" distinct from attention. We tested this hypothesis with hemodynamic measurements (intrinsic-signal optical imaging) from monkey primary visual cortex (V1) while the animals' engagement in a periodic fixation task over several hours was varied through reward size and as animals took breaks. With higher rewards, animals appeared more task-engaged; task-related responses were more temporally precise at the task period (approximately 10-20 seconds) and modestly stronger. The 2-5 minute blocks of high-reward trials led to ramp-like decreases in mean local blood volume; these reversed with ramp-like increases during low reward. The blood volume increased even more sharply when the animal shut his eyes and disengaged completely from the task (5-10 minutes). We propose a mechanism that controls vascular tone, likely along with local neural responses in a manner that reflects task engagement over the full range of timescales tested.

The relationship between heart rate and functional connectivity of brain regions involved in autonomic control

The peripheral autonomic nervous system (ANS) adjusts the heart rate (HR) to intrinsic and extrinsic demands. It is controlled by a group of functionally connected brain regions assembling the so-called central autonomic network (CAN). More specifically, forebrain cortical regions, limbic and brainstem structures within the CAN have been identified as important components of circuits involved in HR regulation. The present study aimed to investigate whether functional connectivity (FC) between these regions varies in subjects with different heart rates. Thus, 84 healthy subjects were separated according to their HR in slow, medium and fast. We observed a direct association between HR and FC in CAN regions, where stronger FC was related to slower HR. This relationship, however, is non-linear, follows an exponential course and is not restricted to CAN areas only. The network-based analysis (NBS) using time series from 262 independent anatomical ROIs revealed significantly increased functional connectivity in subjects with slow HR compared to subjects with fast HR mainly in regions being part of the salience network, but also of the default-mode network. We additionally simulated the effect of aliasing on the functional connectivity using several TRs and heart rates to exclude the possibility that FC differences might be due to different aliasing effects in the data. The result of the simulation indicated that aliasing cannot explain our findings. Thus, present results imply a functionally meaningful coupling between FC and HR that need to be accounted for in future studies. Moreover, given the established link between HR and emotional, cognitive and social processes, present findings may also be considered to explain individual differences in brain activation or connectivity when using corresponding paradigms in the MR scanner to investigate such processes.

Global signal regression strengthens association between resting-state functional connectivity and behavior

Global signal regression (GSR) is one of the most debated preprocessing strategies for resting-state functional MRI. GSR effectively removes global artifacts driven by motion and respiration, but also discards globally distributed neural information and introduces negative correlations between certain brain regions. The vast majority of previous studies have focused on the effectiveness of GSR in removing imaging artifacts, as well as its potential biases. Given the growing interest in functional connectivity fingerprinting, here we considered the utilitarian question of whether GSR strengthens or weakens associations between resting-state functional connectivity (RSFC) and multiple behavioral measures across cognition, personality and emotion. By applying the variance component model to the Brain Genomics Superstruct Project (GSP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 47% across 23 behavioral measures after GSR. In the Human Connectome Project (HCP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 40% across 58 behavioral measures, when GSR was applied after ICA-FIX de-noising. To ensure generalizability, we repeated our analyses using kernel regression. GSR improved behavioral prediction accuracies by an average of 64% and 12% in the GSP and HCP datasets respectively. Importantly, the results were consistent across methods. A behavioral measure with greater RSFC-explained variance (using the variance component model) also exhibited greater prediction accuracy (using kernel regression). A behavioral measure with greater improvement in behavioral variance explained after GSR (using the variance component model) also enjoyed greater improvement in prediction accuracy after GSR (using kernel regression). Furthermore, GSR appeared to benefit task performance measures more than self-reported measures. Since GSR was more effective at removing motion-related and respiratory-related artifacts, GSR-related increases in variance explained and prediction accuracies were unlikely the result of motion-related or respiratory-related artifacts. However, it is worth emphasizing that the current study focused on whole-brain RSFC, so it remains unclear whether GSR improves RSFC-behavioral associations for specific connections or networks. Overall, our results suggest that at least in the case for young healthy adults, GSR strengthens the associations between RSFC and most (although not all) behavioral measures. Code for the variance component model and ridge regression can be found here: https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/preprocessing/Li2019_GSR.

Instability of brain connectivity during nonrapid eye movement sleep reflects altered properties of information integration

Nonrapid eye movement (NREM) sleep is associated with fading consciousness in humans. Recent neuroimaging studies have demonstrated the spatiotemporal alterations of the brain functional connectivity (FC) in NREM sleep, suggesting the changes of information integration in the sleeping brain. However, the common stationarity assumption in FC does not satisfactorily explain the dynamic process of information integration during sleep. The dynamic FC (dFC) across brain networks is speculated to better reflect the time-varying information propagation during sleep. Accordingly, we conducted simultaneous EEG-fMRI recordings involving 12 healthy men during sleep and observed dFC across sleep stages using the sliding-window approach. We divided dFC into two aspects: mean dFC (dFC_mean ) and variance dFC (dFC_var ). A high dFC_mean indicates stable brain network integrity, whereas a high dFC_var indicates instability of information transfer within and between functional networks. For the network-based dFC, the dFC_var were negatively correlated with the dFC_mean across the waking and three NREM sleep stages. As sleep deepened, the dFC_mean decreased (N0~N1 > N2 > N3), whereas the dFC_var peaked during the N2 stage (N0~N1 < N3 < N2). The highest dFC_var during the N2 stage indicated the unstable synchronizations across the entire brain. In the N3 stage, the overall disrupted network integration was observed through the lowest dFC_mean and elevated dFC_var, compared with N0 and N1. Conclusively, when the network specificity (dFC_mean ) breaks down, the consciousness dissipates with increasing variability of information exchange (dFC_var ).