Validation of real-time fMRI neurofeedback procedure for cognitive training using counterbalanced active-sham study design

Investigation of neural mechanisms of real-time functional MRI neurofeedback (rtfMRI-nf) training requires an efficient study control approach. A common rtfMRI-nf study design involves an experimental group, receiving active rtfMRI-nf, and a control group, provided with sham rtfMRI-nf. We report the first study in which rtfMRI-nf procedure is controlled through counterbalancing training runs with active and sham rtfMRI-nf for each participant. Healthy volunteers (n = 18) used rtfMRI-nf to upregulate fMRI activity of an individually defined target region in the left dorsolateral prefrontal cortex (DLPFC) while performing tasks that involved mental generation of a random numerical sequence and serial summation of numbers in the sequence. Sham rtfMRI-nf was provided based on fMRI activity of a different brain region, not involved in these tasks. The experimental procedure included two training runs with the active rtfMRI-nf and two runs with the sham rtfMRI-nf, in a randomized order. The participants achieved significantly higher fMRI activation of the left DLPFC target region during the active rtfMRI-nf conditions compared to the sham rtfMRI-nf conditions. fMRI functional connectivity of the left DLPFC target region with the nodes of the central executive network was significantly enhanced during the active rtfMRI-nf conditions relative to the sham conditions. fMRI connectivity of the target region with the nodes of the default mode network was similarly enhanced. fMRI connectivity changes between the active and sham conditions exhibited meaningful associations with individual performance measures on the Working Memory Multimodal Attention Task, the Approach-Avoidance Task, and the Trail Making Test. Our results demonstrate that the counterbalanced active-sham study design can be efficiently used to investigate mechanisms of active rtfMRI-nf in direct comparison to those of sham rtfMRI-nf. Further studies with larger group sizes are needed to confirm the reported findings and evaluate clinical utility of this study control approach.

Dynamic Functional Connectivity in Pediatric Mild Traumatic Brain Injury

Resting-state fMRI can be used to identify recurrent oscillatory patterns of functional connectivity within the human brain, also known as dynamic brain states. Alterations in dynamic brain states are highly likely to occur following pediatric mild traumatic brain injury (pmTBI) due to the active developmental changes. The current study used resting-state fMRI to investigate dynamic brain states in 200 patients with pmTBI (ages 8-18 years, median = 14 years) at the subacute (∼1-week post-injury) and early chronic (∼ 4 months post-injury) stages, and in 179 age- and sex-matched healthy controls (HC). A k-means clustering analysis was applied to the dominant time-varying phase coherence patterns to obtain dynamic brain states. In addition, correlations between brain signals were computed as measures of static functional connectivity. Dynamic connectivity analyses showed that patients with pmTBI spend less time in a frontotemporal default mode/limbic brain state, with no evidence of change as a function of recovery post-injury. Consistent with models showing traumatic strain convergence in deep grey matter and midline regions, static interhemispheric connectivity was affected between the left and right precuneus and thalamus, and between the right supplementary motor area and contralateral cerebellum. Changes in static or dynamic connectivity were not related to symptom burden or injury severity measures, such as loss of consciousness and post-traumatic amnesia. In aggregate, our study shows that brain dynamics are altered up to 4 months after pmTBI, in brain areas that are known to be vulnerable to TBI. Future longitudinal studies are warranted to examine the significance of our findings in terms of long-term neurodevelopment.

Multifaceted neural and vascular pathologies after pediatric mild traumatic brain injury

Dynamic changes in neurodevelopment and cognitive functioning occur during adolescence, including a switch from reactive to more proactive forms of cognitive control, including response inhibition. Pediatric mild traumatic brain injury (pmTBI) affects these cognitions immediately post-injury, but the role of vascular versus neural injury in cognitive dysfunction remains debated. This study consecutively recruited 214 sub-acute pmTBI (8-18 years) and age/sex-matched healthy controls (HC; N = 186), with high retention rates (>80%) at four months post-injury. Multimodal imaging (functional MRI during response inhibition, cerebral blood flow and cerebrovascular reactivity) assessed for pathologies within the neurovascular unit. Patients exhibited increased errors of commission and hypoactivation of motor circuitry during processing of probes. Evidence of increased/delayed cerebrovascular reactivity within motor circuitry during hypercapnia was present along with normal perfusion. Neither age-at-injury nor post-concussive symptom load were strongly associated with imaging abnormalities. Collectively, mild cognitive impairments and clinical symptoms may continue up to four months post-injury. Prolonged dysfunction within the neurovascular unit was observed during proactive response inhibition, with preliminary evidence that neural and pure vascular trauma are statistically independent. These findings suggest pmTBI is characterized by multifaceted pathologies during the sub-acute injury stage that persist several months post-injury.

Neural correlates of cognitive control deficits in pediatric mild traumatic brain injury

There is a growing body of research showing that cerebral pathophysiological processes triggered by pediatric mild traumatic brain injury (pmTBI) may extend beyond the usual clinical recovery timeline. It is paramount to further unravel these processes, because the possible long-term cognitive effects resulting from ongoing secondary injury in the developing brain are not known. In the current fMRI study, neural processes related to cognitive control were studied in 181 patients with pmTBI at sub-acute (SA; ~1 week) and early chronic (EC; ~4 months) stages post-injury. Additionally, a group of 162 age- and sex-matched healthy controls (HC) were recruited at equivalent time points. Proactive (post-cue) and reactive (post-probe) cognitive control were examined using a multimodal attention fMRI paradigm for either congruent or incongruent stimuli. To study brain network function, the triple-network model was used, consisting of the executive and salience networks (collectively known as the cognitive control network), and the default mode network. Additionally, whole-brain voxel-wise analyses were performed. Decreased deactivation was found within the default mode network at the EC stage following pmTBI during both proactive and reactive control. Voxel-wise analyses revealed sub-acute hypoactivation of a frontal area of the cognitive control network (left pre-supplementary motor area) during proactive control, with a reversed effect at the EC stage after pmTBI. Similar effects were observed in areas outside of the triple-network during reactive control. Group differences in activation during proactive control were limited to the visual domain, whereas for reactive control findings were more pronounced during the attendance of auditory stimuli. No significant correlations were present between task-related activations and (persistent) post-concussive symptoms. In aggregate, current results show alterations in neural functioning during cognitive control in pmTBI up to 4 months post-injury, regardless of clinical recovery. We propose that subacute decreases in activity reflect a general state of hypo-excitability due to the injury, while early chronic hyperactivation represents a compensatory mechanism to prevent default mode interference and to retain cognitive control.

The propensity to sign-track is associated with externalizing behavior and distinct patterns of reward-related brain activation in youth

Externalizing behaviors in childhood often predict impulse control disorders in adulthood; however, the underlying bio-behavioral risk factors are incompletely understood. In animals, the propensity to sign-track, or the degree to which incentive motivational value is attributed to reward cues, is associated with externalizing-type behaviors and deficits in executive control. Using a Pavlovian conditioned approach paradigm, we quantified sign-tracking in 40 healthy 9-12-year-olds. We also measured parent-reported externalizing behaviors and anticipatory neural activations to outcome-predicting cues using the monetary incentive delay fMRI task. Sign-tracking was associated with attentional and inhibitory control deficits and the degree of amygdala, but not cortical, activation during reward anticipation. These findings support the hypothesis that youth with a propensity to sign-track are prone to externalizing tendencies, with an over-reliance on subcortical cue-reactive brain systems. This research highlights sign-tracking as a promising experimental approach delineating the behavioral and neural circuitry of individuals at risk for externalizing disorders.

Prospective Study of Gray Matter Atrophy Following Pediatric Mild Traumatic Brain Injury

BACKGROUND AND OBJECTIVES

The clinical and physiologic time course for recovery following pediatric mild traumatic brain injury (pmTBI) remains actively debated. The primary objective of the current study was to prospectively examine structural brain changes (cortical thickness and subcortical volumes) and age-at-injury effects. A priori study hypotheses predicted reduced cortical thickness and hippocampal volumes up to 4 months postinjury, which would be inversely associated with age at injury.

METHODS

Prospective cohort study design with consecutive recruitment. Study inclusion adapted from American Congress of Rehabilitation Medicine (upper threshold) and Zurich Concussion in Sport Group (minimal threshold) and diagnosed by Emergency Department and Urgent Care clinicians. Major neurologic, psychiatric, or developmental disorders were exclusionary. Clinical (Common Data Element) and structural (3 T MRI) evaluations within 11 days (subacute visit [SA]) and at 4 months (early chronic visit [EC]) postinjury. Age- and sex-matched healthy controls (HC) to control for repeat testing/neurodevelopment. Clinical outcomes based on self-report and cognitive testing. Structural images quantified with FreeSurfer (version 7.1.1).

RESULTS

A total of 208 patients with pmTBI (age = 14.4 ± 2.9; 40.4% female) and 176 HC (age = 14.2 ± 2.9; 42.0% female) were included in the final analyses (>80% retention). Reduced cortical thickness (right rostral middle frontal gyrus; d = -0.49) and hippocampal volumes (d = -0.24) observed for pmTBI, but not associated with age at injury. Hippocampal volume recovery was mediated by loss of consciousness/posttraumatic amnesia. Significantly greater postconcussive symptoms and cognitive deficits were observed at SA and EC visits, but were not associated with the structural abnormalities. Structural abnormalities slightly improved balanced classification accuracy above and beyond clinical gold standards (∆+3.9%), with a greater increase in specificity (∆+7.5%) relative to sensitivity (∆+0.3%).

DISCUSSION

Current findings indicate that structural brain abnormalities may persist up to 4 months post-pmTBI and are partially mediated by initial markers of injury severity. These results contribute to a growing body of evidence suggesting prolonged physiologic recovery post-pmTBI. In contrast, there was no evidence for age-at-injury effects or physiologic correlates of persistent symptoms in our sample.

Machine Learning Evidence for Sex Differences Consistently Influences Resting-State Functional Magnetic Resonance Imaging Fluctuations Across Multiple Independently Acquired Data Sets

Background/Introduction: Sex classification using functional connectivity from resting-state functional magnetic resonance imaging (rs-fMRI) has shown promising results. This suggested that sex difference might also be embedded in the blood-oxygen-level-dependent properties such as the amplitude of low-frequency fluctuation (ALFF) and the fraction of ALFF (fALFF). This study comprehensively investigates sex differences using a reliable and explainable machine learning (ML) pipeline. Five independent cohorts of rs-fMRI with over than 5500 samples were used to assess sex classification performance and map the spatial distribution of the important brain regions. Methods: Five rs-fMRI samples were used to extract ALFF and fALFF features from predefined brain parcellations and then were fed into an unbiased and explainable ML pipeline with a wide range of methods. The pipeline comprehensively assessed unbiased performance for within-sample and across-sample validation. In addition, the parcellation effect, classifier selection, scanning length, spatial distribution, reproducibility, and feature importance were analyzed and evaluated thoroughly in the study. Results: The results demonstrated high sex classification accuracies from healthy adults (area under the curve >0.89), while degrading for nonhealthy subjects. Sex classification showed moderate to good intraclass correlation coefficient based on parcellation. Linear classifiers outperform nonlinear classifiers. Sex differences could be detected even with a short rs-fMRI scan (e.g., 2 min). The spatial distribution of important features overlaps with previous results from studies. Discussion: Sex differences are consistent in rs-fMRI and should be considered seriously in any study design, analysis, or interpretation. Features that discriminate males and females were found to be distributed across several different brain regions, suggesting a complex mosaic for sex differences in rs-fMRI. Impact statement The presented study unraveled that sex differences are embedded in the blood-oxygen-level dependent (BOLD) and can be predicted using unbiased and explainable machine learning pipeline. The study revealed that psychiatric disorders and demographics might influence the BOLD signal and interact with the classification of sex. The spatial distribution of the important features presented here supports the notion that the brain is a mosaic of male and female features. The findings emphasize the importance of controlling for sex when conducting brain imaging analysis. In addition, the presented framework can be adapted to classify other variables from resting-state BOLD signals.

Canonical EEG microstates transitions reflect switching among BOLD resting state networks and predict fMRI signal

Objective.Electroencephalography (EEG) microstates (MSs), which reflect a large topographical representation of coherent electrophysiological brain activity, are widely adopted to study cognitive processes mechanisms and aberrant alterations in brain disorders. MS topographies are quasi-stable lasting between 60-120 ms. Some evidence suggests that MS are the electrophysiological signature of resting-state networks (RSNs). However, the spatial and functional interpretation of MS and their association with functional magnetic resonance imaging (fMRI) remains unclear.Approach. In a cohort of healthy subjects (n= 52), we conducted several statistical and machine learning (ML) approaches analyses on the association among MS spatio-temporal dynamics and the blood-oxygenation-level dependent (BOLD) simultaneous EEG-fMRI data using statistical and ML approaches.Main results.Our results using a generalized linear model showed that MS transitions were largely and negatively associated with BOLD signals in the somatomotor, visual, dorsal attention, and ventral attention fMRI networks with limited association within the default mode network. Additionally, a novel recurrent neural network (RNN) confirmed the association between MS transitioning and fMRI signal while revealing that MS dynamics can model BOLD signals and vice versa.Significance.Results suggest that MS transitions may represent the deactivation of fMRI RSNs and provide evidence that both modalities measure common aspects of undergoing brain neuronal activities. These results may help to better understand the electrophysiological interpretation of MS.

Automated pipeline for EEG artifact reduction (APPEAR) recorded during fMRI

Objective.Simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) recordings offer a high spatiotemporal resolution approach to study human brain and understand the underlying mechanisms mediating cognitive and behavioral processes. However, the high susceptibility of EEG to MRI-induced artifacts hinders a broad adaptation of this approach. More specifically, EEG data collected during fMRI acquisition are contaminated with MRI gradients and ballistocardiogram artifacts, in addition to artifacts of physiological origin. There have been several attempts for reducing these artifacts with manual and time-consuming pre-processing, which may result in biasing EEG data due to variations in selecting steps order, parameters, and classification of artifactual independent components. Thus, there is a strong urge to develop a fully automatic and comprehensive pipeline for reducing all major EEG artifacts. In this work, we introduced an open-access toolbox with a fully automatic pipeline for reducing artifacts from EEG data collected simultaneously with fMRI (refer to APPEAR).Approach.The pipeline integrates average template subtraction and independent component analysis to suppress both MRI-related and physiological artifacts. To validate our results, we tested APPEAR on EEG data recorded from healthy control subjects during resting-state (n= 48) and task-based (i.e. event-related-potentials (ERPs);n= 8) paradigms. The chosen gold standard is an expert manual review of the EEG database.Main results.We compared manually and automated corrected EEG data during resting-state using frequency analysis and continuous wavelet transformation and found no significant differences between the two corrections. A comparison between ERP data recorded during a so-called stop-signal task (e.g. amplitude measures and signal-to-noise ratio) also showed no differences between the manually and fully automatic fMRI-EEG-corrected data.Significance.APPEAR offers the first comprehensive open-source toolbox that can speed up advancement of EEG analysis and enhance replication by avoiding experimenters' preferences while allowing for processing large EEG-fMRI cohorts composed of hundreds of subjects with manageable researcher time and effort.

Hippocampal volume recovery with real-time functional MRI amygdala neurofeedback emotional training for posttraumatic stress disorder

BACKGROUND

Small hippocampal volume is a prevalent neurostructural abnormality in posttraumatic stress disorder (PTSD). However, whether the hippocampal atrophy is the cause of disease symptoms or a pre-existing risk factor and whether it is a reversible alteration or a permanent trait are unclear. The trait- or state-dependent alteration could also differ among the hippocampal subfields.

METHODS

The study examined the longitudinal hippocampal volume changes due to positive emotional training with left amygdala (LA) real-time fMRI neurofeedback (rtfMRI-nf) in combat veterans with PTSD. The participants were trained to increase the neurofeedback signal from LA (experimental group, N = 20) or brain region not involved in emotion processing (control group, N = 9) by recalling a positive autobiographical memory. The pre- and post-training structural MRI brain images were processed with FreeSurfer to evaluate the hippocampal subfield volumes. Hippocampal volumes for healthy controls (N = 43) were also examined to evaluate the baseline abnormality in PTSD.

RESULTS

A significant group difference in volume change was found in the left CA1 head region. This region had the most significant volume reduction at the baseline in PTSD. The experimental group showed a significant volume increase, while the control group showed a significant volume decrease in this region. The volume change in the control group negatively correlated with interval days between the scans.

LIMITATIONS

A cognitive improvement due to the hippocampal volume increase could not be found with symptom scales.

CONCLUSIONS

RtfMRI-nf positive emotional training increased the hippocampus volume among people with PTSD, suggesting that hippocampal atrophy in PTSD is modifiable.

Real-time fMRI neurofeedback amygdala training may influence kynurenine pathway metabolism in major depressive disorder

•

rtfMRI-nf LA emotional training reduces depressive symptoms.

•

rtfMRI-nf LA training increases KynA/3-HK, a neuroprotective index.

•

Baseline KynA/QA is associated with the ability to upregulate the LA.

•

In partial responder group LA upregulation positively correlates with KynA/QA.

•

In partial responder group LA upregulation inversely correlates with MADRS.

•

Modulation of the KP may drive rtfMRI-nf-induced changes in neuroplasticity.

•

Non-specific effects cannot be ruled out due to the lack of a sham control.

Real-time fMRI neurofeedback (rtfMRI-nf) left amygdala (LA) training is a promising intervention for major depressive disorder (MDD). We have previously proposed that rtfMRI-nf LA training may reverse depression-associated regional impairments in neuroplasticity and restore information flow within emotion-regulating neural circuits. Inflammatory cytokines as well as the neuroactive metabolites of an immunoregulatory pathway, i.e. the kynurenine pathway (KP), have previously been implicated in neuroplasticity. Therefore, in this proof-of-principle study, we investigated the association between rtfMRI-nf LA training and circulating inflammatory mediators and KP metabolites. Based on our previous work, the primary variable of interest was the ratio of the NMDA-receptor antagonist, kynurenic acid to the NMDA receptor agonist, quinolinic acid (KynA/QA), a putative neuroprotective index. We tested two main hypotheses. i. Whether rtfMRI-nf acutely modulates KynA/QA, and ii. whether baseline KynA/QA predicts response to rtfMRI-nf.

Twenty-nine unmedicated participants who met DSM-5 criteria for MDD based on the Mini-International Neuropsychiatric Interview and had current depressive symptoms (Montgomery-Åsberg Depression Rating Scale (MADRS) score > 6) completed two rtfMRI-nf sessions to upregulate LA activity (Visit1 and 2), as well as a follow-up (Visit3) without rtfMRI-nf. All visits occurred at two-week intervals. At all three visits, the MADRS was administered to participants and serum samples for the quantification of inflammatory cytokines and KP metabolites were obtained. First, the longitudinal changes in the MADRS score and immune markers were tested by linear mixed effect model analysis. Further, utilizing a linear regression model, we investigated the relationship between rtfMRI-nf performance and immune markers. After two sessions of rtfMRI-nf, MADRS scores were significantly reduced (t[58] = −4.07, p = 0.009, d = 0.56). Thirteen participants showed a ≥ 25% reduction in the MADRS score (the partial responder group). There was a significant effect of visit (F[2,58] = 3.17, p = 0.05) for the neuroprotective index, KynA to 3-hydroxykynurenine (3-HK), that was driven by a significant increase in KynA/3-HK between Visit1 and Visit3 (t[58] = 2.50, p = 0.03, d = 0.38). A higher baseline level of KynA/QA (β = 5.23, p = 0.06; rho = 0.49, p = 0.02) was associated with greater ability to upregulate the LA. Finally, for exploratory purposes correlation analyses were performed between the partial responder and the non-responder groups as well as in the whole sample including all KP metabolites and cytokines. In the partial responder group, greater ability to upregulate the LA was correlated with an increase in KynA/QA after rtfMRI-nf (rho = 0.75, p = 0.03). The results are consistent with the possibility that rtfMRI-nf decreases metabolism down the so-called neurotoxic branch of the KP. Nevertheless, non-specific effects cannot be ruled out due to the lack of a sham control. Future, controlled studies are needed to determine whether the increase in KynA/3HK and KynA/QA is specific to rtfMRI-nf or whether it is a non-specific correlate of the resolution of depressive symptoms. Similarly, replication studies are needed to determine whether KynA/QA has clinical utility as a treatment response biomarker.

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

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

Effects of simultaneous real-time fMRI and EEG neurofeedback in major depressive disorder evaluated with brain electromagnetic tomography

Recently, we reported an emotion self-regulation study (Zotev et al., 2020), in which patients with major depressive disorder (MDD) used simultaneous real-time fMRI and EEG neurofeedback (rtfMRI-EEG-nf) to upregulate two fMRI and two EEG activity measures, relevant to MDD. The target measures included fMRI activities of the left amygdala and left rostral anterior cingulate cortex, and frontal EEG asymmetries in the alpha band (FAA) and high-beta band (FBA). Here we apply the exact low resolution brain electromagnetic tomography (eLORETA) to investigate EEG source activities during the rtfMRI-EEG-nf procedure. The exploratory analyses reveal significant changes in hemispheric lateralities of upper alpha and high-beta current source densities in the prefrontal regions, consistent with upregulation of the FAA and FBA during the rtfMRI-EEG-nf task. Similar laterality changes are observed for current source densities in the amygdala. Prefrontal upper alpha current density changes show significant negative correlations with anhedonia severity. Changes in prefrontal high-beta current density are consistent with reduction in comorbid anxiety. Comparisons with results of previous LORETA studies suggest that the rtfMRI-EEG-nf training is beneficial to MDD patients, and may have the ability to correct functional deficiencies associated with anhedonia and comorbid anxiety in MDD.

Emotion self-regulation training in major depressive disorder using simultaneous real-time fMRI and EEG neurofeedback

Simultaneous real-time fMRI and EEG neurofeedback (rtfMRI-EEG-nf) is an emerging neuromodulation approach, that enables simultaneous volitional regulation of both hemodynamic (BOLD fMRI) and electrophysiological (EEG) brain activities. Here we report the first application of rtfMRI-EEG-nf for emotion self-regulation training in patients with major depressive disorder (MDD). In this proof-of-concept study, MDD patients in the experimental group (n = 16) used rtfMRI-EEG-nf during a happy emotion induction task to simultaneously upregulate two fMRI and two EEG activity measures relevant to MDD. The target measures included BOLD activities of the left amygdala (LA) and left rostral anterior cingulate cortex (rACC), and frontal EEG asymmetries in the alpha band (FAA, [7.5-12.5] Hz) and high-beta band (FBA, [21-30] Hz). MDD patients in the control group (n = 8) were provided with sham feedback signals. An advanced procedure for improved real-time EEG-fMRI artifact correction was implemented. The experimental group participants demonstrated significant upregulation of the LA BOLD activity, FAA, and FBA during the rtfMRI-EEG-nf task, as well as significant enhancement in fMRI connectivity between the LA and left rACC. Average individual FAA changes during the rtfMRI-EEG-nf task positively correlated with depression and anhedonia severities, and negatively correlated with after-vs-before changes in depressed mood ratings. Temporal correlations between the FAA and FBA time courses and the LA BOLD activity were significantly enhanced during the rtfMRI-EEG-nf task. The experimental group participants reported significant mood improvements after the training. Our results suggest that the rtfMRI-EEG-nf may have potential for treatment of MDD.

Consensus on the reporting and experimental design of clinical and cognitive-behavioural neurofeedback studies (CRED-nf checklist)

Neurofeedback has begun to attract the attention and scrutiny of the scientific and medical mainstream. Here, neurofeedback researchers present a consensus-derived checklist that aims to improve the reporting and experimental design standards in the field.

Self-regulation of ventromedial prefrontal cortex activation using real-time fMRI neurofeedback-Influence of default mode network

The ventromedial prefrontal cortex (vmPFC) is involved in regulation of negative emotion and decision-making, emotional and behavioral control, and active resilient coping. This pilot study examined the feasibility of training healthy subjects (n = 27) to self-regulate the vmPFC activity using a real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf). Participants in the experimental group (EG, n = 18) were provided with an ongoing vmPFC hemodynamic activity (rtfMRI-nf signal represented as variable-height bar). Individuals were instructed to raise the bar by self-relevant value-based thinking. Participants in the control group (CG, n = 9) performed the same task; however, they were provided with computer-generated sham neurofeedback signal. Results demonstrate that (a) both the CG and the EG show a higher vmPFC fMRI signal at the baseline than during neurofeedback training; (b) no significant positive training effect was seen in the vmPFC across neurofeedback runs; however, the medial prefrontal cortex, middle temporal gyri, inferior frontal gyri, and precuneus showed significant decreasing trends across the training runs only for the EG; (c) the vmPFC rtfMRI-nf signal associated with the fMRI signal across the default mode network (DMN). These findings suggest that it may be difficult to modulate a single DMN region without affecting other DMN regions. Observed decreased vmPFC activity during the neurofeedback task could be due to interference from the fMRI signal within other DMN network regions, as well as interaction with task-positive networks. Even though participants in the EG did not show significant positive increase in the vmPFC activity among neurofeedback runs, they were able to learn to accommodate the demand of self-regulation task to maintain the vmPFC activity with the help of a neurofeedback signal.

Brain activity mediators of PTSD symptom reduction during real-time fMRI amygdala neurofeedback emotional training

Self-regulation of brain activation with real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf) is emerging as a promising treatment for psychiatric disorders. The association between the regulation and symptom reduction, however, has not been consistent, and the mechanisms underlying the symptom reduction remain poorly understood. The present study investigated brain activity mediators of the amygdala rtfMRI-nf training effect on combat veterans' PTSD symptom reduction. The training was designed to increase a neurofeedback signal either from the left amygdala (experimental group; EG) or from a control region not implicated in emotion regulation (control group; CG) during positive autobiographical memory recall. We employed a structural equation model mapping analysis to identify brain regions that mediated the effects of the rtfMRI-nf training on PTSD symptoms. Symptom reduction was mediated by low activation in the dorsomedial prefrontal cortex (DMPFC) and the middle cingulate cortex. There was a trend toward less activation in these regions for the EG compared to the CG. Low activation in the precuneus, the right superior parietal, the right insula, and the right cerebellum also mediated symptom reduction while their effects were moderated by the neurofeedback signal; a higher signal was linked to less effect on symptom reduction. This moderation was not specific to the EG. MDD comorbidity was associated with high DMPFC activation, which resulted in less effective regulation of the feedback signal. These results indicated that symptom reduction due to the neurofeedback training was not specifically mediated by the neurofeedback target activity, but broad regions were involved in the process.

EEG Microstates Temporal Dynamics Differentiate Individuals with Mood and Anxiety Disorders From Healthy Subjects

Electroencephalography (EEG) measures the brain’s electrophysiological spatio-temporal activities with high temporal resolution. Multichannel and broadband analysis of EEG signals is referred to as EEG microstates (EEG-ms) and can characterize such dynamic neuronal activity. EEG-ms have gained much attention due to the increasing evidence of their association with mental activities and large-scale brain networks identified by functional magnetic resonance imaging (fMRI). Spatially independent EEG-ms are quasi-stationary topographies (e.g., stable, lasting a few dozen milliseconds) typically classified into four canonical classes (microstates A through D). They can be identified by clustering EEG signals around EEG global field power (GFP) maxima points. We examined the EEG-ms properties and the dynamics of cohorts of mood and anxiety (MA) disorders subjects (n = 61) and healthy controls (HCs; n = 52). In both groups, we found four distinct classes of EEG-ms (A through D), which did not differ among cohorts. This suggests a lack of significant structural cortical abnormalities among cohorts, which would otherwise affect the EEG-ms topographies. However, both cohorts’ brain network dynamics significantly varied, as reflected in EEG-ms properties. Compared to HC, the MA cohort features a lower transition probability between EEG-ms B and D and higher transition probability from A to D and from B to C, with a trend towards significance in the average duration of microstate C. Furthermore, we harnessed a recently introduced theoretical approach to analyze the temporal dependencies in EEG-ms. The results revealed that the transition matrices of MA group exhibit higher symmetrical and stationarity properties as compared to HC ones. In addition, we found an elevation in the temporal dependencies among microstates, especially in microstate B for the MA group. The determined alteration in EEG-ms temporal dependencies among the cohorts suggests that brain abnormalities in mood and anxiety disorders reflect aberrant neural dynamics and a temporal dwelling among ceratin brain states (i.e., mood and anxiety disorders subjects have a less dynamicity in switching between different brain states).

Real-time fMRI amygdala neurofeedback positive emotional training normalized resting-state functional connectivity in combat veterans with and without PTSD: a connectome-wide investigation

Self-regulation of brain activation using real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf) is an emerging approach for treating mood and anxiety disorders. The effect of neurofeedback training on resting-state functional connectivity warrants investigation as changes in spontaneous brain activation could reflect the association between sustained symptom relief and brain alteration. We investigated the effect of amygdala-focused rtfMRI-nf training on resting-state functional connectivity in combat veterans with and without posttraumatic stress disorder (PTSD) who were trained to increase a feedback signal reflecting left amygdala activity while recalling positive autobiographical memories (Zotev et al., 2018). The analysis was performed in three stages: i) first, we investigated the connectivity in the left amygdala region; ii) next, we focused on the abnormal resting-state functional connectivity identified in our previous analysis of this data (Misaki et al., 2018); and iii) finally, we performed a novel data-driven longitudinal connectome-wide analysis. We introduced a longitudinal multivariate distance matrix regression (MDMR) analysis to comprehensively examine neurofeedback training effects beyond those associated with abnormal baseline connectivity.

These comprehensive exploratory analyses suggested that abnormal resting-state connectivity for combat veterans with PTSD was partly normalized after the training. This included hypoconnectivities between the left amygdala and the left ventrolateral prefrontal cortex (vlPFC) and between the supplementary motor area (SMA) and the dorsal anterior cingulate cortex (dACC). The increase of SMA-dACC connectivity was associated with PTSD symptom reduction. Longitudinal MDMR analysis found a connectivity change between the precuneus and the left superior frontal cortex. The connectivity increase was associated with a decrease in hyperarousal symptoms. The abnormal connectivity for combat veterans without PTSD - such as hypoconnectivity in the precuneus with a superior frontal region and hyperconnectivity in the posterior insula with several regions - could also be normalized after the training. These results suggested that the rtfMRI-nf training effect was not limited to a feedback target region and symptom relief could be mediated by brain modulation in several regions other than in a feedback target area. While further confirmatory research is needed, the results may provide valuable insight into treatment effects on the whole brain resting-state connectivity.

•

fMRI neurofeedback training effect on resting-state connectivity was examined

•

Left amygdala activity was trained to increase with positive memory

•

Neurofeedback normalized altered connectivity in veterans with and without PTSD

•

PTSD symptom reductions were significant but not specific to group (exp/ctrl)

•

Connectivity-symptom association was seen in mPFC and precuneus

Amygdala real-time functional magnetic resonance imaging neurofeedback for major depressive disorder: A review

Advances in imaging technologies have allowed for the analysis of functional magnetic resonance imaging data in real-time (rtfMRI), leading to the development of neurofeedback (nf) training. This rtfMRI-nf training utilizes functional magnetic resonance imaging (fMRI) tomographic localization capacity to allow a person to see and regulate the localized hemodynamic signal from his or her own brain. In this review, we summarize the results of several studies that have developed and applied neurofeedback training to healthy and depressed individuals with the amygdala as the neurofeedback target and the goal to increase the hemodynamic response during positive autobiographical memory recall. We review these studies and highlight some of the challenges and advances in developing an rtfMRI-nf paradigm for broader use in psychiatric populations. The work described focuses on our line of research aiming to develop the rtfMRI-nf into an intervention, and includes a discussion of the selection of a region of interest for feedback, selecting a control condition, behavioral and cognitive effects of training, and predicting which participants are most likely to respond well to training. While the results of these studies are encouraging and suggest the clinical potential of amygdala rtfMRI-nf in alleviating symptoms of major depressive disorder, larger studies are warranted to confirm its efficacy.

Tracking resting state connectivity dynamics in veterans with PTSD

Posttraumatic stress disorder (PTSD) is a trauma- and stressor-related disorder that may emerge following a traumatic event. Neuroimaging studies have shown evidence of functional abnormality in many brain regions and systems affected by PTSD. Exaggerated threat detection associated with abnormalities in the salience network, as well as abnormalities in executive functions involved in emotions regulations, self-referencing and context evaluation processing are broadly reported in PTSD. Here we aimed to investigate the behavior and dynamic properties of fMRI resting state networks in combat-related PTSD, using a novel, multimodal imaging approach. Simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) was employed to measure neurobiological brain activity among 36 veterans with combat-related PTSD and 20 combat-exposed veterans without PTSD. Based on the recently established method of measuring temporal-independent EEG microstates, we developed a novel strategy to integrate EEG and fMRI by quantifying the fast temporal dynamics associated with the resting state networks. We found distinctive occurrence rates of microstates associated with the dorsal default mode network and salience networks in the PTSD group as compared with control. Furthermore, the occurrence rate of the microstate for the dorsal default mode network was positively correlated with PTSD severity, whereas the occurrence rate of the microstate for the anterior salience network was negatively correlated with hedonic tone reported by participants with PTSD. Our findings reveal a novel aspect of abnormal network dynamics in combat-related PTSD and contribute to a better understanding of the pathophysiology of the disorder. Simultaneous EEG and fMRI will be a valuable tool in continuing to study the neurobiology underlying PTSD.

•

Concurrent EEG-fMRI study of resting brain activity in combat related PTSD.

•

EEG-microstates were associated with fMRI resting state networks in PTSD.

•

PTSD associated with alterations in dorsal default mode and salience networks.

•

Occurrence rates of EEG-microstates were related to PTSD symptoms.

Real-time fMRI neurofeedback training of the amygdala activity with simultaneous EEG in veterans with combat-related PTSD

Posttraumatic stress disorder (PTSD) is a chronic and disabling neuropsychiatric disorder characterized by insufficient top-down modulation of the amygdala activity by the prefrontal cortex. Real-time fMRI neurofeedback (rtfMRI-nf) is an emerging method with potential for modifying the amygdala-prefrontal interactions. We report the first controlled emotion self-regulation study in veterans with combat-related PTSD utilizing rtfMRI-nf of the amygdala activity. PTSD patients in the experimental group (EG, n = 20) learned to upregulate blood‑oxygenation-level-dependent (BOLD) activity of the left amygdala (LA) using the rtfMRI-nf during a happy emotion induction task. PTSD patients in the control group (CG, n = 11) were provided with a sham rtfMRI-nf. The study included three rtfMRI-nf training sessions, and EEG recordings were performed simultaneously with fMRI. PTSD severity was assessed before and after the training using the Clinician-Administered PTSD Scale (CAPS). The EG participants who completed the study showed a significant reduction in total CAPS ratings, including significant reductions in avoidance and hyperarousal symptoms. They also exhibited a significant reduction in comorbid depression severity. Overall, 80% of the EG participants demonstrated clinically meaningful reductions in CAPS ratings, compared to 38% in the CG. No significant difference in the CAPS rating changes was observed between the groups. During the first rtfMRI-nf session, functional connectivity of the LA with the orbitofrontal cortex (OFC) and the dorsolateral prefrontal cortex (DLPFC) was progressively enhanced, and this enhancement significantly and positively correlated with the initial CAPS ratings. Left-lateralized enhancement in upper alpha EEG coherence also exhibited a significant positive correlation with the initial CAPS. Reduction in PTSD severity between the first and last rtfMRI-nf sessions significantly correlated with enhancement in functional connectivity between the LA and the left DLPFC. Our results demonstrate that the rtfMRI-nf of the amygdala activity has the potential to correct the amygdala-prefrontal functional connectivity deficiencies specific to PTSD.

Altered task-based and resting-state amygdala functional connectivity following real-time fMRI amygdala neurofeedback training in major depressive disorder

Background

We have previously shown that in participants with major depressive disorder (MDD) trained to upregulate their amygdala hemodynamic response during positive autobiographical memory (AM) recall with real-time fMRI neurofeedback (rtfMRI-nf) training, depressive symptoms diminish. Here, we assessed the effect of rtfMRI-nf on amygdala functional connectivity during both positive AM recall and rest.

Method

The current manuscript consists of a secondary analysis on data from our published clinical trial of neurofeedback. Patients with MDD completed two rtfMRI-nf sessions (18 received amygdala rtfMRI-nf, 16 received control parietal rtfMRI-nf). One-week prior-to and following training participants also completed a resting-state fMRI scan. A GLM-based functional connectivity analysis was applied using a seed ROI in the left amygdala. We compared amygdala functional connectivity changes while recalling positive AMs from the baseline run to the final transfer run during rtfMRI-nf training, as well during rest from the baseline to the one-week follow-up visit. Finally, we assessed the correlation between change in depression scores and change in amygdala connectivity, as well as correlations between amygdala regulation success and connectivity changes.

Results

Following training, amygdala connectivity during positive AM recall increased with widespread regions in the frontal and limbic network. During rest, amygdala connectivity increased following training within the fronto-temporal-limbic network. During both task and resting-state analyses, amygdala-temporal pole connectivity decreased. We identified increased amygdala-precuneus and amygdala-inferior frontal gyrus connectivity during positive memory recall and increased amygdala-precuneus and amygdala-thalamus connectivity during rest as functional connectivity changes that explained significant variance in symptom improvement. Amygdala-precuneus connectivity changes also explain a significant amount of variance in neurofeedback regulation success.

Conclusions

Neurofeedback training to increase amygdala hemodynamic activity during positive AM recall increased amygdala connectivity with regions involved in self-referential, salience, and reward processing. Results suggest future targets for neurofeedback interventions, particularly interventions involving the precuneus.

•

Changes in amygdala functional connectivity following neurofeedback were examined.

•

Amygdala rtfMRI-nf training alters functional connectivity with prefrontal regions.

•

Increased amygdala-precuneus connectivity may underlie clinical improvements.

Real-time fMRI neurofeedback of the mediodorsal and anterior thalamus enhances correlation between thalamic BOLD activity and alpha EEG rhythm

Real-time fMRI neurofeedback (rtfMRI-nf) with simultaneous EEG allows volitional modulation of BOLD activity of target brain regions and investigation of related electrophysiological activity. We applied this approach to study correlations between thalamic BOLD activity and alpha EEG rhythm. Healthy volunteers in the experimental group (EG, n = 15) learned to upregulate BOLD activity of the target region consisting of the mediodorsal (MD) and anterior (AN) thalamic nuclei using rtfMRI-nf during retrieval of happy autobiographical memories. Healthy subjects in the control group (CG, n = 14) were provided with a sham feedback. The EG participants were able to significantly increase BOLD activities of the MD and AN. Functional connectivity between the MD and the inferior precuneus was significantly enhanced during the rtfMRI-nf task. Average individual changes in the occipital alpha EEG power significantly correlated with the average MD BOLD activity levels for the EG. Temporal correlations between the occipital alpha EEG power and BOLD activities of the MD and AN were significantly enhanced, during the rtfMRI-nf task, for the EG compared to the CG. Temporal correlations with the alpha power were also significantly enhanced for the posterior nodes of the default mode network, including the precuneus/posterior cingulate, and for the dorsal striatum. Our findings suggest that the temporal correlation between the MD BOLD activity and posterior alpha EEG power is modulated by the interaction between the MD and the inferior precuneus, reflected in their functional connectivity. Our results demonstrate the potential of the rtfMRI-nf with simultaneous EEG for noninvasive neuromodulation studies of human brain function.

Connectome-wide investigation of altered resting-state functional connectivity in war veterans with and without posttraumatic stress disorder

Altered resting-state functional connectivity in posttraumatic stress disorder (PTSD) suggests neuropathology of the disorder. While seed-based fMRI connectivity analysis is often used for the studies, such analysis requires defining a seed location a priori, which restricts search scope and could bias findings toward presupposed areas. Recently, a comprehensive exploratory voxel-wise connectivity analysis, the connectome-wide association approach, has been introduced using multivariate distance matrix regression (MDMR) for resting-state functional connectivity analysis. The current study performed a connectome-wide investigation of resting-state functional connectivity for war veterans with and without PTSD compared to non-trauma-exposed healthy controls using MDMR. Thirty-five male combat veterans with PTSD (unmedicated), 18 male combat veterans without PTSD (veterans control, VC), and 28 age-matched non-trauma-exposed healthy males (NC) participated in a resting-state fMRI scan. MDMR analysis was used to identify between-groups differences in regions with altered connectivity. The identified regions were used as a seed for post-hoc functional connectivity analysis. The analysis revealed that PTSD patients had hypoconnectivity between the left lateral prefrontal regions and the salience network regions as well as hypoconnectivity between the parahippocampal gyrus and the visual cortex areas. Connectivity between the ventromedial prefrontal cortex and the middle frontal gyrus and between the parahippocampal gyrus and the anterior insula were negatively correlated with PTSD symptom severity. VC subjects also had altered functional connectivity compared to NC, including increased connectivity between the posterior insula and several brain regions and decreased connectivity between the precuneus region and several other brain areas. The decreased connectivity between the lateral prefrontal regions and the salience network regions in PTSD was consistent with previous reports that indicated lowered emotion-regulation function in these regions. The decreased connectivity between the parahippocampal gyrus and visual cortex supported the dual representation theory of PTSD, which suggests dissociation between sensory and contextual memory representations in PTSD. The theory also supposes that the precuneus is a region that triggers retrieval of sensory memory of traumatic events. The decreased connectivity at the precuneus for VC might be associated with suppressing such a process.

Real-Time Functional Magnetic Resonance Imaging Amygdala Neurofeedback Changes Positive Information Processing in Major Depressive Disorder

BACKGROUND

In participants with major depressive disorder who are trained to upregulate their amygdalar hemodynamic responses during positive autobiographical memory recall with real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf) training, depressive symptoms diminish. This study tested whether amygdalar rtfMRI-nf also changes emotional processing of positive and negative stimuli in a variety of behavioral and imaging tasks.

METHODS

Patients with major depressive disorder completed two rtfMRI-nf sessions (18 received amygdalar rtfMRI-nf, 16 received control parietal rtfMRI-nf). One week before and following rtfMRI-nf training, participants performed tasks measuring responses to emotionally valenced stimuli including a backward-masking task, which measures the amygdalar hemodynamic response to emotional faces presented for traditionally subliminal duration and followed by a mask, and the Emotional Test Battery in which reaction times and performance accuracy are measured during tasks involving emotional faces and words.

RESULTS

During the backward-masking task, amygdalar responses increased while viewing masked happy faces but decreased to masked sad faces in the experimental versus control group following rtfMRI-nf. During the Emotional Test Battery, reaction times decreased to identification of positive faces and during self-identification with positive words and vigilance scores increased to positive faces and decreased to negative faces during the faces dot-probe task in the experimental versus control group following rtfMRI-nf.

CONCLUSIONS

rtfMRI-nf training to increase the amygdalar hemodynamic response to positive memories was associated with changes in amygdalar responses to happy and sad faces and improved processing of positive stimuli during performance of the Emotional Test Battery. These results may suggest that amygdalar rtfMRI-nf training alters responses to emotional stimuli in a manner similar to antidepressant pharmacotherapy.

Randomized Clinical Trial of Real-Time fMRI Amygdala Neurofeedback for Major Depressive Disorder: Effects on Symptoms and Autobiographical Memory Recall

OBJECTIVE

Patients with depression show blunted amygdala hemodynamic activity to positive stimuli, including autobiographical memories. The authors examined the therapeutic efficacy of real-time functional MRI neurofeedback (rtfMRI-nf) training aimed at increasing the amygdala's hemodynamic response to positive memories in patients with depression.

METHOD

In a double-blind, placebo-controlled, randomized clinical trial, unmedicated adults with depression (N=36) were randomly assigned to receive two sessions of rtfMRI-nf either from the amygdala (N=19) or from a parietal control region not involved in emotional processing (N=17). Clinical scores and autobiographical memory performance were assessed at baseline and 1 week after the final rtfMRI-nf session. The primary outcome measure was change in score on the Montgomery-Åsberg Depression Rating Scale (MADRS), and the main analytic approach consisted of a linear mixed-model analysis.

RESULTS

In participants in the experimental group, the hemodynamic response in the amygdala increased relative to their own baseline and to the control group. Twelve participants in the amygdala rtfMRI-nf group, compared with only two in the control group, had a >50% decrease in MADRS score. Six participants in the experimental group, compared with one in the control group, met conventional criteria for remission at study end, resulting in a number needed to treat of 4. In participants receiving amygdala rtfMRI-nf, the percent of positive specific memories recalled increased relative to baseline and to the control group.

CONCLUSIONS

rtfMRI-nf training to increase the amygdala hemodynamic response to positive memories significantly decreased depressive symptoms and increased the percent of specific memories recalled on an autobiographical memory test. These data support a role of the amygdala in recovery from depression.

Correlation between amygdala BOLD activity and frontal EEG asymmetry during real-time fMRI neurofeedback training in patients with depression

Real-time fMRI neurofeedback (rtfMRI-nf) is an emerging approach for studies and novel treatments of major depressive disorder (MDD). EEG performed simultaneously with an rtfMRI-nf procedure allows an independent evaluation of rtfMRI-nf brain modulation effects. Frontal EEG asymmetry in the alpha band is a widely used measure of emotion and motivation that shows profound changes in depression. However, it has never been directly related to simultaneously acquired fMRI data. We report the first study investigating electrophysiological correlates of the rtfMRI-nf procedure, by combining the rtfMRI-nf with simultaneous and passive EEG recordings. In this pilot study, MDD patients in the experimental group (n = 13) learned to upregulate BOLD activity of the left amygdala using an rtfMRI-nf during a happy emotion induction task. MDD patients in the control group (n = 11) were provided with a sham rtfMRI-nf. Correlations between frontal EEG asymmetry in the upper alpha band and BOLD activity across the brain were examined. Average individual changes in frontal EEG asymmetry during the rtfMRI-nf task for the experimental group showed a significant positive correlation with the MDD patients' depression severity ratings, consistent with an inverse correlation between the depression severity and frontal EEG asymmetry at rest. The average asymmetry changes also significantly correlated with the amygdala BOLD laterality. Temporal correlations between frontal EEG asymmetry and BOLD activity were significantly enhanced, during the rtfMRI-nf task, for the amygdala and many regions associated with emotion regulation. Our findings demonstrate an important link between amygdala BOLD activity and frontal EEG asymmetry during emotion regulation. Our EEG asymmetry results indicate that the rtfMRI-nf training targeting the amygdala is beneficial to MDD patients. They further suggest that EEG-nf based on frontal EEG asymmetry in the alpha band would be compatible with the amygdala-based rtfMRI-nf. Combination of the two could enhance emotion regulation training and benefit MDD patients.

Automatic EEG-assisted retrospective motion correction for fMRI (aE-REMCOR)

Head motions during functional magnetic resonance imaging (fMRI) impair fMRI data quality and introduce systematic artifacts that can affect interpretation of fMRI results. Electroencephalography (EEG) recordings performed simultaneously with fMRI provide high-temporal-resolution information about ongoing brain activity as well as head movements. Recently, an EEG-assisted retrospective motion correction (E-REMCOR) method was introduced. E-REMCOR utilizes EEG motion artifacts to correct the effects of head movements in simultaneously acquired fMRI data on a slice-by-slice basis. While E-REMCOR is an efficient motion correction approach, it involves an independent component analysis (ICA) of the EEG data and identification of motion-related ICs. Here we report an automated implementation of E-REMCOR, referred to as aE-REMCOR, which we developed to facilitate the application of E-REMCOR in large-scale EEG-fMRI studies. The aE-REMCOR algorithm, implemented in MATLAB, enables an automated preprocessing of the EEG data, an ICA decomposition, and, importantly, an automatic identification of motion-related ICs. aE-REMCOR has been used to perform retrospective motion correction for 305 fMRI datasets from 16 subjects, who participated in EEG-fMRI experiments conducted on a 3T MRI scanner. Performance of aE-REMCOR has been evaluated based on improvement in temporal signal-to-noise ratio (TSNR) of the fMRI data, as well as correction efficiency defined in terms of spike reduction in fMRI motion parameters. The results show that aE-REMCOR is capable of substantially reducing head motion artifacts in fMRI data. In particular, when there are significant rapid head movements during the scan, a large TSNR improvement and high correction efficiency can be achieved. Depending on a subject's motion, an average TSNR improvement over the brain upon the application of aE-REMCOR can be as high as 27%, with top ten percent of the TSNR improvement values exceeding 55%. The average correction efficiency over the 305 fMRI scans is 18% and the largest achieved efficiency is 71%. The utility of aE-REMCOR on the resting state fMRI connectivity of the default mode network is also examined. The motion-induced position-dependent error in the DMN connectivity analysis is shown to be reduced when aE-REMCOR is utilized. These results demonstrate that aE-REMCOR can be conveniently and efficiently used to improve fMRI motion correction in large clinical EEG-fMRI studies.

Reconstructing Large-Scale Brain Resting-State Networks from High-Resolution EEG: Spatial and Temporal Comparisons with fMRI

Functional magnetic resonance imaging (fMRI) studies utilizing measures of hemodynamic signal, such as the blood oxygenation level-dependent (BOLD) signal, have discovered that resting-state brain activities are organized into multiple large-scale functional networks, coined as resting-state networks (RSNs). However, an important limitation of the available fMRI studies is that hemodynamic signals only provide an indirect measure of the neuronal activity. In contrast, electroencephalography (EEG) directly measures electrophysiological activity of the brain. However, little is known about the brain-wide organization of such spontaneous neuronal population signals at the resting state. It is not entirely clear if or how the network structure built upon slowly fluctuating hemodynamic signals is represented in terms of fast, dynamic, and spontaneous neuronal activity. In this study, we investigated the electrophysiological representation of RSNs from simultaneously acquired EEG and fMRI data in the resting human brain. We developed a data-driven analysis approach that reconstructed multiple large-scale electrophysiological networks from high-resolution EEG data alone. The networks derived from EEG were then compared with RSNs independently derived from simultaneously acquired fMRI in their spatial structures as well as temporal dynamics. Results reveal spatially and temporally specific electrophysiological correlates for the fMRI-RSNs. Findings suggest that the spontaneous activity of various large-scale cortical networks is reflected in macroscopic EEG potentials.

Resting-state functional connectivity modulation and sustained changes after real-time functional magnetic resonance imaging neurofeedback training in depression

Amygdala hemodynamic responses to positive stimuli are attenuated in major depressive disorder (MDD) and normalize with remission. Real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf) training with the goal of upregulating amygdala activity during recall of happy autobiographical memories (AMs) has been suggested, and recently explored, as a novel therapeutic approach that resulted in improvement in self-reported mood in depressed subjects. In this study, we assessed the possibility of sustained brain changes as well as the neuromodulatory effects of rtfMRI-nf training of the amygdala during recall of positive AMs in MDD and matched healthy subjects. MDD and healthy subjects went through one visit of rtfMRI-nf training. Subjects were assigned to receive active neurofeedback from the left amygdale (LA) or from a control region putatively not modulated by AM recall or emotion regulation, that is, the left horizontal segment of the intraparietal sulcus. To assess lasting effects of neurofeedback in MDD, the resting-state functional connectivity before and after rtfMRI-nf in 27 depressed subjects, as well as in 27 matched healthy subjects before rtfMRI-nf was measured. Results show that abnormal hypo-connectivity with LA in MDD is reversed after rtfMRI-nf training by recalling positive AMs. Although such neuromodulatory changes are observed in both MDD groups receiving feedback from respective active and control brain regions, only in the active group are larger decreases of depression severity associated with larger increases of amygdala connectivity and a significant, positive correlation is found between the connectivity changes and the days after neurofeedback. In addition, active neurofeedback training of the amygdala enhances connectivity with temporal cortical regions, including the hippocampus. These results demonstrate lasting brain changes induced by amygdala rtfMRI-nf training and suggest the importance of reinforcement learning in rehabilitating emotion regulation in depression.

Contrast enhancement by combining T1- and T2-weighted structural brain MR Images

PURPOSE

In order to more precisely differentiate cerebral structures in neuroimaging studies, a novel technique for enhancing the tissue contrast based on a combination of T1-weighted (T1w) and T2-weighted (T2w) MRI images was developed.

METHODS

The combined image (CI) was calculated as CI = (T1w - sT2w)/(T1w + sT2w), where sT2w is the scaled T2-weighted image. The scaling factor was calculated to adjust the gray- matter (GM) voxel intensities in the T2w image so that their median value equaled that of the GM voxel intensities in the T1w image. The image intensity homogeneity within a tissue and the discriminability between tissues in the CI versus the separate T1w and T2w images were evaluated using the segmentation by the FMRIB Software Library (FSL) and FreeSurfer (Athinoula A. Martinos Center for Biomedical Imaging at Massachusetts General Hospital, Boston, MA) software.

RESULTS

The combined image significantly improved homogeneity in the white matter (WM) and GM compared to the T1w images alone. The discriminability between WM and GM also improved significantly by applying the CI approach. Significant enhancements to the homogeneity and discriminability also were achieved in most subcortical nuclei tested, with the exception of the amygdala and the thalamus.

CONCLUSION

The tissue discriminability enhancement offered by the CI potentially enables more accurate neuromorphometric analyses of brain structures.

Real-time FMRI neurofeedback training of amygdala activity in patients with major depressive disorder

BACKGROUND

Amygdala hemodynamic responses to positive stimuli are attenuated in major depressive disorder (MDD), and normalize with remission. Real-time functional MRI neurofeedback (rtfMRI-nf) offers a non-invasive method to modulate this regional activity. We examined whether depressed participants can use rtfMRI-nf to enhance amygdala responses to positive autobiographical memories, and whether this ability alters symptom severity.

METHODS

Unmedicated MDD subjects were assigned to receive rtfMRI-nf from either left amygdala (LA; experimental group, n = 14) or the horizontal segment of the intraparietal sulcus (HIPS; control group, n = 7) and instructed to contemplate happy autobiographical memories (AMs) to raise the level of a bar representing the hemodynamic signal from the target region to a target level. This 40s Happy condition alternated with 40s blocks of rest and counting backwards. A final Transfer run without neurofeedback information was included.

RESULTS

Participants in the experimental group upregulated their amygdala responses during positive AM recall. Significant pre-post scan decreases in anxiety ratings and increases in happiness ratings were evident in the experimental versus control group. A whole brain analysis showed that during the transfer run, participants in the experimental group had increased activity compared to the control group in left superior temporal gyrus and temporal polar cortex, and right thalamus.

CONCLUSIONS

Using rtfMRI-nf from the left amygdala during recall of positive AMs, depressed subjects were able to self-regulate their amygdala response, resulting in improved mood. Results from this proof-of-concept study suggest that rtfMRI-nf training with positive AM recall holds potential as a novel therapeutic approach in the treatment of depression.

Correlated slow fluctuations in respiration, EEG, and BOLD fMRI

Low-frequency temporal fluctuations of physiological signals (<0.1 Hz), such as the respiration and cardiac pulse rate, occur naturally during rest and have been shown to be correlated with blood-oxygenation-level-dependent (BOLD) signal fluctuation. Such physiological signal modulations have been considered as sources of noise and their effects on BOLD signal are commonly removed in functional magnetic resonance imaging (fMRI) studies. However, possible neural correlates of the physiological fluctuations have not been considered nor examined in detail. In the present study we investigated this possibility by simultaneously acquiring electroencephalogram (EEG) with BOLD fMRI data, respiratory and cardiac waveforms in healthy human subjects at eyes-closed and eyes-open resting. We quantified the concurrent changes of the EEG power in the alpha frequency band, the respiration volume, and the cardiac pulse rate, then assessed the temporal correlations between alpha EEG power and physiological signal fluctuations. In addition, time-shifted time courses of alpha EEG power or physiological data were included as regressors to examine their correlations with the whole-brain BOLD fMRI signals. We observed a significant correlation between alpha EEG global field power and respiration, particularly at eyes-closed resting condition. Similar spatial patterns were observed between the correlation maps of BOLD with alpha EEG power and respiration, with negative correlations coinciding in the visual cortex, superior/middle temporal gyrus, inferior frontal gyrus, and inferior parietal lobule and positive correlations in the thalamus and caudate. Regressing out the physiological variations in the BOLD signal resulted in reduced correlation between BOLD and alpha EEG power. These results suggest a mutual link of neuronal origin between alpha EEG power, respiration, and BOLD signals.

Spatiotemporal dynamics of the brain at rest--exploring EEG microstates as electrophysiological signatures of BOLD resting state networks

Neuroimaging research suggests that the resting cerebral physiology is characterized by complex patterns of neuronal activity in widely distributed functional networks. As studied using functional magnetic resonance imaging (fMRI) of the blood-oxygenation-level dependent (BOLD) signal, the resting brain activity is associated with slowly fluctuating hemodynamic signals (~10s). More recently, multimodal functional imaging studies involving simultaneous acquisition of BOLD-fMRI and electroencephalography (EEG) data have suggested that the relatively slow hemodynamic fluctuations of some resting state networks (RSNs) evinced in the BOLD data are related to much faster (~100 ms) transient brain states reflected in EEG signals, that are referred to as "microstates". To further elucidate the relationship between microstates and RSNs, we developed a fully data-driven approach that combines information from simultaneously recorded, high-density EEG and BOLD-fMRI data. Using independent component analysis (ICA) of the combined EEG and fMRI data, we identified thirteen microstates and ten RSNs that are organized independently in their temporal and spatial characteristics, respectively. We hypothesized that the intrinsic brain networks that are active at rest would be reflected in both the EEG data and the fMRI data. To test this hypothesis, the rapid fluctuations associated with each microstate were correlated with the BOLD-fMRI signal associated with each RSN. We found that each RSN was characterized further by a specific electrophysiological signature involving from one to a combination of several microstates. Moreover, by comparing the time course of EEG microstates to that of the whole-brain BOLD signal, on a multi-subject group level, we unraveled for the first time a set of microstate-associated networks that correspond to a range of previously described RSNs, including visual, sensorimotor, auditory, attention, frontal, visceromotor and default mode networks. These results extend our understanding of the electrophysiological signature of BOLD RSNs and demonstrate the intrinsic connection between the fast neuronal activity and slow hemodynamic fluctuations.

Self-regulation of amygdala activation using real-time FMRI neurofeedback

Real-time functional magnetic resonance imaging (rtfMRI) with neurofeedback allows investigation of human brain neuroplastic changes that arise as subjects learn to modulate neurophysiological function using real-time feedback regarding their own hemodynamic responses to stimuli. We investigated the feasibility of training healthy humans to self-regulate the hemodynamic activity of the amygdala, which plays major roles in emotional processing. Participants in the experimental group were provided with ongoing information about the blood oxygen level dependent (BOLD) activity in the left amygdala (LA) and were instructed to raise the BOLD rtfMRI signal by contemplating positive autobiographical memories. A control group was assigned the same task but was instead provided with sham feedback from the left horizontal segment of the intraparietal sulcus (HIPS) region. In the LA, we found a significant BOLD signal increase due to rtfMRI neurofeedback training in the experimental group versus the control group. This effect persisted during the Transfer run without neurofeedback. For the individual subjects in the experimental group the training effect on the LA BOLD activity correlated inversely with scores on the Difficulty Identifying Feelings subscale of the Toronto Alexithymia Scale. The whole brain data analysis revealed significant differences for Happy Memories versus Rest condition between the experimental and control groups. Functional connectivity analysis of the amygdala network revealed significant widespread correlations in a fronto-temporo-limbic network. Additionally, we identified six regions — right medial frontal polar cortex, bilateral dorsomedial prefrontal cortex, left anterior cingulate cortex, and bilateral superior frontal gyrus — where the functional connectivity with the LA increased significantly across the rtfMRI neurofeedback runs and the Transfer run. The findings demonstrate that healthy subjects can learn to regulate their amygdala activation using rtfMRI neurofeedback, suggesting possible applications of rtfMRI neurofeedback training in the treatment of patients with neuropsychiatric disorders.

Magnetic Resonance Relaxometry at Low and Ultra low Fields

Nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) are ubiquitous tools in science and medicine. NMR provides powerful probes of local and macromolecular chemical structure and dynamics. Recently it has become possible and practical to perform MR at very low fields (from 1 μT to 1 mT), the so-called ultra-low field (ULF) regime. Pulsed pre-polarizing fields greatly enhance the signal strength and allow flexibility in signal acquisition sequences. Improvements in SQUID sensor technology allow ultra-sensitive detection in a pulsed field environment.In this regime the proton Larmor frequencies (1 Hz - 100 kHz) of ULF MR overlap (on a time scale of 10 μs to 100 ms) with "slow" molecular dynamic processes such as diffusion, intra-molecular motion, chemical reactions, and biological processes such as protein folding, catalysis and ligand binding. The frequency dependence of relaxation at ultra-low fields may provide a probe for biomolecular dynamics on the millisecond timescale (protein folding and aggregation, conformational motions of enzymes, binding and structural fluctuations of coupled domains in allosteric mechanisms) relevant to host-pathogen interactions, biofuels, and biomediation. Also this resonance-enhanced coupling at ULF can greatly enhance contrast in medical applications of ULF-MRI resulting in better diagnostic techniques.We have developed a number of instruments and techniques to study relaxation vs. frequency at the ULF regime. Details of the techniques and results are presented.Ultra-low field methods are already being applied at LANL in brain imaging, and detection of liquid explosives at airports. However, the potential power of ultra-low field MR remains to be fully exploited.