Cerebrovascular reactivity mapping in patients with low grade gliomas undergoing presurgical sensorimotor mapping with BOLD fMRI

Recommendations for quantitative cerebral perfusion MRI using multi-timepoint arterial spin labeling: Acquisition, quantification, and clinical applications

Accurate assessment of cerebral perfusion is vital for understanding the hemodynamic processes involved in various neurological disorders and guiding clinical decision-making. This guidelines article provides a comprehensive overview of quantitative perfusion imaging of the brain using multi-timepoint arterial spin labeling (ASL), along with recommendations for its acquisition and quantification. A major benefit of acquiring ASL data with multiple label durations and/or post-labeling delays (PLDs) is being able to account for the effect of variable arterial transit time (ATT) on quantitative perfusion values and additionally visualize the spatial pattern of ATT itself, providing valuable clinical insights. Although multi-timepoint data can be acquired in the same scan time as single-PLD data with comparable perfusion measurement precision, its acquisition and postprocessing presents challenges beyond single-PLD ASL, impeding widespread adoption. Building upon the 2015 ASL consensus article, this work highlights the protocol distinctions specific to multi-timepoint ASL and provides robust recommendations for acquiring high-quality data. Additionally, we propose an extended quantification model based on the 2015 consensus model and discuss relevant postprocessing options to enhance the analysis of multi-timepoint ASL data. Furthermore, we review the potential clinical applications where multi-timepoint ASL is expected to offer significant benefits. This article is part of a series published by the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group, aiming to guide and inspire the advancement and utilization of ASL beyond the scope of the 2015 consensus article.

Attenuation of the BOLD fMRI Signal and Changes in Functional Connectivity Affecting the Whole Brain in Presence of Brain Metastasis

To date, there are almost no investigations addressing functional connectivity (FC) in patients with brain metastases (BM). In this retrospective study, we investigate the influence of BM on hemodynamic brain signals derived from functional magnetic resonance imaging (fMRI) and FC. Motor-fMRI data of 29 patients with BM and 29 matched healthy controls were analyzed to assess percent signal changes (PSC) in the ROIs motor cortex, premotor cortex, and supplementary motor cortex and FC in the sensorimotor, default mode, and salience networks using Statistical Parametric Mapping (SPM12) and marsbar and CONN toolboxes. In the PSC analysis, an attenuation of the BOLD signal in the metastases-affected hemisphere compared to the contralateral hemisphere was significant only in the supplementary motor cortex during hand movement. In the FC analysis, we found alterations in patients' FC compared to controls in all examined networks, also in the hemisphere contralateral to the metastasis. This indicates a qualitative attenuation of the BOLD signal in the affected hemisphere and also that FC is altered by the presence of BM, similarly to what is known for primary brain tumors. This transformation is not only visible in the infiltrated hemisphere, but also in the contralateral one, suggesting an influence of BM beyond local damage.

Transient deoxyhemoglobin formation as a contrast for perfusion MRI studies in patients with brain tumors: a feasibility study

Background: Transient hypoxia-induced deoxyhemoglobin (dOHb) has recently been shown to represent a comparable contrast to gadolinium-based contrast agents for generating resting perfusion measures in healthy subjects. Here, we investigate the feasibility of translating this non-invasive approach to patients with brain tumors. Methods: A computer-controlled gas blender was used to induce transient precise isocapnic lung hypoxia and thereby transient arterial dOHb during echo-planar-imaging acquisition in a cohort of patients with different types of brain tumors (n = 9). We calculated relative cerebral blood volume (rCBV), cerebral blood flow (rCBF), and mean transit time (MTT) using a standard model-based analysis. The transient hypoxia induced-dOHb MRI perfusion maps were compared to available clinical DSC-MRI. Results: Transient hypoxia induced-dOHb based maps of resting perfusion displayed perfusion patterns consistent with underlying tumor histology and showed high spatial coherence to gadolinium-based DSC MR perfusion maps. Conclusion: Non-invasive transient hypoxia induced-dOHb was well-tolerated in patients with different types of brain tumors, and the generated rCBV, rCBF and MTT maps appear in good agreement with perfusion maps generated with gadolinium-based DSC MR perfusion.

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.

Breath-hold BOLD fMRI without CO2 sampling enables estimation of venous cerebral blood volume: potential use in normalization of stimulus-evoked BOLD fMRI data

BOLD fMRI signal has been used in conjunction with vasodilatory stimulation as a marker of cerebrovascular reactivity (CVR): the relative change in cerebral blood flow (CBF) arising from a unit change in the vasodilatory stimulus. Using numerical simulations, we demonstrate that the variability in the relative BOLD signal change induced by vasodilation is strongly influenced by the variability in deoxyhemoglobin-containing cerebral blood volume (CBV), as this source of variability is likely to be more prominent than that of CVR. It may, therefore, be more appropriate to describe the relative BOLD signal change induced by an isometabolic vasodilation as a proxy of deoxygenated CBV (CBVdHb) rather than CVR. With this in mind, a new method was implemented to map a marker of CBVdHb, termed BOLD-CBV, based on the normalization of voxel-wise BOLD signal variation by an estimate of the intravascular venous BOLD signal from voxels filled with venous blood. The intravascular venous BOLD signal variation, recorded during repeated breath-holding, was extracted from the superior sagittal sinus in a cohort of 27 healthy volunteers and used as a regressor across the whole brain, yielding maps of BOLD-CBV. In the same cohort, we demonstrated the potential use of BOLD-CBV for the normalization of stimulus-evoked BOLD fMRI by comparing group-level BOLD fMRI responses to a visuomotor learning task with and without the inclusion of voxel-wise vascular covariates of BOLD-CBV and the BOLD signal change per mmHg variation in end-tidal carbon dioxide (BOLD-CVR). The empirical measure of BOLD-CBV accounted for more between-subject variability in the motor task-induced BOLD responses than BOLD-CVR estimated from end-tidal carbon dioxide recordings. The new method can potentially increase the power of group fMRI studies by including a measure of vascular characteristics and has the strong practical advantage of not requiring experimental measurement of end-tidal carbon dioxide, unlike traditional methods to estimate BOLD-CVR. It also more closely represents a specific physiological characteristic of brain vasculature than BOLD-CVR, namely blood volume.

Modern preoperative imaging and functional mapping in patients with intracranial glioma

Magnetic resonance imaging (MRI) in therapy-naïve intracranial glioma is paramount for neuro-oncological diagnostics, and it provides images that are helpful for surgery planning and intraoperative guidance during tumor resection, including assessment of the involvement of functionally eloquent brain structures. This study reviews emerging MRI techniques to depict structural information, diffusion characteristics, perfusion alterations, and metabolism changes for advanced neuro-oncological imaging. In addition, it reflects current methods to map brain function close to a tumor, including functional MRI and navigated transcranial magnetic stimulation with derived function-based tractography of subcortical white matter pathways. We conclude that modern preoperative MRI in neuro-oncology offers a multitude of possibilities tailored to clinical needs, and advancements in scanner technology (e. g., parallel imaging for acceleration of acquisitions) make multi-sequence protocols increasingly feasible. Specifically, advanced MRI using a multi-sequence protocol enables noninvasive, image-based tumor grading and phenotyping in patients with glioma. Furthermore, the add-on use of preoperatively acquired MRI data in combination with functional mapping and tractography facilitates risk stratification and helps to avoid perioperative functional decline by providing individual information about the spatial location of functionally eloquent tissue in relation to the tumor mass. KEY POINTS::   · Advanced preoperative MRI allows for image-based tumor grading and phenotyping in glioma.. · Multi-sequence MRI protocols nowadays make it possible to assess various tumor characteristics (incl. perfusion, diffusion, and metabolism).. · Presurgical MRI in glioma is increasingly combined with functional mapping to identify and enclose individual functional areas.. · Advancements in scanner technology (e. g., parallel imaging) facilitate increasing application of dedicated multi-sequence imaging protocols.. CITATION FORMAT: · Sollmann N, Zhang H, Kloth C et al. Modern preoperative imaging and functional mapping in patients with intracranial glioma. Fortschr Röntgenstr 2023; 195: 989 - 1000.

Heterogeneous motor BOLD-fMRI responses in brain areas exhibiting negative BOLD cerebrovascular reactivity indicate that steal phenomenon does not always result from exhausted cerebrovascular reserve capacity

INTRODUCTION

Brain areas exhibiting negative blood oxygenation-level dependent cerebrovascular reactivity (BOLD-CVR) responses to carbon dioxide (CO2) are thought to suffer from a completely exhausted autoregulatory cerebrovascular reserve capacity and exhibit vascular steal phenomenon. If this assumption is correct, the presence of vascular steal phenomenon should subsequently result in an equal negative fMRI signal response during a motor-task based BOLD-fMRI study (increase in metabolism without an increase in cerebral blood flow due to exhausted reserve capacity) in otherwise functional brain tissue. To investigate this premise, the aim of this study was to further investigate motor-task based BOLD-fMRI signal responses in brain areas exhibiting negative BOLD-CVR.

MATERIAL AND METHODS

Seventy-one datasets of patients with cerebrovascular steno-occlusive disease without motor defects, who underwent a CO2-calibrated motor task-based BOLD-fMRI study with a fingertapping paradigm and a subsequent BOLD-CVR study with a precisely controlled CO2-challenge during the same MRI examination, were included. We compared BOLD-fMRI signal responses in the bilateral pre- and postcentral gyri - i.e. Region of Interest (ROI) with the corresponding BOLD-CVR in this ROI. The ROI was determined using a second level group analysis of the BOLD-fMRI task study of 42 healthy individuals undergoing the same study protocol.

RESULTS

An overall decrease in BOLD-CVR was associated with a decrease in BOLD-fMRI signal response within the ROI. For patients exhibiting negative BOLD-CVR, we found both positive and negative motor-task based BOLD-fMRI signal responses.

CONCLUSION

We show that the presence of negative BOLD-CVR responses to CO2 is associated with heterogeneous motor task-based BOLD-fMRI signal responses, where some patients show -more presumed- negative BOLD-fMRI signal responses, while other patient showed positive BOLD-fMRI signal responses. This finding may indicate that the autoregulatory vasodilatory reserve capacity does not always need to be completely exhausted for vascular steal phenomenon to occur.

Detection and Mitigation of Neurovascular Uncoupling in Brain Gliomas

Functional magnetic resonance imaging (fMRI) with blood oxygen level-dependent (BOLD) technique is useful for preoperative mapping of brain functional networks in tumor patients, providing reliable in vivo detection of eloquent cortex to help reduce the risk of postsurgical morbidity. BOLD task-based fMRI (tb-fMRI) is the most often used noninvasive method that can reliably map cortical networks, including those associated with sensorimotor, language, and visual functions. BOLD resting-state fMRI (rs-fMRI) is emerging as a promising ancillary tool for visualization of diverse functional networks. Although fMRI is a powerful tool that can be used as an adjunct for brain tumor surgery planning, it has some constraints that should be taken into consideration for proper clinical interpretation. BOLD fMRI interpretation may be limited by neurovascular uncoupling (NVU) induced by brain tumors. Cerebrovascular reactivity (CVR) mapping obtained using breath-hold methods is an effective method for evaluating NVU potential.

Dependence of resting-state-based cerebrovascular reactivity (CVR) mapping on spatial resolution

Cerebrovascular reactivity (CVR) is typically assessed with a carbon dioxide (CO2) stimulus combined with BOLD fMRI. Recently, resting-state (RS) BOLD fMRI has been shown capable of generating CVR maps, providing a potential for broader CVR applications in neuroimaging studies. However, prior RS-CVR studies have primarily been performed at a spatial resolution of 3-4 mm voxel sizes. It remains unknown whether RS-CVR can also be obtained at high-resolution without major degradation in image quality. In this study, we investigated RS-CVR mapping based on resting-state BOLD MRI across a range of spatial resolutions in a group of healthy subjects, in an effort to examine the feasibility of RS-CVR measurement at high resolution. Comparing the results of RS-CVR with the maps obtained by the conventional CO2-inhalation method, our results suggested that good CVR map quality can be obtained at a voxel size as small as 2 mm isotropic. Our results also showed that, RS-CVR maps revealed resolution-dependent sensitivity. However, even at a high resolution of 2 mm isotropic voxel size, the voxel-wise sensitivity is still greater than that of typical task-evoked fMRI. Scan duration affected the sensitivity of RS-CVR mapping, but had no significant effect on its accuracy. These findings suggest that RS-CVR mapping can be applied at a similar resolution as state-of-the-art fMRI studies, which will broaden the use of CVR mapping in basic science and clinical applications including retrospective analysis of previously collected fMRI data.

Beyond Broca's and Wernicke's: Functional Mapping of Ancillary Language Centers Prior to Brain Tumor Surgery

Functional MRI is a well-established tool used for pre-surgical planning to help the neurosurgeon have a roadmap of critical functional areas that should be avoided, if possible, during surgery to minimize morbidity for patients with brain tumors (though this also has applications for surgical resection of epileptogenic tissue and vascular lesions). This article reviews the locations of secondary language centers within the brain along with imaging findings to help improve our confidence in our knowledge on language lateralization. Brief overviews of these language centers and their contributions to the language networks will be discussed. These language centers include primary language centers of "Broca's Area" and "Wernicke's Area". However, there are multiple secondary language centers such as the dorsal lateral prefrontal cortex (DLPFC), frontal eye fields, pre- supplemental motor area (pre-SMA), Basal Temporal Language Area (BTLA), along with other areas of activation. Knowing these foci helps to increase self-assurance when discussing the nature of laterality with the neurosurgeon. By knowing secondary language centers for language lateralization, via fMRI, one can feel confident on providing neurosurgeon colleagues with appropriate information on the laterality of language in preparation for surgery.

Deep-learning-enabled brain hemodynamic mapping using resting-state fMRI

Cerebrovascular disease is a leading cause of death globally. Prevention and early intervention are known to be the most effective forms of its management. Non-invasive imaging methods hold great promises for early stratification, but at present lack the sensitivity for personalized prognosis. Resting-state functional magnetic resonance imaging (rs-fMRI), a powerful tool previously used for mapping neural activity, is available in most hospitals. Here we show that rs-fMRI can be used to map cerebral hemodynamic function and delineate impairment. By exploiting time variations in breathing pattern during rs-fMRI, deep learning enables reproducible mapping of cerebrovascular reactivity (CVR) and bolus arrival time (BAT) of the human brain using resting-state CO2 fluctuations as a natural "contrast media". The deep-learning network is trained with CVR and BAT maps obtained with a reference method of CO2-inhalation MRI, which includes data from young and older healthy subjects and patients with Moyamoya disease and brain tumors. We demonstrate the performance of deep-learning cerebrovascular mapping in the detection of vascular abnormalities, evaluation of revascularization effects, and vascular alterations in normal aging. In addition, cerebrovascular maps obtained with the proposed method exhibit excellent reproducibility in both healthy volunteers and stroke patients. Deep-learning resting-state vascular imaging has the potential to become a useful tool in clinical cerebrovascular imaging.

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.

Pre-Operative Functional Mapping in Patients with Brain Tumors by fMRI and MEG: Advantages and Disadvantages in the Use of One Technique over the Other

Pre-operative mapping of brain functions is crucial to plan neurosurgery and investigate potential plasticity processes. Due to its availability, functional magnetic resonance imaging (fMRI) is widely used for this purpose; on the other hand, the demanding cost and maintenance limit the use of magnetoencephalography (MEG), despite several studies reporting its accuracy in localizing brain functions of interest in patient populations. In this review paper, we discuss the strengths and weaknesses of both techniques, from a methodological perspective first; then, we scrutinized and commented on the findings from 16 studies, identified by a database search, that made pre-operative assessments using both techniques in patients with brain tumors. We commented on the results by accounting for study limitations associated with small sample sizes and variability in the used tasks. Overall, we found that, although some studies reported the superiority for MEG, the majority of them underlined the complementary use of these techniques and suggested assessment using both. Indeed, both fMRI and MEG present some disadvantages, although the development of novel devices and processing procedures has enabled ever more accurate assessments. In particular, the development of new, more feasible MEG devices will allow widespread availability of this technique and its routinely combined use with fMRI.

Task-Based and Resting-State Functional MRI in Observing Eloquent Cerebral Areas Personalized for Epilepsy and Surgical Oncology Patients: A Review of the Current Evidence

Functional magnetic resonance imaging (fMRI) is among the newest techniques of advanced neuroimaging that offer the opportunity for neuroradiologists, neurophysiologists, neuro-oncologists, and neurosurgeons to pre-operatively plan and manage different types of brain lesions. Furthermore, it plays a fundamental role in the personalized evaluation of patients with brain tumors or patients with an epileptic focus for preoperative planning. While the implementation of task-based fMRI has increased in recent years, the existing resources and evidence related to this technique are limited. We have, therefore, conducted a comprehensive review of the available resources to compile a detailed resource for physicians who specialize in managing patients with brain tumors and seizure disorders. This review contributes to the existing literature because it highlights the lack of studies on fMRI and its precise role and applicability in observing eloquent cerebral areas in surgical oncology and epilepsy patients, which we believe is underreported. Taking these considerations into account would help to better understand the role of this advanced neuroimaging technique and, ultimately, improve patient life expectancy and quality of life.

Magnetic Resonance Imaging of Primary Adult Brain Tumors: State of the Art and Future Perspectives

MRI is undoubtedly the cornerstone of brain tumor imaging, playing a key role in all phases of patient management, starting from diagnosis, through therapy planning, to treatment response and/or recurrence assessment. Currently, neuroimaging can describe morphologic and non-morphologic (functional, hemodynamic, metabolic, cellular, microstructural, and sometimes even genetic) characteristics of brain tumors, greatly contributing to diagnosis and follow-up. Knowing the technical aspects, strength and limits of each MR technique is crucial to correctly interpret MR brain studies and to address clinicians to the best treatment strategy. This article aimed to provide an overview of neuroimaging in the assessment of adult primary brain tumors. We started from the basilar role of conventional/morphological MR sequences, then analyzed, one by one, the non-morphological techniques, and finally highlighted future perspectives, such as radiomics and artificial intelligence.

Cerebrovascular Reactivity Mapping Using Resting-State Functional MRI in Patients With Gliomas

BACKGROUND

Recently, a data-driven regression analysis method was developed to utilize the resting-state (rs) blood oxygenation level-dependent signal for cerebrovascular reactivity (CVR) mapping (rs-CVR), which was previously optimized by comparing with the CO₂ inhalation-based method in health subjects and patients with neurovascular diseases.

PURPOSE

To investigate the agreement of rs-CVR and the CVR mapping with breath-hold MRI (bh-CVR) in patients with gliomas.

STUDY TYPE

Retrospective.

POPULATION

Twenty-five patients (12 males, 13 females; mean age ± SD, 48 ± 13 years) with gliomas.

FIELD STRENGTH/SEQUENCE

Dynamic T2*-weighted gradient-echo echo-planar imaging during a breath-hold paradigm and during the rs on a 3-T scanner.

ASSESSMENT

rs-CVR with various frequency ranges and resting-state fluctuation amplitude (RSFA) were assessed. The agreement between each rs-based CVR measurement and bh-CVR was determined by voxel-wise correlation and Dice coefficient in the whole brain, gray matter, and the lesion region of interest (ROI).

STATISTICAL TESTS

Voxel-wise Pearson correlation, Dice coefficient, Fisher Z-transformation, repeated-measure analysis of variance and post hoc test with Bonferroni correction, and nonparametric repeated-measure Friedman test and post hoc test with Bonferroni correction were used. Significance was set at P < 0.05.

RESULTS

Compared with bh-CVR, the highest correlations were found at the frequency bands of 0.04-0.08 Hz and 0.02-0.04 Hz for rs-CVR in both whole brain and the lesion ROI. RSFA had significantly lower correlations than did rs-CVR of 0.02-0.04 Hz and a wider frequency range (0-0.1164 Hz). Significantly higher correlations and Dice coefficient were found in normal tissues than in the lesion ROI for all three methods.

DATA CONCLUSION

The optimal frequency ranges for rs-CVR are determined by comparing with bh-CVR in patients with gliomas. The rs-CVR method outperformed the RSFA. Significantly higher correlation and Dice coefficient between rs- and bh-CVR were found in normal tissue than in the lesion.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 2.

CVR-MRICloud: An online processing tool for CO2-inhalation and resting-state cerebrovascular reactivity (CVR) MRI data

Cerebrovascular Reactivity (CVR) provides an assessment of the brain’s vascular reserve and has been postulated to be a sensitive marker in cerebrovascular diseases. MRI-based CVR measurement typically employs alterations in arterial carbon dioxide (CO2) level while continuously acquiring Blood-Oxygenation-Level-Dependent (BOLD) images. CO2-inhalation and resting-state methods are two commonly used approaches for CVR MRI. However, processing of CVR MRI data often requires special expertise and may become an obstacle in broad utilization of this promising technique. The aim of this work was to develop CVR-MRICloud, a cloud-based CVR processing pipeline, to enable automated processing of CVR MRI data. The CVR-MRICloud consists of several major steps including extraction of end-tidal CO2 (EtCO2) curve from raw CO2 recording, alignment of EtCO2 curve with BOLD time course, computation of CVR value on a whole-brain, regional, and voxel-wise basis. The pipeline also includes standard BOLD image processing steps such as motion correction, registration between functional and anatomic images, and transformation of the CVR images to canonical space. This paper describes these algorithms and demonstrates the performance of the CVR-MRICloud in lifespan healthy subjects and patients with clinical conditions such as stroke, brain tumor, and Moyamoya disease. CVR-MRICloud has potential to be used as a data processing tool for a variety of basic science and clinical applications.

Hemispheric Asymmetry of the Hand Motor Representations in Patients with Highly Malignant Brain Tumors: Implications for Surgery and Clinical Practice

We addressed both brain pre-surgical functional and neurophysiological aspects of the hand representation in 18 right-handed patients harboring a highly malignant brain tumor in the sensorimotor (SM) cortex (10 in the left hemisphere, LH, and 8 in the right hemisphere, RH) and 10 healthy controls, who performed an fMRI hand-clenching task with both hands alternatively. We extracted the main ROI in the SM cortex and compared ROI values and volumes between hemispheres and groups, in addition to their motor neurophysiological measures. Hemispheric asymmetry in the fMRI signal was observed for healthy controls, namely higher signal for the left-hand movements, but not for either patients’ groups. ROI values, although altered in patients vs. controls, did not differ significantly between groups. ROI volumes associated with right-hand movement were lower for both patients’ groups vs. controls, and those associated with left-hand movement were lower in the RH group vs. all groups. These results are relevant to interpret potential preoperative plasticity and make inferences about postoperative plasticity and can be integrated in the surgical planning to increase surgery success and postoperative prognosis and quality of life.

Cerebrovascular Dysregulation in Patients with Glioma Assessed with Time-shifted BOLD fMRI

Background Microscopic vascular events, such as neovascularization and neurovascular uncoupling, are common in cerebral glioma. Mapping the cerebrovascular network remodeling at the macroscopic level may provide an alternative approach to assess hemodynamic dysregulation in patients with glioma. Purpose To investigate cerebrovascular dynamics and their relevance to tumor aggressiveness by using time-shift analysis (TSA) of the systemic low-frequency oscillation (sLFO) of the resting-state blood oxygenation level-dependent signal and a decision tree model. Materials and Methods In this retrospective study, 96 patients with histologically confirmed cerebral glioma were consecutively included (March 2012 to February 2017). TSA was performed to quantify the temporal properties of sLFO signals. Alteration in the time-shift properties was assessed in the tumor region and the contralesional hemisphere relative to the brains of healthy controls by using the Mann-Whitney U test. A decision tree model based on time-shift features was developed to predict the World Health Organization (WHO) glioma grade. Results A total of 88 patients with glioma (WHO grade II, 45; grade III, 21; grade IV, 22; mean age, 42 years; age range, 20-73 years; 51 men) and 40 healthy individuals from the 1000 Functional Connectomes Project (mean age, 32 years; age range, 24-49 years; 19 men) were included. The sLFO of the brain tissues was characterized by increased time shift in the tumor region and enhanced correlation with the global reference signal in the contralesional hemisphere compared with healthy brains. The proportion of tumor voxels with negative correlation to the reference signal significantly increased with the glioma malignancy grade. The decision tree model achieved an accuracy of 91% (80 of 88 patients) in predicting the glioma malignancy grade at the individual level (P = .004) based on the time-shift features. Conclusion Gliomas induced grade-specific cerebrovascular dysregulation in the entire brain, with altered time-shift features of systemic low-frequency oscillation signals. © RSNA, 2022 Online supplemental material is available for this article.

Association between tissue hypoxia, perfusion restrictions, and microvascular architecture alterations with lesion-induced impairment of neurovascular coupling

Functional magnetic resonance imaging (fMRI) has been mainly utilized for the preoperative localization of eloquent cortical areas. However, lesion-induced impairment of neurovascular coupling (NVC) in the lesion border zone may lead to false-negative fMRI results. The purpose of this study was to determine physiological factors impacting the NVC. Twenty patients suffering from brain lesions were preoperatively examined using multimodal neuroimaging including fMRI, magnetoencephalography (MEG) during language or sensorimotor tasks (depending on lesion location), and a novel physiologic MRI approach for the combined quantification of oxygen metabolism, perfusion state, and microvascular architecture. Congruence of brain activity patterns between fMRI and MEG were found in 13 patients. In contrast, we observed missing fMRI activity in perilesional cortex that demonstrated MEG activity in seven patients, which was interpreted as lesion-induced impairment of NVC. In these brain regions with impaired NVC, physiologic MRI revealed significant brain tissue hypoxia, as well as significantly decreased macro- and microvascular perfusion and microvascular architecture. We demonstrated that perilesional hypoxia with reduced vascular perfusion and architecture is associated with lesion-induced impairment of NVC. Our physiologic MRI approach is a clinically applicable method for preoperative risk assessment for the presence of false-negative fMRI results and may prevent severe postoperative functional deficits.

The Effect of Right Temporal Lobe Gliomas on Left and Right Hemisphere Neural Processing During Speech Perception and Production Tasks

Using fMRI, we investigated how right temporal lobe gliomas affecting the posterior superior temporal sulcus alter neural processing observed during speech perception and production tasks. Behavioural language testing showed that three pre-operative neurosurgical patients with grade 2, grade 3 or grade 4 tumours had the same pattern of mild language impairment in the domains of object naming and written word comprehension. When matching heard words for semantic relatedness (a speech perception task), these patients showed under-activation in the tumour infiltrated right superior temporal lobe compared to 61 neurotypical participants and 16 patients with tumours that preserved the right postero-superior temporal lobe, with enhanced activation within the (tumour-free) contralateral left superior temporal lobe. In contrast, when correctly naming objects (a speech production task), the patients with right postero-superior temporal lobe tumours showed higher activation than both control groups in the same right postero-superior temporal lobe region that was under-activated during auditory semantic matching. The task dependent pattern of under-activation during the auditory speech task and over-activation during object naming was also observed in eight stroke patients with right hemisphere infarcts that affected the right postero-superior temporal lobe compared to eight stroke patients with right hemisphere infarcts that spared it. These task-specific and site-specific cross-pathology effects highlight the importance of the right temporal lobe for language processing and motivate further study of how right temporal lobe tumours affect language performance and neural reorganisation. These findings may have important implications for surgical management of these patients, as knowledge of the regions showing functional reorganisation may help to avoid their inadvertent damage during neurosurgery.

Multi-vendor and multisite evaluation of cerebrovascular reactivity mapping using hypercapnia challenge

Cerebrovascular reactivity (CVR), which measures the ability of cerebral blood vessels to dilate or constrict in response to vasoactive stimuli such as CO2 inhalation, is an important index of the brain's vascular health. Quantification of CVR using BOLD MRI with hypercapnia challenge has shown great promises in research and clinical studies. However, in order for it to be used as a potential imaging biomarker in large-scale and multi-site studies, the reliability of CO2-CVR quantification across different MRI acquisition platforms and researchers/raters must be examined. The goal of this report from the MarkVCID small vessel disease biomarkers consortium is to evaluate the reliability of CO2-CVR quantification in three studies. First, the inter-rater reliability of CO2-CVR data processing was evaluated by having raters from 5 MarkVCID sites process the same 30 CVR datasets using a cloud-based CVR data processing pipeline. Second, the inter-scanner reproducibility of CO2-CVR quantification was assessed in 10 young subjects across two scanners of different vendors. Third, test-retest repeatability was evaluated in 20 elderly subjects from 4 sites with a scan interval of less than 2 weeks. In all studies, the CO2 CVR measurements were performed using the fixed inspiration method, where the subjects wore a nose clip and a mouthpiece and breathed room air and 5% CO2 air contained in a Douglas bag alternatively through their mouth. The results showed that the inter-rater CoV of CVR processing was 0.08 ± 0.08% for whole-brain CVR values and ranged from 0.16% to 0.88% in major brain regions, with ICC of absolute agreement above 0.9959 for all brain regions. Inter-scanner CoV was found to be 6.90 ± 5.08% for whole-brain CVR values, and ranged from 4.69% to 12.71% in major brain regions, which are comparable to intra-session CoVs obtained from the same scanners on the same day. ICC of consistency between the two scanners was 0.8498 for whole-brain CVR and ranged from 0.8052 to 0.9185 across major brain regions. In the test-retest evaluation, test-retest CoV across different days was found to be 18.29 ± 17.12% for whole-brain CVR values, and ranged from 16.58% to 19.52% in major brain regions, with ICC of absolute agreement ranged from 0.6480 to 0.7785. These results demonstrated good inter-rater, inter-scanner, and test-retest reliability in healthy volunteers, and suggested that CO2-CVR has suitable instrumental properties for use as an imaging biomarker of cerebrovascular function in multi-site and longitudinal observational studies and clinical trials.

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.

BOLD Asynchrony Elucidates Tumor Burden in IDH-Mutated Gliomas

BACKGROUND

Gliomas comprise the most common type of primary brain tumor, are highly invasive, and often fatal. IDH-mutated gliomas are particularly challenging to image and there is currently no clinically accepted method for identifying the extent of tumor burden in these neoplasms. This uncertainty poses a challenge to clinicians who must balance the need to treat the tumor while sparing healthy brain from iatrogenic damage. The purpose of this study was to investigate the feasibility of using resting-state blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) to detect glioma-related asynchrony in vascular dynamics for distinguishing tumor from healthy brain.

METHODS

Twenty-four stereotactically localized biopsies were obtained during open surgical resection from ten treatment-naïve patients with IDH-mutated gliomas who received standard of care preoperative imaging as well as echo-planar resting-state BOLD fMRI. Signal intensity for BOLD asynchrony and standard of care imaging was compared to cell counts of total cellularity (H&E), tumor density (IDH1 & Sox2), cellular proliferation (Ki67), and neuronal density (NeuN), for each corresponding sample.

RESULTS

BOLD asynchrony was directly related to total cellularity (H&E, p = 4 x 10 -5), tumor density (IDH1, p = 4 x 10 -5; Sox2, p = 3 x 10 -5), cellular proliferation (Ki67, p = 0.002), and as well as inversely related to neuronal density (NeuN, p = 1 x 10 -4).

CONCLUSIONS

Asynchrony in vascular dynamics, as measured by resting-state BOLD fMRI, correlates with tumor burden and provides a radiographic delineation of tumor boundaries in IDH-mutated gliomas.

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 CO₂ 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 CO₂ pressure: a breath-hold task to induce hypercapnia (CO₂ increase) and a cued deep breathing task to induce hypocapnia (CO₂ decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO₂ by systematically shifting the CO₂ 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 CO₂ 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.

fMRI Retinotopic Mapping in Patients with Brain Tumors and Space-Occupying Brain Lesions in the Area of the Occipital Lobe

Simple Summary

Functional magnetic resonance imaging (fMRI) in patients with brain tumors enables the visualization of eloquent cortical areas and can be used for planning surgical interventions and assessing the risk of postoperative functional deficits. While preoperative fMRI paradigms used to determine the localization of speech-critical or motor areas dominate the literature, there are hardly any studies that investigate the retinotopic organization of the visual field in patients with occipital lesions or tumors. The aim of this study was to evaluate the effect of a brain tumor or space-occupying brain lesions on the retinotopic organization of the occipital cortex, the activation of and the functional connectivity between cortical areas involved in visual processing. We found a high degree of similarity in the activation profiles of patients and healthy controls, indicating that the retinotopic organization of the visual cortex can reliably be described by fMRI retinotopic mapping as part of the preoperative examination of patients with tumors and space-occupying brain lesions.

Abstract

Functional magnetic resonance imaging (fMRI) is a valuable tool in the clinical routine of neurosurgery when planning surgical interventions and assessing the risk of postoperative functional deficits. Here, we examined how the presence of a brain tumor or lesion in the area of the occipital lobe affects the results of fMRI retinotopic mapping. fMRI data were evaluated on a retrospectively selected sample of 12 patients with occipital brain tumors, 7 patients with brain lesions and 19 control subjects. Analyses of the cortical activation, percent signal change, cluster size of the activated voxels and functional connectivity were carried out using Statistical Parametric Mapping (SPM12) and the CONN and Marsbar toolboxes. We found similar but reduced patterns of cortical activation and functional connectivity between the two patient groups compared to a healthy control group. Here, we found that retinotopic organization was well-preserved in the patients and was comparable to that of the age-matched controls. The results also showed that, compared to the tumor patients, the lesion patients showed higher percent signal changes but lower values in the cluster sizes of the activated voxels in the calcarine fissure region. Our results suggest that the lesion patients exhibited results that were more similar to those of the control subjects in terms of the BOLD signal, whereas the extent of the activation was comparable to that of the tumor patients.

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.

Cerebrovascular Reactivity Measurement Using Magnetic Resonance Imaging: A Systematic Review

Cerebrovascular reactivity (CVR) magnetic resonance imaging (MRI) probes cerebral haemodynamic changes in response to a vasodilatory stimulus. CVR closely relates to the health of the vasculature and is therefore a key parameter for studying cerebrovascular diseases such as stroke, small vessel disease and dementias. MRI allows in vivo measurement of CVR but several different methods have been presented in the literature, differing in pulse sequence, hardware requirements, stimulus and image processing technique. We systematically reviewed publications measuring CVR using MRI up to June 2020, identifying 235 relevant papers. We summarised the acquisition methods, experimental parameters, hardware and CVR quantification approaches used, clinical populations investigated, and corresponding summary CVR measures. CVR was investigated in many pathologies such as steno-occlusive diseases, dementia and small vessel disease and is generally lower in patients than in healthy controls. Blood oxygen level dependent (BOLD) acquisitions with fixed inspired CO₂ gas or end-tidal CO₂ forcing stimulus are the most commonly used methods. General linear modelling of the MRI signal with end-tidal CO₂ as the regressor is the most frequently used method to compute CVR. Our survey of CVR measurement approaches and applications will help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Perinidal Angiogenesis Is a Predictor for Neurovascular Uncoupling in the Periphery of Brain Arteriovenous Malformations: A Task-Based and Resting-State fMRI Study

BACKGROUND

Potential neurovascular uncoupling (NVU) related to perinidal angiogenesis (PA) of brain arteriovenous malformations (AVMs) may cause inappropriate presurgical mapping using functional magnetic resonance imaging (fMRI), resulting in overconfident resection and postoperative morbidity.

PURPOSE

To evaluate the potential impact of PA upon fMRI blood oxygen level-dependent signal in the periphery of AVMs.

STUDY TYPE

Prospective.

POPULATION

Twenty-one patients with AVMs located in the primary sensorimotor cortex (SM1) undergoing task-based fMRI (hand motor), and 19 patients with supratentorial AVMs undergoing resting-state fMRI.

FIELD STRENGTH/SEQUENCE

3.0T, echo-planar, time-of-flight, and magnetization-prepared rapid gradient-echo.

ASSESSMENT

The presence of PA was determined by three observers (Y.C., J.W., and X.C.) according to digital subtraction angiography and MR angiography. Interhemispheric asymmetry of fMRI activations contralateral to hand movements was evaluated with the interhemispheric ratio of the average t-value within ipsilesional SM1 to contralesional SM1. Regional homogeneity (ReHo) and fractional amplitude of low-frequency fluctuations (fALFF) were extracted from ring-shaped perinidal regions and contralateral homologous regions, and the corresponding interhemispheric ratios were calculated. The effect of PA on the interhemispheric asymmetry of motor activations, ReHo, and fALFF was estimated.

STATISTICAL TESTS

Pearson analysis, paired and independent t-test, multiple linear regression, Friedman test, and factorial analysis of variance were used.

RESULTS

Motor activations were significantly reduced in ipsilesional SM1 compared to contralesional SM1 (P < 0.05). The presence of PA was the independent predictor of activation loss in ipsilateral SM1(P < 0.05). Furthermore, perinidal regions exhibited reduced ReHo compared to the homologous regions (P < 0.05). PA was significantly associated with the decline of ReHo and fALFF in perinidal regions (P < 0.05, for both).

DATA CONCLUSION

The presence of PA can predict perinidal NVU that may confound the interpretation of both task-based and resting-state fMRI, highlighting the importance of alternative approaches of brain functional localization in improving treatment of AVMs.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

Relationships of Language Lateralization with Diffusion Tensor Imaging Metrics of Corpus Callosum, Tumor Grade, and Tumors Distance to Language-Eloquent Areas in Glial Neoplasms

PURPOSE

The aim of the study was to search relationships between language lateralization and corpus callosum (CC) connectivity, tumor grade, and tumors distance to language-eloquent areas in glial neoplasms.

MATERIALS AND METHODS

The functional magnetic resonance imaging and CC diffusion tensor imaging (DTI) metrics of 42 patients with glial neoplasm were evaluated for relationships of language lateralization (left, right, and bilateral) with CC DTI metrics (tracts number, voxel, volume, length, fractional anisotropy [FA], and apparent diffusion coefficient), tumor grade, and tumors distance to language-eloquent areas and relationships of CC DTI metrics with tumor grade. Kruskal-Wallis, Mann-Whitney U, and χ tests were used. All were repeated in 26 patients with left hemispheric masses.

RESULTS

In glial masses, language bilateralism was more common than normal population and more pronounced in low grade than high grade. In right lateralism and bilateralism, tumor settlement nearby language-eloquent areas was more common. In the left lateralism, highest CC tract number, higher tumor grade, and more remote tumor settlements were noted. There was no relationship between CC DTI metrics and tumor grade but increase in CC tracts number and FA with increasing tumor grade.

CONCLUSIONS

Increased bilateralism in glial masses than normal population and in low grade tumors than high grade and increased nearby tumor settlement in right lateralism and bilateralism support interhemispheric reorganization and plasticity. This is more pronounced in low grade because of higher life span. Highest CC tract number, higher tumor grade, and more remote tumor settlement in left lateralized group suggest intact CC integrity with limited hemispheric destruction. Increasing CC tracts number and FA with increasing tumor grade support preserved CC integrity in the shorter life span of high-grade tumors.

The Problem of Neurovascular Uncoupling

Neurovascular uncoupling (NVU) is one of the most important confounds of blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMR imaging) in the setting of focal brain lesions such as brain tumors. This article reviews the assessment of NVU related to focal brain lesions with emphasis on the use of cerebrovascular reactivity mapping measurement methods and resting state BOLD fMR imaging metrics in the detection of NVU, as well as the use of amplitude of low-frequency fluctuation metrics to mitigate the effects of NVU on clinical fMR imaging activation.

Emerging Techniques in Imaging of Glioma Microenvironment

Magnetic resonance imaging (MRI) has been the cornerstone of imaging of brain tumors in the past 4 decades. Conventional MRI remains the workhorse for neuro-oncologic imaging, not only for basic information such as location, extent, and navigation but also able to provide information regarding proliferation and infiltration, angiogenesis, hemorrhage, and more. More sophisticated MRI sequences have extended the ability to assess and quantify these features; for example, permeability and perfusion acquisitions can assess blood-brain barrier disruption and angiogenesis, diffusion techniques can assess cellularity and infiltration, and spectroscopy can address metabolism. Techniques such as fMRI and diffusion fiber tracking can be helpful in diagnostic planning for resection and radiation therapy, and more sophisticated iterations of these techniques can extend our understanding of neurocognitive effects of these tumors and associated treatment responses and effects. More recently, MRI has been used to go beyond such morphological, physiological, and functional characteristics to assess the tumor microenvironment. The current review highlights multiple recent and emerging approaches in MRI to characterize the tumor microenvironment.

A Novel Sequence: ZOOMit-Blood Oxygen Level-Dependent for Motor-Cortex Localization

BACKGROUND

Use of conventional blood oxygen level-dependent functional magnetic resonance imaging (conventional-BOLD-fMRI) presents challenges in accurately identifying the hand-motor cortex when a glioma involves the ipsilateral hand-knob. Zoomed imaging technique with parallel transmission (ZOOMit)-BOLD is a novel sequence allowing high spatial resolution with a relatively small field of view that may solve this problem.

OBJECTIVE

To compare the accuracy of ZOOMit-BOLD and conventional-BOLD in hand-motor cortex identification.

METHODS

A total of 20 patients with gliomas involving the sensorimotor cortex were recruited to identify the hand-motor cortex by both ZOOMit-BOLD and conventional-BOLD. Based on whether the entire or partial glioma directly invaded (was located within) the hand-knob or indirectly affected it by proximity, patients were placed into the involved or uninvolved groups, respectively. Direct cortical stimulation was applied intraoperatively to verify the location of the hand-motor cortex. Overlap indices were used to evaluate the accuracy of the hand-motor cortex identification. An overlap index equal to 0, indicating lack of overlap, was classified as inaccurate classification.

RESULTS

The accuracy of motor-cortex identification with ZOOMit-BOLD was 100% compared to only 65% with conventional-BOLD. The average overlap index yielded by ZOOMit-BOLD was higher than that of conventional-BOLD, regardless of whether gliomas directly invaded the hand-knob (P = .008) or not (P = .004). The overlap index in the involved group was significantly lower than that in the uninvolved group with both ZOOMit-BOLD (P = .002) and conventional-BOLD (P < .001).

CONCLUSION

ZOOMit-BOLD may potentially replace conventional-BOLD to identify the hand-motor cortex, particularly in cases in which gliomas directly invade the hand-knob.

Glioma-Induced Alterations in Neuronal Activity and Neurovascular Coupling during Disease Progression

Diffusely infiltrating gliomas are known to cause alterations in cortical function, vascular disruption, and seizures. These neurological complications present major clinical challenges, yet their underlying mechanisms and causal relationships to disease progression are poorly characterized. Here, we follow glioma progression in awake Thy1-GCaMP6f mice using in vivo wide-field optical mapping to monitor alterations in both neuronal activity and functional hemodynamics. The bilateral synchrony of spontaneous neuronal activity gradually decreases in glioma-infiltrated cortical regions, while neurovascular coupling becomes progressively disrupted compared to uninvolved cortex. Over time, mice develop diverse patterns of high amplitude discharges and eventually generalized seizures that appear to originate at the tumors' infiltrative margins. Interictal and seizure events exhibit positive neurovascular coupling in uninfiltrated cortex; however, glioma-infiltrated regions exhibit disrupted hemodynamic responses driving seizure-evoked hypoxia. These results reveal a landscape of complex physiological interactions occurring during glioma progression and present new opportunities for exploring novel biomarkers and therapeutic targets.

Cerebrovascular reactivity mapping using intermittent breath modulation

Cerebrovascular reactivity (CVR), an index of brain vessel's dilatory capacity, is typically measured using hypercapnic gas inhalation or breath-holding as a vasoactive challenge. However, these methods require considerable subject cooperation and could be challenging in clinical studies. More recently, there have been attempts to use resting-state BOLD data to map CVR by utilizing spontaneous changes in breathing pattern. However, in subjects who have small fluctuations in their spontaneous breathing pattern, the CVR results could be noisy and unreliable. In this study, we aim to develop a new method for CVR mapping that does not require gas-inhalation yet provides substantially higher sensitivity than resting-state CVR mapping. This new method is largely based on resting-state scan, but introduces intermittent modulation of breathing pattern in the subject to enhance fluctuations in their end-tidal CO2 (EtCO2) level. Here we examined the comfort level, sensitivity, and accuracy of this method in two studies. First, in 8 healthy young subjects, we developed the intermittent breath-modulation method using two different modulation frequencies, 6 ​s per breath and 12 ​s per breath, respectively, and compared the results to three existing CVR methods, specifically hypercapnic gas inhalation, breath-holding, and resting-state. Our results showed that the comfort level of the 6-s breath-modulation method was significantly higher than breath-holding (p ​= ​0.007) and CO2-inhalation (p ​= ​0.015) methods, while not different from the resting-state, i.e. free breathing method (p ​= ​0.52). When comparing the sensitivity of CVR methods, the breath-modulation methods revealed higher Z-statistics compared to the resting-state scan (p ​< ​0.008) and was comparable to breath-holding results. Next, we tested the feasibility of breath-modulation CVR mapping (6 ​s per breath) in 21 cognitively normal elderly participants and compared quantitative CVR values to that obtained with the CO2-inhalation method. Whole-brain CVR was found to be 0.150 ​± ​0.055 and 0.154 ​± ​0.032 %ΔBOLD/mmHg for the breath-modulation and CO2-inhalation method, respectively, with a significant correlation between them (y ​= ​0.97x, p ​= ​0.007). CVR mapping with intermittent breath modulation may be a useful method that combines the advantages of resting-state and CO2-inhalation based approaches.

Decreased Hand Motor Resting-State Functional Connectivity in Patients with Glioma: Analysis of Factors including Neurovascular Uncoupling

Background Resting-state functional MRI holds substantial potential for clinical application, but limitations exist in current understanding of how tumors exert local effects on resting-state functional MRI readings. Purpose To investigate the association between tumors, tumor characteristics, and changes in resting-state connectivity, to explore neurovascular uncoupling as a mechanism underlying these changes, and to evaluate seeding methodologies as a clinical tool. Materials and Methods Institutional review board approval was obtained for this HIPAA-compliant observational retrospective study of patients with glioma who underwent MRI and resting-state functional MRI between January 2016 and July 2017. Interhemispheric symmetry of connectivity was assessed in the hand motor region, incorporating tumor position, perfusion, grade, and connectivity generated from seed-based correlation. Statistical analysis was performed by using one-tailed t tests, Wilcoxon rank sum tests, one-way analysis of variance, Pearson correlation, and Spearman rank correlation, with significance at P < .05. Results Data in a total of 45 patients with glioma (mean age, 51.3 years ± 14.3 [standard deviation]) were compared with those in 10 healthy control subjects (mean age, 50.3 years ± 17.2). Patients showed loss of symmetry in measures of hand motor resting-state connectivity compared with control subjects (P < .05). Tumor distance from the ipsilateral hand motor (IHM) region correlated with the degree (R = 0.38, P = .01) and strength (R = 0.33, P = .03) of resting-state connectivity. In patients with World Health Organization grade IV glioblastomas 40 mm or less from the IHM region, loss of symmetry in strength of resting-state connectivity was correlated with tumor perfusion (R = 0.74, P < .01). In patients with gliomas 40 mm or less from the IHM region, seeding the nontumor hemisphere yielded less asymmetric hand motor resting-state connectivity than seeding the tumor hemisphere (connectivity seeded:contralateral = 1.34 nontumor vs 1.38 tumor hemisphere seeded; P = .03, false discovery rate threshold = 0.01). Conclusion Hand motor resting-state connectivity was less symmetrical in a tumor distance-dependent manner in patients with glioma. Differences in resting-state connectivity may be false-negative results driven by a neurovascular uncoupling mechanism. Seeding from the nontumor hemisphere may attenuate asymmetry in patients with tumors near ipsilateral hand motor cortices. © RSNA, 2020 Online supplemental material is available for this article.

Neurovascular Unit: Basic and Clinical Imaging with Emphasis on Advantages of Ferumoxytol

Physiological and pathological processes that increase or decrease the central nervous system's need for nutrients and oxygen via changes in local blood supply act primarily at the level of the neurovascular unit (NVU). The NVU consists of endothelial cells, associated blood-brain barrier tight junctions, basal lamina, pericytes, and parenchymal cells, including astrocytes, neurons, and interneurons. Knowledge of the NVU is essential for interpretation of central nervous system physiology and pathology as revealed by conventional and advanced imaging techniques. This article reviews current strategies for interrogating the NVU, focusing on vascular permeability, blood volume, and functional imaging, as assessed by ferumoxytol an iron oxide nanoparticle.

Functional Connectivity in Adult Brain Tumor Patients: A Systematic Review

Brain tumor (BT) patients often experience reduced cognitive abilities and disrupted adaptive functioning before and after treatment. An innovative approach to understanding the underlying brain networks associated with these outcomes has been to study the brain's functional connectivity (FC), the spatially distributed and temporally correlated activity throughout the brain, and how it can be affected by a tumor. The present review synthesized the extant BT FC literature that utilizes functional magnetic resonance imaging to study FC strength of commonly observed networks during rest and task. A systematic review of English articles using PubMed was conducted. Search terms included brain tumor OR glioma AND functional connectivity, independent component analysis, ICA, psychophysiological interaction, OR PPI. Studies in which participants were diagnosed with BTs as adults that evaluated specific networks of interest using independent component analysis or seed-based component analysis were included. Twenty-five studies met inclusion criteria. BT patients often presented with decreases in FC strength within well-established networks and increases in atypical FC patterns. Network differences were tumor adjacent and distal, and left hemisphere tumors generally had a greater impact on FC. FC alterations often correlated with behavioral or cognitive outcomes when assessed. Overall, BTs appear to lead to various alterations in FC across different functional networks, and the most common change is a decrease in expected FC strength. More longitudinal studies are needed to determine the time course of network alterations across treatment and recovery, the role of medical treatments in BT survivors' FC, and the potential of FC patterns as biomarkers of cognitive outcomes.

Paradoxical Activation in the Cerebellum During Language fMRI in Patients with Brain Tumors: Possible Explanations Based on Neurovascular Uncoupling and Functional Reorganization

The cerebellum is known for its crossed activation pattern with the contralateral cerebral hemisphere during language functional magnetic resonance imaging (fMRI) tasks in healthy patients. Crossed cerebro-cerebellar activation has been previously shown to occur in patients with brain tumors not affecting the activation areas. However, the presence of a tumor in left Broca's area in the inferior frontal gyrus is known to disrupt cerebral activation during language tasks. This study investigated if crossed cerebro-cerebellar activation patterns for language tasks would still occur in such patients. A total of 43 right-handed patients with a glioma affecting left Broca's area were examined for their cerebral and cerebellar activation during an fMRI language task. Only 13 of the 43 patients exhibited crossed cerebro-cerebellar activation patterns. Statistically significant differences of atypical cerebro-cerebellar activation patterns were found between cerebral right-dominant (RD) and cerebral co-dominant (CD) (p < 0.001) as well as cerebral RD and cerebral left-dominant (LD) patients (p < 0.01), while no differences were found when patients were divided based on cerebellar dominance (p > 0.75) or tumor grade (p > 0.5). No relation was found between the cerebellar and cerebral laterality index (LI) values (ρ = - 0.20; p = 0.21). Atypical activation patterns are suspected to have been caused by the tumor, perhaps a result of contralateral reorganization in some cases and false negative activation in left Broca's area from neurovascular uncoupling (NVU) in others. Cerebellar activation may also potentially indicate cerebral false negative behavior and future cerebral contralateral reorganization.

A vascular-task response dependency and its application in functional imaging of brain tumors

PURPOSE

Preoperative functional MRI (fMRI) is limited by a muted BOLD response caused by abnormal vasoreactivity and resultant neurovascular uncoupling adjacent to malignant brain tumors. We propose to overcome this limitation and more accurately identify eloquent areas adjacent to brain tumors by independently assessing vasoreactivity using breath-holding and incorporating these data into the fMRI analysis.

METHODS

Local vasoreactivity using a breath-holding paradigm with the same timing as the functional motor and language tasks was determined in 16 patients (9 glioblastomas, 1 anaplastic astrocytoma, 5 low grade astrocytomas, and 1 metastasis) and 6 healthy control subjects. We derived an fMRI model based on an observed vaso-task response dependency that takes into account the altered hemodynamics adjacent to brain tumors.

RESULTS

In both healthy controls and brain tumor subjects, we found a statistical dependency between breath-hold and task BOLD response. In tumor subjects, activation maps that take into account this vaso-task dependency demonstrated clinically meaningful areas of activation that were not seen using the task-only analysis in about half of the cases studied. This included localization of language areas adjacent to brain tumors.

CONCLUSIONS

The present preliminary results demonstrate that neurovascular uncoupling known to affect the accuracy of BOLD fMRI adjacent to brain tumors may be, at least partially, overcome by incorporating an observed vaso-task dependency in the BOLD signal analysis.

Evaluation of cerebrovascular reserve in patients with cerebrovascular diseases using resting-state MRI: A feasibility study

PURPOSE

To demonstrate the feasibility of mapping cerebrovascular reactivity (CVR) using resting-state functional MRI (fMRI) data without gas or other challenges in patients with cerebrovascular diseases and to show that brain regions affected by the diseases have diminished vascular reactivity.

MATERIALS AND METHODS

Two sub-studies were performed on patients with stroke and Moyamoya disease. In Study 1, 20 stroke patients (56.3 ± 9.7 years, 7 females) were enrolled and resting-state blood‑oxygenation-level-dependent (rs-BOLD) fMRI data were collected, from which CVR maps were computed. CVR values were compared across lesion, perilesional and control ROIs defined on anatomic images. Reproducibility of the CVR measurement was tested in 6 patients with follow-up scans. In Study 2, rs-BOLD fMRI and dynamic susceptibility contrast (DSC) MRI scans were collected in 5 patients with Moyamoya disease (32.4 ± 8.2 years, 4 females). Cerebral blood flow (CBF), cerebral blood volume (CBV), and time-to-peak (TTP) maps were obtained from the DSC MRI data. CVR values were compared between stenotic brain regions and control regions perfused by non-stenotic arteries.

RESULTS

In stroke patients, lesion CVR (0.250 ± 0.055 relative unit (r.u.)) was lower than control CVR (0.731 ± 0.088 r.u., p = 0.0002). CVR was also lower in the perilesional regions in a graded manner (perilesion 1 CVR = 0.422 ± 0.051 r.u., perilesion 2 CVR = 0.492 ± 0.046 r.u.), relative to that in the control regions (p = 0.005 and 0.036, respectively). In the repeatability analysis, a strong correlation was observed between lesion CVR (r² = 0.91, p = 0.006) measured at two time points, as well as between control CVR (r² = 0.79, p = 0.036) at two time points. In Moyamoya patients, CVR in the perfusion deficit regions delineated by DSC TTP maps (0.178 ± 0.189 r.u.) was lower than that in the control regions (0.868 ± 0.214 r.u., p = 0.013). Furthermore, the extent of reduction in CVR was significantly correlated with the extent of lengthening in TTP (r² = 0.91, p = 0.033).

CONCLUSION

Our findings suggested that rs-BOLD data can be used to reproducibly evaluate CVR in patients with cerebrovascular diseases without the use of any vasoactive challenges.

Functional Magnetic Resonance Imaging Activation Optimization in the Setting of Brain Tumor-Induced Neurovascular Uncoupling Using Resting-State Blood Oxygen Level-Dependent Amplitude of Low Frequency Fluctuations

The goal of this study was to demonstrate that a novel resting state BOLD ALFF (amplitude of low frequency fluctuations)-based correction method can substantially enhance the detectability of motor task activation in the presence of tumor-induced neurovascular uncoupling (NVU). Twelve de novo brain tumor patients who underwent comprehensive clinical BOLD fMRI exams including task fMRI and resting state fMRI (rsfMRI) were evaluated. Each patient displayed decreased/absent task fMRI activation in the ipsilesional primary motor cortex in the absence of corresponding motor deficit or suboptimal task performance, consistent with NVU. Z-score maps for the motor tasks were obtained from general linear model (GLM) analysis (reflecting motor activation vs. rest). ALFF maps were calculated from rsfMRI data. Precentral and postcentral gyri in contralesional (CL) and ipsilesional (IL) hemispheres were parcellated using an Automated Anatomical Labeling (AAL) template for each patient. A novel ALFF-based correction method was used to identify the NVU affected voxels in the ipsilesional primary motor cortex (PMC), and a correction factor was applied to normalize the baseline Z-scores for these voxels. In all cases, substantially greater activation was seen on post-ALFF correction motor activation maps within the ipsilesional precentral gyri than in the pre-ALFF correction activation maps. We have demonstrated the feasibility of a new resting state ALFF-based technique for effective correction of brain tumor-related NVU in the primary motor cortex.

Optimized preoperative motor cortex mapping in brain tumors using advanced processing of transcranial magnetic stimulation data

Background and objective

Transcranial magnetic stimulation (TMS) is a useful technique to help localize motor function prior to neurosurgical procedures. Adequate modelling of the effect of TMS on the brain is a prerequisite to obtain reliable data.

Methods

Twelve patients were included with perirolandic tumors to undergo TMS-based motor mapping. Several models were developed to analyze the mapping data, from a projection to the nearest brain surface to motor evoked potential (MEP) amplitude informed weighted average of the induced electric fields over a multilayer detailed individual head model. The probability maps were compared with direct cortical stimulation (DCS) data in all patients for the hand and in three for the foot. The gold standard was defined as the results of the DCS sampling (with on average 8 DCS-points per surgery) extrapolated over the exposed cortex (of the tailored craniotomy), and the outcome parameters were based on the similarity of the probability maps with this gold standard.

Results

All models accurately gauge the location of the motor cortex, with point-cloud based mapping algorithms having an accuracy of 83–86%, with similarly high specificity. To delineate the whole area of the motor cortex representation, the model based on the weighted average of the induced electric fields calculated with a realistic head model performs best. The optimal single threshold to visualize the field based maps is 40% of the maximal value for the anisotropic model and 50% for the isotropic model, but dynamic thresholding adds information for clinical practice.

Conclusions

The method with which TMS mapping data are analyzed clearly affects the predicted area of the primary motor cortex representation. Realistic electric field based modelling is feasible in clinical practice and improves delineation of the motor cortex representation compared to more simple point-cloud based methods.

•

Probability maps of the motor cortex representation were created from a TMS mapping.

•

The MEP-weighted averaged tissue specific induced fields based map performed best.

•

This map can gauge both motor cortex outline and hotspot, by varying the threshold.

Limitations of Resting-State Functional MR Imaging in the Setting of Focal Brain Lesions

Methods of image acquisition and analysis for resting-state functional MR imaging (rsfMR imaging) are still evolving. Neurovascular uncoupling and susceptibility artifact are important confounds of rsfMR imaging in the setting of focal brain lesions such as brain tumors. This article reviews the detection of these confounds using rsfMR imaging metrics in the setting of focal brain lesions. In the near future, with the wide range of ongoing research in rsfMR imaging, these issues likely will be overcome and will open new windows into brain function and connectivity.

Functional magnetic resonance imaging in glioma patients: from clinical applications to future perspectives

Functional magnetic resonance imaging (fMRI) allows the non-invasive assessment of human brain activity in vivo. In glioma patients, fMRI is frequently used to determine the individual functional anatomy of the motor and language network in a presurgical setting to optimize surgical procedures and prevent extensive damage to functionally eloquent areas. Novel developments based on resting-state fMRI may help to improve presurgical planning for patients which are unable to perform structured tasks and might extend presurgical mapping to include additional functional networks. Recent advances indicate a promising potential for future applications of fMRI in glioma patients which might help to identify neoplastic tissue or predict the long-term functional outcome of individual patients.