Daniel A. Handwerker

The impact of global signal regression on resting state correlations: Are anti-correlated networks introduced?

Low-frequency fluctuations in fMRI signal have been used to map several consistent resting state networks in the brain. Using the posterior cingulate cortex as a seed region, functional connectivity analyses have found not only positive correlations in the default mode network but negative correlations in another resting state network related to attentional processes. The interpretation is that the human brain is intrinsically organized into dynamic, anti-correlated functional networks. Global variations of the BOLD signal are often considered nuisance effects and are commonly removed using a general linear model (GLM) technique. This global signal regression method has been shown to introduce negative activation measures in standard fMRI analyses. The topic of this paper is whether such a correction technique could be the cause of anti-correlated resting state networks in functional connectivity analyses. Here we show that, after global signal regression, correlation values to a seed voxel must sum to a negative value. Simulations also show that small phase differences between regions can lead to spurious negative correlation values. A combination breath holding and visual task demonstrates that the relative phase of global and local signals can affect connectivity measures and that, experimentally, global signal regression leads to bell-shaped correlation value distributions, centred on zero. Finally, analyses of negatively correlated networks in resting state data show that global signal regression is most likely the cause of anti-correlations. These results call into question the interpretation of negatively correlated regions in the brain when using global signal regression as an initial processing step.

Dynamic functional connectivity: promise, issues, and interpretations.

The brain must dynamically integrate, coordinate, and respond to internal and external stimuli across multiple time scales. Non-invasive measurements of brain activity with fMRI have greatly advanced our understanding of the large-scale functional organization supporting these fundamental features of brain function. Conclusions from previous resting-state fMRI investigations were based upon static descriptions of functional connectivity (FC), and only recently studies have begun to capitalize on the wealth of information contained within the temporal features of spontaneous BOLD FC. Emerging evidence suggests that dynamic FC metrics may index changes in macroscopic neural activity patterns underlying critical aspects of cognition and behavior, though limitations with regard to analysis and interpretation remain. Here, we review recent findings, methodological considerations, neural and behavioral correlates, and future directions in the emerging field of dynamic FC investigations.

Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses

Estimates of hemodynamic response functions (HRF) are often integral parts of event-related fMRI analyses. Although HRFs vary across individuals and brain regions, few studies have investigated how variations affect the results of statistical analyses using the general linear model (GLM). In this study, we empirically estimated HRFs from primary motor and visual cortices and frontal and supplementary eye fields (SEF) in 20 subjects. We observed more variability across subjects than regions and correlated variation of time-to-peak values across several pairs of regions. Simulations examined the effects of observed variability on statistical results and ways different experimental designs and statistical models can limit these effects. Widely spaced and rapid event-related experimental designs with two sampling rates were tested. Statistical models compared an empirically derived HRF to a canonical HRF and included the first derivative of the HRF in the GLM. Small differences between the estimated and true HRFs did not cause false negatives, but larger differences within an observed range of variation, such as a 2.5-s time-to-onset misestimate, led to false negatives. Although small errors minimally affected detection of activity, time-to-onset misestimates as small as 1 s influenced model parameter estimation and therefore random effects analyses across subjects. Experiment and analysis design methods such as decreasing the sampling rate or including the HRFs temporal derivative in the GLM improved results, but did not eliminate errors caused by HRF misestimates. These results highlight the benefits of determining the best possible HRF estimate and potential negative consequences of assuming HRF consistency across subjects or brain regions.

Whole-brain, time-locked activation with simple tasks revealed using massive averaging and model-free analysis

The brain is the bodys largest energy consumer, even in the absence of demanding tasks. Electrophysiologists report on-going neuronal firing during stimulation or task in regions beyond those of primary relationship to the perturbation. Although the biological origin of consciousness remains elusive, it is argued that it emerges from complex, continuous whole-brain neuronal collaboration. Despite converging evidence suggesting the whole brain is continuously working and adapting to anticipate and actuate in response to the environment, over the last 20 y, task-based functional MRI (fMRI) have emphasized a localizationist view of brain function, with fMRI showing only a handful of activated regions in response to task/stimulation. Here, we challenge that view with evidence that under optimal noise conditions, fMRI activations extend well beyond areas of primary relationship to the task; and blood-oxygen level-dependent signal changes correlated with task-timing appear in over 95% of the brain for a simple visual stimulation plus attention control task. Moreover, we show that response shape varies substantially across regions, and that whole-brain parcellations based on those differences produce distributed clusters that are anatomically and functionally meaningful, symmetrical across hemispheres, and reproducible across subjects. These findings highlight the exquisite detail lying in fMRI signals beyond what is normally examined, and emphasize both the pervasiveness of false negatives, and how the sparseness of fMRI maps is not a result of localized brain function, but a consequence of high noise and overly strict predictive response models.

fMRI in the presence of task-correlated breathing variations

Variations in the subjects heart rate and breathing pattern have been shown to result in significant fMRI signal changes, mediated in part by non-neuronal physiological mechanisms such as global changes in levels of arterial CO(2). When these physiological changes are correlated with a task, as may happen in response to emotional stimuli or tasks that change levels of arousal, a concern arises that non-neuronal physiologically-induced signal changes may be misinterpreted as reflecting task-related neuronal activation. The purpose of this study is to provide information that can help in determining whether task activation maps are influenced by task-correlated physiological noise, particularly task-correlated breathing changes. We also compare different strategies to reduce the influence of physiological noise. Two paradigms are investigated--1) a lexical decision task where some subjects showed task-related breathing changes, and 2) a task where subjects were instructed to hold their breath during the presentation of contrast-reversing checkerboard, an extreme case of task-correlated physiological noise. Consistent with previous literature, we find that MRI signal changes correlated with variations in breathing depth and rate have a characteristic spatial and temporal profile that is different from the typical activation-induced BOLD response. The delineation of activation in the presence of task correlated breathing changes was improved either by independent component analysis, or by including specific nuisance regressors in a regression analysis. The difference in the spatial and temporal characteristics of physiological-induced and neuronal-induced fluctuations exploited by these strategies suggests that activation can be studied even in the presence of task-correlated physiological changes.

Tracking ongoing cognition in individuals using brief, whole-brain functional connectivity patterns

Significance Recently, it was shown that functional connectivity patterns exhibit complex spatiotemporal dynamics at the scale of tens of seconds. Of particular interest is the observation of a limited set of quasi-stable, whole-brain, recurring configurations—commonly referred to as functional connectivity states (FC states)—hypothesized to reflect the continuous flux of cognitive processes. Here, to test this hypothesis, subjects were continuously scanned as they engaged in and transitioned between mental states dictated by tasks. We demonstrate that there is a strong relationship between FC states and ongoing cognition that permits accurate tracking of mental states in individual subjects. We also demonstrate how informative changes in connectivity are not restricted solely to those regions with sustained elevations in activity during task performance. Functional connectivity (FC) patterns in functional MRI exhibit dynamic behavior on the scale of seconds, with rich spatiotemporal structure and limited sets of whole-brain, quasi-stable FC configurations (FC states) recurring across time and subjects. Based on previous evidence linking various aspects of cognition to group-level, minute-to-minute FC changes in localized connections, we hypothesized that whole-brain FC states may reflect the global, orchestrated dynamics of cognitive processing on the scale of seconds. To test this hypothesis, subjects were continuously scanned as they engaged in and transitioned between mental states dictated by tasks. FC states computed within windows as short as 22.5 s permitted robust tracking of cognition in single subjects with near perfect accuracy. Accuracy dropped markedly for subjects with the lowest task performance. Spatially restricting FC information decreased accuracy at short time scales, emphasizing the distributed nature of whole-brain FC dynamics, beyond univariate magnitude changes, as valuable markers of cognition.

Reducing vascular variability of fMRI data across aging populations using a breathholding task.

The magnitude and shape of blood oxygen level‐dependent (BOLD) responses in functional MRI (fMRI) studies vary across brain regions, subjects, and populations. This variability may be secondary to neural activity or vasculature differences, thus complicating interpretations of BOLD signal changes in fMRI experiments. We compare the BOLD responses to neural activity and a vascular challenge and test a method to dissociate these influences in 26 younger subjects (ages 18–36) and 24 older subjects (ages 51–78). Each subject performed a visuomotor saccade task (a vascular response to neural activity) and a breathholding task (vascular dilation induced by hypercapnia) during separate runs in the same scanning session. For the saccade task, signal magnitude showed a significant decrease with aging in FEF, SEF, and V1, and a delayed time‐to‐peak with aging in V1. The signal magnitudes from the saccade and hypercapnia tasks showed significant linear regressions within subjects and across individuals and populations. These two tasks had weaker, but sometimes significant linear regressions for time‐to‐peak and coherence phase measures. The significant magnitude decrease with aging in V1 remained after dividing the saccade task magnitude by the hypercapnia task magnitude, implying that the signal decrease is neural in origin. These findings may lead to a method to identify vascular reactivity‐induced differences in the BOLD response across populations and the development of methods to account for the influence of these vasculature differences in a simple, noninvasive manner. Hum Brain Mapp 2006.

Long-term neural and physiological phenotyping of a single human

Psychiatric disorders are characterized by major fluctuations in psychological function over the course of weeks and months, but the dynamic characteristics of brain function over this timescale in healthy individuals are unknown. Here, as a proof of concept to address this question, we present the MyConnectome project. An intensive phenome-wide assessment of a single human was performed over a period of 18 months, including functional and structural brain connectivity using magnetic resonance imaging, psychological function and physical health, gene expression and metabolomics. A reproducible analysis workflow is provided, along with open access to the data and an online browser for results. We demonstrate dynamic changes in brain connectivity over the timescales of days to months, and relations between brain connectivity, gene expression and metabolites. This resource can serve as a testbed to study the joint dynamics of human brain and metabolic function over time, an approach that is critical for the development of precision medicine strategies for brain disorders.

The continuing challenge of understanding and modeling hemodynamic variation in fMRI

Interpretation of fMRI data depends on our ability to understand or model the shape of the hemodynamic response (HR) to a neural event. Although the HR has been studied almost since the beginning of fMRI, we are still far from having robust methods to account for the full range of known HR variation in typical fMRI analyses. This paper reviews how the authors and others contributed to our understanding of HR variation. We present an overview of studies that describe HR variation across voxels, healthy volunteers, populations, and dietary or pharmaceutical modulations. We also describe efforts to minimize the effects of HR variation in intrasubject, group, population, and connectivity analyses and the limits of these methods.

Connectivity trajectory across lifespan differentiates the precuneus from the default network

The default network of the human brain has drawn much attention due to its relevance to various brain disorders, cognition, and behavior. However, its functional components and boundaries have not been precisely defined. There is no consensus as to whether the precuneus, a hub in the functional connectome, acts as part of the default network. This discrepancy is more critical for brain development and aging studies: it is not clear whether age has a stronger impact on the default network or precuneus, or both. We used Generalized Ranking and Averaging Independent Component Analysis by Reproducibility (gRAICAR) to investigate the lifespan trajectories of intrinsic functional networks. By estimating individual-specific spatial components and aligning them across subjects, gRAICAR measures the spatial variation of component maps across a population without constraining the same components to appear in every subject. In a cross-lifespan fMRI dataset (N=126, 7-85years old), we observed stronger age dependence in the spatial pattern of a precuneus-dorsal posterior cingulate cortex network compared to the default network, despite the fact that the two networks exhibit considerable spatial overlap and temporal correlation. These results remained even when analyses were restricted to a subpopulation with very similar head motion across age. Our analyses further showed that the two networks tend to merge with increasing age. Post-hoc analyses of functional connectivity confirmed the distinguishable cross-lifespan trajectories between the two networks. Based on these observations, we proposed a dynamic model of cross-lifespan functional segregation and integration between the two networks, suggesting that the precuneus network may have a different functional role than the default network, which declines with age. These findings have implications for understanding the functional roles of the default network, gaining insight into its dynamics throughout life, and guiding interpretation of alterations in brain disorders.