fMRI BOLD Correlates of EEG Independent Components: Spatial Correspondence With the Default Mode Network

Mindfulness-based Neurofeedback: A Systematic Review of EEG and fMRI studies

Neurofeedback concurrent with mindfulness meditation may reveal meditation effects on the brain and facilitate improved mental health outcomes. Here, we systematically reviewed EEG and fMRI studies of mindfulness meditation with neurofeedback (mbNF) and followed PRISMA guidelines. We identified 10 fMRI reports, consisting of 177 unique participants, and 9 EEG reports, consisting of 242 participants. Studies of fMRI focused primarily on downregulating the default-mode network (DMN). Although studies found decreases in DMN activations during neurofeedback, there is a lack of evidence for transfer effects, and the majority of studies did not employ adequate controls, e.g. sham neurofeedback. Accordingly, DMN decreases may have been confounded by general task-related deactivation. EEG studies typically examined alpha, gamma, and theta frequency bands, with the most robust evidence supporting the modulation of theta band activity. Both EEG and fMRI mbNF have been implemented with high fidelity in clinical populations. However, the mental health benefits of mbNF have not been established. In general, mbNF studies would benefit from sham-controlled RCTs, as well as clear reporting (e.g. CRED-NF).

Longitudinal network-based brain grey matter MRI measures are clinically relevant and sensitive to treatment effects in multiple sclerosis

In multiple sclerosis clinical trials, MRI outcome measures are typically extracted at a whole-brain level, but pathology is not homogeneous across the brain and so whole-brain measures may overlook regional treatment effects. Data-driven methods, such as independent component analysis, have shown promise in identifying regional disease effects but can only be computed at a group level and cannot be applied prospectively. The aim of this work was to develop a technique to extract longitudinal independent component analysis network-based measures of co-varying grey matter volumes, derived from T1-weighted volumetric MRI, in individual study participants, and assess their association with disability progression and treatment effects in clinical trials. We used longitudinal MRI and clinical data from 5089 participants (22 045 visits) with multiple sclerosis from eight clinical trials. We included people with relapsing-remitting, primary and secondary progressive multiple sclerosis. We used data from five negative clinical trials (2764 participants, 13 222 visits) to extract the independent component analysis-based measures. We then trained and cross-validated a least absolute shrinkage and selection operator regression model (which can be applied prospectively to previously unseen data) to predict the independent component analysis measures from the same regional MRI volume measures and applied it to data from three positive clinical trials (2325 participants, 8823 visits). We used nested mixed-effect models to determine how networks differ across multiple sclerosis phenotypes are associated with disability progression and to test sensitivity to treatment effects. We found 17 consistent patterns of co-varying regional volumes. In the training cohort, volume loss was faster in four networks in people with secondary progressive compared with relapsing-remitting multiple sclerosis and three networks with primary progressive multiple sclerosis. Volume changes were faster in secondary compared with primary progressive multiple sclerosis in four networks. In the combined positive trials cohort, eight independent component analysis networks and whole-brain grey matter volume measures showed treatment effects, and the magnitude of treatment-placebo differences in the network-based measures was consistently greater than with whole-brain grey matter volume measures. Longitudinal network-based analysis of grey matter volume changes is feasible using clinical trial data, showing differences cross-sectionally and longitudinally between multiple sclerosis phenotypes, associated with disability progression, and treatment effects. Future work is required to understand the pathological mechanisms underlying these regional changes.

Electrocortical theta activity may reflect sensory prediction errors during adaptation to a gradual gait perturbation

Locomotor adaptation to abrupt and gradual perturbations are likely driven by fundamentally different neural processes. The aim of this study was to quantify brain dynamics associated with gait adaptation to a gradually introduced gait perturbation, which typically results in smaller behavioral errors relative to an abrupt perturbation. Loss of balance during standing and walking elicits transient increases in midfrontal theta oscillations that have been shown to scale with perturbation intensity. We hypothesized there would be no significant change in anterior cingulate theta power (4-7 Hz) with respect to pre-adaptation when a gait perturbation is introduced gradually because the gradual perturbation acceleration and stepping kinematic errors are small relative to an abrupt perturbation. Using mobile electroencephalography (EEG), we measured gait-related spectral changes near the anterior cingulate, posterior cingulate, sensorimotor, and posterior parietal cortices as young, neurotypical adults (n = 30) adapted their gait to an incremental split-belt treadmill perturbation. Most cortical clusters we examined (>70%) did not exhibit changes in electrocortical activity between 2-50 Hz. However, we did observe gait-related theta synchronization near the left anterior cingulate cortex during strides with the largest errors, as measured by step length asymmetry. These results suggest gradual adaptation with small gait asymmetry and perturbation magnitude may not require significant cortical resources beyond normal treadmill walking. Nevertheless, the anterior cingulate may remain actively engaged in error monitoring, transmitting sensory prediction error information via theta oscillations.

Towards an electroencephalographic measure of awareness based on the reactivity of oscillatory macrostates to hearing a subject's own name

We proposed that the brain's electrical activity is composed of a sequence of alternating states with repeating topographic spectral distributions on scalp electroencephalogram (EEG), referred to as oscillatory macrostates. The macrostate showing the largest decrease in the probability of occurrence, measured as a percentage (reactivity), during sensory stimulation was labelled as the default EEG macrostate (DEM). This study aimed to assess the influence of awareness on DEM reactivity (DER). We included 11 middle cerebral artery ischaemic stroke patients with impaired awareness having a median Glasgow Coma Scale (GCS) of 6/15 and a group of 11 matched healthy controls. EEG recordings were carried out during auditory 1 min stimulation epochs repeating either the subject's own name (SON) or the SON in reverse (rSON). The DEM was identified across three SON epochs alternating with three rSON epochs. Compared with the patients, the DEM of controls contained more posterior theta activity reflecting source dipoles that could be mapped in the posterior cingulate cortex. The DER was measured from the 1 min quiet baseline preceding each stimulation epoch. The difference in mean DER between the SON and rSON epochs was measured by the salient EEG reactivity (SER) theoretically ranging from -100% to 100%. The SER was 12.4 ± 2.7% (Mean ± standard error of the mean) in controls and only 1.3 ± 1.9% in the patient group (P < 0.01). The patient SER decreased with the Glasgow Coma Scale. Our data suggest that awareness increases DER to SON as measured by SER.

Conflict- and error-related theta activities are coupled to BOLD signals in different brain regions

Both conflict and error processing have been linked to the midfrontal theta power (4-8 Hz) increase as indicated by EEG studies and greater hemodynamic activity in the anterior midcingulate cortex (aMCC) as indicated by fMRI studies. Conveniently, the source of the midfrontal theta power was estimated in or nearby aMCC. However, previous studies using concurrent EEG and fMRI recordings in resting-state or other cognitive tasks observed only a negative relationship between theta power and BOLD signal in the brain regions typically showing task-related deactivations. In this study, we used a simultaneous EEG-fMRI technique to investigate a trial-by-trial coupling between theta power and hemodynamic activity during the performance of two conflict tasks. Independent component analysis (ICA) was applied to denoise the EEG signal and select individual midfrontal EEG components, whereas group ICA was applied to fMRI data to obtain a functional parcellation of the frontal cortex. Using a linear mixed-effect model, theta power was coupled with the peak of hemodynamic responses from various frontal, cingulate, and insular cortical sites to unravel the potential brain sources that contribute to conflict- and error-related theta variability. Although several brain regions exhibited conflict-related increases in hemodynamic activity, the conflict pre-response theta showed only a negative correlation to BOLD signal in the midline area 9 (MA9), a region exhibiting conflict-sensitive deactivation. Conversely, and more expectedly, error-related theta showed a positive relationship to activity in the aMCC. Our results provide novel evidence suggesting that the amplitude of pre-response theta reflects the process of active inhibition that suppresses the MA9 activity. This process is affected independently by the stimulus congruency, reaction times variance, and is susceptible to the time-on-task effect. Finally, it predicts the commitment of an omission error. Together, our findings highlight that conflict- and error-related theta oscillations represent fundamentally different processes.

EEG and fMRI coupling and decoupling based on joint independent component analysis (jICA)

BACKGROUND

Meaningful integration of functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) requires knowing whether these measurements reflect the activity of the same neural sources, i.e., estimating the degree of coupling and decoupling between the neuroimaging modalities.

NEW METHOD

This paper proposes a method to quantify the coupling and decoupling of fMRI and EEG signals based on the mixing matrix produced by joint independent component analysis (jICA). The method is termed fMRI/EEG-jICA.

RESULTS

fMRI and EEG acquired during a syllable detection task with variable syllable presentation rates (0.25-3 Hz) were separated with jICA into two spatiotemporally distinct components, a primary component that increased nonlinearly in amplitude with syllable presentation rate, putatively reflecting an obligatory auditory response, and a secondary component that declined nonlinearly with syllable presentation rate, putatively reflecting an auditory attention orienting response. The two EEG subcomponents were of similar amplitude, but the secondary fMRI subcomponent was ten folds smaller than the primary one.

COMPARISON TO EXISTING METHOD

FMRI multiple regression analysis yielded a map more consistent with the primary than secondary fMRI subcomponent of jICA, as determined by a greater area under the curve (0.5 versus 0.38) in a sensitivity and specificity analysis of spatial overlap.

CONCLUSION

fMRI/EEG-jICA revealed spatiotemporally distinct brain networks with greater sensitivity than fMRI multiple regression analysis, demonstrating how this method can be used for leveraging EEG signals to inform the detection and functional characterization of fMRI signals. fMRI/EEG-jICA may be useful for studying neurovascular coupling at a macro-level, e.g., in neurovascular disorders.

A method to assess the default EEG macrostate and its reactivity to stimulation

OBJECTIVE

The default mode network (DMN) is deactivated by stimulation. We aimed to assess the DMN reactivity impairment by routine EEG recordings in stroke patients with impaired consciousness.

METHODS

Binocular light flashes were delivered at 1 Hz in 1-minute epochs, following a 1-minute baseline (PRE). The EEG was decomposed in a series of binary oscillatory macrostates by topographic spectral clustering. The most deactivated macrostate was labeled the default EEG macrostate (DEM). Its reactivity (DER) was quantified as the decrease in DEM occurrence probability during stimulation. A normalized DER index (DERI) was calculated as DER/PRE. The measures were compared between 14 healthy controls and 32 comatose patients under EEG monitoring following an acute stroke.

RESULTS

The DEM was mapped to the posterior DMN hubs. In the patients, these DEM source dipoles were 3-4 times less frequent and were associated with an increased theta activity. Even in a reduced 6-channel montage, a DER below 6.26% corresponding to a DERI below 0.25 could discriminate the patients with sensitivity and specificity well above 80%.

CONCLUSION

The method detected the DMN impairment in post-stroke coma patients.

SIGNIFICANCE

The DEM and its reactivity to stimulation could be useful to monitor the DMN function at bedside.

Stationary EEG pattern relates to large-scale resting state networks - An EEG-fMRI study connecting brain networks across time-scales

Relating brain dynamics acting on time scales that differ by at least an order of magnitude is a fundamental issue in brain research. The same is true for the observation of stable dynamical structures in otherwise highly non-stationary signals. The present study addresses both problems by the analysis of simultaneous resting state EEG-fMRI recordings of 53 patients with epilepsy. Confirming previous findings, we observe a generic and temporally stable average correlation pattern in EEG recordings. We design a predictor for the General Linear Model describing fluctuations around the stationary EEG correlation pattern and detect resting state networks in fMRI data. The acquired statistical maps are contrasted to several surrogate tests and compared with maps derived by spatial Independent Component Analysis of the fMRI data. By means of the proposed EEG-predictor we observe core nodes of known fMRI resting state networks with high specificity in the default mode, the executive control and the salience network. Our results suggest that both, the stationary EEG pattern as well as resting state fMRI networks are different expressions of the same brain activity. This activity is interpreted as the dynamics on (or close to) a stable attractor in phase space that is necessary to maintain the brain in an efficient operational mode. We discuss that this interpretation is congruent with the theoretical framework of complex systems as well as with the brain's energy balance.

Predicting disability progression and cognitive worsening in multiple sclerosis using patterns of grey matter volumes

OBJECTIVE

In multiple sclerosis (MS), MRI measures at the whole brain or regional level are only modestly associated with disability, while network-based measures are emerging as promising prognostic markers. We sought to demonstrate whether data-driven patterns of covarying regional grey matter (GM) volumes predict future disability in secondary progressive MS (SPMS).

METHODS

We used cross-sectional structural MRI, and baseline and longitudinal data of Expanded Disability Status Scale, Nine-Hole Peg Test (9HPT) and Symbol Digit Modalities Test (SDMT), from a clinical trial in 988 people with SPMS. We processed T1-weighted scans to obtain GM probability maps and applied spatial independent component analysis (ICA). We repeated ICA on 400 healthy controls. We used survival models to determine whether baseline patterns of covarying GM volume measures predict cognitive and motor worsening.

RESULTS

We identified 15 patterns of regionally covarying GM features. Compared with whole brain GM, deep GM and lesion volumes, some ICA components correlated more closely with clinical outcomes. A mainly basal ganglia component had the highest correlations at baseline with the SDMT and was associated with cognitive worsening (HR=1.29, 95% CI 1.09 to 1.52, p<0.005). Two ICA components were associated with 9HPT worsening (HR=1.30, 95% CI 1.06 to 1.60, p<0.01 and HR=1.21, 95% CI 1.01 to 1.45, p<0.05). ICA measures could better predict SDMT and 9HPT worsening (C-index=0.69-0.71) compared with models including only whole and regional MRI measures (C-index=0.65-0.69, p value for all comparison <0.05).

CONCLUSIONS

The disability progression was better predicted by some of the covarying GM regions patterns, than by single regional or whole-brain measures. ICA, which may represent structural brain networks, can be applied to clinical trials and may play a role in stratifying participants who have the most potential to show a treatment effect.

Brain variability in dynamic resting-state networks identified by fuzzy entropy: a scalp EEG study

Objective.Exploring the temporal variability in spatial topology during the resting state attracts growing interest and becomes increasingly useful to tackle the cognitive process of brain networks. In particular, the temporal brain dynamics during the resting state may be delineated and quantified aligning with cognitive performance, but few studies investigated the temporal variability in the electroencephalogram (EEG) network as well as its relationship with cognitive performance.Approach.In this study, we proposed an EEG-based protocol to measure the nonlinear complexity of the dynamic resting-state network by applying the fuzzy entropy. To further validate its applicability, the fuzzy entropy was applied into simulated and two independent datasets (i.e. decision-making and P300).Main results.The simulation study first proved that compared to the existing methods, this approach could not only exactly capture the pattern dynamics in time series but also overcame the magnitude effect of time series. Concerning the two EEG datasets, the flexible and robust network architectures of the brain cortex at rest were identified and distributed at the bilateral temporal lobe and frontal/occipital lobe, respectively, whose variability metrics were found to accurately classify different groups. Moreover, the temporal variability of resting-state network property was also either positively or negatively related to individual cognitive performance.Significance.This outcome suggested the potential of fuzzy entropy for evaluating the temporal variability of the dynamic resting-state brain networks, and the fuzzy entropy is also helpful for uncovering the fluctuating network variability that accounts for the individual decision differences.

Reduced Dynamic Interactions Within Intrinsic Functional Brain Networks in Early Blind Patients

Neuroimaging studies in early blind (EB) patients have shown altered connections or brain networks. However, it remains unclear how the causal relationships are disrupted within intrinsic brain networks. In our study, we used spectral dynamic causal modeling (DCM) to estimate the causal interactions using resting-state data in a group of 20 EB patients and 20 healthy controls (HC). Coupling parameters in specific regions were estimated, including the medial prefrontal cortex (mPFC), posterior cingulate cortex (PCC), and inferior parietal lobule (IPC) in the default mode network (DMN); dorsal anterior cingulate cortex (dACC) and bilateral anterior insulae (AI) in the salience network (SN), and bilateral frontal eye fields (FEF) and superior parietal lobes (SPL) within the dorsal attention network (DAN). Statistical analyses found that all endogenous connections and the connections from the mPFC to bilateral IPCs in EB patients were significantly reduced within the DMN, and the effective connectivity from the PCC and lIPC to the mPFC, and from the mPFC to the PCC were enhanced. For the SN, all significant connections in EB patients were significantly decreased, except the intrinsic right AI connections. Within the DAN, more significant effective connections were observed to be reduced between the EB and HC groups, while only the connections from the right SPL to the left SPL and the intrinsic connection in the left SPL were significantly enhanced. Furthermore, discovery of more decreased effective connections in the EB subjects suggested that the disrupted causal interactions between specific regions are responsive to the compensatory brain plasticity in early deprivation.