Differences in the resting-state fMRI global signal amplitude between the eyes open and eyes closed states are related to changes in EEG vigilance

Revisiting alpha resting state dynamics underlying hallucinatory vulnerability: Insights from Hidden semi-Markov Modeling

BACKGROUND

Resting state (RS) brain activity is inherently non-stationary. Hidden semi-Markov Models (HsMM) can characterize continuous RS data as a sequence of recurring and distinct brain states along with their spatio-temporal dynamics.

NEW METHOD

Recent explorations suggest that HsMM state dynamics in the alpha frequency band link to auditory hallucination proneness (HP) in non-clinical individuals. The present study aimed to replicate these findings to elucidate robust neural correlates of hallucinatory vulnerability. Specifically, we aimed to investigate the reproducibility of HsMM states across different data sets and within-data set variants as well as the replicability of the association between alpha brain state dynamics and HP.

RESULTS

We found that most brain states are reproducible in different data sets, confirming that the HsMM characterized robust and generalizable EEG RS dynamics on a sub-second timescale. Brain state topographies and temporal dynamics of different within-data set variants showed substantial similarities and were robust against reduced data length and number of electrodes. However, the association with HP was not directly reproducible across data sets.

COMPARISON WITH EXISTING METHODS

The HsMM optimally leverages the high temporal resolution of EEG data and overcomes time-domain restrictions of other state allocation methods.

CONCLUSION

The results indicate that the sensitivity of brain state dynamics to capture individual variability in HP may depend on the data recording characteristics and individual variability in RS cognition, such as mind wandering. Future studies should consider that the order in which eyes-open and eyes-closed RS data are acquired directly influences an individual's attentional state and generation of spontaneous thoughts, and thereby might mediate the link to hallucinatory vulnerability.

The neural correlates of arousal: Ventral posterolateral nucleus-global transient co-activation

Arousal and awareness are two components of consciousness whose neural mechanisms remain unclear. Spontaneous peaks of global (brain-wide) blood-oxygenation-level-dependent (BOLD) signal have been found to be sensitive to changes in arousal. By contrasting BOLD signals at different arousal levels, we find decreased activation of the ventral posterolateral nucleus (VPL) during transient peaks in the global signal in low arousal and awareness states (non-rapid eye movement sleep and anesthesia) compared to wakefulness and in eyes-closed compared to eyes-open conditions in healthy awake individuals. Intriguingly, VPL-global co-activation remains high in patients with unresponsive wakefulness syndrome (UWS), who exhibit high arousal without awareness, while it reduces in rapid eye movement sleep, a state characterized by low arousal but high awareness. Furthermore, lower co-activation is found in individuals during N3 sleep compared to patients with UWS. These results demonstrate that co-activation of VPL and global activity is critical to arousal but not to awareness.

Modular brain network in volitional eyes closing: enhanced integration with a marked impact on hubs

Volitional eyes closing would shift brain's information processing modes from the "exteroceptive" to "interoceptive" state. This transition induced by the eyes closing is underpinned by a large-scale reconfiguration of brain network, which is still not fully comprehended. Here, we investigated the eyes-closing-relevant network reconfiguration by examining the functional integration among intrinsic modules. Our investigation utilized a publicly available dataset with 48 subjects being scanned in both eyes closed and eyes open conditions. It was found that the modular integration was significantly enhanced during the eyes closing, including lower modularity index, higher participation coefficient, less provincial hubs, and more connector hubs. Moreover, the eyes-closing-enhanced integration was particularly noticeable in the hubs of network, mainly located in the default-mode network. Finally, the hub-dominant modular enhancement was positively correlated to the eyes-closing-reduced entropy of BOLD signal, suggesting a close connection to the diminished consciousness of individuals. Collectively, our findings strongly suggested that the enhanced modular integration with substantially reorganized hubs characterized the large-scale cortical underpinning of the volitional eyes closing.

Neurobiological Underpinnings of Hyperarousal in Depression: A Comprehensive Review

Patients with major depressive disorder (MDD) exhibit an abnormal physiological arousal pattern known as hyperarousal, which may contribute to their depressive symptoms. However, the neurobiological mechanisms linking this abnormal arousal to depressive symptoms are not yet fully understood. In this review, we summarize the physiological and neural features of arousal, and review the literature indicating abnormal arousal in depressed patients. Evidence suggests that a hyperarousal state in depression is characterized by abnormalities in sleep behavior, physiological (e.g., heart rate, skin conductance, pupil diameter) and electroencephalography (EEG) features, and altered activity in subcortical (e.g., hypothalamus and locus coeruleus) and cortical regions. While recent studies highlight the importance of subcortical-cortical interactions in arousal, few have explored the relationship between subcortical-cortical interactions and hyperarousal in depressed patients. This gap limits our understanding of the neural mechanism through which hyperarousal affects depressive symptoms, which involves various cognitive processes and the cerebral cortex. Based on the current literature, we propose that the hyperconnectivity in the thalamocortical circuit may contribute to both the hyperarousal pattern and depressive symptoms. Future research should investigate the relationship between thalamocortical connections and abnormal arousal in depression, and explore its implications for non-invasive treatments for depression.

Brain Connectivity Correlates of Breathing and Cardiac Irregularities in SUDEP: A Resting-State fMRI Study

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of premature mortality among people with epilepsy. Evidence from witnessed and monitored SUDEP cases indicate seizure-induced cardiovascular and respiratory failures; yet, the underlying mechanisms remain obscure. SUDEP occurs often during the night and early morning hours, suggesting that sleep or circadian rhythm-induced changes in physiology contribute to the fatal event. Resting-state fMRI studies have found altered functional connectivity between brain structures involved in cardiorespiratory regulation in later SUDEP cases and in individuals at high-risk of SUDEP. However, those connectivity findings have not been related to changes in cardiovascular or respiratory patterns. Here, we compared fMRI patterns of brain connectivity associated with regular and irregular cardiorespiratory rhythms in SUDEP cases with those of living epilepsy patients of varying SUDEP risk, and healthy controls. We analysed resting-state fMRI data from 98 patients with epilepsy (9 who subsequently succumbed to SUDEP, 43 categorized as low SUDEP risk (no tonic-clonic seizures (TCS) in the year preceding the fMRI scan), and 46 as high SUDEP risk (>3 TCS in the year preceding the scan)) and 25 healthy controls. The global signal amplitude (GSA), defined as the moving standard deviation of the fMRI global signal, was used to identify periods with regular ('low state') and irregular ('high state') cardiorespiratory rhythms. Correlation maps were derived from seeds in twelve regions with a key role in autonomic or respiratory regulation, for the low and high states. Following principal component analysis, component weights were compared between the groups. We found widespread alterations in connectivity of precuneus/posterior cingulate cortex in epilepsy compared to controls, in the low state (regular cardiorespiratory activity). In the low state, and to a lesser degree in the high state, reduced anterior insula connectivity (mainly with anterior and posterior cingulate cortex) in epilepsy appeared, relative to healthy controls. For SUDEP cases, the insula connectivity differences were inversely related to the interval between the fMRI scan and death. The findings suggest that anterior insula connectivity measures may provide a biomarker of SUDEP risk. The neural correlates of autonomic brain structures associated with different cardiorespiratory rhythms may shed light on the mechanisms underlying terminal apnea observed in SUDEP.

Brain State Relays Self-Processing and Heartbeat-Evoked Cortical Responses

The self has been proposed to be grounded in interoceptive processing, with heartbeat-evoked cortical activity as a neurophysiological marker of this processing. However, inconsistent findings have been reported on the relationship between heartbeat-evoked cortical responses and self-processing (including exteroceptive- and mental-self-processing). In this review, we examine previous research on the association between self-processing and heartbeat-evoked cortical responses and highlight the divergent temporal-spatial characteristics and brain regions involved. We propose that the brain state relays the interaction between self-processing and heartbeat-evoked cortical responses and thus accounts for the inconsistency. The brain state, spontaneous brain activity which highly and continuously changes in a nonrandom way, serves as the foundation upon which the brain functions and was proposed as a point in an extremely high-dimensional space. To elucidate our assumption, we provide reviews on the interactions between dimensions of brain state with both self-processing and heartbeat-evoked cortical responses. These interactions suggest the relay of self-processing and heartbeat-evoked cortical responses by brain state. Finally, we discuss possible approaches to investigate whether and how the brain state impacts the self-heart interaction.

Altered global signal topography in Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is a neurodegenerative disease associated with widespread disruptions in intrinsic local specialization and global integration in the functional system of the brain. These changes in integration may further disrupt the global signal (GS) distribution, which might represent the local relative contribution to global activity in functional magnetic resonance imaging (fMRI).

METHODS

fMRI scans from a discovery dataset (n = 809) and a validated dataset (n = 542) were used in the analysis. We investigated the alteration of GS topography using the GS correlation (GSCORR) in patients with mild cognitive impairment (MCI) and AD. The association between GS alterations and functional network properties was also investigated based on network theory. The underlying mechanism of GSCORR alterations was elucidated using imaging-transcriptomics.

FINDINGS

Significantly increased GS topography in the frontal lobe and decreased GS topography in the hippocampus, cingulate gyrus, caudate, and middle temporal gyrus were observed in patients with AD (P_adj < 0.05). Notably, topographical GS changes in these regions correlated with cognitive ability (P < 0.05). The changes in GS topography also correlated with the changes in functional network segregation (ρ = 0.5). Moreover, the genes identified based on GS topographical changes were enriched in pathways associated with AD and neurodegenerative diseases.

INTERPRETATION

Our findings revealed significant changes in GS topography and its molecular basis, confirming the informative role of GS in AD and further contributing to the understanding of the relationship between global and local neuronal activities in patients with AD.

FUNDING

Beijing Natural Science Funds for Distinguished Young Scholars, China; Fundamental Research Funds for the Central Universities, China; National Natural Science Foundation, China.

Gestational age-related changes in the fetal functional connectome: in utero evidence for the global signal

The human brain begins to develop in the third gestational week and rapidly grows and matures over the course of pregnancy. Compared to fetal structural neurodevelopment, less is known about emerging functional connectivity in utero. Here, we investigated gestational age (GA)-associated in vivo changes in functional brain connectivity during the second and third trimesters in a large dataset of 110 resting-state functional magnetic resonance imaging scans from a cohort of 95 healthy fetuses. Using representational similarity analysis, a multivariate analytical technique that reveals pair-wise similarity in high-order space, we showed that intersubject similarity of fetal functional connectome patterns was strongly related to between-subject GA differences (r = 0.28, P < 0.01) and that GA sensitivity of functional connectome was lateralized, especially at the frontal area. Our analysis also revealed a subnetwork of connections that were critical for predicting age (mean absolute error = 2.72 weeks); functional connectome patterns of individual fetuses reliably predicted their GA (r = 0.51, P < 0.001). Lastly, we identified the primary principal brain network that tracked fetal brain maturity. The main network showed a global synchronization pattern resembling global signal in the adult brain.

Predicting age from resting-state scalp EEG signals with deep convolutional neural networks on TD-brain dataset

Introduction

Brain age prediction has been shown to be clinically relevant, with errors in its prediction associated with various psychiatric and neurological conditions. While the prediction from structural and functional magnetic resonance imaging data has been feasible with high accuracy, whether the same results can be achieved with electroencephalography is unclear.

Methods

The current study aimed to create a new deep learning solution for brain age prediction using raw resting-state scalp EEG. To this end, we utilized the TD-BRAIN dataset, including 1,274 subjects (both healthy controls and individuals with various psychiatric disorders, with a total of 1,335 recording sessions). To achieve the best age prediction, we used data augmentation techniques to increase the diversity of the training set and developed a deep convolutional neural network model.

Results

The model's training was done with 10-fold cross-subject cross-validation, with the EEG recordings of the subjects used for training not considered to test the model. In training, using the relative rather than the absolute loss function led to a better mean absolute error of 5.96 years in cross-validation. We found that the best performance could be achieved when both eyes-open and eyes-closed states are used simultaneously. The frontocentral electrodes played the most important role in age prediction.

Discussion

The architecture and training method of the proposed deep convolutional neural networks (DCNN) improve state-of-the-art metrics in the age prediction task using raw resting-state EEG data by 13%. Given that brain age prediction might be a potential biomarker of numerous brain diseases, inexpensive and precise EEG-based estimation of brain age will be in demand for clinical practice.

An orderly sequence of autonomic and neural events at transient arousal changes

Resting-state functional magnetic resonance imaging (rsfMRI) allows the study of functional brain connectivity based on spatially structured variations in neuronal activity. Proper evaluation of connectivity requires removal of non-neural contributions to the fMRI signal, in particular hemodynamic changes associated with autonomic variability. Regression analysis based on autonomic indicator signals has been used for this purpose, but may be inadequate if neuronal and autonomic activities covary. To investigate this potential co-variation, we performed rsfMRI experiments while concurrently acquiring electroencephalography (EEG) and autonomic indicator signals, including heart rate, respiratory depth, and peripheral vascular tone. We identified a recurrent and systematic spatiotemporal pattern of fMRI (named as fMRI cascade), which features brief signal reductions in salience and default-mode networks and the thalamus, followed by a biphasic global change with a sensory-motor dominance. This fMRI cascade, which was mostly observed during eyes-closed condition, was accompanied by large EEG and autonomic changes indicative of arousal modulations. Importantly, the removal of the fMRI cascade dynamics from rsfMRI diminished its correlations with various signals. These results suggest that the rsfMRI correlations with various physiological and neural signals are not independent but arise, at least partly, from the fMRI cascades and associated neural and physiological changes at arousal modulations.

A reusable benchmark of brain-age prediction from M/EEG resting-state signals

Population-level modeling can define quantitative measures of individual aging by applying machine learning to large volumes of brain images. These measures of brain age, obtained from the general population, helped characterize disease severity in neurological populations, improving estimates of diagnosis or prognosis. Magnetoencephalography (MEG) and Electroencephalography (EEG) have the potential to further generalize this approach towards prevention and public health by enabling assessments of brain health at large scales in socioeconomically diverse environments. However, more research is needed to define methods that can handle the complexity and diversity of M/EEG signals across diverse real-world contexts. To catalyse this effort, here we propose reusable benchmarks of competing machine learning approaches for brain age modeling. We benchmarked popular classical machine learning pipelines and deep learning architectures previously used for pathology decoding or brain age estimation in 4 international M/EEG cohorts from diverse countries and cultural contexts, including recordings from more than 2500 participants. Our benchmarks were built on top of the M/EEG adaptations of the BIDS standard, providing tools that can be applied with minimal modification on any M/EEG dataset provided in the BIDS format. Our results suggest that, regardless of whether classical machine learning or deep learning was used, the highest performance was reached by pipelines and architectures involving spatially aware representations of the M/EEG signals, leading to R^2 scores between 0.60-0.71. Hand-crafted features paired with random forest regression provided robust benchmarks even in situations in which other approaches failed. Taken together, this set of benchmarks, accompanied by open-source software and high-level Python scripts, can serve as a starting point and quantitative reference for future efforts at developing M/EEG-based measures of brain aging. The generality of the approach renders this benchmark reusable for other related objectives such as modeling specific cognitive variables or clinical endpoints.

Aperiodic measures of neural excitability are associated with anticorrelated hemodynamic networks at rest: A combined EEG-fMRI study

The hallmark of resting EEG spectra are distinct rhythms emerging from a broadband, aperiodic background. This aperiodic neural signature accounts for most of total EEG power, although its significance and relation to functional neuroanatomy remains obscure. We hypothesized that aperiodic EEG reflects a significant metabolic expenditure and therefore might be associated with the default mode network while at rest. During eyes-open, resting-state recordings of simultaneous EEG-fMRI, we find that aperiodic and periodic components of EEG power are only minimally associated with activity in the default mode network. However, a whole-brain analysis identifies increases in aperiodic power correlated with hemodynamic activity in an auditory-salience-cerebellar network, and decreases in aperiodic power are correlated with hemodynamic activity in prefrontal regions. Desynchronization in residual alpha and beta power is associated with visual and sensorimotor hemodynamic activity, respectively. These findings suggest that resting-state EEG signals acquired in an fMRI scanner reflect a balance of top-down and bottom-up stimulus processing, even in the absence of an explicit task.

A multi-measure approach for assessing the performance of fMRI preprocessing strategies in resting-state functional connectivity

It is well established that head motion and physiological processes (e.g. cardiac and breathing activity) should be taken into consideration when analyzing and interpreting results in fMRI studies. However, even though recent studies aimed to evaluate the performance of different preprocessing pipelines there is still no consensus on the optimal strategy. This is partly due to the fact that the quality control (QC) metrics used to evaluate differences in performance across pipelines have often yielded contradictory results. Furthermore, preprocessing techniques based on physiological recordings or data decomposition techniques (e.g. aCompCor) have not been comprehensively examined. Here, to address the aforementioned issues, we propose a framework that summarizes the scores from eight previously proposed and novel QC metrics to a reduced set of two QC metrics that reflect the signal-to-noise ratio and the reduction in motion artifacts and biases in the preprocessed fMRI data. Using this framework, we evaluate the performance of three commonly used practices on the quality of data: 1) Removal of nuisance regressors from fMRI data, 2) discarding motion-contaminated volumes (i.e., scrubbing) before regression, and 3) low-pass filtering the data and the nuisance regressors before their removal. Using resting-state fMRI data from the Human Connectome Project, we show that the scores of the examined QC metrics improve the most when the global signal (GS) and about 17% of principal components from white matter (WM) are removed from the data. Finally, we observe a small further improvement with low-pass filtering at 0.20 Hz and milder variants of WM denoising, but not with scrubbing.

Frequency-specific altered global signal topography in drug-naïve first-episode patients with adolescent-onset schizophrenia

Adolescent-onset schizophrenia (AOS) is a severe neuropsychiatric disease associated with frequency-specific abnormalities across distributed neural systems in a slow rhythm. Recently, functional magnetic resonance imaging (fMRI) studies have determined that the global signal. (GS) is an important source of the local neuronal activity in 0.01-0.1 Hz frequency band. However, it remains unknown whether the effects follow a specific spatially preferential pattern in different frequency bands in schizophrenia. To address this issue, resting-state fMRI data from 39 drug-naïve AOS patients and 31 healthy controls (HCs) were used to assess the changes in GS topography patterns in the slow-4 (0.027-0.073 Hz) and slow-5 bands (0.01-0.027 Hz). Results revealed that GS mainly affects the default mode network (DMN) in slow-4 and sensory regions in the slow-5 band respectively, and GS has a stronger driving effect in the slow-5 band. Moreover, significant frequency-by-group interaction was observed in the frontoparietal network. Compared with HCs, patients with AOS exhibited altered GS topography mainly located in the DMN. Our findings demonstrated that the influence of the GS on brain networks altered in a frequency-specific way in schizophrenia.

Brain Activity Fluctuations Propagate as Waves Traversing the Cortical Hierarchy

The brain exhibits highly organized patterns of spontaneous activity as measured by resting-state functional magnetic resonance imaging (fMRI) fluctuations that are being widely used to assess the brain's functional connectivity. Some evidence suggests that spatiotemporally coherent waves are a core feature of spontaneous activity that shapes functional connectivity, although this has been difficult to establish using fMRI given the temporal constraints of the hemodynamic signal. Here, we investigated the structure of spontaneous waves in human fMRI and monkey electrocorticography. In both species, we found clear, repeatable, and directionally constrained activity waves coursed along a spatial axis approximately representing cortical hierarchical organization. These cortical propagations were closely associated with activity changes in distinct subcortical structures, particularly those related to arousal regulation, and modulated across different states of vigilance. The findings demonstrate a neural origin of spatiotemporal fMRI wave propagation at rest and link it to the principal gradient of resting-state fMRI connectivity.

Reduced coupling between cerebrospinal fluid flow and global brain activity is linked to Alzheimer disease-related pathology

The glymphatic system plays an important role in clearing the amyloid-β (Aβ) and tau proteins that are closely linked to Alzheimer disease (AD) pathology. Glymphatic clearance, as well as Aβ accumulation, is highly dependent on sleep, but the sleep-dependent driving forces behind cerebrospinal fluid (CSF) movements essential to the glymphatic flux remain largely unclear. Recent studies have reported that widespread, high-amplitude spontaneous brain activations in the drowsy state and during sleep, which are shown as large global signal peaks in resting-state functional magnetic resonance imaging (rsfMRI), are coupled with CSF movements, suggesting their potential link to glymphatic flux and metabolite clearance. By analyzing multimodal data from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) project, here we showed that the coupling between the global fMRI signal and CSF influx is correlated with AD-related pathology, including various risk factors for AD, the severity of AD-related diseases, the cortical Aβ level, and cognitive decline over a 2-year follow-up. These results provide critical initial evidence for involvement of sleep-dependent global brain activity, as well as the associated physiological modulations, in the clearance of AD-related brain waste.

This study reveals strong coupling between the global fMRI signal and cerebrospinal fluid influx, finding that this is correlated with Alzheimer’s disease-related pathology, disease severity, and cognitive decline. This supports a link between spontaneous low-frequency brain dynamics and Alzheimer’s disease pathology, presumably due to their role in glymphatic clearance.

Respiratory, cardiac, EEG, BOLD signals and functional connectivity over multiple microsleep episodes

Falling asleep is common in fMRI studies. By using long eyelid closures to detect microsleep onset, we showed that the onset and termination of short sleep episodes invokes a systematic sequence of BOLD signal changes that are large, widespread, and consistent across different microsleep durations. The signal changes are intimately intertwined with shifts in respiration and heart rate, indicating that autonomic contributions are integral to the brain physiology evaluated using fMRI and cannot be simply treated as nuisance signals. Additionally, resting state functional connectivity (RSFC) was altered in accord with the frequency of falling asleep and in a manner that global signal regression does not eliminate. Our findings point to the need to develop a consensus among neuroscientists using fMRI on how to deal with microsleep intrusions. SIGNIFICANCE STATEMENT: Sleep, breathing and cardiac action are influenced by common brainstem nuclei. We show that falling asleep and awakening are associated with a sequence of BOLD signal changes that are large, widespread and consistent across varied durations of sleep onset and awakening. These signal changes follow closely those associated with deceleration and acceleration of respiration and heart rate, calling into question the separation of the latter signals as 'noise' when the frequency of falling asleep, which is commonplace in RSFC studies, correlates with the extent of RSFC perturbation. Autonomic and central nervous system contributions to BOLD signal have to be jointly considered when interpreting fMRI and RSFC studies.

Global Signal Topography of the Human Brain: A Novel Framework of Functional Connectivity for Psychological and Pathological Investigations

The global signal (GS), which was once regarded as a nuisance of functional magnetic resonance imaging, has been proven to convey valuable neural information. This raised the following question: what is a GS represented in local brain regions? In order to answer this question, the GS topography was developed to measure the correlation between global and local signals. It was observed that the GS topography has an intrinsic structure characterized by higher GS correlation in sensory cortices and lower GS correlation in higher-order cortices. The GS topography could be modulated by individual factors, attention-demanding tasks, and conscious states. Furthermore, abnormal GS topography has been uncovered in patients with schizophrenia, major depressive disorder, bipolar disorder, and epilepsy. These findings provide a novel insight into understanding how the GS and local brain signals coactivate to organize information in the human brain under various brain states. Future directions were further discussed, including the local-global confusion embedded in the GS correlation, the integration of spatial information conveyed by the GS, and temporal information recruited by the connection analysis. Overall, a unified psychopathological framework is needed for understanding the GS topography.

Beyond MRI: on the scientific value of combining non-human primate neuroimaging with metadata

Sharing and pooling large amounts of non-human primate neuroimaging data offer new exciting opportunities to understand the primate brain. The potential of big data in non-human primate neuroimaging could however be tremendously enhanced by combining such neuroimaging data with other types of information. Here we describe metadata that have been identified as particularly valuable by the non-human primate neuroimaging community, including behavioural, genetic, physiological and phylogenetic data.

Distinct alterations in resting-state electroencephalogram during eyes closed and eyes open and between morning and evening are present in first-episode psychosis patients

Abnormalities in resting-state electroencephalogram (rs-EEG) activity have been previously reported in schizophrenia. While most rs-EEG recordings were performed in patients with chronic schizophrenia during eyes closed (EC), only a handful of studies have investigated rs-EEG activity during both EC and eyes open (EO) conditions. It is also unknown whether EC and EO rs-EEG alterations are present at illness onset, and whether they change during the day. Here, we performed EC and EO rs-EEG recordings in the morning (AM) and evening (PM) in twenty-six first-episode psychosis (FEP) patients and seventeen matched healthy controls (HC). In AM/EC rs-EEG, a widespread reduction was found in low alpha power in FEP relative to HC. In PM/EC, the FEP group demonstrated a trend toward decreased theta power in parietal regions, while decreased high alpha power in frontal and left parietal regions was present during PM/EO. Moreover, reduced low alpha power during AM/EC was associated with worse positive symptoms. Altogether, those findings indicate that rs-EEG alterations are present in FEP patients at illness onset, that they are linked to the severity of their psychosis, and that distinct RS abnormalities can be detected in different conditions of visual alertness and time of the day. Future work should therefore account for those factors, which will help reduce variability of rs-EEG findings across studies and may serve as monitoring biomarkers of illness severity in schizophrenia and related psychotic disorders.

Amplitude of fNIRS Resting-State Global Signal Is Related to EEG Vigilance Measures: A Simultaneous fNIRS and EEG Study

Recently, functional near-infrared spectroscopy (fNIRS) has been utilized to image the hemodynamic activities and connectivity in the human brain. With the advantage of economic efficiency, portability, and fewer physical constraints, fNIRS enables studying of the human brain at versatile environment and various body positions, including at bed side and during exercise, which complements the use of functional magnetic resonance imaging (fMRI). However, like fMRI, fNIRS imaging can be influenced by the presence of a strong global component. Yet, the nature of the global signal in fNIRS has not been established. In this study, we investigated the relationship between fNIRS global signal and electroencephalogram (EEG) vigilance using simultaneous recordings in resting healthy subjects in high-density and whole-head montage. In Experiment 1, data were acquired at supine, sitting, and standing positions. Results found that the factor of body positions significantly affected the amplitude of the resting-state fNIRS global signal, prominently in the frequency range of 0.05-0.1 Hz but not in the very low frequency range of less than 0.05 Hz. As a control, the task-induced fNIRS or EEG responses to auditory stimuli did not differ across body positions. However, EEG vigilance plays a modulatory role in the fNIRS signals in the frequency range of less than 0.05 Hz: resting-state sessions of low EEG vigilance measures are associated with high amplitudes of fNIRS global signals. Moreover, in Experiment 2, we further examined the epoch-to-epoch fluctuations in concurrent fNIRS and EEG data acquired from a separate group of subjects and found a negative temporal correlation between EEG vigilance measures and fNIRS global signal amplitudes. Our study for the first time revealed that vigilance as a neurophysiological factor modulates the resting-state dynamics of fNIRS, which have important implications for understanding and processing the noises in fNIRS signals.

Prevalent and sex-biased breathing patterns modify functional connectivity MRI in young adults

Resting state functional connectivity magnetic resonance imaging (fMRI) is a tool for investigating human brain organization. Here we identify, visually and algorithmically, two prevalent influences on fMRI signals during 440 h of resting state scans in 440 healthy young adults, both caused by deviations from normal breathing which we term deep breaths and bursts. The two respiratory patterns have distinct influences on fMRI signals and signal covariance, distinct timescales, distinct cardiovascular correlates, and distinct tendencies to manifest by sex. Deep breaths are not sex-biased. Bursts, which are serial taperings of respiratory depth typically spanning minutes at a time, are more common in males. Bursts share features of chemoreflex-driven clinical breathing patterns that also occur primarily in males, with notable neurological, psychiatric, medical, and lifespan associations. These results identify common breathing patterns in healthy young adults with distinct influences on functional connectivity and an ability to differentially influence resting state fMRI studies.

Functional connectivity measured from fMRI data is widely used in neuroscience. Here the authors report an association between two types of breathing signature and obtained BOLD data, and associated sex differences.

Evidence for modulation of EEG microstate sequence by vigilance level

The momentary global functional state of the brain is reflected in its electric field configuration and cluster analytical approaches have consistently shown four configurations, referred to as EEG microstate classes A to D. Changes in microstate parameters are associated with a number of neuropsychiatric disorders, task performance, and mental state establishing their relevance for cognition. However, the common practice to use eye-closed resting state data to assess the temporal dynamics of microstate parameters might induce systematic confounds related to vigilance levels. Here, we studied the dynamics of microstate parameters in two independent data sets and showed that the parameters of microstates are strongly associated with vigilance level assessed both by EEG power analysis and fMRI global signal. We found that the duration and contribution of microstate class C, as well as transition probabilities towards microstate class C were positively associated with vigilance, whereas the sign was reversed for microstate classes A and B. Furthermore, in looking for the origins of the correspondence between microstates and vigilance level, we found Granger-causal effects of vigilance levels on microstate sequence parameters. Collectively, our findings suggest that duration and occurrence of microstates have a different origin and possibly reflect different physiological processes. Finally, our findings indicate the need for taking vigilance levels into consideration in resting-sate EEG investigations.

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.

Loss of frontal regulator of vigilance during sleep inertia: A simultaneous EEG-fMRI study

Sleep inertia refers to a distinct physiological state of waking up from sleep accompanied by performance impairments and sleepiness. The neural substrates of sleep inertia are unknown, but growing evidence suggests that this inertia state maintains certain sleep features. To investigate the neurophysiological mechanisms of sleep inertia, a comparison of pre-sleep and post-sleep wakefulness with eyes-open resting-state was performed using simultaneous EEG-fMRI, which has the potential to reveal the dynamic details of neuroelectric and hemodynamic responses with high temporal resolution. Our data suggested sleep-like features of slow EEG power and decreased BOLD activity were persistent during sleep inertia. In the pre-sleep phase, participants with stronger EEG vigilance showed stronger activity in the fronto-parietal network (FPN), but this phenomenon disappeared during sleep inertia. A time course analysis confirmed a decreased correlation between EEG vigilance and the FPN activity during sleep inertia. This simultaneous EEG-fMRI study advanced our understanding of sleep inertia and revealed the importance of the FPN in maintaining awareness. This is the first study to reveal the dynamic brain network changes from multi-modalities perspective during sleep inertia.

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.

Opening or closing eyes at rest modulates the functional connectivity of V1 with default and salience networks

Current evidence suggests that volitional opening or closing of the eyes modulates brain activity and connectivity. However, how the eye state influences the functional connectivity of the primary visual cortex has been poorly investigated. Using the same scanner, fMRI data from two groups of participants similar in age, sex and educational level were acquired. One group (n = 105) performed a resting state with eyes closed, and the other group (n = 63) performed a resting state with eyes open. Seed-based voxel-wise functional connectivity whole-brain analyses were performed to study differences in the connectivity of the primary visual cortex. This region showed higher connectivity with the default mode and sensorimotor networks in the eyes closed group, but higher connectivity with the salience network in the eyes open group. All these findings were replicated using an open source shared dataset. These results suggest that opening or closing the eyes may set brain functional connectivity in an interoceptive or exteroceptive state.

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.

Evaluating global brain connectivity as an imaging marker for depression: influence of preprocessing strategies and placebo-controlled ketamine treatment

Major depressive disorder (MDD) is associated with altered global brain connectivity (GBC), as assessed via resting-state functional magnetic resonance imaging (rsfMRI). Previous studies found that antidepressant treatment with ketamine normalized aberrant GBC changes in the prefrontal and cingulate cortices, warranting further investigations of GBC as a putative imaging marker. These results were obtained via global signal regression (GSR). This study is an independent replication of that analysis using a separate dataset. GBC was analyzed in 28 individuals with MDD and 22 healthy controls (HCs) at baseline, post-placebo, and post-ketamine. To investigate the effects of preprocessing, three distinct pipelines were used: (1) regression of white matter (WM)/cerebrospinal fluid (CSF) signals only (BASE); (2) WM/CSF + GSR (GSR); and (3) WM/CSF + physiological parameter regression (PHYSIO). Reduced GBC was observed in individuals with MDD only at baseline in the anterior and medial cingulate cortices, as well as in the prefrontal cortex only after regressing the global signal. Ketamine had no effect compared to baseline or placebo in either group in any pipeline. PHYSIO did not resemble GBC preprocessed with GSR. These results concur with several studies that used GSR to study GBC. Further investigations are warranted into disease-specific components of global fMRI signals that may drive these results and of GBCr as a potential imaging marker in MDD.

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.

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.

Assessing Brain Networks by Resting-State Dynamic Functional Connectivity: An fNIRS-EEG Study

The coordination of brain activity between disparate neural populations is highly dynamic. Investigations into intrinsic brain organization by evaluating dynamic resting-state functional connectivity (dRSFC) have attracted great attention in recent years. However, there are few dRSFC studies based on functional near-infrared spectroscopy (fNIRS) even though it has some advantages for studying the temporal evolution of brain function. In this research, we recruited 20 young adults and measured their resting-state brain fluctuations in several areas of the frontal, parietal, temporal, and occipital lobes using fNIRS-electroencephalography (EEG) simultaneous recording. Based on a sliding-window approach, we found that the variability of the dRSFC within any region of interest was significantly lower than the connections between region of interests but noticeably greater than the correlation between the channels with a short interoptode distance, which mainly consist of physiological fluctuations occurring in the superficial layers. Furthermore, based on a time-resolved k-means clustering analysis, the temporal evolution was extracted for three dominant functional networks. These networks were roughly consistent between different subject subgroups and in varying sliding time window lengths of 20, 30, and 60 s. Between these three functional networks, there were obvious time-varied and system-specific synchronous relationships. In addition, the oscillation of the frontal-parietal-temporal network showed significant correlation with the switching of one EEG microstate, a finding which is consistent with a previous functional MRI-EEG study. All this evidence implies the functional significance of fNIRS-dRSFC and demonstrates the feasibility of fNIRS for extracting the dominant functional networks based on RSFC dynamics.

Percent amplitude of fluctuation: A simple measure for resting-state fMRI signal at single voxel level

The amplitude of low-frequency fluctuation (ALFF) measures resting-state functional magnetic resonance imaging (RS-fMRI) signal of each voxel. However, the unit of blood oxygenation level-dependent (BOLD) signal is arbitrary and hence ALFF is sensitive to the scale of raw signal. A well-accepted standardization procedure is to divide each voxel's ALFF by the global mean ALFF, named mALFF. Although fractional ALFF (fALFF), a ratio of the ALFF to the total amplitude within the full frequency band, offers possible solution of the standardization, it actually mixes with the fluctuation power within the full frequency band and thus cannot reveal the true amplitude characteristics of a given frequency band. The current study borrowed the percent signal change in task fMRI studies and proposed percent amplitude of fluctuation (PerAF) for RS-fMRI. We firstly applied PerAF and mPerAF (i.e., divided by global mean PerAF) to eyes open (EO) vs. eyes closed (EC) RS-fMRI data. PerAF and mPerAF yielded prominently difference between EO and EC, being well consistent with previous studies. We secondly performed test-retest reliability analysis and found that (PerAF ≈ mPerAF ≈ mALFF) > (fALFF ≈ mfALFF). Head motion regression (Friston-24) increased the reliability of PerAF, but decreased all other metrics (e.g. mPerAF, mALFF, fALFF, and mfALFF). The above results suggest that mPerAF is a valid, more reliable, more straightforward, and hence a promising metric for voxel-level RS-fMRI studies. Future study could use both PerAF and mPerAF metrics. For prompting future application of PerAF, we implemented PerAF in a new version of REST package named RESTplus.

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.

Unstable wakefulness during resting-state fMRI and its associations with network connectivity and affective psychopathology in young adults

BACKGROUND

Drifts between wakefulness and sleep are common during resting state functional MRI (rsfMRI). Among healthy adults, within-scanner sleep can impact functional connectivity of default mode (DMN), task-positive (TPN), and thalamo-cortical networks. Because dysfunctional arousal states (i.e., sleepiness, sleep disturbance) are common in affective disorders, individuals with affective psychopathology may be more prone to unstable wakefulness during rsfMRI, hampering the estimation of clinically meaningful functional connectivity biomarkers.

METHODS

A transdiagnostic sample of 150 young adults (68 psychologically distressed; 82 psychiatrically healthy) completed rsfMRI and reported whether they experienced within-scanner sleep. Symptom scales were reduced into depression/anxiety and mania proneness dimensions using principal component analysis. We evaluated associations between within-scanner sleep, clinical status, and functional connectivity of the DMN, TPN, and thalamus.

RESULTS

Within-scanner sleep during rsfMRI was reported by 44% of participants (n = 66) but was unrelated to psychiatric diagnoses or mood symptom severity (p-values > 0.05). Across all participants, self-reported within-scanner sleep was associated with connectivity signatures akin to objectively-assessed sleep, including lower within-DMN connectivity, lower DMN-TPN anti-correlation, and altered thalamo-cortical connectivity (p < 0.05, corrected). Among participants reporting sustained wakefulness (n = 84), depression/anxiety severity positively associated with averaged DMN-TPN connectivity and mania proneness negatively associated with averaged thalamus-DMN connectivity (p-values < 0.05). Both relationships were attenuated and became non-significant when participants reporting within-scanner sleep were included (p-values > 0.05).

LIMITATIONS

Subjective report of within-scanner sleep.

CONCLUSIONS

Findings implicate within-scanner sleep as a source of variance in network connectivity; careful monitoring and correction for within-scanner sleep may enhance our ability to characterize network signatures underlying affective psychopathology.

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.

Decreased static and increased dynamic global signal topography in major depressive disorder

Major depressive disorder (MDD) has been linked to imbalanced communication among large-scale brain networks. However, the details of altered large-scale coordination of MDD remains unknown. To explore the altered large-scale functional organization in MDD. We used static and dynamic global signal (GS) topography, which are data-driven methods to explore altered relationship between global and local neuronal activities in MDD. Sixty three MDD patients and matched 63 healthy controls (HCs) were recruited in current study. Patients with MDD presented decreased static GS topography in bilateral parahippocampal gyrus and hippocampus gyrus. Meanwhile, patients with MDD presented increased variability of dynamic GS topography in the right ventromedial prefrontal cortex. This result may reflect the decreased and unstable whole brain functional coherence in MDD. The decreased static GS topography in the right parahippocampal gyrus was correlated with psychomotor retardation in patients with MDD. Our results presented that the altered static and dynamic GS topography can provide distinct evidence on the physiological mechanisms of MDD.

Resting State fMRI: Going Through the Motions

Resting state functional magnetic resonance imaging (rs-fMRI) has become an indispensable tool in neuroscience research. Despite this, rs-fMRI signals are easily contaminated by artifacts arising from movement of the head during data collection. The artifacts can be problematic even for motions on the millimeter scale, with complex spatiotemporal properties that can lead to substantial errors in functional connectivity estimates. Effective correction methods must be employed, therefore, to distinguish true functional networks from motion-related noise. Research over the last three decades has produced numerous correction methods, many of which must be applied in combination to achieve satisfactory data quality. Subject instruction, training, and mild restraints are helpful at the outset, but usually insufficient. Improvements come from applying multiple motion correction algorithms retrospectively after rs-fMRI data are collected, although residual artifacts can still remain in cases of elevated motion, which are especially prevalent in patient populations. Although not commonly adopted at present, “real-time” correction methods are emerging that can be combined with retrospective methods and that promise better correction and increased rs-fMRI signal sensitivity. While the search for the ideal motion correction protocol continues, rs-fMRI research will benefit from good disclosure practices, such as: (1) reporting motion-related quality control metrics to provide better comparison between studies; and (2) including motion covariates in group-level analyses to limit the extent of motion-related confounds when studying group differences.

Nuisance effects in inter-scan functional connectivity estimates before and after nuisance regression

In resting-state functional MRI, the correlation between blood-oxygenation-level-dependent (BOLD) signals across brain regions is used to estimate the functional connectivity (FC) of the brain. FC estimates are prone to the influence of nuisance factors including scanner-related artifacts and physiological modulations of the BOLD signal. Nuisance regression is widely performed to reduce the effect of nuisance factors on FC estimates on a per-scan basis. However, a dedicated analysis of nuisance effects on the variability of FC metrics across a collection of scans has been lacking. This work investigates the effects of nuisance factors on the variability of FC estimates across a collection of scans both before and after nuisance regression. Inter-scan variations in FC estimates are shown to be significantly correlated with the geometric norms of various nuisance terms, including head motion measurements, signals derived from white-matter and cerebrospinal regions, and the whole-brain global signal (GS) both before and after nuisance regression. In addition, it is shown that GS regression (GSR) can introduce GS norm-related fluctuations that are negatively correlated with inter-scan FC estimates. The empirical results are shown to be largely consistent with the predictions of a theoretical framework previously developed for the characterization of dynamic FC measures. This work shows that caution must be exercised when interpreting inter-scan FC measures across scans both before and after nuisance regression.

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.

The Not-So-Global Blood Oxygen Level-Dependent Signal

Global signal regression is a controversial processing step for resting-state functional magnetic resonance imaging, partly because the source of the global blood oxygen level-dependent (BOLD) signal remains unclear. On the one hand, nuisance factors such as motion can readily introduce coherent BOLD changes across the whole brain. On the other hand, the global signal has been linked to neural activity and vigilance levels, suggesting that it contains important neurophysiological information and should not be discarded. Any widespread pattern of coordinated activity is likely to contribute appreciably to the global signal. Such patterns may include large-scale quasiperiodic spatiotemporal patterns, known also to be tied to performance on vigilance tasks. This uncertainty surrounding the separability of the global BOLD signal from concurrent neurological processes motivated an examination of the global BOLD signal's spatial distribution. The results clarify that although the global signal collects information from all tissue classes, a diverse subset of the BOLD signal's independent components contribute the most to the global signal. Further, the timing of each network's contribution to the global signal is not consistent across volunteers, confirming the independence of a constituent process that comprises the global signal.

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.

Spectral entropy indicates electrophysiological and hemodynamic changes in drug-resistant epilepsy - A multimodal MREG study

Objective

Epilepsy causes measurable irregularity over a range of brain signal frequencies, as well as autonomic nervous system functions that modulate heart and respiratory rate variability. Imaging dynamic neuronal signals utilizing simultaneously acquired ultra-fast 10 Hz magnetic resonance encephalography (MREG), direct current electroencephalography (DC-EEG), and near-infrared spectroscopy (NIRS) can provide a more comprehensive picture of human brain function. Spectral entropy (SE) is a nonlinear method to summarize signal power irregularity over measured frequencies. SE was used as a joint measure to study whether spectral signal irregularity over a range of brain signal frequencies based on synchronous multimodal brain signals could provide new insights in the neural underpinnings of epileptiform activity.

Methods

Ten patients with focal drug-resistant epilepsy (DRE) and ten healthy controls (HC) were scanned with 10 Hz MREG sequence in combination with EEG, NIRS (measuring oxygenated, deoxygenated, and total hemoglobin: HbO, Hb, and HbT, respectively), and cardiorespiratory signals. After pre-processing, voxelwise SE_MREG was estimated from MREG data. Different neurophysiological and physiological subfrequency band signals were further estimated from MREG, DC-EEG, and NIRS: fullband (0–5 Hz, FB), near FB (0.08–5 Hz, NFB), brain pulsations in very-low (0.009–0.08 Hz, VLFP), respiratory (0.12–0.4 Hz, RFP), and cardiac (0.7–1.6 Hz, CFP) frequency bands. Global dynamic fluctuations in MREG and NIRS were analyzed in windows of 2 min with 50% overlap.

Results

Right thalamus, cingulate gyrus, inferior frontal gyrus, and frontal pole showed significantly higher SE_MREG in DRE patients compared to HC. In DRE patients, SE of cortical Hb was significantly reduced in FB (p = .045), NFB (p = .017), and CFP (p = .038), while both HbO and HbT were significantly reduced in RFP (p = .038, p = .045, respectively). Dynamic SE of HbT was reduced in DRE patients in RFP during minutes 2 to 6. Fitting to the frontal MREG and NIRS results, DRE patients showed a significant increase in SE_EEG in FB in fronto-central and parieto-occipital regions, in VLFP in parieto-central region, accompanied with a significant decrease in RFP in frontal pole and parietal and occipital (O2, Oz) regions.

Conclusion

This is the first study to show altered spectral entropy from synchronous MREG, EEG, and NIRS in DRE patients. Higher SE_MREG in DRE patients in anterior cingulate gyrus together with SE_EEG and SE_NIRS results in 0.12–0.4 Hz can be linked to altered parasympathetic function and respiratory pulsations in the brain. Higher SE_MREG in thalamus in DRE patients is connected to disturbances in anatomical and functional connections in epilepsy. Findings suggest that spectral irregularity of both electrophysiological and hemodynamic signals are altered in specific way depending on the physiological frequency range.

•

Simultaneous imaging methods indicate spectral irregularity in neurovascular and electrophysiological brain pulsations in DRE.

•

Altered spectral entropy in EEG, NIRS and BOLD indicate dysfunctional brain pulsations in respiratory frequency in epilepsy.

•

Spectral irregularity (0-5 Hz) of BOLD in right thalamus supports previous structural and functional findings in epilepsy.

Global fluctuations of cerebral blood flow indicate a global brain network independent of systemic factors

Global synchronization across specialized brain networks is a common feature of network models and in-vivo electrical measurements. Although the imaging of specialized brain networks with blood oxygenation sensitive resting state functional magnetic resonance imaging (rsfMRI) has enabled detailed study of regional networks, the study of globally correlated fluctuations with rsfMRI is confounded by spurious contributions to the global signal from systemic physiologic factors and other noise sources. Here we use an alternative rsfMRI method, arterial spin labeled perfusion MRI, to characterize global correlations and their relationship to correlations and anti-correlations between regional networks. Global fluctuations that cannot be explained by systemic factors dominate the fluctuations in cerebral blood flow. Power spectra of these fluctuations are band limited to below 0.05 Hz, similar to prior measurements of regional network fluctuations in the brain. Removal of these global fluctuations prior to measurement of regional networks reduces all regional network fluctuation amplitudes to below the global fluctuation amplitude and changes the strength and sign of inter network correlations. Our findings support large amplitude, globally synchronized activity across networks that require a reassessment of regional network amplitude and correlation measures.

Sleep as a model to understand neuroplasticity and recovery after stroke: Observational, perturbational and interventional approaches

Our own experiences with disturbances to sleep demonstrate its crucial role in the recovery of cognitive functions. This importance is likely enhanced in the recovery from stroke; both in terms of its physiology and cognitive abilities. Decades of experimental research have highlighted which aspects and mechanisms of sleep are likely to underlie these forms of recovery. Conversely, damage to certain areas of the brain, as well as the indirect effects of stroke, may disrupt sleep. However, only limited research has been conducted which seeks to directly explore this bidirectional link between both the macro and micro-architecture of sleep and stroke. Here we describe a series of semi-independent approaches that aim to establish this link through observational, perturbational, and interventional experiments. Our primary aim is to describe the methodology for future clinical and translational research needed to delineate competing accounts of the current data. At the observational level we suggest the use of high-density EEG recording, combined analysis of macro and micro-architecture of sleep, detailed analysis of the stroke lesion, and sensitive measures of functional recovery. The perturbational approach attempts to find the causal links between sleep and stroke. We promote the use of transcranial magnetic stimulation combined with EEG to examine the cortical dynamics of the peri-infarct stroke area. Translational research should take this a step further using optogenetic techniques targeting more specific cell populations. The interventional approach focuses on how the same clinical and translational perturbational techniques can be adapted to influence long-term recovery of function.

Variance of the global signal as a pretreatment predictor of antidepressant treatment response in drug-naïve major depressive disorder

Several behavioral and neuroimaging markers could be used to predict eventual antidepressant medication (ADM) outcomes in patients with major depressive disorder (MDD). However, these predictors are either subjective or complex, which has limited their clinical use. Thus, we aimed to identify an objective and easy-to-get marker to predict early therapeutic efficacy. Forty-seven drug-naïve patients with MDD and 47 age-, gender- and education-matched healthy controls underwent resting-state functional magnetic resonance imaging (fMRI) scans. We calculated the variable coefficient (VC) of the global signal for each subject. Baseline Hamilton Rating Scale for Depression (HRSD) score and that after 2 weeks of ADM were assessed for patients. Although there was no difference in VC between patients with MDD and healthy controls, we found a significant positive correlation between the VC and the decline rate of HRSD scores in the patients. Compared with the non-responding depression (NRD) group, the treatment-responsive depression (TRD) group had a higher VC. Receiver operator characteristic curve analysis revealed that the VC exhibited a good ability to differentiate TRD from NRD. In addition, the linear and logistic regression analyses showed that the VC was a significant predictor of the decline rate of HRSD scores and the antidepressant treatment response. These findings suggest that variance of the global signal may serve as a useful marker to help clinicians find an appropriate drug for individuals with MDD at the earliest opportunity and then further to facilitate personalized therapy.

Eyes-Open and Eyes-Closed Resting States With Opposite Brain Activity in Sensorimotor and Occipital Regions: Multidimensional Evidences From Machine Learning Perspective

Studies have demonstrated that there are widespread significant differences in spontaneous brain activity between eyes-open (EO) and eyes-closed (EC) resting states. However, it remains largely unclear whether spontaneous brain activity is effectively related to EO and EC resting states. The amplitude, local functional concordance, inter-hemisphere functional synchronization, and network centrality of spontaneous brain activity were measured by the fraction amplitude of low frequency fluctuation (fALFF), regional homogeneity (ReHo), voxel-mirrored homotopic connectivity (VMHC) and degree centrality (DC), respectively. Using the public Eyes-open/Eyes-closed dataset, we employed the support vector machine (SVM) and bootstrap technique to establish linking models for the fALFF, ReHo, VMHC and DC dimensions. The classification accuracies of linking models are 0.72 (0.59, 0.82), 0.88 (0.79, 0.97), 0.82 (0.74, 0.91) and 0.70 (0.62, 0.79), respectively. Specifically, we observed that brain activity in the EO condition is significantly greater in attentional system areas, including the fusiform gyrus, occipital and parietal cortex, but significantly lower in sensorimotor system areas, including the precentral/postcentral gyrus, paracentral lobule (PCL) and temporal cortex compared to the EC condition from the four dimensions. The results consistently indicated that spontaneous brain activity is effectively related to EO and EC resting states, and the two resting states are of opposite brain activity in sensorimotor and occipital regions. It may provide new insight into the neural substrate of the resting state and help computational neuroscientists or neuropsychologists to choose an appropriate resting state condition to investigate various mental disorders from the resting state functional magnetic resonance imaging (fMRI) technique.