PHYCAA+: an optimized, adaptive procedure for measuring and controlling physiological noise in BOLD fMRI

Advances in functional and structural imaging of the brainstem: implications for disease

PURPOSE OF REVIEW

The brainstem's complex anatomy and relatively small size means that structural and functional assessment of this structure is done less frequently compared to other brain areas. However, recent years have seen substantial progress in brainstem imaging, enabling more detailed investigations into its structure and function, as well as its role in neuropathology.

RECENT FINDINGS

Advancements in ultrahigh field MRI technology have allowed for unprecedented spatial resolution in brainstem imaging, facilitating the new creation of detailed brainstem-specific atlases. Methodological improvements have significantly enhanced the accuracy of physiological (cardiac and respiratory) noise correction within brainstem imaging studies. These technological and methodological advancements have allowed for in-depth analyses of the brainstem's anatomy, including quantitative assessments and examinations of structural connectivity within both gray and white matter. Furthermore, functional studies, including assessments of activation patterns and functional connectivity, have revealed the brainstem's roles in both specialized functions and broader neural integration. Notably, these investigations have identified alterations in brainstem structure and function associated with various neurological disorders.

SUMMARY

The aforementioned developments have allowed for a greater appreciation of the importance of the brainstem in the wider context of neuroscience and clinical neurology.

Large Individual Differences in Functional Connectivity in the Context of Major Depression and Antidepressant Pharmacotherapy

Clinical studies of major depression (MD) generally focus on group effects, yet interindividual differences in brain function are increasingly recognized as important and may even impact effect sizes related to group effects. Here, we examine the magnitude of individual differences in relation to group differences that are commonly investigated (e.g., related to MD diagnosis and treatment response). Functional MRI data from 107 participants (63 female, 44 male) were collected at baseline, 2, and 8 weeks during which patients received pharmacotherapy (escitalopram, N = 68) and controls (N = 39) received no intervention. The unique contributions of different sources of variation were examined by calculating how much variance in functional connectivity was shared across all participants and sessions, within/across groups (patients vs controls, responders vs nonresponders, female vs male participants), recording sessions, and individuals. Individual differences and common connectivity across groups, sessions, and participants contributed most to the explained variance (>95% across analyses). Group differences related to MD diagnosis, treatment response, and biological sex made significant but small contributions (0.3-1.2%). High individual variation was present in cognitive control and attention areas, while low individual variation characterized primary sensorimotor regions. Group differences were much smaller than individual differences in the context of MD and its treatment. These results could be linked to the variable findings and difficulty translating research on MD to clinical practice. Future research should examine brain features with low and high individual variation in relation to psychiatric symptoms and treatment trajectories to explore the clinical relevance of the individual differences identified here.

Functional connectivity MRI quality control procedures in CONN

Quality control (QC) for functional connectivity magnetic resonance imaging (FC-MRI) is critical to ensure the validity of neuroimaging studies. Noise confounds are common in MRI data and, if not accounted for, may introduce biases in functional measures affecting the validity, replicability, and interpretation of FC-MRI study results. Although FC-MRI analysis rests on the assumption of adequate data processing, QC is underutilized and not systematically reported. Here, we describe a quality control pipeline for the visual and automated evaluation of MRI data implemented as part of the CONN toolbox. We analyzed publicly available resting state MRI data (N = 139 from 7 MRI sites) from the FMRI Open QC Project. Preprocessing steps included realignment, unwarp, normalization, segmentation, outlier identification, and smoothing. Data denoising was performed based on the combination of scrubbing, motion regression, and aCompCor - a principal component characterization of noise from minimally eroded masks of white matter and of cerebrospinal fluid tissues. Participant-level QC procedures included visual inspection of raw-level data and of representative images after each preprocessing step for each run, as well as the computation of automated descriptive QC measures such as average framewise displacement, average global signal change, prevalence of outlier scans, MNI to anatomical and functional overlap, anatomical to functional overlap, residual BOLD timeseries variability, effective degrees of freedom, and global correlation strength. Dataset-level QC procedures included the evaluation of inter-subject variability in the distributions of edge connectivity in a 1,000-node graph (FC distribution displays), and the estimation of residual associations across participants between functional connectivity strength and potential noise indicators such as participant's head motion and prevalence of outlier scans (QC-FC analyses). QC procedures are demonstrated on the reference dataset with an emphasis on visualization, and general recommendations for best practices are discussed in the context of functional connectivity and other fMRI analysis. We hope this work contributes toward the dissemination and standardization of QC testing performance reporting among peers and in scientific journals.

Altered connectivity of default mode and executive control networks among female patients with persistent post-concussion symptoms

OBJECTIVE

To examine the roles of the default mode network (DMN) and executive control network (ECN) in prolonged recovery after mild traumatic brain injury (mTBI), and relationships with indices of white matter microstructural injury.

METHODS

Seventeen mTBI patients with persistent symptoms were imaged an average of 21.5 months post-injury, along with 23 healthy controls. Resting-state functional magnetic resonance imaging (rs-fMRI) was used to evaluate functional connectivity (FC) of the DMN and ECN. Diffusion tensor imaging (DTI) quantified fractional anisotropy, along with mean, axial and radial diffusivity of white matter tracts.

RESULTS

Compared to controls, patients with mTBI had increased functional connectivity of the DMN and ECN to brain regions implicated in salience and frontoparietal networks, and increased white matter diffusivity within the cerebrum and brainstem. Among the patients, FC was correlated with better neurocognitive test scores, while diffusivity was correlated with more severe self-reported symptoms. The FC and diffusivity values within abnormal brain regions were not significantly correlated.

CONCLUSION

For female mTBI patients with prolonged symptoms, hyper-connectivity may represent a compensatory response that helps to mitigate the effects of mTBI on cognition. These effects are unrelated to indices of microstructural injury, which are correlated with symptom severity, suggesting that rs-fMRI and DTI may capture distinct aspects of pathophysiology.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR)

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

WHOCARES: WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions

Objective. To spatio-temporally resolve cardiac signals in functional magnetic resonance imaging (fMRI) time-series of the human brain using neither external physiological measurements nor ad hoc modelling assumptions.Approach. Cardiac pulsation is a physiological confound of fMRI time-series that introduces spurious signal fluctuations in proximity to blood vessels. fMRI alone is not sufficiently fast to resolve cardiac pulsation. Depending on the ratio between the instantaneous heart-rate and the acquisition sampling frequency (1/TR, with TR being the repetition time), the cardiac signal may alias into the frequency band of neural activation so that its removal through spectral filtering techniques is generally not possible. In this paper, we show that it is feasible to temporally and spatially resolve cardiac signals throughout the brain even when cardiac aliasing occurs by combining fMRI hyper-sampling with simultaneous multislice (SMS) imaging. The technique, which we name WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions (WHOCARES), was developed on 695 healthy subjects selected from the Human Connectome Project and its performance validated against the RETROICOR, HAPPY and the pulse oxymeter signal regression methods.Main results.WHOCARES is capable of retrieving voxel-wise cardiac signal regressors. This is achieved without employing external physiological recordings nor through ad hoc modelling assumptions. The performance of WHOCARES was, on average, superior to RETROICOR, HAPPY and the pulse oxymeter regression methods.Significance.WHOCARES holds basis for the reliable mapping of cardiac activity in fMRI time-series. WHOCARES can be employed for the retrospective removal of cardiac noise in publicly available fMRI datasets where physiological recordings are not available. WHOCARES is freely available athttps://github.com/gferrazzi/WHOCARES.

Bilingualism contributes to reserve and working memory efficiency: Evidence from structural and functional neuroimaging

This study compared brain and behavioral outcomes for monolingual and bilingual older adults who reported no cognitive or memory problems on three types of memory that typically decline in older age, namely, working memory (measured by n-back), item, and associative recognition. The results showed that bilinguals were faster on the two-back working memory task than monolinguals but used a set of frontostriatal regions less than monolinguals. There was no group difference on an item/associative recognition task. In brain structure, gray matter volume and white matter integrity (fractional anisotropy) were generally lower in bilinguals than in monolinguals, but bilinguals had better white matter integrity than monolinguals in the bilateral superior corona radiata and better gray matter density in the left inferior temporal gyrus. These regions may help preserve bilinguals' executive functions despite generally more significant atrophy throughout the brain than monolinguals in that these structures contribute to efficient communication between executive frontal regions and subcortical motor regions, and perceptual pathways. Reliable negative correlations between brain structure and age were only observed in bilinguals, and to the extent that bilinguals (but not monolinguals) had better brain structure, their performance was enhanced. Collectively, the findings provide evidence for reserve in bilingual older adults.

Baseline Functional Connectivity in Resting State Networks Associated with Depression and Remission Status after 16 Weeks of Pharmacotherapy: A CAN-BIND Report

Understanding the neural underpinnings of major depressive disorder (MDD) and its treatment could improve treatment outcomes. So far, findings are variable and large sample replications scarce. We aimed to replicate and extend altered functional connectivity associated with MDD and pharmacotherapy outcomes in a large, multisite sample. Resting-state fMRI data were collected from 129 patients and 99 controls through the Canadian Biomarker Integration Network in Depression. Symptoms were assessed with the Montgomery-Åsberg Depression Rating Scale (MADRS). Connectivity was measured as correlations between four seeds (anterior and posterior cingulate cortex, insula and dorsolateral prefrontal cortex) and all other brain voxels. Partial least squares was used to compare connectivity prior to treatment between patients and controls, and between patients reaching remission (MADRS ≤ 10) early (within 8 weeks), late (within 16 weeks), or not at all. We replicated previous findings of altered connectivity in patients. In addition, baseline connectivity of the anterior/posterior cingulate and insula seeds differentiated patients with different treatment outcomes. The stability of these differences was established in the largest single-site subsample. Our replication and extension of altered connectivity highlighted previously reported and new differences between patients and controls, and revealed features that might predict remission prior to pharmacotherapy. Trial registration:  ClinicalTrials.gov: NCT01655706.

Mapping dependencies of BOLD signal change to end-tidal CO2: linear and nonlinear modeling, and effect of physiological noise correction

BACKGROUND

Disentangling physiological noise and signal of interest is a major issue when evaluating BOLD-signal changes in response to breath holding. Currently-adopted approaches for retrospective noise correction are general-purpose, and have non-negligible effects in studies on hypercapnic challenges.

NEW METHOD

We provide a novel approach to the analysis of specific and non-specific BOLD-signal changes related to end-tidal CO₂ (P_ETCO₂) in breath-hold fMRI studies. Multiple-order nonlinear predictors for P_ETCO₂ model a region-dependent nonlinear input-output relationship hypothesized in literature and possibly playing a crucial role in disentangling noise. We explore Retrospective Image-based Correction (RETROICOR) effects on the estimated BOLD response, applying our analysis both with and without RETROICOR and analyzing the linear and non-linear correlation between P_ETCO₂ and RETROICOR regressors.

RESULTS

The RETROICOR model of noise related to respiratory activity correlated with P_ETCO₂ both linearly and non-linearly. The correction affected the shape of the estimated BOLD response to hypercapnia but allowed to discard spurious activity in ventricles and white matter. Activation clusters were best detected using non-linear components in the BOLD response model.

COMPARISON WITH EXISTING METHOD

We evaluated the side-effects of standard physiological noise correction procedure, tailoring our analysis on challenging understudied brainstem and subcortical regions. Our novel approach allowed to characterize delays and non-linearities in BOLD response.

CONCLUSIONS

RETROICOR successfully avoided false positives, still broadly affecting the estimated non-linear BOLD responses. Non-linearities in the model better explained CO₂-related BOLD signal fluctuations. The necessity to modify the standard procedure for physiological-noise correction in breath-hold studies was addressed, stating its crucial importance.

Exploring brain connectivity changes in major depressive disorder using functional-structural data fusion: A CAN-BIND-1 study

There is a growing interest in examining the wealth of data generated by fusing functional and structural imaging information sources. These approaches may have clinical utility in identifying disruptions in the brain networks that underlie major depressive disorder (MDD). We combined an existing software toolbox with a mathematically dense statistical method to produce a novel processing pipeline for the fast and easy implementation of data fusion analysis (FATCAT‐awFC). The novel FATCAT‐awFC pipeline was then utilized to identify connectivity (conventional functional, conventional structural and anatomically weighted functional connectivy) changes in MDD patients compared to healthy comparison participants (HC). Data were acquired from the Canadian Biomarker Integration Network for Depression (CAN‐BIND‐1) study. Large‐scale resting‐state networks were assessed. We found statistically significant anatomically‐weighted functional connectivity (awFC) group differences in the default mode network and the ventral attention network, with a modest effect size (d < 0.4). Functional and structural connectivity seemed to overlap in significance between one region‐pair within the default mode network. By combining structural and functional data, awFC served to heighten or reduce the magnitude of connectivity differences in various regions distinguishing MDD from HC. This method can help us more fully understand the interconnected nature of structural and functional connectivity as it relates to depression.

Insular Connectivity Is Associated With Self-Appraisal of Cognitive Function After a Concussion

Concussion is associated with acute cognitive impairments, with declines in processing speed and reaction time being common. In the clinical setting, these issues are identified via symptom assessments and neurocognitive test (NCT) batteries. Practice guidelines recommend integrating both symptoms and NCTs into clinical decision-making, but correlations between these measures are often poor. This suggests that many patients experience difficulties in the self-appraisal of cognitive issues. It is presently unclear what neural mechanisms give rise to appraisal mismatch after a concussion. One promising target is the insula, which regulates aspects of cognition, particularly interoception and self-monitoring. The present study tested the hypothesis that appraisal mismatch is due to altered functional connectivity of the insula to frontal and midline structures, with hypo-connectivity leading to under-reporting of cognitive issues and hyper-connectivity leading to over-reporting. Data were collected from 59 acutely concussed individuals and 136 normative controls, including symptom assessments, NCTs and magnetic resonance imaging (MRI) data. Analysis of resting-state functional MRI supported the hypothesis, identifying insular networks that were associated with appraisal mismatch in concussed athletes that included frontal, sensorimotor, and cingulate connections. Subsequent analysis of diffusion tensor imaging also determined that symptom over-reporting was associated with reduced fractional anisotropy and increased mean diffusivity of posterior white matter. These findings provide new insights into the mechanisms of cognitive appraisal mismatch after a concussion. They are of particular interest given the central role of symptom assessments in the diagnosis and clinical management of concussion.

Resting state fMRI scanner instabilities revealed by longitudinal phantom scans in a multi-center study

Quality assurance (QA) is crucial in longitudinal and/or multi-site studies, which involve the collection of data from a group of subjects over time and/or at different locations. It is important to regularly monitor the performance of the scanners over time and at different locations to detect and control for intrinsic differences (e.g., due to manufacturers) and changes in scanner performance (e.g., due to gradual component aging, software and/or hardware upgrades, etc.). As part of the Ontario Neurodegenerative Disease Research Initiative (ONDRI) and the Canadian Biomarker Integration Network in Depression (CAN-BIND), QA phantom scans were conducted approximately monthly for three to four years at 13 sites across Canada with 3T research MRI scanners. QA parameters were calculated for each scan using the functional Biomarker Imaging Research Network's (fBIRN) QA phantom and pipeline to capture between- and within-scanner variability. We also describe a QA protocol to measure the full-width-at-half-maximum (FWHM) of slice-wise point spread functions (PSF), used in conjunction with the fBIRN QA parameters. Variations in image resolution measured by the FWHM are a primary source of variance over time for many sites, as well as between sites and between manufacturers. We also identify an unexpected range of instabilities affecting individual slices in a number of scanners, which may amount to a substantial contribution of unexplained signal variance to their data. Finally, we identify a preliminary preprocessing approach to reduce this variance and/or alleviate the slice anomalies, and in a small human data set show that this change in preprocessing can have a significant impact on seed-based connectivity measurements for some individual subjects. We expect that other fMRI centres will find this approach to identifying and controlling scanner instabilities useful in similar studies.

Vascular origins of low-frequency oscillations in the cerebrospinal fluid signal in resting-state fMRI: Interpretation using photoplethysmography

In vivo mapping of cerebrovascular oscillations in the 0.05–0.15 Hz remains difficult. Oscillations in the cerebrospinal fluid (CSF) represent a possible avenue for noninvasively tracking these oscillations using resting‐state functional MRI (rs‐fMRI), and have been used to correct for vascular oscillations in rs‐fMRI functional connectivity. However, the relationship between low‐frequency CSF and vascular oscillations remains unclear. In this study, we investigate this relationship using fast simultaneous rs‐fMRI and photoplethysmogram (PPG), examining the 0.1 Hz PPG signal, heart‐rate variability (HRV), pulse‐intensity ratio (PIR), and the second derivative of the PPG (SDPPG). The main findings of this study are: (a) signals in different CSF regions are not equivalent in their associations with vascular and tissue rs‐fMRI signals; (b) the PPG signal is maximally coherent with the arterial and CSF signals at the cardiac frequency, but coherent with brain tissue at ~0.2 Hz; (c) PIR is maximally coherent with the CSF signal near 0.03 Hz; and (d) PPG‐related vascular oscillations only contribute to ~15% of the CSF (and arterial) signal in rs‐fMRI. These findings caution against averaging all CSF regions when extracting physiological nuisance regressors in rs‐fMRI applications, and indicate the drivers of the CSF signal are more than simply cardiac. Our study is an initial attempt at the refinement and standardization of how the CSF signal in rs‐fMRI can be used and interpreted. It also paves the way for using rs‐fMRI in the CSF as a potential tool for tracking cerebrovascular health through, for instance, the potential relationship between PIR and the CSF signal.

A retrospective physiological noise correction method for oscillating steady-state imaging

PURPOSE

Oscillating steady-state imaging (OSSI) is an SNR-efficient steady-state sequence with T 2 ∗ sensitivity suitable for FMRI. Due to the frequency sensitivity of the signal, respiration- and drift-induced field changes can create unwanted signal fluctuations. This study aims to address this issue by developing retrospective signal correction methods that utilize OSSI signal properties to denoise task-based OSSI FMRI experiments.

METHODS

A retrospective denoising approach was developed that leverages the unique signal properties of OSSI to perform denoising without a manually specified noise region of interest and works with both voxel timecourses (oscillating steady-state correction [OSSCOR]) or FID timecourses (F-OSSCOR). Simulations were performed to estimate the number of principal components optimal for denoising. In vivo experiments at 3 T field strength were conducted to compare the performance of proposed methods against a standard principal component analysis-based method, measured using mean t score within an region of interest, number of activations, and mean temporal SNR.

RESULTS

Correction using OSSCOR was significantly better than the standard method in all metrics. Correction using F-OSSCOR was not significantly different from the standard method using an equal number of principal components. Increasing the number of OSSCOR principal components decreased activation strength and increased the number of suspected false positives. However, increasing the number of principal components in F-OSSCOR increased activation strength with little to no increase in false activation.

CONCLUSION

Both OSSCOR and F-OSSCOR substantially reduce physiological noise components and increase temporal SNR, improving the functional results of task-based OSSI functional experiments. F-OSSCOR demonstrates a proof of concept utilization of coil-localized FID signal information for physiological noise correction.

Functional magnetic resonance imaging of the trail-making test in older adults

The trail-making test (TMT) is a popular neuropsychological test, which is used extensively to measure cognitive impairment associated with neurodegenerative disorders in older adults. Behavioural performance on the TMT has been investigated in older populations, but there is limited research on task-related brain activity in older adults. The current study administered a naturalistic version of the TMT to a healthy older-aged population in an MRI environment using a novel, MRI-compatible tablet. Functional MRI was conducted during task completion, allowing characterization of the brain activity associated with the TMT. Performance on the TMT was evaluated using number of errors and seconds per completion of each link. Results are reported for 36 cognitively healthy older adults between the ages of 52 and 85. Task-related activation was observed in extensive regions of the bilateral frontal, parietal, temporal and occipital lobes as well as key motor areas. Increased age was associated with reduced brain activity and worse task performance. Specifically, older age was correlated with decreased task-related activity in the bilateral occipital, temporal and parietal lobes. These results suggest that healthy older aging significantly affects brain function during the TMT, which consequently may result in performance decrements. The current study reveals the brain activation patterns underlying TMT performance in a healthy older aging population, which functions as an important, clinically-relevant control to compare to pathological aging in future investigations.

Reliability of a functional magnetic resonance imaging task of emotional conflict in healthy participants

Task-based functional neuroimaging methods are increasingly being used to identify biomarkers of treatment response in psychiatric disorders. To facilitate meaningful interpretation of neural correlates of tasks and their potential changes with treatment over time, understanding the reliability of the blood-oxygen-level dependent (BOLD) signal of such tasks is essential. We assessed test-retest reliability of an emotional conflict task in healthy participants collected as part of the Canadian Biomarker Integration Network in Depression. Data for 36 participants, scanned at three time points (weeks 0, 2, and 8) were analyzed, and intra-class correlation coefficients (ICC) were used to quantify reliability. We observed moderate reliability (median ICC values between 0.5 and 0.6), within occipital, parietal, and temporal regions, specifically for conditions of lower cognitive complexity, that is, face, congruent or incongruent trials. For these conditions, activation was also observed within frontal and sub-cortical regions, however, their reliability was poor (median ICC < 0.2). Clinically relevant prognostic markers based on task-based fMRI require high predictive accuracy at an individual level. For this to be achieved, reliability of BOLD responses needs to be high. We have shown that reliability of the BOLD response to an emotional conflict task in healthy individuals is moderate. Implications of these findings to further inform studies of treatment effects and biomarker discovery are discussed.

The Contribution of Functional Magnetic Resonance Imaging to the Understanding of the Effects of Acute Physical Exercise on Cognition

The fact that a single bout of acute physical exercise has a positive impact on cognition is well-established in the literature, but the neural correlates that underlie these cognitive improvements are not well understood. Here, the use of neuroimaging techniques, such as functional magnetic resonance imaging (fMRI), offers great potential, which is just starting to be recognized. This review aims at providing an overview of those studies that used fMRI to investigate the effects of acute physical exercises on cerebral hemodynamics and cognition. To this end, a systematic literature survey was conducted by two independent reviewers across five electronic databases. The search returned 668 studies, of which 14 studies met the inclusion criteria and were analyzed in this systematic review. Although the findings of the reviewed studies suggest that acute physical exercise (e.g., cycling) leads to profound changes in functional brain activation, the small number of available studies and the great variability in the study protocols limits the conclusions that can be drawn with certainty. In order to overcome these limitations, new, more well-designed trials are needed that (i) use a more rigorous study design, (ii) apply more sophisticated filter methods in fMRI data analysis, (iii) describe the applied processing steps of fMRI data analysis in more detail, and (iv) provide a more precise exercise prescription.

Escitalopram ameliorates differences in neural activity between healthy comparison and major depressive disorder groups on an fMRI Emotional conflict task: A CAN-BIND-1 study

BACKGROUND

Identifying objective biomarkers can assist in predicting remission/non-remission to treatment, improving remission rates, and reducing illness burden in major depressive disorder (MDD).

METHODS

Sixteen MDD 8-week remitters (MDD-8), twelve 16-week remitters (MDD-16), 14 non-remitters (MDD-NR) and 30 healthy comparison participants (HC) completed a functional magnetic resonance imaging emotional conflict task at baseline, prior to treatment with escitalopram, and 8 weeks after treatment initiation. Patients were followed 16 weeks to assess remitter status.

RESULTS

All groups demonstrated emotional Stroop in reaction time (RT) at baseline and Week 8. There were no baseline differences between HC and MDD-8, MDD-16, or MDD-NR in RT or accuracy. By Week 8, MDD-8 demonstrated poorer accuracy compared to HC. Compared to HC, the baseline blood-oxygen level dependent (BOLD) signal was decreased in MDD-8 in brain-stem and thalamus; in MDD-16 in lateral occipital cortex, middle temporal gyrus, and cuneal cortex; in MDD-NR in lingual and occipital fusiform gyri, thalamus, putamen, caudate, cingulate gyrus, insula, cuneal cortex, and middle temporal gyrus. By Week 8, there were no BOLD activity differences between MDD groups and HC.

LIMITATIONS

The Emotional Conflict Task lacks a neutral (non-emotional) condition, restricting interpretation of how mood may influence perception of non-emotionally valenced stimuli.

CONCLUSIONS

The Emotional Conflict Task is not an objective biomarker for remission trajectory in patients with MDD receiving escitalopram treatment. Escitalopram may have influenced emotion recognition in MDD groups in terms of augmented accuracy and BOLD signal in response to an Emotional Conflict Task, following 8 weeks of escitalopram treatment.

Meta-analytic and functional connectivity evidence from functional magnetic resonance imaging for an anterior to posterior gradient of function along the hippocampal axis

There is considerable evidence from non-human animal studies that the anterior and posterior regions of the hippocampus have different anatomical connections and support different behavioural functions. Although there are some recent human studies using functional magnetic resonance imaging (fMRI) that have addressed this idea directly in the memory and spatial processing domains and provided support for it, there has been no broader meta-analysis of the fMRI literature to determine if there is consistent evidence for functional dissociations in anterior and posterior hippocampus across all of the different cognitive domains in which the hippocampus participates. The purpose of this review is to address this gap in our knowledge using three approaches. One approach involved PubMed searches to identify relevant fMRI papers reporting hippocampal activation during episodic encoding and retrieval, semantic retrieval, working memory, spatial navigation, simulation/scene construction, transitive inference, and social cognition tasks. The second was to use a large meta-analytic database (neurosynth) to find text terms and coactivation maps associated with the anterior and posterior hippocampal regions identified in the literature search. The third approach was to contrast the resting-state functional connectivity of the anterior and posterior hippocampal regions using a publicly available database that includes a large sample of adults. These three approaches provided converging evidence that not only are cognitive processes differently distributed along the hippocampal axis, but there also are distinct areas coactivated and functionally connected with the anterior and posterior segments. This anterior/posterior distinction involving multiple cognitive domains is consistent with the animal literature and provides strong support from fMRI for the idea of functional dissociations across the long axis of the hippocampus.

Model-based physiological noise removal in fast fMRI

Recent improvements in the speed and sensitivity of fMRI acquisition techniques suggest that fast fMRI can be used to detect and precisely localize sub-second neural dynamics. This enhanced temporal resolution has enormous potential for neuroscientists. However, physiological noise poses a major challenge for the analysis of fast fMRI data. Physiological noise scales with sensitivity, and its autocorrelation structure is altered in rapidly sampled data, suggesting that new approaches are needed for physiological noise removal in fast fMRI. Existing strategies either rely on external physiological recordings, which can be noisy or difficult to collect, or employ data-driven approaches which make assumptions that may not hold true in fast fMRI. We created a statistical model of harmonic regression with autoregressive noise (HRAN) to estimate and remove cardiac and respiratory noise from the fMRI signal directly. This technique exploits the fact that cardiac and respiratory noise signals are fully sampled (rather than aliasing) when imaging at fast rates, allowing us to track and model physiology over time without requiring external physiological measurements. We then created a joint model of neural hemodynamics, and physiological and autocorrelated noise to more accurately remove noise. We first verified that HRAN accurately estimates cardiac and respiratory dynamics and that our model demonstrates goodness-of-fit in fast fMRI data. In task-driven data, we then demonstrated that HRAN is able to remove physiological noise while leaving the neural signal intact, thereby increasing detection of task-driven voxels. Finally, we established that in both simulations and fast fMRI data HRAN is able to improve statistical inferences as compared with gold-standard physiological noise removal techniques. In conclusion, we created a tool that harnesses the novel information in fast fMRI to remove physiological noise, enabling broader use of the technology to study human brain function.

Reduced accuracy accompanied by reduced neural activity during the performance of an emotional conflict task by unmedicated patients with major depression: A CAN-BIND fMRI study

METHODS

We studied 48 MDD and 30 HC who performed an emotional conflict task in a functional magnetic resonance imaging (fMRI) scanner.

RESULTS

On the emotional conflict task, MDD and HC demonstrated a robust emotional Stroop effect in reaction time and accuracy. Overall, accuracy was lower in MDD compared to HC with no significant reaction time differences. The fMRI data indicated lower BOLD activation in MDD compared to HC on comparisons of all trials, congruent, incongruent, and incongruent > congruent trials in regions including right inferior temporal gyrus, lateral occipital cortex, and occipital fusiform gyrus. Behavioural and neuroimaging data indicated no group differences in fearful versus happy face processing.

LIMITATIONS

Inclusion of a neutral condition may have provided a valuable contrast to how MDD and HC process stimuli without emotional valence compared to stimuli with a strong emotional valence.

CONCLUSIONS

MDD and HC demonstrated a robust emotional Stroop effect. Compared to HC, MDD demonstrated an overall reduced accuracy on the emotional conflict task and reduced BOLD activity in regions important for face perception and emotion information processing, with no differences in responding to fearful versus happy faces. These findings provide support for the theory of emotion context insensitivity in individuals with depression.

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

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

Social motivational processing and interpersonal function in aging cocaine smokers

Illicit drug use among aging cohorts is increasing, yet little is known about functional impairments in older drug users. Given the importance of social integration for aging and documented social decrements in cocaine users, we examined social function and its neurocognitive substrates in aging cocaine users relative to carefully matched non-cocaine users. Regular (≥twice/week), long-term (≥15 years) cocaine smokers 50-60 years old (COCs; n = 22; four women) and controls (CTRLs; n = 19; four women) underwent standardized probes of social reward and threat processing during functional magnetic resonance imaging and a behavioral facial affect recognition task. Self-report and peer-report of daily interpersonal function were also collected. COCs, and CTRLs reporting current marijuana or alcohol use, were tested after four drug-free inpatient days. COCs had pronounced problems in daily social function relative to CTRLs indicated by both their own and their peers' reports. Compared with CTRLs, COCs had stronger amygdala responses to social threat versus control stimuli, with no other differences in social processing or cognition. Aging cocaine users appear to have marked, generalized difficulties in 'real-world' interpersonal function but largely intact social processing on laboratory-based measures when compared with appropriately matched controls and tested under well-controlled conditions. Daily social difficulties may be related to transient factors such as acute/residual drug effects or cocaine-related changes in health behaviors (e.g. disrupted sleep and poor diet). These data suggest that interpersonal function may be a valid intervention target for aging cocaine users and warrants further study in older drug users.

Functional MRI of Letter Cancellation Task Performance in Older Adults

The Letter Cancellation Task (LCT) is a widely used pen-and-paper probe of attention in clinical and research settings. Despite its popularity, the neural correlates of the task are not well understood. The present study uses functional magnetic resonance imaging (fMRI) and specialized tablet technology to identify the neural correlates of the LCT in 32 healthy older adults between 50-85 years of age, and further investigates the effect of healthy aging on performance. Subjects performed the LCT in its standard pen-and-paper administration and with the tablet during fMRI. Performance on the tablet was significantly slower than on pen-and-paper, with both response modes showing slower performance as a function of age. Across all ages, bilateral brain activation was observed in the cerebellum, superior temporal lobe, precentral gyrus, frontal gyrus, and occipital and parietal areas. Increasing age correlated with reduced brain activity in the supplementary motor area, middle occipital gyrus, medial and inferior frontal gyrus, cerebellum and putamen. Better LCT performance was correlated with increased activity in the middle frontal gyrus, and reduced activity in the cerebellum. The brain regions activated are associated with visuospatial attention and motor control, and are consistent with the neural correlates of LCT performance previously identified in lesion studies.

BOLD signal variability and complexity in children and adolescents with and without autism spectrum disorder

•

Resting-state BOLD signal variability and complexity were examined.

•

No significant group differences were observed in youth with and without autism.

•

A continuum of brain-behavior relationships was observed across diagnostic groups.

•

Positive correlations were found between brain measures, age and global efficiency.

•

Negative correlations were found between the brain measures and behavioral severity.

Variability of neural signaling is an important index of healthy brain functioning, as is signal complexity, which relates to information processing capacity. Alterations in variability and complexity may underlie certain brain dysfunctions. Here, resting-state fMRI was used to examine variability and complexity in children and adolescents with and without autism spectrum disorder (ASD). Variability was measured using the mean square successive difference (MSSD) of the time series, and complexity was assessed using sample entropy. A categorical approach was implemented to determine if the brain measures differed between diagnostic groups (ASD and controls). A dimensional approach was used to examine the continuum of relationships between each brain measure and behavioral severity, age, IQ, and the global efficiency (GE) of each participant’s structural connectome, which reflects the structural capacity for information processing. Using the categorical approach, no significant group differences were found for neither MSSD nor entropy. The dimensional approach revealed significant positive correlations between each brain measure, GE, and age. Negative correlations were observed between each brain measure and the severity of ASD behaviors across all participants. These results reveal the nature of variability and complexity of BOLD signals in children and adolescents with and without ASD.

The Neural Correlates of the Clock-Drawing Test in Healthy Aging

Importance: The clock-drawing test (CDT) is an important neurocognitive assessment tool, widely used as a screening test for dementia. Behavioral performance on the test has been studied extensively, but there is scant literature on the underlying neural correlates.

Purpose: To administer the CDT naturalistically to a healthy older aging population in an MRI environment, and characterize the brain activity associated with test completion.

Main Outcome and Measure: Blood-oxygen-level dependent (BOLD) functional MRI was conducted as participants completed the CDT using novel tablet technology. Brain activity during CDT performance was contrasted to rest periods of visual fixation. Performance on the CDT was evaluated using a standardized scoring system (Rouleau score) and time to test completion. To assess convergent validity, performance during fMRI was compared to performance on a standard paper version of the task, administered in a psychometric testing room.

Results: Study findings are reported for 33 cognitively healthy older participants aged 52–85. Activation was observed in the bilateral frontal, occipital and parietal lobes as well as the supplementary motor area and precentral gyri. Increased age was significantly correlated with Rouleau scores on the clock number drawing (R2) component (rho = -0.55, p < 0.001); the clock hand drawing (R3) component (rho = -0.50, p < 0.005); and the total clock (rho = -0.62, p < 0.001). Increased age was also associated with decreased activity in the bilateral parietal and occipital lobes as well as the right temporal lobe and right motor areas.

Conclusion and Relevance: This imaging study characterizes the brain activity underlying performance of the CDT in a healthy older aging population using the most naturalistic version of the task to date. The results suggest that the functions of the occipital and parietal lobe are significantly altered by the normal aging process, which may lead to performance decrements.

Brain connectivity tracks effects of chemotherapy separately from behavioral measures

Several studies in cancer research have suggested that cognitive dysfunction following chemotherapy, referred to in lay terms as "chemobrain", is a serious problem. At present, the changes in integrative brain function that underlie such dysfunction remain poorly understood. Recent developments in neuroimaging suggest that patterns of functional connectivity can provide a broadly applicable neuromarker of cognitive performance and other psychometric measures. The current study used multivariate analysis methods to identify patterns of disruption in resting state functional connectivity of the brain due to chemotherapy and the degree to which the disruptions can be linked to behavioral measures of distress and cognitive performance. Sixty two women (22 healthy control, 18 patients treated with adjuvant chemotherapy, and 22 treated without chemotherapy) were evaluated with neurocognitive measures followed by self-report questionnaires and open eyes resting-state fMRI scanning at three time points: diagnosis (M0, pre-adjuvant treatment), 1 month (M1), and 7 months (M7) after treatment. The results indicated deficits in cognitive health of breast cancer patients immediately after chemotherapy that improved over time. This psychological trajectory was paralleled by a disruption and later recovery of resting-state functional connectivity, mostly in the parietal and frontal brain regions. Mediation analysis showed that the functional connectivity alteration pattern is a separable treatment symptom from the decreased cognitive health. Current study indicates that more targeted support for patients should be developed to ameliorate these multi-faceted side effects of chemotherapy treatment on neural functioning and cognitive health.

Different neural activities support auditory working memory in musicians and bilinguals

Musical training and bilingualism benefit executive functioning and working memory (WM)-however, the brain networks supporting this advantage are not well specified. Here, we used functional magnetic resonance imaging and the n-back task to assess WM for spatial (sound location) and nonspatial (sound category) auditory information in musician monolingual (musicians), nonmusician bilinguals (bilinguals), and nonmusician monolinguals (controls). Musicians outperformed bilinguals and controls on the nonspatial WM task. Overall, spatial and nonspatial WM were associated with greater activity in dorsal and ventral brain regions, respectively. Increasing WM load yielded similar recruitment of the anterior-posterior attention network in all three groups. In both tasks and both levels of difficulty, musicians showed lower brain activity than controls in superior prefrontal frontal gyrus and dorsolateral prefrontal cortex (DLPFC) bilaterally, a finding that may reflect improved and more efficient use of neural resources. Bilinguals showed enhanced activity in language-related areas (i.e., left DLPFC and left supramarginal gyrus) relative to musicians and controls, which could be associated with the need to suppress interference associated with competing semantic activations from multiple languages. These findings indicate that the auditory WM advantage in musicians and bilinguals is mediated by different neural networks specific to each life experience.

Age differences in brain signal variability are robust to multiple vascular controls

A host of studies support that younger, better performing adults express greater moment-to-moment blood oxygen level-dependent (BOLD) signal variability (SD_BOLD) in various cortical regions, supporting an emerging view that the aging brain may undergo a generalized reduction in dynamic range. However, the exact physiological nature of age differences in SD_BOLD remains understudied. In a sample of 29 younger and 45 older adults, we examined the contribution of vascular factors to age group differences in fixation-based SD_BOLD using (1) a dual-echo BOLD/pseudo-continuous arterial spin labeling (pCASL) sequence, and (2) hypercapnia via a computer-controlled gas delivery system. We tested the hypothesis that, although SD_BOLD may relate to individual differences in absolute cerebral blood flow (CBF), BOLD cerebrovascular reactivity (CVR), or maximum BOLD signal change (M), robust age differences in SD_BOLD would remain after multiple statistical controls for these vascular factors. As expected, our results demonstrated that brain regions in which younger adults expressed higher SD_BOLD persisted after comprehensive control of vascular effects. Our findings thus further establish BOLD signal variability as an important marker of the aging brain.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

An information theory framework for dynamic functional domain connectivity

BACKGROUND

Dynamic functional network connectivity (dFNC) analyzes time evolution of coherent activity in the brain. In this technique dynamic changes are considered for the whole brain. This paper proposes an information theory framework to measure information flowing among subsets of functional networks call functional domains.

NEW METHOD

Our method aims at estimating bits of information contained and shared among domains. The succession of dynamic functional states is estimated at the domain level. Information quantity is based on the probabilities of observing each dynamic state. Mutual information measurement is then obtained from probabilities across domains. Thus, we named this value the cross domain mutual information (CDMI).

RESULTS

Strong CDMIs were observed in relation to the subcortical domain. Domains related to sensorial input, motor control and cerebellum form another CDMI cluster. Information flow among other domains was seldom found.

COMPARISON WITH EXISTING METHODS

Other methods of dynamic connectivity focus on whole brain dFNC matrices. In the current framework, information theory is applied to states estimated from pairs of multi-network functional domains. In this context, we apply information theory to measure information flow across functional domains.

CONCLUSION

Identified CDMI clusters point to known information pathways in the basal ganglia and also among areas of sensorial input, patterns found in static functional connectivity. In contrast, CDMI across brain areas of higher level cognitive processing follow a different pattern that indicates scarce information sharing. These findings show that employing information theory to formally measured information flow through brain domains reveals additional features of functional connectivity.

Optimizing fMRI preprocessing pipelines for block-design tasks as a function of age

Functional Magnetic Resonance Imaging (fMRI) is a powerful neuroimaging tool, which is often hampered by significant noise confounds. There is evidence that our ability to detect activations in task fMRI is highly dependent on the preprocessing steps used to control noise and artifact. However, the vast majority of studies examining preprocessing pipelines in fMRI have focused on young adults. Given the widespread use of fMRI for characterizing the neurobiology of aging, it is critical to examine how the impact of preprocessing choices varies as a function of age. In this study, we employ the NPAIRS cross-validation framework, which optimizes pipelines based on metrics of prediction accuracy (P) and spatial reproducibility (R), to compare the effects of pipeline optimization between young (21-33 years) and older (61-82 years) cohorts, for three different block-design contrasts. Motion is shown to be a greater issue in the older cohort, and we introduce new statistical approaches to control for potential biases due to head motion during pipeline optimization. In comparison, data-driven methods of physiological noise correction show comparable benefits for both young and old cohorts. Using our optimization framework, we demonstrate that the optimal pipelines tend to be highly similar across age cohorts. In addition, there is a comparable, significant benefit of pipeline optimization across age cohorts, for (P, R) metrics and independent validation measures of activation overlap (both between-subject, within-session and within-subject, between-session). The choice of task contrast consistently shows a greater impact than the age cohort, for (P, R) metrics and activation overlap. Finally, adaptive pipeline optimization per task run shows improved sensitivity to age-related changes in brain activity, particularly for weaker, more complex cognitive contrasts. The current study provides the first detailed examination of preprocessing pipelines across age cohorts, demonstrating a significant benefit of adaptive pipeline optimization across age groups.

Quantifying functional connectivity in multi-subject fMRI data using component models

Functional magnetic resonance imaging (fMRI) is increasingly used to characterize functional connectivity between brain regions. Given the vast number of between-voxel interactions in high-dimensional fMRI data, it is an ongoing challenge to detect stable and generalizable functional connectivity in the brain among groups of subjects. Component models can be used to define subspace representations of functional connectivity that are more interpretable. It is, however, unclear which component model provides the optimal representation of functional networks for multi-subject fMRI datasets. A flexible cross-validation approach that assesses the ability of the models to predict voxel-wise covariance in new data, using three different measures of generalization was proposed. This framework is used to compare a range of component models with varying degrees of flexibility in their representation of functional connectivity, evaluated on both simulated and experimental resting-state fMRI data. It was demonstrated that highly flexible subject-specific component subspaces, as well as very constrained average models, are poor predictors of whole-brain functional connectivity, whereas the best-generalizing models account for subject variability within a common spatial subspace. Within this set of models, spatial Independent Component Analysis (sICA) on concatenated data provides more interpretable brain patterns, whereas a consistent-covariance model that accounts for subject-specific network scaling (PARAFAC2) provides greater stability in functional connectivity relationships between components and their spatial representations. The proposed evaluation framework is a promising quantitative approach to evaluating component models, and reveals important differences between subspace models in terms of predictability, robustness, characterization of subject variability, and interpretability of the model parameters. Hum Brain Mapp 38:882-899, 2017. © 2016 Wiley Periodicals, Inc.

The Functional Segregation and Integration Model: Mixture Model Representations of Consistent and Variable Group-Level Connectivity in fMRI

The brain consists of specialized cortical regions that exchange information between each other, reflecting a combination of segregated (local) and integrated (distributed) processes that define brain function. Functional magnetic resonance imaging (fMRI) is widely used to characterize these functional relationships, although it is an ongoing challenge to develop robust, interpretable models for high-dimensional fMRI data. Gaussian mixture models (GMMs) are a powerful tool for parcellating the brain, based on the similarity of voxel time series. However, conventional GMMs have limited parametric flexibility: they only estimate segregated structure and do not model interregional functional connectivity, nor do they account for network variability across voxels or between subjects. To address these issues, this letter develops the functional segregation and integration model (FSIM). This extension of the GMM framework simultaneously estimates spatial clustering and the most consistent group functional connectivity structure. It also explicitly models network variability, based on voxel- and subject-specific network scaling profiles. We compared the FSIM to standard GMM in a predictive cross-validation framework and examined the importance of different model parameters, using both simulated and experimental resting-state data. The reliability of parcellations is not significantly altered by flexibility of the FSIM, whereas voxel- and subject-specific network scaling profiles significantly improve the ability to predict functional connectivity in independent test data. Moreover, the FSIM provides a set of interpretable parameters to characterize both consistent and variable aspects functional connectivity structure. As an example of its utility, we use subject-specific network profiles to identify brain regions where network expression predicts subject age in the experimental data. Thus, the FSIM is effective at summarizing functional connectivity structure in group-level fMRI, with applications in modeling the relationships between network variability and behavioral/demographic variables.

The suppression of scale-free fMRI brain dynamics across three different sources of effort: aging, task novelty and task difficulty

There is growing evidence that fluctuations in brain activity may exhibit scale-free (“fractal”) dynamics. Scale-free signals follow a spectral-power curve of the form P(f ) ∝ f^−β, where spectral power decreases in a power-law fashion with increasing frequency. In this study, we demonstrated that fractal scaling of BOLD fMRI signal is consistently suppressed for different sources of cognitive effort. Decreases in the Hurst exponent (H), which quantifies scale-free signal, was related to three different sources of cognitive effort/task engagement: 1) task difficulty, 2) task novelty, and 3) aging effects. These results were consistently observed across multiple datasets and task paradigms. We also demonstrated that estimates of H are robust across a range of time-window sizes. H was also compared to alternative metrics of BOLD variability (SD_BOLD) and global connectivity (Gconn), with effort-related decreases in H producing similar decreases in SD_BOLD and Gconn. These results indicate a potential global brain phenomenon that unites research from different fields and indicates that fractal scaling may be a highly sensitive metric for indexing cognitive effort/task engagement.

Building a Science of Individual Differences from fMRI

To date, fMRI research has been concerned primarily with evincing generic principles of brain function through averaging data from multiple subjects. Given rapid developments in both hardware and analysis tools, the field is now poised to study fMRI-derived measures in individual subjects, and to relate these to psychological traits or genetic variations. We discuss issues of validity, reliability and statistical assessment that arise when the focus shifts to individual subjects and that are applicable also to other imaging modalities. We emphasize that individual assessment of neural function with fMRI presents specific challenges and necessitates careful consideration of anatomical and vascular between-subject variability as well as sources of within-subject variability.

Test-retest reliability of the default mode network in a multi-centric fMRI study of healthy elderly: Effects of data-driven physiological noise correction techniques

Understanding how to reduce the influence of physiological noise in resting state fMRI data is important for the interpretation of functional brain connectivity. Limited data is currently available to assess the performance of physiological noise correction techniques, in particular when evaluating longitudinal changes in the default mode network (DMN) of healthy elderly participants. In this 3T harmonized multisite fMRI study, we investigated how different retrospective physiological noise correction (rPNC) methods influence the within-site test-retest reliability and the across-site reproducibility consistency of DMN-derived measurements across 13 MRI sites. Elderly participants were scanned twice at least a week apart (five participants per site). The rPNC methods were: none (NPC), Tissue-based regression, PESTICA and FSL-FIX. The DMN at the single subject level was robustly identified using ICA methods in all rPNC conditions. The methods significantly affected the mean z-scores and, albeit less markedly, the cluster-size in the DMN; in particular, FSL-FIX tended to increase the DMN z-scores compared to others. Within-site test-retest reliability was consistent across sites, with no differences across rPNC methods. The absolute percent errors were in the range of 5-11% for DMN z-scores and cluster-size reliability. DMN pattern overlap was in the range 60-65%. In particular, no rPNC method showed a significant reliability improvement relative to NPC. However, FSL-FIX and Tissue-based physiological correction methods showed both similar and significant improvements of reproducibility consistency across the consortium (ICC = 0.67) for the DMN z-scores relative to NPC. Overall these findings support the use of rPNC methods like tissue-based or FSL-FIX to characterize multisite longitudinal changes of intrinsic functional connectivity. Hum Brain Mapp 37:2114-2132, 2016. © 2016 Wiley Periodicals, Inc.

Reliability of Task-Based fMRI for Preoperative Planning: A Test-Retest Study in Brain Tumor Patients and Healthy Controls

Background

Functional magnetic resonance imaging (fMRI) continues to develop as a clinical tool for patients with brain cancer, offering data that may directly influence surgical decisions. Unfortunately, routine integration of preoperative fMRI has been limited by concerns about reliability. Many pertinent studies have been undertaken involving healthy controls, but work involving brain tumor patients has been limited. To develop fMRI fully as a clinical tool, it will be critical to examine these reliability issues among patients with brain tumors. The present work is the first to extensively characterize differences in activation map quality between brain tumor patients and healthy controls, including the effects of tumor grade and the chosen behavioral testing paradigm on reliability outcomes.

Method

Test-retest data were collected for a group of low-grade (n = 6) and high-grade glioma (n = 6) patients, and for matched healthy controls (n = 12), who performed motor and language tasks during a single fMRI session. Reliability was characterized by the spatial overlap and displacement of brain activity clusters, BOLD signal stability, and the laterality index. Significance testing was performed to assess differences in reliability between the patients and controls, and low-grade and high-grade patients; as well as between different fMRI testing paradigms.

Results

There were few significant differences in fMRI reliability measures between patients and controls. Reliability was significantly lower when comparing high-grade tumor patients to controls, or to low-grade tumor patients. The motor task produced more reliable activation patterns than the language tasks, as did the rhyming task in comparison to the phonemic fluency task.

Conclusion

In low-grade glioma patients, fMRI data are as reliable as healthy control subjects. For high-grade glioma patients, further investigation is required to determine the underlying causes of reduced reliability. To maximize reliability outcomes, testing paradigms should be carefully selected to generate robust activation patterns.

An Automated, Adaptive Framework for Optimizing Preprocessing Pipelines in Task-Based Functional MRI

BOLD fMRI is sensitive to blood-oxygenation changes correlated with brain function; however, it is limited by relatively weak signal and significant noise confounds. Many preprocessing algorithms have been developed to control noise and improve signal detection in fMRI. Although the chosen set of preprocessing and analysis steps (the “pipeline”) significantly affects signal detection, pipelines are rarely quantitatively validated in the neuroimaging literature, due to complex preprocessing interactions. This paper outlines and validates an adaptive resampling framework for evaluating and optimizing preprocessing choices by optimizing data-driven metrics of task prediction and spatial reproducibility. Compared to standard “fixed” preprocessing pipelines, this optimization approach significantly improves independent validation measures of within-subject test-retest, and between-subject activation overlap, and behavioural prediction accuracy. We demonstrate that preprocessing choices function as implicit model regularizers, and that improvements due to pipeline optimization generalize across a range of simple to complex experimental tasks and analysis models. Results are shown for brief scanning sessions (<3 minutes each), demonstrating that with pipeline optimization, it is possible to obtain reliable results and brain-behaviour correlations in relatively small datasets.

Amphetamine modulates brain signal variability and working memory in younger and older adults

Better-performing younger adults typically express greater brain signal variability relative to older, poorer performers. Mechanisms for age and performance-graded differences in brain dynamics have, however, not yet been uncovered. Given the age-related decline of the dopamine (DA) system in normal cognitive aging, DA neuromodulation is one plausible mechanism. Hence, agents that boost systemic DA [such as d-amphetamine (AMPH)] may help to restore deficient signal variability levels. Furthermore, despite the standard practice of counterbalancing drug session order (AMPH first vs. placebo first), it remains understudied how AMPH may interact with practice effects, possibly influencing whether DA up-regulation is functional. We examined the effects of AMPH on functional-MRI-based blood oxygen level-dependent (BOLD) signal variability (SD(BOLD)) in younger and older adults during a working memory task (letter n-back). Older adults expressed lower brain signal variability at placebo, but met or exceeded young adult SD(BOLD) levels in the presence of AMPH. Drug session order greatly moderated change-change relations between AMPH-driven SD(BOLD) and reaction time means (RT(mean)) and SDs (RT(SD)). Older adults who received AMPH in the first session tended to improve in RT(mean) and RT(SD) when SD(BOLD) was boosted on AMPH, whereas younger and older adults who received AMPH in the second session showed either a performance improvement when SD(BOLD) decreased (for RT(mean)) or no effect at all (for RT(SD)). The present findings support the hypothesis that age differences in brain signal variability reflect aging-induced changes in dopaminergic neuromodulation. The observed interactions among AMPH, age, and session order highlight the state- and practice-dependent neurochemical basis of human brain dynamics.

A Novel Brain Networks Enhancement Model (BNEM) for BOLD fMRI Data Analysis With Highly Spatial Reproducibility

Independent component analysis aiming at detecting the functional connectivity among discrete cortical brain regions has been extensively used to explore the functional magnetic resonance imaging data. Although the independent components (ICs) were with relatively high quality, the noise embedding in ICs has a great impact on the true active/inactive region inference and the reproducibility, in postprocessing stage, e.g., the extraction of statistical parametrical maps (SPMs). In this paper, a novel brain network enhancement model (BNEM) is proposed, which mainly consists of two key techniques: 1) 3-D wavelet noise filter (3DWNF) for the meaningful ICs, which greatly suppresses noise and enforces the real activation inference of SPMs; and 2) a spatial reproducibility enhancement algorithm (SREA), aiming to improve the reproducibility of SPMs. The simulated experiment demonstrated that the postfiltering signals by 3DWNF were with higher correlation and less normalized mean square error to the ground truths than the prefiltering ones; SREA could further enhance the quality of most postfiltering ones, preserving the consistency with 3DWNF. The real data experiments also revealed that 1) 3DWNF could lead to more accurate preservation of the true positive voxels by correctly identifying the high proportionally misclassified voxels of the nonenhanced SPMs; 2) SREA could further improve the classification accuracy of the active/inactive voxels of SPMs corresponding to the 3DWNF denoised ICs; and 3) both 3DWNF and SREA contribute to the reproducibility enhancement of the reproduced SPMs by BNEM. Thus, BNEM is expected to have wide applicability in the neuroscience and clinical domain.

Graph analysis of functional brain networks: practical issues in translational neuroscience

The brain can be regarded as a network: a connected system where nodes, or units, represent different specialized regions and links, or connections, represent communication pathways. From a functional perspective, communication is coded by temporal dependence between the activities of different brain areas. In the last decade, the abstract representation of the brain as a graph has allowed to visualize functional brain networks and describe their non-trivial topological properties in a compact and objective way. Nowadays, the use of graph analysis in translational neuroscience has become essential to quantify brain dysfunctions in terms of aberrant reconfiguration of functional brain networks. Despite its evident impact, graph analysis of functional brain networks is not a simple toolbox that can be blindly applied to brain signals. On the one hand, it requires the know-how of all the methodological steps of the pipeline that manipulate the input brain signals and extract the functional network properties. On the other hand, knowledge of the neural phenomenon under study is required to perform physiologically relevant analysis. The aim of this review is to provide practical indications to make sense of brain network analysis and contrast counterproductive attitudes.

An open science resource for establishing reliability and reproducibility in functional connectomics

Efforts to identify meaningful functional imaging-based biomarkers are limited by the ability to reliably characterize inter-individual differences in human brain function. Although a growing number of connectomics-based measures are reported to have moderate to high test-retest reliability, the variability in data acquisition, experimental designs, and analytic methods precludes the ability to generalize results. The Consortium for Reliability and Reproducibility (CoRR) is working to address this challenge and establish test-retest reliability as a minimum standard for methods development in functional connectomics. Specifically, CoRR has aggregated 1,629 typical individuals’ resting state fMRI (rfMRI) data (5,093 rfMRI scans) from 18 international sites, and is openly sharing them via the International Data-sharing Neuroimaging Initiative (INDI). To allow researchers to generate various estimates of reliability and reproducibility, a variety of data acquisition procedures and experimental designs are included. Similarly, to enable users to assess the impact of commonly encountered artifacts (for example, motion) on characterizations of inter-individual variation, datasets of varying quality are included.

Tracking cerebral blood flow in BOLD fMRI using recursively generated regressors

BOLD functional MRI (fMRI) data are dominated by low frequency signals, many of them of unclear origin. We have recently shown that some portions of the low frequency oscillations found in BOLD fMRI are systemic signals closely related to the blood circulation (Tong et al. [2013]: NeuroImage 76:202-215). They are commonly treated as physiological noise in fMRI studies. In this study, we propose and test a novel data-driven analytical method that uses these systemic low frequency oscillations in the BOLD signal as a tracer to follow cerebral blood flow dynamically. Our findings demonstrate that: (1) systemic oscillations pervade the BOLD signal; (2) the temporal traces evolve as the blood propagates though the brain; and, (3) they can be effectively extracted via a recursive procedure and used to derive the cerebral circulation map. Moreover, this method is independent from functional analyses, and thus allows simultaneous and independent assessment of information about cerebral blood flow to be conducted in parallel with the functional studies. In this study, the method was applied to data from the resting state scans, acquired using a multiband EPI sequence (fMRI scan with much shorter TRs), of seven healthy participants. Dynamic maps with consistent features resembling cerebral blood circulation were derived, confirming the robustness and repeatability of the method.

Preoperative mapping of the sensorimotor cortex: comparative assessment of task-based and resting-state FMRI

Resting state fMRI (rs-fMRI) has recently been considered as a possible complement or alternative to task-based fMRI (tb-fMRI) for presurgical mapping. However, evidence of its usefulness remains scant, because existing studies have investigated relatively small samples and focused primarily on qualitative evaluation. The aim of this study is to investigate the clinical usefulness of rs-fMRI in the context of presurgical mapping of motor functions, and in particular to determine the degree of correspondence with tb-fMRI which, while not a gold-standard, is commonly used in preoperative setting. A group of 13 patients with lesions close to the sensorimotor cortex underwent rs-fMRI and tb-fMRI to localize the hand, foot and mouth motor areas. We assessed quantitatively the degree of correspondence between multiple rs-fMRI analyses (independent-component and seed-based analyses) and tb-fMRI, with reference to sensitivity and specificity of rs-fMRI with respect to tb-fMRI, and centre-of-mass distances. Agreement with electro-cortical stimulation (ECS) was also investigated, and a traditional map thresholding approach based on agreement between two experienced operators was compared to an automatic threshold determination method. Rs-fMRI can localize the sensorimotor cortex successfully, providing anatomical specificity for hand, foot and mouth motor subregions, in particular with seed-based analyses. Agreement with tb-fMRI was only partial and rs-fMRI tended to provide larger patterns of correlated activity. With respect to the ECS data available, rs-fMRI and tb-fMRI performed comparably, even though the shortest distance to stimulation points was observed for the latter. Notably, the results of both were on the whole robust to thresholding procedure. Localization performed by rs-fMRI is not equivalent to tb-fMRI, hence rs-fMRI cannot be considered as an outright replacement for tb-fMRI. Nevertheless, since there is significant agreement between the two techniques, rs-fMRI can be considered with caution as a potential alternative to tb-fMRI when patients are unable to perform the task.