Expansion of a frontostriatal salience network in individuals with depression.

Hundreds of neuroimaging studies spanning two decades have revealed differences in brain structure and functional connectivity in depression, but with modest effect sizes, complicating efforts to derive mechanistic pathophysiologic insights or develop biomarkers. 1 Furthermore, although depression is a fundamentally episodic condition, few neuroimaging studies have taken a longitudinal approach, which is critical for understanding cause and effect and delineating mechanisms that drive mood state transitions over time. The emerging field of precision functional mapping using densely-sampled longitudinal neuroimaging data has revealed unexpected, functionally meaningful individual differences in brain network topology in healthy individuals, 2-5 but these approaches have never been applied to individuals with depression. Here, using precision functional mapping techniques and 11 datasets comprising n=187 repeatedly sampled individuals and >21,000 minutes of fMRI data, we show that the frontostriatal salience network is expanded two-fold in most individuals with depression. This effect was replicable in multiple samples, including large-scale, group-average data (N=1,231 subjects), and caused primarily by network border shifts affecting specific functional systems, with three distinct modes of encroachment occurring in different individuals. Salience network expansion was unexpectedly stable over time, unaffected by changes in mood state, and detectable in children before the subsequent onset of depressive symptoms in adolescence. Longitudinal analyses of individuals scanned up to 62 times over 1.5 years identified connectivity changes in specific frontostriatal circuits that tracked fluctuations in specific symptom domains and predicted future anhedonia symptoms before they emerged. Together, these findings identify a stable trait-like brain network topology that may confer risk for depression and mood-state dependent connectivity changes in frontostriatal circuits that predict the emergence and remission of depressive symptoms over time.

Denoising task-correlated head motion from motor-task fMRI data with multi-echo ICA.

Motor-task functional magnetic resonance imaging (fMRI) is crucial in the study of several clinical conditions, including stroke and Parkinson's disease. However, motor-task fMRI is complicated by task-correlated head motion, which can be magnified in clinical populations and confounds motor activation results. One method that may mitigate this issue is multi-echo independent component analysis (ME-ICA), which has been shown to separate the effects of head motion from the desired BOLD signal but has not been tested in motor-task datasets with high amounts of motion. In this study, we collected an fMRI dataset from a healthy population who performed a hand grasp task with and without task-correlated amplified head motion to simulate a motor-impaired population. We analyzed these data using three models: single-echo (SE), multi-echo optimally combined (ME-OC), and ME-ICA. We compared the models' performance in mitigating the effects of head motion on the subject level and group level. On the subject level, ME-ICA better dissociated the effects of head motion from the BOLD signal and reduced noise. Both ME models led to increased t-statistics in brain motor regions. In scans with high levels of motion, ME-ICA additionally mitigated artifacts and increased stability of beta coefficient estimates, compared to SE. On the group level, all three models produced activation clusters in expected motor areas in scans with both low and high motion, indicating that group-level averaging may also sufficiently resolve motion artifacts that vary by subject. These findings demonstrate that ME-ICA is a useful tool for subject-level analysis of motor-task data with high levels of task-correlated head motion. The improvements afforded by ME-ICA are critical to improve reliability of subject-level activation maps for clinical populations in which group-level analysis may not be feasible or appropriate, for example in a chronic stroke cohort with varying stroke location and degree of tissue damage.

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI.

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO2 (PETCO2) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable PETCO2 recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with PETCO2. We define a threshold to determine whether a PETCO2 recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient PETCO2 and 6 datasets with insufficient PETCO2. When PETCO2 quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from PETCO2 (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to PETCO2 with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when PETCO2 quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient PETCO2 recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow.

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function.

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2-3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO2 pressure: a breath-hold task to induce hypercapnia (CO2 increase) and a cued deep breathing task to induce hypocapnia (CO2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO2 by systematically shifting the CO2 regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

Brainhack: Developing a culture of open, inclusive, community-driven neuroscience.

Brainhack is an innovative meeting format that promotes scientific collaboration and education in an open, inclusive environment. This NeuroView describes the myriad benefits for participants and the research community and how Brainhacks complement conventional formats to augment scientific progress.

ICA-based denoising strategies in breath-hold induced cerebrovascular reactivity mapping with multi echo BOLD fMRI.

Performing a BOLD functional MRI (fMRI) acquisition during breath-hold (BH) tasks is a non-invasive, robust method to estimate cerebrovascular reactivity (CVR). However, movement and breathing-related artefacts caused by the BH can substantially hinder CVR estimates due to their high temporal collinearity with the effect of interest, and attention has to be paid when choosing which analysis model should be applied to the data. In this study, we evaluate the performance of multiple analysis strategies based on lagged general linear models applied on multi-echo BOLD fMRI data, acquired in ten subjects performing a BH task during ten sessions, to obtain subject-specific CVR and haemodynamic lag estimates. The evaluated approaches range from conventional regression models, i.e. including drifts and motion timecourses as nuisance regressors, applied on single-echo or optimally-combined data, to more complex models including regressors obtained from multi-echo independent component analysis with different grades of orthogonalization in order to preserve the effect of interest, i.e. the CVR. We compare these models in terms of their ability to make signal intensity changes independent from motion, as well as the reliability as measured by voxelwise intraclass correlation coefficients of both CVR and lag maps over time. Our results reveal that a conservative independent component analysis model applied on the optimally-combined multi-echo fMRI signal offers the largest reduction of motion-related effects in the signal, while yielding reliable CVR amplitude and lag estimates, although a conventional regression model applied on the optimally-combined data results in similar estimates. This work demonstrates the usefulness of multi-echo based fMRI acquisitions and independent component analysis denoising for precision mapping of CVR in single subjects based on BH paradigms, fostering its potential as a clinically-viable neuroimaging tool for individual patients. It also proves that the way in which data-driven regressors should be incorporated in the analysis model is not straight-forward due to their complex interaction with the BH-induced BOLD response.

BIDScoin: A User-Friendly Application to Convert Source Data to Brain Imaging Data Structure.

Analyses of brain function and anatomy using shared neuroimaging data is an important development, and have acquired the potential to be scaled up with the specification of a new Brain Imaging Data Structure (BIDS) standard. To date, a variety of software tools help researchers in converting their source data to BIDS but often require programming skills or are tailored to specific institutes, data sets, or data formats. In this paper, we introduce BIDScoin, a cross-platform, flexible, and user-friendly converter that provides a graphical user interface (GUI) to help users finding their way in BIDS standard. BIDScoin does not require programming skills to be set up and used and supports plugins to extend their functionality. In this paper, we show its design and demonstrate how it can be applied to a downloadable tutorial data set. BIDScoin is distributed as free and open-source software to foster the community-driven effort to promote and facilitate the use of BIDS standard.

Voxelwise optimization of hemodynamic lags to improve regional CVR estimates in breath-hold fMRI.

Cerebrovascular Reactivity (CVR), the responsiveness of blood vessels to a vasodilatory stimulus, is an important indicator of cerebrovascular health. Assessing CVR with fMRI, we can measure the change in the Blood Oxygen Level Dependent (BOLD) response induced by a change in CO2 pressure (%BOLD/mmHg). However, there exists a temporal offset between the recorded CO2 pressure and the local BOLD response, due to both measurement and physiological delays. If this offset is not corrected for, voxel-wise CVR values will not be accurate. In this paper, we propose a framework for mapping hemodynamic lag in breath-hold fMRI data. As breath-hold tasks drive task-correlated head motion artifacts in BOLD fMRI data, our framework for lag estimation fits a model that includes polynomial terms and head motion parameters, as well as a shifted variant of the CO2 regressor (±9 s in 0.3 s increments), and the hemodynamic lag at each voxel is the shift producing the maximum total model R2 within physiological constraints. This approach is evaluated in 8 subjects with multi-echo fMRI data, resulting in robust maps of hemodynamic delay that show consistent regional variation across subjects, and improved contrast-to-noise compared to methods where motion regression is ignored or performed earlier in preprocessing.Clinical Relevance- We map hemodynamic lag using breathhold fMRI, providing insight into vascular transit times and improving the regional accuracy of cerebrovascular reactivity measurements.

A deconvolution algorithm for multi-echo functional MRI: Multi-echo Sparse Paradigm Free Mapping.

This work introduces a novel algorithm for deconvolution of the BOLD signal in multi-echo fMRI data: Multi-echo Sparse Paradigm Free Mapping (ME-SPFM). Assuming a linear dependence of the BOLD percent signal change on the echo time (TE) and using sparsity-promoting regularized least squares estimation, ME-SPFM yields voxelwise time-varying estimates of the changes in the apparent transverse relaxation (ΔR2⁎) without prior knowledge of the timings of individual BOLD events. Our results in multi-echo fMRI data collected during a multi-task event-related paradigm at 3 Tesla demonstrate that the maps of R2⁎ changes obtained with ME-SPFM at the times of the stimulus trials show high spatial and temporal concordance with the activation maps and BOLD signals obtained with standard model-based analysis. This method yields estimates of ΔR2⁎ having physiologically plausible values. Owing to its ability to blindly detect events, ME-SPFM also enables us to map ΔR2⁎ associated with spontaneous, transient BOLD responses occurring between trials. This framework is a step towards deciphering the dynamic nature of brain activity in naturalistic paradigms, resting-state or experimental paradigms with unknown timing of the BOLD events.