The Role of Cerebrovascular-Reactivity Mapping in Functional MRI: Calibrated fMRI and Resting-State fMRI

The significance and limited influence of cerebrovascular reactivity on age and sex effects in task- and resting-state brain activity

Functional MRI measures the blood-oxygen-level dependent signals, which provide an indirect measure of neural activity mediated by neurovascular responses. Cerebrovascular reactivity affects both task-induced and resting-state blood-oxygen-level dependent activity and may confound inter-individual effects, such as those related to aging and biological sex. We examined a large dataset containing breath-holding, checkerboard, and resting-state tasks. We used the breath-holding task to measure cerebrovascular reactivity, used the checkerboard task to obtain task-based activations, and quantified resting-state activity with amplitude of low-frequency fluctuations and regional homogeneity. We hypothesized that cerebrovascular reactivity would be correlated with blood-oxygen-level dependent measures and that accounting for these correlations would result in better estimates of age and sex effects. We found that cerebrovascular reactivity was correlated with checkerboard task activations in the visual cortex and with amplitude of low-frequency fluctuations and regional homogeneity in widespread fronto-parietal regions, as well as regions with large vessels. We also found significant age and sex effects in cerebrovascular reactivity, some of which overlapped with those observed in amplitude of low-frequency fluctuations and regional homogeneity. However, correcting for the effects of cerebrovascular reactivity had very limited influence on the estimates of age and sex. Our results highlight the limitations of accounting for cerebrovascular reactivity with the current breath-holding task.

Multifaceted neural and vascular pathologies after pediatric mild traumatic brain injury

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

Functional MRI for Acute Covert Consciousness: Emerging Data and Implementation Case Series

Although research studies have begun to demonstrate relationships between disorders of consciousness and brain network biomarkers, there are limited data on the practical aspects of obtaining such network biomarkers to potentially guide care. As the state of knowledge continues to evolve, guidelines from professional societies such as the American and European Academies of Neurology and many experts have advocated that the risk-benefit ratio for the assessment of network biomarkers has begun to favor their application toward potentially detecting covert consciousness. Given the lack of detailed operationalization guidance and the context of the ethical implications, herein we offer a roadmap based on local institutional experience with the implementation of functional MRI in the neonatal, pediatric, and adult intensive care units of our local government-supported health system. We provide a case-based demonstrative approach intended to review the current literature and to assist with the initiation of such services at other facilities.

Acquisition and processing methods of whole-brain layer-fMRI VASO and BOLD: The Kenshu dataset

Cortical depth-dependent functional magnetic resonance image (fMRI), also known as layer-fMRI, has the potential to capture directional neural information flow of brain computations within and across large-scale cortical brain networks. E.g., layer-fMRI can differentiate feedforward and feedback cortical input in hierarchically organized brain networks. Recent advancements in 3D-EPI sampling approaches and MR contrast generation strategies have allowed proof-of-principle studies showing that layer-fMRI can provide sufficient data quality for capturing laminar changes in functional connectivity. These studies have however not shown how reliable the signal is and how repeatable the respective results are. It is especially unclear whether whole-brain layer-fMRI functional connectivity protocols are widely applicable across common neuroscience-driven analysis approaches. Moreover, there are no established preprocessing fMRI methods that are optimized to work for whole-brain layer-fMRI datasets. In this work, we aimed to serve the field of layer-fMRI and build tools for future routine whole-brain layer-fMRI in application-based neuroscience research. We have developed publicly available sequences, acquisition protocols, and processing pipelines for whole-brain layer-fMRI. These protocols are validated across 60 hours of scanning in nine participants. Specifically, we identified and exploited methodological advancements for maximizing tSNR efficiency and test-retest reliability. We are sharing an extensive multi-modal whole-brain layer-fMRI dataset (20 scan hours of movie-watching in a single participant) for the purpose of benchmarking future method developments: The Kenshu dataset. With this dataset, we are also exemplifying the usefulness of whole brain layer-fMRI for commonly applied analysis approaches in modern cognitive neuroscience fMRI studies. This includes connectivity analyses, representational similarity matrix estimations, general linear model analyses, principal component analysis clustering, etc. We believe that this work paves the road for future routine measurements of directional functional connectivity across the entire brain.

The Significance and Limited Influence of Cerebrovascular Reactivity on Age and Sex Effects in Task- and Resting-State Brain Activity

Functional MRI (fMRI) measures the blood-oxygen-level dependent (BOLD) signals, which provide an indirect measure of neural activity mediated by neurovascular responses. Cerebrovascular reactivity affects both task-induced and resting-state BOLD activity and may confound inter-individual effects observed in BOLD-based measures, such as those related to aging and biological sex. To investigate this, we examined a large open-access fMRI dataset containing a breath-holding task, checkerboard task, and resting-state scans. We used the breath-holding task to measure cerebrovascular reactivity, used the checkerboard task to obtain task-based activations, and from the resting-state data, we quantified the resting-state amplitude of low-frequency fluctuations (ALFF), and resting-state regional homogeneity (ReHo). We hypothesized that cerebrovascular reactivity would be correlated with BOLD measures and that accounting for these correlations would result in better estimates of age and sex effects. Our analysis showed that cerebrovascular reactivity was correlated with checkerboard task activations in the visual cortex and with ALFF and ReHo in widespread fronto-parietal regions, as well as regions with large vessels. We also found significant age and sex effects in cerebrovascular reactivity, some of which overlapped with those observed in ALFF and ReHo scores. Finally, we demonstrated that correcting for the effects of cerebrovascular reactivity had very limited influence on the estimates of age and sex. Our results highlight the limitations of accounting for cerebrovascular reactivity with the current breath-holding task.

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.

Agreement in cerebrovascular reactivity assessed with diffuse correlation spectroscopy across experimental paradigms improves with short separation regression

Significance

Cerebrovascular reactivity (CVR), i.e., the ability of cerebral vasculature to dilate or constrict in response to vasoactive stimuli, is a biomarker of vascular health. Exogenous administration of inhaled carbon dioxide, i.e., hypercapnia (HC), remains the "gold-standard" intervention to assess CVR. More tolerable paradigms that enable CVR quantification when HC is difficult/contraindicated have been proposed. However, because these paradigms feature mechanistic differences in action, an assessment of agreement of these more tolerable paradigms to HC is needed.

Aim

We aim to determine the agreement of CVR assessed during HC, breath-hold (BH), and resting state (RS) paradigms.

Approach

Healthy adults were subject to HC, BH, and RS paradigms. End tidal carbon dioxide (EtCO2) and cerebral blood flow (CBF, assessed with diffuse correlation spectroscopy) were monitored continuously. CVR (%/mmHg) was quantified via linear regression of CBF versus EtCO2 or via a general linear model (GLM) that was used to minimize the influence of systemic and extracerebral signal contributions.

Results

Strong agreement ( CCC ≥ 0.69 ; R ≥ 0.76 ) among CVR paradigms was demonstrated when utilizing a GLM to regress out systemic/extracerebral signal contributions. Linear regression alone showed poor agreement across paradigms ( CCC ≤ 0.35 ; R ≤ 0.45 ).

Conclusions

More tolerable experimental paradigms coupled with regression of systemic/extracerebral signal contributions may offer a viable alternative to HC for assessing CVR.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

Generating dynamic carbon-dioxide traces from respiration-belt recordings: Feasibility using neural networks and application in functional magnetic resonance imaging

INTRODUCTION

In the context of functional magnetic resonance imaging (fMRI), carbon dioxide (CO2) is a well-known vasodilator that has been widely used to monitor and interrogate vascular physiology. Moreover, spontaneous fluctuations in end-tidal carbon dioxide (PETCO2) reflects changes in arterial CO2 and has been demonstrated as the largest physiological noise source for denoising the low-frequency range of the resting-state fMRI (rs-fMRI) signal. However, the majority of rs-fMRI studies do not involve CO2 recordings, and most often only heart rate and respiration are recorded. While the intrinsic link between these latter metrics and CO2 led to suggested possible analytical models, they have not been widely applied.

METHODS

In this proof-of-concept study, we propose a deep-learning (DL) approach to reconstruct CO2 and PETCO2 data from respiration waveforms in the resting state.

RESULTS

We demonstrate that the one-to-one mapping between respiration and CO2 recordings can be well predicted using fully convolutional networks (FCNs), achieving a Pearson correlation coefficient (r) of 0.946 ± 0.056 with the ground truth CO2. Moreover, dynamic PETCO2 can be successfully derived from the predicted CO2, achieving r of 0.512 ± 0.269 with the ground truth. Importantly, the FCN-based methods outperform previously proposed analytical methods. In addition, we provide guidelines for quality assurance of respiration recordings for the purposes of CO2 prediction.

DISCUSSION

Our results demonstrate that dynamic CO2 can be obtained from respiration-volume using neural networks, complementing the still few reports in DL of physiological fMRI signals, and paving the way for further research in DL based bio-signal processing.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

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.

Hemodynamic Imaging in Cerebral Diffuse Glioma-Part A: Concept, Differential Diagnosis and Tumor Grading

Diffuse gliomas are the most common primary malignant intracranial neoplasms. Aside from the challenges pertaining to their treatment-glioblastomas, in particular, have a dismal prognosis and are currently incurable-their pre-operative assessment using standard neuroimaging has several drawbacks, including broad differentials diagnosis, imprecise characterization of tumor subtype and definition of its infiltration in the surrounding brain parenchyma for accurate resection planning. As the pathophysiological alterations of tumor tissue are tightly linked to an aberrant vascularization, advanced hemodynamic imaging, in addition to other innovative approaches, has attracted considerable interest as a means to improve diffuse glioma characterization. In the present part A of our two-review series, the fundamental concepts, techniques and parameters of hemodynamic imaging are discussed in conjunction with their potential role in the differential diagnosis and grading of diffuse gliomas. In particular, recent evidence on dynamic susceptibility contrast, dynamic contrast-enhanced and arterial spin labeling magnetic resonance imaging are reviewed together with perfusion-computed tomography. While these techniques have provided encouraging results in terms of their sensitivity and specificity, the limitations deriving from a lack of standardized acquisition and processing have prevented their widespread clinical adoption, with current efforts aimed at overcoming the existing barriers.

Metabolism modulates network synchrony in the aging brain

Brain aging is associated with hypometabolism and global changes in functional connectivity. Using functional MRI (fMRI), we show that network synchrony, a collective property of brain activity, decreases with age. Applying quantitative methods from statistical physics, we provide a generative (Ising) model for these changes as a function of the average communication strength between brain regions. We find that older brains are closer to a critical point of this communication strength, in which even small changes in metabolism lead to abrupt changes in network synchrony. Finally, by experimentally modulating metabolic activity in younger adults, we show how metabolism alone-independent of other changes associated with aging-can provide a plausible candidate mechanism for marked reorganization of brain network topology.

Effects of Thresholding on Voxel-Wise Correspondence of Breath-Hold and Resting-State Maps of Cerebrovascular Reactivity

Functional magnetic resonance imaging for presurgical brain mapping enables neurosurgeons to identify viable tissue near a site of operable pathology which might be at risk of surgery-induced damage. However, focal brain pathology (e.g., tumors) may selectively disrupt neurovascular coupling while leaving the underlying neurons functionally intact. Such neurovascular uncoupling can result in false negatives on brain activation maps thereby compromising their use for surgical planning. One way to detect potential neurovascular uncoupling is to map cerebrovascular reactivity using either an active breath-hold challenge or a passive resting-state scan. The equivalence of these two methods has yet to be fully established, especially at a voxel level of resolution. To quantitatively compare breath-hold and resting-state maps of cerebrovascular reactivity, we first identified threshold settings that optimized coverage of gray matter while minimizing false responses in white matter. When so optimized, the resting-state metric had moderately better gray matter coverage and specificity. We then assessed the spatial correspondence between the two metrics within cortical gray matter, again, across a wide range of thresholds. Optimal spatial correspondence was strongly dependent on threshold settings which if improperly set tended to produce statistically biased maps. When optimized, the two CVR maps did have moderately good correspondence with each other (mean accuracy of 73.6%). Our results show that while the breath-hold and resting-state maps may appear qualitatively similar they are not quantitatively identical at a voxel level of resolution.