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

Macrovascular blood flow and microvascular cerebrovascular reactivity are regionally coupled in adolescence

Cerebrovascular imaging assessments are particularly challenging in adolescent cohorts, where not all modalities are appropriate, and rapid brain maturation alters hemodynamics at both macro- and microvascular scales. In a preliminary sample of healthy adolescents (n=12, 8-25 years), we investigated relationships between 4D flow MRI-derived blood velocity and blood flow in bilateral anterior, middle, and posterior cerebral arteries and BOLD cerebrovascular reactivity in associated vascular territories. As hypothesized, higher velocities in large arteries are associated with an earlier response to a vasodilatory stimulus (cerebrovascular reactivity delay) in the downstream territory. Higher blood flow through these arteries is associated with a larger BOLD response to a vasodilatory stimulus (cerebrovascular reactivity amplitude) in the associated territory. These trends are consistent in a case study of adult moyamoya disease. In our small adolescent cohort, macrovascular-microvascular relationships for velocity/delay and flow/CVR change with age, though underlying mechanisms are unclear. Our work emphasizes the need to better characterize this key stage of human brain development, when cerebrovascular hemodynamics are changing, and standard imaging methods offer limited insight into these processes. We provide important normative data for future comparisons in pathology, where combining macro- and microvascular assessments may better help us prevent, stratify, and treat cerebrovascular disease.

MRI mapping of hemodynamics in the human spinal cord

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

Understanding recovery of language after stroke: insights from neurovascular MRI studies

Stroke causes a disruption in blood flow to the brain that can lead to profound language impairments. Understanding the mechanisms of language recovery after stroke is crucial for the prognosis and effective rehabilitation of people with aphasia. While the role of injured brain structures and disruptions in functional connectivity have been extensively explored, the relationship between neurovascular measures and language recovery in both early and later stages has not received sufficient attention in the field. Fully functioning healthy brain tissue requires oxygen and nutrients to be delivered promptly via its blood supply. Persistent decreases in blood flow after a stroke to the remaining non-lesioned tissue have been shown to contribute to poor language recovery. The goal of the current paper is to critically examine stroke studies looking at the relationship between different neurovascular measures and language deficits and mechanisms of language recovery via changes in neurovascular metrics. Measures of perfusion or cerebral blood flow (CBF) and cerebrovascular reactivity (CVR) provide complementary approaches to understanding neurovascular mechanisms post stroke by capturing both cerebral metabolic demands and mechanical vascular properties. While CBF measures indicate the amount of blood delivered to a certain region and serve as a proxy for metabolic demands of that area, CVR indices reflect the ability of the vasculature to recruit blood flow in response to a shortage of oxygen, such as when one is holding their breath. Increases in CBF during recovery beyond the site of the lesion have been shown to promote language gains. Similarly, CVR changes, when collateral vessels are recruited to help reorganize the flow of blood in hypoperfused regions, have been related to functional recovery post stroke. In the current review, we highlight the main findings in the literature investigating neurovascular changes in stroke recovery with a particular emphasis on how language abilities can be affected by changes in CBF and CVR. We conclude by summarizing existing methodological challenges and knowledge gaps that need to be addressed in future work in this area, outlining a promising avenue of research.

Correspondence between BOLD fMRI task response and cerebrovascular reactivity across the cerebral cortex

BOLD sensitivity to baseline perfusion and blood volume is a well-acknowledged fMRI confound. Vascular correction techniques based on cerebrovascular reactivity (CVR) might reduce variance due to baseline cerebral blood volume, however this is predicated on an invariant linear relationship between CVR and BOLD signal magnitude. Cognitive paradigms have relatively low signal, high variance and involve spatially heterogenous cortical regions; it is therefore unclear whether the BOLD response magnitude to complex paradigms can be predicted by CVR. The feasibility of predicting BOLD signal magnitude from CVR was explored in the present work across two experiments using different CVR approaches. The first utilized a large database containing breath-hold BOLD responses and 3 different cognitive tasks. The second experiment, in an independent sample, calculated CVR using the delivery of a fixed concentration of carbon dioxide and a different cognitive task. An atlas-based regression approach was implemented for both experiments to evaluate the shared variance between task-invoked BOLD responses and CVR across the cerebral cortex. Both experiments found significant relationships between CVR and task-based BOLD magnitude, with activation in the right cuneus (R ² = 0.64) and paracentral gyrus (R ² = 0.71), and the left pars opercularis (R ² = 0.67), superior frontal gyrus (R ² = 0.62) and inferior parietal cortex (R ² = 0.63) strongly predicted by CVR. The parietal regions bilaterally were highly consistent, with linear regressions significant in these regions for all four tasks. Group analyses showed that CVR correction increased BOLD sensitivity. Overall, this work suggests that BOLD signal response magnitudes to cognitive tasks are predicted by CVR across different regions of the cerebral cortex, providing support for the use of correction based on baseline vascular physiology.

Hemodynamic timing in resting-state and breathing-task BOLD fMRI

The blood flow response to a vasoactive stimulus demonstrates regional heterogeneity across both the healthy brain and in cerebrovascular pathology. The timing of a regional hemodynamic response is emerging as an important biomarker of cerebrovascular dysfunction, as well as a confound within fMRI analyses. Previous research demonstrated that hemodynamic timing is more robustly characterized when a larger systemic vascular response is evoked by a breathing challenge, compared to when only spontaneous fluctuations in vascular physiology are present (i.e., in resting-state data). However, it is not clear whether hemodynamic delays in these two conditions are physiologically interchangeable, and how methodological signal-to-noise factors may limit their agreement. To address this, we generated whole-brain maps of hemodynamic delays in nine healthy adults. We assessed the agreement of voxel-wise gray matter (GM) hemodynamic delays between two conditions: resting-state and breath-holding. We found that delay values demonstrated poor agreement when considering all GM voxels, but increasingly greater agreement when limiting analyses to voxels showing strong correlation with the GM mean time-series. Voxels showing the strongest agreement with the GM mean time-series were primarily located near large venous vessels, however these voxels explain some, but not all, of the observed agreement in timing. Increasing the degree of spatial smoothing of the fMRI data enhanced the correlation between individual voxel time-series and the GM mean time-series. These results suggest that signal-to-noise factors may be limiting the accuracy of voxel-wise timing estimates and hence their agreement between the two data segments. In conclusion, caution must be taken when using voxel-wise delay estimates from resting-state and breathing-task data interchangeably, and additional work is needed to evaluate their relative sensitivity and specificity to aspects of vascular physiology and pathology.

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 CO₂ (P_ETCO₂) 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 P_ETCO₂ 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 P_ETCO₂. We define a threshold to determine whether a P_ETCO₂ recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient P_ETCO₂ and 6 datasets with insufficient P_ETCO₂. When P_ETCO₂ quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from P_ETCO₂ (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 P_ETCO₂ with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when P_ETCO₂ 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 P_ETCO₂ 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.

Presurgical Magnetic Resonance Imaging Indicators of Revascularization Response in Adults With Moyamoya Vasculopathy

BACKGROUND

Moyamoya is a progressive intracranial vasculopathy, primarily affecting distal segments of the internal carotid and middle cerebral arteries. Treatment may comprise angiogenesis-inducing surgical revascularization; however, lack of randomized trials often results in subjective treatment decisions.

HYPOTHESIS

Compensatory presurgical posterior vertebrobasilar artery (VBA) flow-territory reactivity, including greater cerebrovascular reactivity (CVR) and reduced vascular delay time, portends greater neoangiogenic response verified on digital subtraction angiography (DSA) at 1-year follow-up.

STUDY TYPE

Prospective intervention cohort.

SUBJECTS

Thirty-one patients with moyamoya (26 females; age = 45 ± 13 years; 41 revascularized hemispheres).

METHODS

Anatomical MRI, hypercapnic CVR MRI, and DSA acquired presurgically in adult moyamoya participants scheduled for clinically indicated surgical revascularization. One-year postsurgery, DSA was repeated to evaluate collateralization.

FIELD STRENGTH

3 T.

SEQUENCE

Hypercapnic T 2 * -weighted gradient-echo blood-oxygenation-level-dependent, T2 -weighted turbo-spin-echo fluid-attenuated-inversion-recovery, T1 -weighted magnetization-prepared-rapid-gradient-echo, and T2 -weighted diffusion-weighted-imaging.

ASSESSMENT

Presurgical maximum CVR and response times were evaluated in VBA flow-territories. Revascularization success was determined using an ordinal scoring system of neoangiogenic collateralization from postsurgical DSA by two cerebrovascular neurosurgeons (R.V.C. with 8 years of experience; M.R.F. with 9 years of experience) and one neuroradiologist (L.T.D. with 8 years of experience). Stroke risk factors (age, sex, race, vasculopathy, and diabetes) were recorded.

STATISTICAL TESTS

Fisher's exact and Wilcoxon rank-sum tests were applied to compare presurgical variables between cohorts with angiographically confirmed good (>1/3 middle cerebral artery [MCA] territory revascularized) vs. poor (<1/3 MCA territory revascularized) outcomes.

SIGNIFICANCE

two-sided P < 0.05. Normalized odds ratios (ORs) were calculated.

RESULTS

Criteria for good collateralization were met in 25 of the 41 revascularized hemispheres. Presurgical normalized VBA flow-territory CVR was significantly higher in those with good (1.12 ± 0.13 unitless) vs. poor (1.04 ± 0.05 unitless) outcomes. Younger (OR = -0.60 ± 0.67) and White (OR = -1.81 ± 1.40) participants had highest revascularization success (good outcomes: age = 42 ± 14 years, race = 84% White; poor outcomes: age = 49 ± 11 years, race = 44% White).

DATA CONCLUSION

Presurgical MRI-measures of VBA flow-territory CVR are highest in moyamoya participants with better angiographic responses to surgical revascularization.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY STAGE: 4.

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.

Future Directions for Chemosensory Connectomes: Best Practices and Specific Challenges

Ecological chemosensory stimuli almost always evoke responses in more than one sensory system. Moreover, any sensory processing takes place along a hierarchy of brain regions. So far, the field of chemosensory neuroimaging is dominated by studies that examine the role of brain regions in isolation. However, to completely understand neural processing of chemosensation, we must also examine interactions between regions. In general, the use of connectivity methods has increased in the neuroimaging field, providing important insights to physical sensory processing, such as vision, audition, and touch. A similar trend has been observed in chemosensory neuroimaging, however, these established techniques have largely not been rigorously applied to imaging studies on the chemical senses, leaving network insights overlooked. In this article, we first highlight some recent work in chemosensory connectomics and we summarize different connectomics techniques. Then, we outline specific challenges for chemosensory connectome neuroimaging studies. Finally, we review best practices from the general connectomics and neuroimaging fields. We recommend future studies to develop or use the following methods we perceive as key to improve chemosensory connectomics: (1) optimized study designs, (2) reporting guidelines, (3) consensus on brain parcellations, (4) consortium research, and (5) data sharing.

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.