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

Value of Frequency Domain Resting-State Functional Magnetic Resonance Imaging Metrics Amplitude of Low-Frequency Fluctuation and Fractional Amplitude of Low-Frequency Fluctuation in the Assessment of Brain Tumor-Induced Neurovascular Uncoupling

The aim of this study was to explore whether the phenomenon of brain tumor-related neurovascular uncoupling (NVU) in resting-state blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) (rsfMRI) may also affect the resting-state fMRI (rsfMRI) frequency domain metrics the amplitude of low-frequency fluctuation (ALFF) and fractional ALFF (fALFF). Twelve de novo brain tumor patients, who underwent clinical fMRI examinations, including task-based fMRI (tbfMRI) and rsfMRI, were included in this Institutional Review Board-approved study. Each patient displayed decreased/absent tbfMRI activation in the primary ipsilesional (IL) sensorimotor cortex in the absence of a corresponding motor deficit or suboptimal task performance, consistent with NVU. Z-score maps for the motor tasks were obtained from general linear model analysis (reflecting motor activation vs. rest). Seed-based correlation analysis (SCA) maps of sensorimotor network, ALFF, and fALFF were calculated from rsfMRI data. Precentral and postcentral gyri in contralesional (CL) and IL hemispheres were parcellated using an automated anatomical labeling template for each patient. Region of interest (ROI) analysis was performed on four maps: tbfMRI, SCA, ALFF, and fALFF. Voxel values in the CL and IL ROIs of each map were divided by the corresponding global mean of ALFF and fALFF in the cortical brain tissue. Group analysis revealed significantly decreased IL ALFF (p = 0.02) and fALFF (p = 0.03) metrics compared with CL ROIs, consistent with similar findings of significantly decreased IL BOLD signal for tbfMRI (p = 0.0005) and SCA maps (p = 0.0004). The frequency domain metrics ALFF and fALFF may be markers of lesion-induced NVU in rsfMRI similar to previously reported alterations in tbfMRI activation and SCA-derived resting-state functional connectivity maps.

Cerebrovascular reactivity (CVR) MRI with CO2 challenge: A technical review

Cerebrovascular reactivity (CVR) is an indicator of cerebrovascular reserve and provides important information about vascular health in a range of brain conditions and diseases. Unlike steady-state vascular parameters, such as cerebral blood flow (CBF) and cerebral blood volume (CBV), CVR measures the ability of cerebral vessels to dilate or constrict in response to challenges or maneuvers. Therefore, CVR mapping requires a physiological challenge while monitoring the corresponding hemodynamic changes in the brain. The present review primarily focuses on methods that use CO2 inhalation as a physiological challenge while monitoring changes in hemodynamic MRI signals. CO2 inhalation has been increasingly used in CVR mapping in recent literature due to its potency in causing vasodilation, rapid onset and cessation of the effect, as well as advances in MRI-compatible gas delivery apparatus. In this review, we first discuss the physiological basis of CVR mapping using CO2 inhalation. We then review the methodological aspects of CVR mapping, including gas delivery apparatus, the timing paradigm of the breathing challenge, the MRI imaging sequence, and data analysis. In addition, we review alternative approaches for CVR mapping that do not require CO2 inhalation.

Key clinical benefits of neuroimaging at 7T

The growing interest in ultra-high field MRI, with more than 35.000 MR examinations already performed at 7T, is related to improved clinical results with regard to morphological as well as functional and metabolic capabilities. Since the signal-to-noise ratio increases with the field strength of the MR scanner, the most evident application at 7T is to gain higher spatial resolution in the brain compared to 3T. Of specific clinical interest for neuro applications is the cerebral cortex at 7T, for the detection of changes in cortical structure, like the visualization of cortical microinfarcts and cortical plaques in Multiple Sclerosis. In imaging of the hippocampus, even subfields of the internal hippocampal anatomy and pathology may be visualized with excellent spatial resolution. Using Susceptibility Weighted Imaging, the plaque-vessel relationship and iron accumulations in Multiple Sclerosis can be visualized, which may provide a prognostic factor of disease. Vascular imaging is a highly promising field for 7T which is dealt with in a separate dedicated article in this special issue. The static and dynamic blood oxygenation level-dependent contrast also increases with the field strength, which significantly improves the accuracy of pre-surgical evaluation of vital brain areas before tumor removal. Improvement in acquisition and hardware technology have also resulted in an increasing number of MR spectroscopic imaging studies in patients at 7T. More recent parallel imaging and short-TR acquisition approaches have overcome the limitations of scan time and spatial resolution, thereby allowing imaging matrix sizes of up to 128×128. The benefits of these acquisition approaches for investigation of brain tumors and Multiple Sclerosis have been shown recently. Together, these possibilities demonstrate the feasibility and advantages of conducting routine diagnostic imaging and clinical research at 7T.

Brain tumors disrupt the resting-state connectome

Brain tumor patients often experience functional deficits that extend beyond the tumor site. While resting-state functional MRI (rsfMRI) has been used to map such functional connectivity changes in brain tumor patients, the interplay between abnormal tumor vasculature and the rsfMRI signal is still not well understood. Therefore, there is an exigent need for new tools to elucidate how the blood‑oxygenation-level-dependent (BOLD) rsfMRI signal is modulated in brain cancer. In this initial study, we explore the utility of a preclinical model for quantifying brain tumor-induced changes on the rsfMRI signal and resting-state brain connectivity. We demonstrate that brain tumors induce brain-wide alterations of resting-state networks that extend to the contralateral hemisphere, accompanied by global attenuation of the rsfMRI signal. Preliminary histology suggests that some of these alterations in brain connectivity may be attributable to tumor-related remodeling of the neurovasculature. Moreover, this work recapitulates clinical rsfMRI findings from brain tumor patients in terms of the effects of tumor size on the neurovascular microenvironment. Collectively, these results lay the foundation of a preclinical platform for exploring the usefulness of rsfMRI as a potential new biomarker in patients with brain cancer.

PreSurgMapp: a MATLAB Toolbox for Presurgical Mapping of Eloquent Functional Areas Based on Task-Related and Resting-State Functional MRI

The main goal of brain tumor surgery is to maximize tumor resection while minimizing the risk of irreversible postoperative functional sequelae. Eloquent functional areas should be delineated preoperatively, particularly for patients with tumors near eloquent areas. Functional magnetic resonance imaging (fMRI) is a noninvasive technique that demonstrates great promise for presurgical planning. However, specialized data processing toolkits for presurgical planning remain lacking. Based on several functions in open-source software such as Statistical Parametric Mapping (SPM), Resting-State fMRI Data Analysis Toolkit (REST), Data Processing Assistant for Resting-State fMRI (DPARSF) and Multiple Independent Component Analysis (MICA), here, we introduce an open-source MATLAB toolbox named PreSurgMapp. This toolbox can reveal eloquent areas using comprehensive methods and various complementary fMRI modalities. For example, PreSurgMapp supports both model-based (general linear model, GLM, and seed correlation) and data-driven (independent component analysis, ICA) methods and processes both task-based and resting-state fMRI data. PreSurgMapp is designed for highly automatic and individualized functional mapping with a user-friendly graphical user interface (GUI) for time-saving pipeline processing. For example, sensorimotor and language-related components can be automatically identified without human input interference using an effective, accurate component identification algorithm using discriminability index. All the results generated can be further evaluated and compared by neuro-radiologists or neurosurgeons. This software has substantial value for clinical neuro-radiology and neuro-oncology, including application to patients with low- and high-grade brain tumors and those with epilepsy foci in the dominant language hemisphere who are planning to undergo a temporal lobectomy.

The Resting-State Functional Magnetic Resonance Imaging Regional Homogeneity Metrics-Kendall's Coefficient of Concordance-Regional Homogeneity and Coherence-Regional Homogeneity-Are Valid Indicators of Tumor-Related Neurovascular Uncoupling

The aim of this study is to determine whether regional homogeneity (ReHo) of resting-state blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (rsfMRI) data based on Kendall's coefficient of concordance (KCC-ReHo) and coherence (Cohe-ReHo) metrics may allow detection of brain tumor-induced neurovascular uncoupling (NVU) in the sensorimotor network similar to findings in standard motor task-based BOLD fMRI (tbfMRI) activation. Twelve de novo brain tumor patients undergoing clinical fMRI exams (tbfMRI and rsfMRI) were included in this Institutional Review Board (IRB)-approved study. Each patient displayed decreased/absent tbfMRI activation in the primary ipsilesional sensorimotor cortex in the absence of corresponding motor deficit or suboptimal task performance, consistent with NVU. Z-score maps for motor tasks were obtained from the general linear model (GLM) analysis (reflecting motor activation vs. rest). KCC-ReHo and Cohe-ReHo 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. Similar region of interest (ROI) analysis was performed on tbfMRI, KCC-ReHo, and Cohe-ReHo maps to allow direct comparison of results. Voxel values in CL and IL ROIs of each map were divided by the corresponding global mean of KCC-ReHo and Cohe-ReHo in bihemispheric cortical brain tissue. Group analysis revealed significantly decreased IL mean KCC-ReHo (p = 0.02) and Cohe-ReHo (p = 0.04) metrics compared with respective values in the CL ROIs, consistent with similar findings of significantly decreased ipsilesional BOLD signal for tbfMRI (p = 0.0005). Ipsilesional abnormalities in ReHo derived from rsfMRI may serve as potential indicators of NVU in patients with brain tumors and other resectable brain lesions; as such, ReHo findings may complement findings on tbfMRI used for presurgical planning.

Local Glioma Cells Are Associated with Vascular Dysregulation

BACKGROUND AND PURPOSE

Malignant glioma is a highly infiltrative malignancy that causes variable disruptions to the structure and function of the cerebrovasculature. While many of these structural disruptions have known correlative histopathologic alterations, the mechanisms underlying vascular dysfunction identified by resting-state blood oxygen level-dependent imaging are not yet known. The purpose of this study was to characterize the alterations that correlate with a blood oxygen level-dependent biomarker of vascular dysregulation.

MATERIALS AND METHODS

Thirty-two stereotactically localized biopsies were obtained from contrast-enhancing (n = 16) and nonenhancing (n = 16) regions during open surgical resection of malignant glioma in 17 patients. Preoperative resting-state blood oxygen level-dependent fMRI was used to evaluate the relationships between radiographic and histopathologic characteristics. Signal intensity for a blood oxygen level-dependent biomarker was compared with scores of tumor infiltration and microvascular proliferation as well as total cell and neuronal density.

RESULTS

Biopsies corresponded to a range of blood oxygen level-dependent signals, ranging from relatively normal (z = -4.79) to markedly abnormal (z = 8.84). Total cell density was directly related to blood oxygen level-dependent signal abnormality (P = .013, R² = 0.19), while the neuronal labeling index was inversely related to blood oxygen level-dependent signal abnormality (P = .016, R² = 0.21). The blood oxygen level-dependent signal abnormality was also related to tumor infiltration (P = .014) and microvascular proliferation (P = .045).

CONCLUSIONS

The relationship between local, neoplastic characteristics and a blood oxygen level-dependent biomarker of vascular function suggests that local effects of glioma cell infiltration contribute to vascular dysregulation.

[Evaluation of vascular reactivity to overcome limitations of neurovascular uncoupling in BOLD fMRI of malignant brain tumors]

The sensitivity of fMRI in identification of eloquent cortical centers in the case of large infiltrative growing tumors and pronounced peritumoral edema may be reduced or significantly limited in some cases. The main cause is an attenuated Blood-Oxygen-Level-Dependent response (BOLD) caused by pathological vascular reactivity and subsequent neurovascular uncoupling of fMRI. In our study, we attempted to overcome these limitations and increase the sensitivity of this technique in identification of eloquent cortical areas adjacent to brain tumors by using vasoreactivity features of a breath-holding test and including these data in the BOLD analysis. Local vasoreactivity using a breath-holding paradigm with the same block design of both motor and speech tests was determined in 5 healthy volunteers and 3 patients in the preoperative period (two patients with high grade gliomas and one patient with single metastasis). A coherence-based model was developed for analysis of BOLD fMRI, which took into account altered hemodynamics in peritumoral zones. Obtained coherence maps demonstrated clinically more significant activation zones that were not seen with standard methods of fMRI processing. Thus, neurovascular uncoupling that is known to affect the accuracy of the BOLD fMRI response adjacent to brain tumors may be partially overcome by including an independent measurement of vasoreactivity using a breath-holding test in the BOLD analysis.

Challenges and techniques for presurgical brain mapping with functional MRI

Functional magnetic resonance imaging (fMRI) is increasingly used for preoperative counseling and planning, and intraoperative guidance for tumor resection in the eloquent cortex. Although there have been improvements in image resolution and artifact correction, there are still limitations of this modality. In this review, we discuss clinical fMRI's applications, limitations and potential solutions. These limitations depend on the following parameters: foundations of fMRI, physiologic effects of the disease, distinctions between clinical and research fMRI, and the design of the fMRI study. We also compare fMRI to other brain mapping modalities which should be considered as alternatives or adjuncts when appropriate, and discuss intraoperative use and validation of fMRI. These concepts direct the clinical application of fMRI in neurosurgical patients.

•

fMRI is increasingly used for presurgical brain mapping for surgical planning.

•

Understanding of the limitations of fMRI is critical for its clinical use.

•

Clinical fMRI's challenges and potential solutions are discussed.

•

Intraoperative use and validation of fMRI are discussed.

Implications of neurovascular uncoupling in functional magnetic resonance imaging (fMRI) of brain tumors

Functional magnetic resonance imaging (fMRI) serves as a critical tool for presurgical mapping of eloquent cortex and changes in neurological function in patients diagnosed with brain tumors. However, the blood-oxygen-level-dependent (BOLD) contrast mechanism underlying fMRI assumes that neurovascular coupling remains intact during brain tumor progression, and that measured changes in cerebral blood flow (CBF) are correlated with neuronal function. Recent preclinical and clinical studies have demonstrated that even low-grade brain tumors can exhibit neurovascular uncoupling (NVU), which can confound interpretation of fMRI data. Therefore, to avoid neurosurgical complications, it is crucial to understand the biophysical basis of NVU and its impact on fMRI. Here we review the physiology of the neurovascular unit, how it is remodeled, and functionally altered by brain cancer cells. We first discuss the latest findings about the components of the neurovascular unit. Next, we synthesize results from preclinical and clinical studies to illustrate how brain tumor induced NVU affects fMRI data interpretation. We examine advances in functional imaging methods that permit the clinical evaluation of brain tumors with NVU. Finally, we discuss how the suppression of anomalous tumor blood vessel formation with antiangiogenic therapies can "normalize" the brain tumor vasculature, and potentially restore neurovascular coupling.

Resting-state Functional Magnetic Resonance Imaging in Presurgical Functional Mapping: Sensorimotor Localization

This article compares resting-state functional magnetic resonance (fMR) imaging with task fMR imaging for presurgical functional mapping of the sensorimotor (SM) region. Before tumor resection, 38 patients were scanned using both methods. The SM area was anatomically defined using 2 different software tools. Overlap of anatomic regions of interest with task activation maps and resting-state networks was measured in the SM region. A paired t-test showed higher overlap between resting-state maps and anatomic references compared with task activation when using a maximal overlap criterion. Resting state-derived maps are more comprehensive than those derived from task fMR imaging.

Anatomic Location of Tumor Predicts the Accuracy of Motor Function Localization in Diffuse Lower-Grade Gliomas Involving the Hand Knob Area

BACKGROUND AND PURPOSE

The accuracy of preoperative blood oxygen level-dependent fMRI remains controversial. This study assessed the association between the anatomic location of a tumor and the accuracy of fMRI-based motor function mapping in diffuse lower-grade gliomas.

MATERIALS AND METHODS

Thirty-five patients with lower-grade gliomas involving motor areas underwent preoperative blood oxygen level-dependent fMRI scans with grasping tasks and received intraoperative direct cortical stimulation. Patients were classified into an overlapping group and a nonoverlapping group, depending on the extent to which blood oxygen level-dependent fMRI and direct cortical stimulation results concurred. Tumor location was quantitatively measured, including the shortest distance from the tumor to the hand knob and the deviation distance of the midpoint of the hand knob in the lesion hemisphere relative to the midline compared with the normal contralateral hemisphere.

RESULTS

A 4-mm shortest distance from the tumor to the hand knob value was identified as optimal for differentiating the overlapping and nonoverlapping group with the receiver operating characteristic curve (sensitivity, 84.6%; specificity, 77.8%). The shortest distances from the tumor to the hand knob of ≤4 mm were associated with inaccurate fMRI-based localizations of the hand motor cortex. The shortest distances from the tumor to the hand knob were larger (P = .002), and the deviation distances for the midpoint of the hand knob in the lesion hemisphere were smaller (P = .003) in the overlapping group than in the nonoverlapping group.

CONCLUSIONS

This study suggests that the shortest distance from the tumor to the hand knob and the deviation distance for the midpoint of the hand knob on the lesion hemisphere are predictive of the accuracy of blood oxygen level-dependent fMRI results. Smaller shortest distances from the tumor to the hand knob and larger deviation distances for the midpoint of hand knob on the lesion hemisphere are associated with less accuracy of motor cortex localization with blood oxygen level-dependent fMRI. Preoperative fMRI data for surgical planning should be used cautiously when the shortest distance from the tumor to the hand knob is ≤4 mm, especially for lower-grade gliomas anterior to the central sulcus.

Functional MRI for Surgery of Gliomas

PURPOSE OF REVIEW

Advanced neuroimaging techniques such as functional MRI (fMRI) and diffusion MR tractography have been increasingly used at every stage of the surgical management of brain gliomas, as a means to improve tumor resection while preserving brain functions. This review provides an overview of the last advancements in the field of functional MRI techniques, with a particular focus on their current clinical use and reliability in the preoperative and intraoperative setting, as well as their future perspectives for personalized multimodal management of patients with gliomas.

RECENT FINDINGS

fMRI and diffusion MR tractography give relevant insights on the anatomo-functional organization of eloquent cortical areas and subcortical connections near or inside a tumor. Task-based fMRI and diffusion tensor imaging (DTI) tractography have proven to be valid and highly sensitive tools for localizing the distinct eloquent cortical and subcortical areas before surgery in glioma patients; they also show good accuracy when compared with intraoperative stimulation mapping data. Resting-state fMRI functional connectivity as well as new advanced HARDI (high angular resolution diffusion imaging) tractography methods are improving and reshaping the role of functional MRI for surgery of gliomas, with potential benefit for personalized treatment strategies. Noninvasive functional MRI techniques may offer the opportunity to perform a multimodal assessment in brain tumors, to be integrated with intraoperative mapping and clinical data for improving surgical management and oncological and functional outcome in patients affected by gliomas.

Robust presurgical functional MRI at 7 T using response consistency

Functional MRI is valuable in presurgical planning due to its non-invasive nature, repeatability, and broad availability. Using ultra-high field MRI increases the specificity and sensitivity, increasing the localization reliability and reducing scan time. Ideally, fMRI analysis for this application should identify unreliable runs and work even if the patient deviates from the prescribed task timing or if there are changes to the hemodynamic response due to pathology. In this study, a model-free analysis method-UNBIASED-based on the consistency of fMRI responses over runs was applied, to ultra-high field fMRI localizations of the hand area. Ten patients with brain tumors and epilepsy underwent 7 Tesla fMRI with multiple runs of a hand motor task in a block design. FMRI data were analyzed with the proposed approach (UNBIASED) and the conventional General Linear Model (GLM) approach. UNBIASED correctly identified and excluded fMRI runs that contained little or no activation. Generally, less motion artifact contamination was present in UNBIASED than in GLM results. Some cortical regions were identified as activated in UNBIASED but not GLM results. These were confirmed to show reproducible delayed or transient activation, which was time-locked to the task. UNBIASED is a robust approach to generating activation maps without the need for assumptions about response timing or shape. In presurgical planning, UNBIASED can complement model-based methods to aid surgeons in making prudent choices about optimal surgical access and resection margins for each patient, even if the hemodynamic response is modified by pathology. Hum Brain Mapp 38:3163-3174, 2017. © 2017 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Presurgical Brain Mapping of the Ventral Somatomotor Network in Patients with Brain Tumors Using Resting-State fMRI

BACKGROUND AND PURPOSE

Resting-state fMRI readily identifies the dorsal but less consistently the ventral somatomotor network. Our aim was to assess the relative utility of resting-state fMRI in the identification of the ventral somatomotor network via comparison with task-based fMRI in patients with brain tumor.

MATERIALS AND METHODS

We identified 26 surgically naïve patients referred for presurgical fMRI brain mapping who had undergone both satisfactory ventral motor activation tasks and resting-state fMRI. Following standard preprocessing for task-based fMRI and resting-state fMRI, general linear model analysis of the ventral motor tasks and independent component analysis of resting-state fMRI were performed with the number of components set to 20, 30, 40, and 50. Visual overlap of task-based fMRI and resting-state fMRI at different component levels was assessed and categorized as full match, partial match, or no match. Rest-versus-task-fMRI concordance was calculated with Dice coefficients across varying fMRI thresholds before and after noise removal. Multithresholded Dice coefficient volume under the surface was calculated.

RESULTS

The ventral somatomotor network was identified in 81% of patients. At the subject level, better matches between resting-state fMRI and task-based fMRI were seen with an increasing order of components (53% of cases for 20 components versus 73% for 50 components). Noise-removed group-mean volume under the surface improved as component numbers increased from 20 to 50, though ANOVA demonstrated no statistically significant difference among the 4 groups.

CONCLUSIONS

In most patients, the ventral somatomotor network can be identified with an increase in the probability of a better match at a higher component number. There is variable concordance of the ventral somatomotor network at the single-subject level between resting-state and task-based fMRI.

Vascular Reactivity Maps in Patients with Gliomas Using Breath-Holding BOLD fMRI

BACKGROUND AND PURPOSE

To evaluate whether breath-holding (BH) blood oxygenation level-dependent (BOLD) fMRI can quantify differences in vascular reactivity (VR), as there is a need for improved contrast mechanisms in gliomas.

METHODS

16 patients (gliomas, grade II = 5, III = 2, IV = 9) were evaluated using the BH paradigm: 4-second single deep breath followed by 16 seconds of BH and 40 seconds of regular breathing for five cycles. VR was defined as the difference in BOLD signal between the minimal signal seen at the end of the deep breath and maximal signal seen at the end of BH (peak-to-trough). VR was measured for every voxel and compared for gray versus white matter and tumor versus normal contralateral brain. VR maps were compared to the areas of enhancement and FLAIR/T2 abnormality.

RESULTS

VR was significantly lower in normal white matter than gray matter (P < .05) and in tumors compared to the normal, contralateral brain (P < 0.002). The area of abnormal VR (1103 ± 659 mm²) was significantly greater (P = .019) than the enhancement (543 ± 530 mm²), but significantly smaller (P = .0011) than the FLAIR abnormality (2363 ± 1232 mm²). However, the variability in the areas of gadolinium contrast enhancement versus VR abnormality indicates that the contrast mechanism elicited by BH (caused by abnormal arteriolar smooth muscles) appears to be fundamentally different from the contrast mechanism of gadolinium enhancement (caused by the presence of "leaky" gap junctions).

CONCLUSIONS

BH maps based on peak-to-trough can be used to characterize VR in brain tumors. VR maps in brain tumor patients appear to be caused by a different mechanism than gadolinium enhancement.

Cerebrovascular reactivity mapping without gas challenges

Cerebrovascular reactivity (CVR), the ability of cerebral vessels to dilate or constrict, has been shown to provide valuable information in the diagnosis and treatment evaluation of patients with various cerebrovascular conditions. CVR mapping is typically performed using hypercapnic gas inhalation as a vasoactive challenge while collecting BOLD images, but the inherent need of gas inhalation and the associated apparatus setup present a practical obstacle in applying it in routine clinical use. Therefore, we aimed to develop a new method to map CVR using resting-state BOLD data without the need of gas inhalation. This approach exploits the natural variation in respiration and measures its influence on BOLD MRI signal. In this work, we first identified a surrogate of the arterial CO₂ fluctuation during spontaneous breathing from the global BOLD signal. Second, we tested the feasibility and reproducibility of the proposed approach to use the above-mentioned surrogate as a regressor to estimate voxel-wise CVR. Third, we validated the "resting-state CVR map" with a conventional CVR map obtained with hypercapnic gas inhalation in healthy volunteers. Finally, we tested the utility of this new approach in detecting abnormal CVR in a group of patients with Moyamoya disease, and again validated the results using the conventional gas inhalation method. Our results showed that global BOLD signal fluctuation in the frequency range of 0.02-0.04Hz contains the most prominent contribution from natural variation in arterial CO₂. The CVR map calculated using this signal as a regressor is reproducible across runs (ICC=0.91±0.06), and manifests a strong spatial correlation with results measured with a conventional hypercapnia-based method in healthy subjects (r=0.88, p<0.001). We also found that resting-state CVR was able to identify vasodilatory deficit in patients with steno-occlusive disease, the spatial pattern of which matches that obtained using the conventional gas method (r=0.71±0.18). These results suggest that CVR obtained with resting-state BOLD may be a useful alternative in detecting vascular deficits in clinical applications when gas challenge is not feasible.

Brain Tumors: The Influence of Tumor Type and Routine MR Imaging Characteristics at BOLD Functional MR Imaging in the Primary Motor Gyrus

Purpose To evaluate the effects of histologic features and anatomic magnetic resonance (MR) imaging characteristics of brain tumors on the functional MR imaging signal in the primary motor cortex (PMC), as false-negative blood oxygen level-dependent (BOLD) functional MR imaging activation can limit the accurate localization of eloquent cortices. Materials and Methods Institutional review board approval was obtained, and informed consent was waived for this HIPAA-compliant retrospective study. It comprised 63 patients referred between 2006 and 2014 for preoperative functional MR imaging localization of the Rolandic cortex. The patients had glioblastoma multiforme (GBM) (n = 20), metastasis (n = 21), or meningioma (n = 22). The volumes of functional MR imaging activation were measured during performance of a bilateral hand motor task. Ratios of functional MR imaging activation were normalized to PMC volume. Statistical analysis was performed for the following: (a) differences between hemispheres within each histologic tumor type (paired Wilcoxon test), (b) differences across tumor types (Kruskal-Wallis and Fisher tests), (c) pairwise tests between tumor types (Mann-Whitney U test), (d) relationships between fast fluid-attenuated inversion recovery (FLAIR) data and enhancement volume with activation (Spearman rank correlation coefficient), and (e) differences in activation volumes by tumor location (Mann-Whitney U test). Results A significant interhemispheric difference was found between the activation volumes in GBMs (mean, 511.43 voxels ± 307.73 [standard deviation] and 330.78 voxels ± 278.95; P < .01) but not in metastases (504.68 voxels ± 220.98 and 460.22 voxels ± 276.83; P = .15) or meningiomas (424.07 voxels ± 247.58 and 415.18 voxels ± 222.36; P = .85). GBMs showed significantly lower activation ratios (median, 0.49; range, 0.04-1.15) than metastases (median, 0.79; range, 0.28-1.66; P = .043) and meningiomas (median, 0.91; range, 0.52-2.05; P < .01). There was a moderate correlation with the volumes of FLAIR abnormality in metastases (ρ = -0.50) and meningiomas (ρ = -0.55). Enhancement volume (ρ = -0.11) and tumor distance from the PMC (median, 0.73 and range, 0.04-2.05 for near and median, 0.82 and range, 0.39-1.66 for far; P = .14) did not influence activation. Conclusion BOLD functional MR imaging activation in the ipsilateral PMC is influenced by tumor type and is significantly reduced in GBMs. FLAIR abnormality correlates moderately with the activation ratios in metastases and meningiomas. ^© RSNA, 2016 Online supplemental material is available for this article.

Presurgical motor, somatosensory and language fMRI: Technical feasibility and limitations in 491 patients over 13 years

OBJECTIVES

To analyse the long-term feasibility and limitations of presurgical fMRI in a cohort of tumour and epilepsy patients with different MR-scanners at 1.5 and 3.0 T.

METHODS

Four hundred and ninety-one consecutive patients undergoing presurgical fMRI between 2000 and 2012 on five different MR-scanners using established paradigms and semi-automated data processing were included. Success rates of task performance and BOLD-activation were determined for motor and somatosensory somatotopic mapping and language localisation. Procedural success, failures and imaging artifacts were analysed. MR-field strengths were compared.

RESULTS

Two thousand three hundred fifteen of 2348 (98.6 %) attempted paradigms (1033 motor, 1220 speech, 95 somatosensory) were successfully performed. 100 paradigms (4.3 %) were repetition runs. 23 speech, 6 motor and 2 sensory paradigms failed for non-compliance and technical issues. Most language paradigm failures were noted in overt sentence generation. Average significant BOLD-activation was higher for motor than language paradigms (95.8 vs. 81.6 %). Most language paradigms showed significantly higher activation rates at 3 T compared to 1.5 T, whereas no significant difference was found for motor paradigms.

CONCLUSIONS

fMRI proved very robust for the presurgical localisation of the different motor and somatosensory body representations, as well as Broca's and Wernicke's language areas across different MR-scanners at 1.5 and 3.0 T over 13 years.

KEY POINTS

• Standardised presurgical motor and language fMRI is robust across various MRI platforms. • Motor fMRI is less dependent on field strength than language fMRI. • fMRI task failures are relatively low and are reduced by paradigm repetition.

Improving the clinical potential of ultra-high field fMRI using a model-free analysis method based on response consistency

Objective

To develop an analysis method that is sensitive to non-model-conform responses often encountered in ultra-high field presurgical planning fMRI. Using the consistency of time courses over a number of experiment repetitions, it should exclude low quality runs and generate activation maps that reflect the reliability of responses.

Materials and methods

7 T fMRI data were acquired from six healthy volunteers: three performing purely motor tasks and three a visuomotor task. These were analysed with the proposed approach (UNBIASED) and the GLM.

Results

UNBIASED results were generally less affected by false positive results than the GLM. Runs that were identified as being of low quality were confirmed to contain little or no activation. In two cases, regions were identified as activated in UNBIASED but not GLM results. Signal changes in these areas were time-locked to the task, but were delayed or transient.

Conclusion

UNBIASED is shown to be a reliable means of identifying consistent task-related signal changes regardless of response timing. In presurgical planning, UNBIASED could be used to rapidly generate reliable maps of the consistency with which eloquent brain regions are activated without recourse to task timing and despite modified hemodynamics.

Electronic supplementary material

The online version of this article (doi:10.1007/s10334-016-0533-8) contains supplementary material, which is available to authorized users.

Assessing Region of Interest Schemes for the Corticospinal Tract in Patients With Brain Tumors

Diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) techniques are widely used for identifying the corticospinal tract (CST) white matter pathways as part of presurgical planning. However, mass effects in patients with brain tumors tend to cause anatomical distortions and compensatory functional reorganization of the cortex, which may lead to inaccurate mapping of white matter tracts. To overcome these problems, we compared different region-of-interest (ROI) selection schemes to track CST fibers in patients with brain tumors. Our study investigated the CSTs of 16 patients with intracranial tumors. The patients were classified into 3 subgroups according to the spatial relationships of the lesion and the primary motor cortex (PMC)/internal capsule. Specifically, we investigated the key factors that cause distorted tractography in patients with tumors. We compared 3 CST tractography methods that used different ROI selection schemes. The results indicate that CST fiber tracking methods based only on anatomical ROIs could possibly lead to distortions near the PMC region and may be unable to effectively localize the PMC. In contrast, the dual ROI method, which uses ROIs that have been selected from both blood oxygen level-dependent functional MRI (BOLD-fMRI) activation and anatomical landmarks, enabled the tracking of fibers to the motor cortex. The results demonstrate that the dual ROI method can localize the entire CST fiber pathway and can accurately describe the spatial relationships of CST fibers relative to the tumor. These results illustrate the reliability of using fMRI-guided DTT in patients with tumors. The combination of fMRI and anatomical information enhances the identification of tracts of interest in brains with anatomical deformations, which provides neurosurgeons with a more accurate approach for visualizing and localizing white matter fiber tracts in patients with brain tumors. This approach enhances surgical performance and perserves brain function.

Demonstration of Brain Tumor-Induced Neurovascular Uncoupling in Resting-State fMRI at Ultrahigh Field

To demonstrate in a small case series for the first time the phenomenon of brain tumor-related neurovascular uncoupling (NVU) in resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) at ultrahigh field (7T). Two de novo (i.e., untreated) brain tumor patients underwent both BOLD resting-state fMRI (rsfMRI) on a 7T MRI system and motor task-based BOLD fMRI at 3T. Ipsilesional (i.e., ipsilateral to tumor or IL) and contralesional (i.e., contralateral to tumor or CL) region of interest (ROI) analysis was performed on both 3T motor task-related general linear model-derived activation maps and on 7T rsfMRI independent component analysis (ICA)-derived sensorimotor network maps for each case. Asymmetry scores (ASs) were computed based on numbers of suprathreshold voxels in the IL and CL ROIs. In each patient, ASs derived from ROI analysis of suprathreshold voxels in IL and CL ROIs in task-related activation maps and rsfMRI ICA-derived sensorimotor component maps indicate greater number of suprathreshold voxels in contralesional than ipsilesional sensorimotor cortex in both maps. In patient 1, an AS of 0.2 was obtained from the suprathreshold Z-score spectrum (voxels with Z-scores >5.0) of the task-based activation map and AS of 1.0 was obtained from the suprathreshold Z-score spectrum (Z-scores >5.0) of the ICA-derived sensorimotor component map. Similarly, in patient 2, an AS of 1.0 was obtained from the suprathreshold Z-score spectrum (Z-scores >5.0) of the task-based activation map and an AS of 1.0 was obtained from the suprathreshold Z-score spectrum (Z-scores >5.0) of the ICA-derived sensorimotor component map. Overall, decreased BOLD signal was noted in IL compared with CL ROIs on both task-based activation maps and ultrahigh field resting-state maps, indicating the presence of NVU. We have demonstrated evidence of NVU on ultrahigh field 7T rsfMRI comparable with the findings on standard 3T motor task-based fMRI in both cases.

Blood oxygen-level dependent functional assessment of cerebrovascular reactivity: Feasibility for intraoperative 3 Tesla MRI

PURPOSE

To assess the feasibility of functional blood oxygen-level dependent (BOLD) MRI to evaluate intraoperative cerebrovascular reactivity (CVR) at 3 Tesla field strength.

METHODS

Ten consecutive neurosurgical subjects scheduled for a clinical intraoperative MRI examination were enrolled in this study. In addition to the clinical protocol a BOLD sequence was implemented with three cycles of 44 s apnea to calculate CVR values on a voxel-by-voxel basis throughout the brain. The CVR range was then color-coded and superimposed on an anatomical volume to create high spatial resolution CVR maps.

RESULTS

Ten subjects (mean age 34.8 ± 13.4; 2 females) uneventfully underwent the intraoperative BOLD protocol, with no complications occurring. Whole-brain CVR for all subjects was (mean ± SD) 0.69 ± 0.42, whereas CVR was markedly higher for tumor subjects as compared to vascular subjects, 0.81 ± 0.44 versus 0.33 ± 0.10, respectively. Furthermore, color-coded functional maps could be robustly interpreted for a whole-brain assessment of CVR.

CONCLUSION

We demonstrate that intraoperative BOLD MRI is feasible in creating functional maps to assess cerebrovascular reactivity throughout the brain in subjects undergoing a neurosurgical procedure. Magn Reson Med 77:806-813, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Neurovascular uncoupling in resting state fMRI demonstrated in patients with primary brain gliomas

BACKGROUND

To demonstrate that the problem of brain tumor-related neurovascular uncoupling (NVU) is a significant issue with respect to resting state blood oxygen level dependent (BOLD) functional MRI (rsfMRI) similar to task-based BOLD fMRI, in which signal detectability can be compromised by breakdown of normal neurovascular coupling.

METHODS

We evaluated seven de novo brain tumor patients who underwent resting state fMRI as part of comprehensive clinical fMRI exams at 3 Tesla. For each of the seven patients who demonstrated evidence of NVU on task-based motor fMRI, we performed both an independent component analysis (ICA) and an atlas-based parcellation-based seed correlation analysis (SCA) of the resting state fMRI data. For each patient, ipsilesional (IL) and contralesional (CL) regions of interest (ROIs) comprising primary motor and somatosensory cortices were used to evaluate BOLD signal changes on Z score maps derived from both ICA and SCA analysis for evidence of NVU. A subsequent two-tailed t-test was performed to determine whether statistically significant differences between the two sides were present that were consistent with NVU.

RESULTS

In seven patients, overall decreased BOLD signal (based on suprathreshold voxels in ICA and SCA-derived Z-score maps) was noted in IL compared with CL ROIs (P < 0.01), consistent with NVU.

CONCLUSION

We have demonstrated that NVU can result in false negative BOLD signal changes on rsfMRI comparable to previously published findings on standard motor task-based fMRI.

BOLD fMRI of cerebrovascular reactivity in the middle cerebral artery territory: A 100 volunteers' study

BACKGROUND AND PURPOSE

The assessment of cerebrovascular reactivity (CVR) has shown promising results for its use in medical diagnosis and prognosis, especially in patients suffering from severe intracranial arterial stenosis. However, its quantification remains uncertain because of a large variability inherent in brain anatomy and in methodological settings. To overcome this variability, we provide lateralization index (LI) values for CVR within the middle cerebral artery territory to detect CVR impairment.

MATERIALS AND METHODS

We assessed CVR in 100 volunteers (41 females; 47.52 ± 21.58 years) without cervico-encephalic arterial stenosis using BOLD-fMRI contrast with a block-design hypercapnic challenge. Averaged end-tidal CO2 was used as a physiological regressor for statistical analyses with a general linear model. We measured %BOLD signal-change in segmented gray matter regions of interest in the middle cerebral artery territory (MCA). We calculated a laterality index according to the following formula: LI=(CVRleft-CVRright)/(CVRleft+CVRright). We tested the effects of methodological settings (i.e. hypercapnic gas, gas administration means, MR acquisition and sex) on %BOLD signal change and LI values with analysis of variance.

RESULTS

No adverse effects of the hypercapnic challenge were reported. LI values were independent of experimental conditions. Mean LI calculated in MCA territories was 0.016 ± 0.031, giving the lower and upper limits of 95% (m ± 2SD) of this population distribution at]-0.05; 0.08[.

CONCLUSION

LI can effectively help us to overcome measurement variabilities. Therefore, it can be used to detect abnormal asymmetries in CVR and identify patients at higher risk of ischemic stroke.

A case study of magnetic resonance imaging of cerebrovascular reactivity: a powerful imaging marker for mild traumatic brain injury

PRIMARY OBJECTIVE

To use breath-hold functional magnetic resonance imaging (fMRI) to localize the brain regions with impaired cerebrovascular reactivity (CVR) in a female patient diagnosed with mild traumatic brain injury (mTBI). The extent of impaired CVR was evaluated 2 months after concussion. Follow-up scan was performed 1 year post-mTBI using the same breath-hold fMRI technique.

RESEARCH DESIGN

Case report.

METHODS AND PROCEDURES

fMRI blood oxygenation dependent level (BOLD) signals were measured under breath-hold challenge in a female mTBI patient 2 months after concussion followed by a second fMRI with breath-hold challenge 1 year later. CVR was expressed as the percentage change of BOLD signals per unit time of breath-hold.

MAIN OUTCOMES

In comparison with CVR measurement of normal control subjects, statistical maps of CVR revealed substantial neurovascular deficits and hemispheric asymmetry within grey and white matter in the initial breath-hold fMRI scan. Follow-up breath-hold fMRI performed 1 year post-mTBI demonstrated normalization of CVR accompanied with symptomatic recovery.

CONCLUSIONS

CVR may serve as an imaging biomarker to detect subtle deficits in both grey and white matter for individual diagnosis of mTBI. The findings encourage further investigation of hypercapnic fMRI as a diagnostic tool for mTBI.

Special considerations/technical limitations of blood-oxygen-level-dependent functional magnetic resonance imaging

In this review, limitations affecting the results of presurgical mapping with blood-oxygen-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) are discussed. There is a great need to standardize fMRI acquisition and analysis methods and establish guidelines to address quality control issues. Several national and international organizations are formulating guidelines and standards for both clinical and research applications of BOLD fMRI. Consensus regarding management of these issues will likely both improve the clinical standard of care and enhance future research applications of fMRI.

Preoperative mapping of the sensorimotor cortex: comparative assessment of task-based and resting-state FMRI

Resting state fMRI (rs-fMRI) has recently been considered as a possible complement or alternative to task-based fMRI (tb-fMRI) for presurgical mapping. However, evidence of its usefulness remains scant, because existing studies have investigated relatively small samples and focused primarily on qualitative evaluation. The aim of this study is to investigate the clinical usefulness of rs-fMRI in the context of presurgical mapping of motor functions, and in particular to determine the degree of correspondence with tb-fMRI which, while not a gold-standard, is commonly used in preoperative setting. A group of 13 patients with lesions close to the sensorimotor cortex underwent rs-fMRI and tb-fMRI to localize the hand, foot and mouth motor areas. We assessed quantitatively the degree of correspondence between multiple rs-fMRI analyses (independent-component and seed-based analyses) and tb-fMRI, with reference to sensitivity and specificity of rs-fMRI with respect to tb-fMRI, and centre-of-mass distances. Agreement with electro-cortical stimulation (ECS) was also investigated, and a traditional map thresholding approach based on agreement between two experienced operators was compared to an automatic threshold determination method. Rs-fMRI can localize the sensorimotor cortex successfully, providing anatomical specificity for hand, foot and mouth motor subregions, in particular with seed-based analyses. Agreement with tb-fMRI was only partial and rs-fMRI tended to provide larger patterns of correlated activity. With respect to the ECS data available, rs-fMRI and tb-fMRI performed comparably, even though the shortest distance to stimulation points was observed for the latter. Notably, the results of both were on the whole robust to thresholding procedure. Localization performed by rs-fMRI is not equivalent to tb-fMRI, hence rs-fMRI cannot be considered as an outright replacement for tb-fMRI. Nevertheless, since there is significant agreement between the two techniques, rs-fMRI can be considered with caution as a potential alternative to tb-fMRI when patients are unable to perform the task.

Motor network plasticity and low-frequency oscillations abnormalities in patients with brain gliomas: a functional MRI study

Brain plasticity is often associated with the process of slow-growing tumor formation, which remodels neural organization and optimizes brain network function. In this study, we aimed to investigate whether motor function plasticity would display deficits in patients with slow-growing brain tumors located in or near motor areas, but who were without motor neurological deficits. We used resting-state functional magnetic resonance imaging to probe motor networks in 15 patients with histopathologically confirmed brain gliomas and 15 age-matched healthy controls. All subjects performed a motor task to help identify individual motor activity in the bilateral primary motor cortex (PMC) and supplementary motor area (SMA). Frequency-based analysis at three different frequencies was then used to investigate possible alterations in the power spectral density (PSD) of low-frequency oscillations. For each group, the average PSD was determined for each brain region and a nonparametric test was performed to determine the difference in power between the two groups. Significantly reduced inter-hemispheric functional connectivity between the left and right PMC was observed in patients compared with controls (P<0.05). We also found significantly decreased PSD in patients compared to that in controls, in all three frequency bands (low: 0.01-0.02 Hz; middle: 0.02-0.06 Hz; and high: 0.06-0.1 Hz), at three key motor regions. These findings suggest that in asymptomatic patients with brain tumors located in eloquent regions, inter-hemispheric connection may be more vulnerable. A comparison of the two approaches indicated that power spectral analysis is more sensitive than functional connectivity analysis for identifying the neurological abnormalities underlying motor function plasticity induced by slow-growing tumors.

Cerebrovascular reactivity mapping: an evolving standard for clinical functional imaging

SUMMARY

This review article explains the methodology of breath-hold cerebrovascular reactivity mapping, both in terms of acquisition and analysis, and reviews applications of this method to presurgical mapping, particularly with respect to blood oxygen level-dependent fMRI. Its main application in clinical fMRI is for the assessment of neurovascular uncoupling potential. Neurovascular uncoupling is potentially a major limitation of clinical fMRI, particularly in the setting of mass lesions in the brain such as brain tumors and intracranial vascular malformations that are associated with alterations in regional hemodynamics on either an acquired or congenital basis. As such, breath-hold cerebrovascular reactivity mapping constitutes an essential component of quality control analysis in clinical fMRI, particularly when performed for presurgical mapping of eloquent cortex. Exogenous carbon dioxide challenges used for cerebrovascular reactivity mapping will also be discussed, and their applications to the evaluation of cerebrovascular reserve and cerebrovascular disease will be described.

GRE T2∗-weighted MRI: principles and clinical applications

The sequence of a multiecho gradient recalled echo (GRE) T2*-weighted imaging (T2*WI) is a relatively new magnetic resonance imaging (MRI) technique. In contrast to T2 relaxation, which acquires a spin echo signal, T2* relaxation acquires a gradient echo signal. The sequence of a GRE T2*WI requires high uniformity of the magnetic field. GRE T2*WI can detect the smallest changes in uniformity in the magnetic field and can improve the rate of small lesion detection. In addition, the T2* value can indirectly reflect changes in tissue biochemical components. Moreover, it can be used for the early diagnosis and quantitative diagnosis of some diseases. This paper reviews the principles and clinical applications as well as the advantages and disadvantages of GRE T2*WI.