Hemodynamic response imaging: a potential tool for the assessment of angiogenesis in brain tumors

Patterns of gray and white matter functional networks involvement in glioblastoma patients: indirect mapping from clinical MRI scans

Background

Resting-state functional-MRI studies identified several cortical gray matter functional networks (GMNs) and white matter functional networks (WMNs) with precise anatomical localization. Here, we aimed at describing the relationships between brain's functional topological organization and glioblastoma (GBM) location. Furthermore, we assessed whether GBM distribution across these networks was associated with overall survival (OS).

Materials and methods

We included patients with histopathological diagnosis of IDH-wildtype GBM, presurgical MRI and survival data. For each patient, we recorded clinical-prognostic variables. GBM core and edema were segmented and normalized to a standard space. Pre-existing functional connectivity-based atlases were used to define network parcellations: 17 GMNs and 12 WMNs were considered in particular. We computed the percentage of lesion overlap with GMNs and WMNs, both for core and edema. Differences between overlap percentages were assessed through descriptive statistics, ANOVA, post-hoc tests, Pearson's correlation tests and canonical correlations. Multiple linear and non-linear regression tests were employed to explore relationships with OS.

Results

99 patients were included (70 males, mean age 62  years). The most involved GMNs included ventral somatomotor, salient ventral attention and default-mode networks; the most involved WMNs were ventral frontoparietal tracts, deep frontal white matter, and superior longitudinal fasciculus system. Superior longitudinal fasciculus system and dorsal frontoparietal tracts were significantly more included in the edema (p < 0.001). 5 main patterns of GBM core distribution across functional networks were found, while edema localization was less classifiable. ANOVA showed significant differences between mean overlap percentages, separately for GMNs and WMNs (p-values<0.0001). Core-N12 overlap predicts higher OS, although its inclusion does not increase the explained OS variance.

Discussion and conclusion

Both GBM core and edema preferentially overlap with specific GMNs and WMNs, especially associative networks, and GBM core follows five main distribution patterns. Some inter-related GMNs and WMNs were co-lesioned by GBM, suggesting that GBM distribution is not independent of the brain's structural and functional organization. Although the involvement of ventral frontoparietal tracts (N12) seems to have some role in predicting survival, network-topology information is overall scarcely informative about OS. fMRI-based approaches may more effectively demonstrate the effects of GBM on brain networks and survival.

Hemodynamic Imaging in Cerebral Diffuse Glioma—Part B: Molecular Correlates, Treatment Effect Monitoring, Prognosis, and Future Directions

Gliomas, and glioblastoma in particular, exhibit an extensive intra- and inter-tumoral molecular heterogeneity which represents complex biological features correlating to the efficacy of treatment response and survival. From a neuroimaging point of view, these specific molecular and histopathological features may be used to yield imaging biomarkers as surrogates for distinct tumor genotypes and phenotypes. The development of comprehensive glioma imaging markers has potential for improved glioma characterization that would assist in the clinical work-up of preoperative treatment planning and treatment effect monitoring. In particular, the differentiation of tumor recurrence or true progression from pseudoprogression, pseudoresponse, and radiation-induced necrosis can still not reliably be made through standard neuroimaging only. Given the abundant vascular and hemodynamic alterations present in diffuse glioma, advanced hemodynamic imaging approaches constitute an attractive area of clinical imaging development. In this context, the inclusion of objective measurable glioma imaging features may have the potential to enhance the individualized care of diffuse glioma patients, better informing of standard-of-care treatment efficacy and of novel therapies, such as the immunotherapies that are currently increasingly investigated. In Part B of this two-review series, we assess the available evidence pertaining to hemodynamic imaging for molecular feature prediction, in particular focusing on isocitrate dehydrogenase (IDH) mutation status, MGMT promoter methylation, 1p19q codeletion, and EGFR alterations. The results for the differentiation of tumor progression/recurrence from treatment effects have also been the focus of active research and are presented together with the prognostic correlations identified by advanced hemodynamic imaging studies. Finally, the state-of-the-art concepts and advancements of hemodynamic imaging modalities are reviewed together with the advantages derived from the implementation of radiomics and machine learning analyses pipelines.

Feasibility of glioblastoma tissue response mapping with physiologic BOLD imaging using precise oxygen and carbon dioxide challenge

OBJECTIVES

Innovative physiologic MRI development focuses on depiction of heterogenous vascular and metabolic features in glioblastoma. For this feasibility study, we employed blood oxygenation level-dependent (BOLD) MRI with standardized and precise carbon dioxide (CO₂) and oxygen (O₂) modulation to investigate specific tumor tissue response patterns in patients with newly diagnosed glioblastoma.

MATERIALS AND METHODS

Seven newly diagnosed untreated patients with suspected glioblastoma were prospectively included to undergo a BOLD study with combined CO₂ and O₂ standardized protocol. %BOLD signal change/mmHg during hypercapnic, hypoxic, and hyperoxic stimulus was calculated in the whole brain, tumor lesion and segmented volumes of interest (VOI) [contrast-enhancing (CE) - tumor, necrosis and edema] to analyze their tissue response patterns.

RESULTS

Quantification of BOLD signal change after gas challenges can be used to identify specific responses to standardized stimuli in glioblastoma patients. Integration of this approach with automatic VOI segmentation grants improved characterization of tumor subzones and edema. Magnitude of BOLD signal change during the 3 stimuli can be visualized at voxel precision through color-coded maps overlayed onto whole brain and identified VOIs.

CONCLUSIONS

Our preliminary investigation shows good feasibility of BOLD with standardized and precise CO₂ and O₂ modulation as an emerging physiologic imaging technique to detail specific glioblastoma characteristics. The unique tissue response patterns generated can be further investigated to better detail glioblastoma lesions and gauge treatment response.

Chronic hypoxia in pregnant mice impairs the placental and fetal vascular response to acute hypercapnia in BOLD-MRI hemodynamic response imaging

INTRODUCTION

Brief hypercapnic challenge causes acute placental hypoperfusion with fetal brain sparing on BOLD-MRI. We hypothesize that this non-invasive imaging strategy can distinguish between normal pregnancy and chronic placental hypoperfusion (using the maternal hypoxia model).

METHODS

Eighteen pregnant female ICR mice were randomized to three groups: normoxia, late-onset hypoxia (12%O₂;E13.5-17.5) and early-onset hypoxia (12%O₂;E10.5-17.5). On E17.5, animals were imaged in a 4.7-T Bruker-Biospec MRI scanner. Fast coronal True-FISP was performed to identify organs of interest (placenta and fetal heart, liver and brain). BOLD-MRI was performed at baseline and during a 4-min hypercapnic challenge (5%CO₂). %-change in placental and fetal signal was analyzed from T2*-weighted gradient echo MR images. Following MRI, fetuses and placentas were harvested, weighed and immuno-stained.

RESULTS

In normoxic mice, hypercapnia caused reduction in BOLD-MRI signal in placenta (-44% ± 7%; p < 0.0001), fetal liver (-32% ± 7%; p < 0.0001) and fetal heart (-54% ± 12%; p < 0.002), with relative fetal brain sparing (-12% ± 5%; p < 0.0001). These changes were markedly attenuated in both hypoxia groups. Baseline fetal brain/placenta SI ratio was highest in normoxic mice (1.14 ± 0.017) and reduced with increasing duration of hypoxia (late-onset hypoxia: 1.00 ± 0.026; early-onset hypoxia: 0.91 ± 0.016; p = 0.02). Both hypoxic groups exhibited fetal growth restriction with prominent placental glycogen-containing cells, particularly in early-onset hypoxia. There was increased fetal neuro- and intestinal-apoptosis in early-onset hypoxia only.

CONCLUSIONS

BOLD-MRI with brief hypercapnic challenge distinguished between normoxia and both hypoxia groups, while fetal neuroapoptosis was only observed after early-onset hypoxia. This suggests that BOLD-MRI with hypercapnic challenge can identify chronic fetal asphyxia before the onset of irreversible brain injury.

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.

Dependency of the blood oxygen level dependent-response to hyperoxic challenges on the order of gas administration in intracranial malignancies

PURPOSE

Literature reports contradicting results on the response of brain tumors to vascular stimuli measured in T2*-weighted MRI. Here, we analyzed the potential dependency of the MRI-response to (hypercapnic) hyperoxia on the order of the gas administration.

METHODS

T2* values were quantified at 3 Tesla in eight consenting patients at rest and during inhalation of hyperoxic/hypercapnic gas mixtures. Patients were randomly divided into two groups undergoing different gas administration protocols (group A: medical air-pure oxygen-carbogen; group B: medical air-carbogen-pure oxygen). Mann-Whitney U test and Wilcoxon signed rank test have been used to proof differences in T2* regarding respiratory challenge or different groups, respectively.

RESULTS

T2* values at rest for gray and white matter were 50.3 ± 2.6 ms and 46.1 ± 2.0 ms, respectively, and slightly increased during challenge. In tumor areas, T2* at rest were: necrosis = 74.1 ± 10.1 ms; edema = 60.3 ± 17.6 ms; contrast-enhancing lesions = 48.6 ± 20.7 ms; and solid T2-hyperintense lesions = 45.0 ± 3.0 ms. Contrast-enhancing lesions strongly responded to oxygen (+ 20.7%) regardless on the gas protocol (p = 0.482). However, the response to carbogen significantly depended on the order of gas administration (group A, + 18.6%; group B, - 6.4%, p = 0.042). In edemas, a different trend between group was found when breathing oxygen (group A, - 9.9%; group B, + 19.5%, p = 0.057).

CONCLUSION

Preliminary results show a dependency of the T2* response of contrast-enhancing brain tumor lesions on the order of the gas administration. The gas administration protocol is an important factor in the interpretation of the T2*-response in areas of abnormal vascular growth.

Increased putamen hypercapnic vasoreactivity in levodopa-induced dyskinesia

In a rodent model of Parkinson's disease (PD), levodopa-induced involuntary movements have been linked to striatal angiogenesis - a process that is difficult to document in living human subjects. Angiogenesis can be accompanied by localized increases in cerebral blood flow (CBF) responses to hypercapnia. We therefore explored the possibility that, in the absence of levodopa, local hypercapnic CBF responses are abnormally increased in PD patients with levodopa-induced dyskinesias (LID) but not in their nondyskinetic (NLID) counterparts. We used H215O PET to scan 24 unmedicated PD subjects (12 LID and 12 NLID) and 12 matched healthy subjects in the rest state under normocapnic and hypercapnic conditions. Hypercapnic CBF responses were compared to corresponding levodopa responses from the same subjects. Group differences in hypercapnic vasoreactivity were significant only in the posterior putamen, with greater CBF responses in LID subjects compared with the other subjects. Hypercapnic and levodopa-mediated CBF responses measured in this region exhibited distinct associations with disease severity: the former correlated with off-state motor disability ratings but not symptom duration, whereas the latter correlated with symptom duration but not motor disability. These are the first in vivo human findings linking LID to microvascular changes in the basal ganglia.

Diffuse gliomas exhibit whole brain impaired cerebrovascular reactivity

PURPOSE

Cerebral diffuse gliomas exhibit perilesional impaired cerebrovascular reactivity (CVR), yet the degree of impairment as well as its full spatial extent in the brain remains unknown. With quantitative fMRI, we studied twelve subjects with untreated brain diffuse glioma and twelve healthy controls to assess CVR impairment and determine its distribution throughout the brain.

METHODS

In a prospective case-control study, quantitative CVR measurements were derived from BOLD fMRI volumes during standardized iso-oxic changes in carbon dioxide. Whole brain CVR was assessed with additional detailed analyses using specific tumor and tissue masks and compared to datasets of healthy controls.

RESULTS

Whole brain CVR was significantly impaired compared to healthy controls (0.11±0.10 versus 0.28±0.8, p<0.01). All diffuse glioma patients exhibited even more severely impaired intralesional CVR (mean 0.01±0.06). Increasing tumor volume significantly correlated with severity of intralesional CVR impairment (p<0.05, R²=0.38), and whole brain CVR impairment (p<0.05, R²=0.55).

CONCLUSION

Patients with brain diffuse glioma exhibit intralesional and whole brain impaired CVR with severity correlating to tumor volume. Quantitative fMRI may be entertained to study antitumor therapy efficacy by tracking CVR changes and may have a complementary role to better interpret BOLD associated neurovascular uncoupling.

Carbogen-induced increases in tumor oxygenation depend on the vascular status of the tumor: A multiparametric MRI study in two rat glioblastoma models

The alleviation of hypoxia in glioblastoma with carbogen to improve treatment has met with limited success. Our hypothesis is that the eventual benefits of carbogen depend on the capacity for vasodilation. We examined, with MRI, changes in fractional cerebral blood volume, blood oxygen saturation, and blood oxygenation level dependent signals in response to carbogen. The analyses were performed in two xenograft models of glioma (U87 and U251) recognized to have different vascular patterns. Carbogen increased fractional cerebral blood volume, blood oxygen saturation, and blood oxygenation level dependent signals in contralateral tissues. In the tumor core and peritumoral regions, changes were dependent on the capacity to vasodilate rather than on resting fractional cerebral blood volume. In the highly vascularised U87 tumor, carbogen induced a greater increase in fractional cerebral blood volume and blood oxygen saturation in comparison to the less vascularized U251 tumor. The blood oxygenation level dependent signal revealed a delayed response in U251 tumors relative to the contralateral tissue. Additionally, we highlight the considerable heterogeneity of fractional cerebral blood volume, blood oxygen saturation, and blood oxygenation level dependent within U251 tumor in which multiple compartments co-exist (tumor core, rim and peritumoral regions). Finally, our study underlines the complexity of the flow/metabolism interactions in different models of glioblastoma. These irregularities should be taken into account in order to palliate intratumoral hypoxia in clinical trials.

Impact of physiological noise correction on detecting blood oxygenation level-dependent contrast in the breast

Physiological fluctuations are expected to be a dominant source of noise in blood oxygenation level-dependent (BOLD) magnetic resonance imaging (MRI) experiments to assess tumour oxygenation and angiogenesis. This work investigates the impact of various physiological noise regressors: retrospective image correction (RETROICOR), heart rate (HR) and respiratory volume per unit time (RVT), on signal variance and the detection of BOLD contrast in the breast in response to a modulated respiratory stimulus. BOLD MRI was performed at 3 T in ten volunteers at rest and during cycles of oxygen and carbogen gas breathing. RETROICOR was optimized using F-tests to determine which cardiac and respiratory phase terms accounted for a significant amount of signal variance. A nested regression analysis was performed to assess the effect of RETROICOR, HR and RVT on the model fit residuals, temporal signal-to-noise ratio, and BOLD activation parameters. The optimized RETROICOR model accounted for the largest amount of signal variance ([Formula: see text]  =  3.3  ±  2.1%) and improved the detection of BOLD activation (P  =  0.002). Inclusion of HR and RVT regressors explained additional signal variance, but had a negative impact on activation parameter estimation (P  <  0.001). Fluctuations in HR and RVT appeared to be correlated with the stimulus and may contribute to apparent BOLD signal reactivity.

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.

Using patient-specific hemodynamic response function in epileptic spike analysis of human epilepsy: a study based on EEG-fNIRS

Functional near-infrared spectroscopy (fNIRS) can be combined with electroencephalography (EEG) to continuously monitor the hemodynamic signal evoked by epileptic events such as seizures or interictal epileptiform discharges (IEDs, aka spikes). As estimation methods assuming a canonical shape of the hemodynamic response function (HRF) might not be optimal, we sought to model patient-specific HRF (sHRF) with a simple deconvolution approach for IED-related analysis with EEG-fNIRS data. Furthermore, a quadratic term was added to the model to account for the nonlinearity in the response when IEDs are frequent. Prior to analyzing clinical data, simulations were carried out to show that the HRF was estimable by the proposed deconvolution methods under proper conditions. EEG-fNIRS data of five patients with refractory focal epilepsy were selected due to the presence of frequent clear IEDs and their unambiguous focus localization. For each patient, both the linear sHRF and the nonlinear sHRF were estimated at each channel. Variability of the estimated sHRFs was seen across brain regions and different patients. Compared with the SPM8 canonical HRF (cHRF), including these sHRFs in the general linear model (GLM) analysis led to hemoglobin activations with higher statistical scores as well as larger spatial extents on all five patients. In particular, for patients with frequent IEDs, nonlinear sHRFs were seen to provide higher sensitivity in activation detection than linear sHRFs. These observations support using sHRFs in the analysis of IEDs with EEG-fNIRS data.

Childhood extracranial neoplasms: the role of imaging in drug development and clinical trials

Cancer is the leading cause of death in children older than 1 year of age and new drugs are necessary to improve outcomes. Imaging is crucial to the drug development process and assessment of therapeutic response. In adults, tumours are often assessed with CT using size criteria. Unfortunately, techniques established in adults are not necessarily applicable in children due to differing pathophysiology, ability to cooperate and increased susceptibility to ionising radiation. MRI, in particular quantitative MRI, has to date not been fully utilised in children with extracranial neoplasms. The specific challenges of imaging in children, the potential for functional imaging techniques to inform upon and their inclusion in clinical trials are discussed.

Increased hippocampal blood volume and normal blood flow in schizophrenia

Neuroimaging studies have provided compelling evidence for abnormal hippocampal activity in schizophrenia. Most studies made inferences about baseline hippocampal activity using a single hemodynamic parameter (e.g., blood volume or blood flow). Here we studied several hemodynamic measures in the same cohort to test the hypothesis of increased hippocampal activity in schizophrenia. We used dynamic susceptibility contrast- (DSC-) magnetic resonance imaging (MRI) to assess blood volume, blood flow, and mean transit time in the hippocampus of 15 patients with chronic schizophrenia and 15 healthy controls. Left and right hippocampal measurements were combined for absolute measures of cerebral blood volume (CBV), cerebral blood flow (CBF), and mean transit time (MTT). We found significantly increased hippocampal CBV, but normal CBF and MTT, in schizophrenia. The uncoupling of CBV and CBF could be due to several factors, including antipsychotic medication, loss of cerebral perfusion pressure, or angiogenesis. Further studies need to incorporate several complementary imaging modalities to better characterize hippocampal dysfunction in schizophrenia.

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.

Temporal assessment of vascular reactivity and functionality using MRI during postischemic proangiogenenic vascular remodeling

Postischemic angiogenesis is an important recovery mechanism. Both arteries and veins are upregulated during angiogenesis, but eventually there are more angiogenic veins than arteries in terms of number and length. It is critical to understand how the veins are modulated after ischemia and then transitioned into angiogenic vessels during the proangiogenic stage to finally serve as a restorative strength to the injured area. Using a rat model of transient focal cerebral ischemia, the hypercapnic blood oxygen level-dependent (BOLD) response was used to evaluate vascular reactivity, while the hyperoxic BOLD and tissue oxygen level-dependent (TOLD) responses were used to evaluate the vascular functionality at 1, 3, and 7days after ischemia. Vessel-like venous signals appeared on R2* maps on days 3 and 7, but not on day 1. The large hypercapnic BOLD responses on days 3 and 7 indicated that these areas have high vascular reactivity. The temporal correlation between vascular reactivity and the immunoreactivity to desmin and VEGF further indicates that the integrity of vascular reactivity is associated with the pericyte coverage as regulated by the VEGF level. Vascular functionality remained low on days 1, 3, and 7, as reflected by the small hyperoxic BOLD and large hyperoxic TOLD responses, indicating the low oxygen consumption of the ischemic tissues. These functional changes in proangiogenic veins may be critical for angiogenesis.

In vivo quantification of cerebral r2*-response to graded hyperoxia at 3 tesla

Objectives:

This study aims to quantify the response of the transverse relaxation rate of the magnetic resonance (MR) signal of the cerebral tissue in healthy volunteers to the administration of air with step-wise increasing percentage of oxygen.

Materials and Methods:

The transverse relaxation rate (R2*) of the MR signal was quantified in seven volunteers under respiratory intake of normobaric gas mixtures containing 21, 50, 75, and 100% oxygen, respectively. End-tidal breath composition, arterial blood saturation (SaO₂), and heart pulse rate were monitored during the challenge. R2* maps were computed from multi-echo, gradient-echo magnetic resonance imaging (MRI) data, acquired at 3.0T. The average values in the segmented white matter (WM) and gray matter (GM) were tested by the analysis of variance (ANOVA), with Bonferroni post-hoc correction. The GM R2*-reactivity to hyperoxia was modeled using the Hill's equation.

Results:

Graded hyperoxia resulted in a progressive and significant (P < 0.05) decrease of the R2* in GM. Under normoxia the GM-R2* was 17.2 ± 1.1 s^-1. At 75% O₂ supply, the R2* had reached a saturation level, with 16.4 ± 0.7 s^-1 (P = 0.02), without a significant further decrease for 100% O₂. The R2*-response of GM correlated positively with CO₂ partial pressure (R = 0.69 ± 0.19) and negatively with SaO₂ (R = -0.74 ± 0.17). The WM showed a similar progressive, but non-significant, decrease in the relaxation rates, with an increase in oxygen intake (P = 0.055). The Hill's model predicted a maximum R2* response of the GM, of 3.5%, with half the maximum at 68% oxygen concentration.

Conclusions:

The GM-R2* responds to hyperoxia in a concentration-dependent manner, suggesting that monitoring and modeling of the R2*-response may provide new oxygenation biomarkers for tumor therapy or assessment of cerebrovascular reactivity in patients.

Monitoring brain tumor vascular heamodynamic following anti-angiogenic therapy with advanced magnetic resonance imaging in mice

Advanced MR imaging methods have an essential role in classification, grading, follow-up and therapeutic management in patients with brain tumors. With the introduction of new therapeutic options, the challenge for better tissue characterization and diagnosis increase, calling for new reliable non-invasive imaging methods. In the current study we evaluated the added value of a combined protocol of blood oxygen level dependent (BOLD) imaging during hyperoxic challenge (termed hemodynamic response imaging (HRI)) in an orthotopic mouse model for glioblastoma under anti-angiogenic treatment with B20-4.1.1, an anti-VEGF antibody. In glioblastoma tumors, the elevated HRI indicated progressive angiogenesis as further confirmed by histology. In the current glioblastoma model, B20-treatment caused delayed tumor progression with no significant changes in HRI yet with slightly reduced tumor vascularity as indicated by histology. Furthermore, fewer apoptotic cells and higher proliferation index were detected in the B20-treated tumors compared to control-treated tumors. In conclusion, HRI provides an easy, safe and contrast agent free method for the assessment of the brain hemodynamic function, an additionally important clinical information.

Photoacoustic imaging of vascular hemodynamics: validation with blood oxygenation level-dependent MR imaging

PURPOSE

To noninvasively assess vascular hemodynamics with photoacoustic imaging (PAI) and blood oxygenation level-dependent (BOLD) magnetic resonance (MR) imaging in phantoms and in an animal model.

MATERIALS AND METHODS

In vivo studies were performed with institutional animal care and use committee approval. In vitro experiments were performed by using a tissue-mimicking phantom in multiple oxygenation conditions (n = 6) to compare PAI measurements and BOLD MR imaging measurements. PAI and T2-weighted spin-echo-based BOLD MR imaging were performed to assess tumor response to carbogen (95% O2, 5% CO2) in mice with head and neck tumors before (n = 11) and after (n = 9) treatment with a vascular disrupting agent (VDA). Two-tailed Pearson correlation analysis was performed to determine the correlation between the parameters measured with PAI and BOLD MR imaging in vitro. Two-tailed paired t tests were used to compare change in tumor hemoglobin oxygen saturation (sO2) levels and BOLD signal in response to carbogen. Changes in PAI and BOLD signal intensity before and after VDA treatment were analyzed for significance by using analysis of variance with repeated measures.

RESULTS

Phantom measurements yielded good correlation between photoacoustically derived sO2 levels and BOLD signal intensity (r = 0.937, P = .005) and partial pressure of oxygen (r = 0.981, P = .005). In vivo hemodynamic response to carbogen was characterized by a significant increase in tumor sO2 levels (P = .003) and BOLD signal (P = .001). When compared with pretreatment estimates, treatment with VDA resulted in a significant reduction in the tumor hemodynamic response to carbogen at PAI (P = .030).

CONCLUSION

Carbogen-based functional imaging with PAI and BOLD MR imaging enables monitoring of early changes in tumor hemodynamics after vascular targeted therapy.

Multimodal imaging of gliomas in the context of evolving cellular and molecular therapies

The vast majority of malignant gliomas relapse after surgery and standard radio-chemotherapy. Novel molecular and cellular therapies are thus being developed, targeting specific aspects of tumor growth. While histopathology remains the gold standard for tumor classification, neuroimaging has over the years taken a central role in the diagnosis and treatment follow up of brain tumors. It is used to detect and localize lesions, define the target area for biopsies, plan surgical and radiation interventions and assess tumor progression and treatment outcome. In recent years the application of novel drugs including anti-angiogenic agents that affect the tumor vasculature, has drastically modulated the outcome of brain tumor imaging. To properly evaluate the effects of emerging experimental therapies and successfully support treatment decisions, neuroimaging will have to evolve. Multi-modal imaging systems with existing and new contrast agents, molecular tracers, technological advances and advanced data analysis can all contribute to the establishment of disease relevant biomarkers that will improve disease management and patient care. In this review, we address the challenges of glioma imaging in the context of novel molecular and cellular therapies, and take a prospective look at emerging experimental and pre-clinical imaging techniques that bear the promise of meeting these challenges.

Systematic protocol for assessment of the validity of BOLD MRI in a rabbit model of inflammatory arthritis at 1.5 tesla

BACKGROUND

Blood-oxygen-level-dependent (BOLD) MRI has the potential to identify regions of early hypoxic and vascular joint changes in inflammatory arthritis. There is no standard protocol for analysis of BOLD MRI measurements in musculoskeletal disorders.

OBJECTIVE

To optimize the following BOLD MRI reading parameters: (1) statistical threshold values (low, r > 0.01 versus high, r > 0.2); (2) summary measures of BOLD contrast (percentage of activated voxels [PT%] versus percentage signal difference between on-and-off signal intensities [diff_on_off]); and (3) direction of BOLD response (positive, negative and positive + negative).

MATERIALS AND METHODS

Using BOLD MRI protocols at 1.5 T, arthritic (n = 21) and contralateral (n = 21) knees of 21 juvenile rabbits were imaged at baseline and on days 1, 14 and 28 after a unilateral intra-articular injection of carrageenan. Nine non-injected rabbits served as external control knees (n = 18). By comparing arthritic to contralateral knees, receiver operating characteristic curves were used to determine diagnostic accuracy.

RESULTS

Using diff_on_off and positive + negative responses, a threshold of r > 0.01 was more accurate than r > 0.2 (P = 0.03 at day 28). Comparison of summary measures yielded no statistically significant difference (P > 0.05). Although positive + negative (AUC = 0.86 at day 28) and negative responses (AUC = 0.90 at day 28) for PT% were the most diagnostically accurate, positive + negative responses for diff_on_off (AUC = 0.78 at day 28) also had acceptable accuracy.

CONCLUSIONS

The most clinically relevant reading parameters included a lower threshold of r > 0.01 and a positive + negative BOLD response. We propose that diff_on_off is a more clinically relevant summary measure of BOLD MRI, while PT% can be used as an ancillary measure.

Quantitative functional MRI biomarkers improved early detection of colorectal liver metastases

PURPOSE

To implement and evaluate the performance of a computerized statistical tool designed for robust and quantitative analysis of hemodynamic response imaging (HRI) -derived maps for the early identification of colorectal liver metastases (CRLM).

MATERIALS AND METHODS

CRLM-bearing mice were scanned during the early stage of tumor growth and subsequently during the advanced-stage. Three experienced radiologists marked various suspected-foci on the early stage anatomical images and classified each as either highly certain or as suspected tumors. The statistical model construction was based on HRI maps (functional-MRI combined with hypercapnia and hyperoxia) using a supervised learning paradigm which was further trained either with the advanced-stage sets (late training; LT) or with the early stage sets (early training; ET). For each group of foci, the classifier results were compared with the ground-truth.

RESULTS

The ET-based classification significantly improved the manual classification of the highly certain foci (P < 0.05) and was superior compared with the LT-based classification (P < 0.05). Additionally, the ET-based classification, offered high sensitivity (57-63%), accompanied with high positive predictive value (>94%) and high specificity (>98%) for suspected-foci.

CONCLUSION

The ET-based classifier can strengthen the radiologist's classification of highly certain foci. Additionally, it can aid in classifying suspected-foci, thus enabling earlier intervention which can often be lifesaving.

Opportunities and challenges for successful use of bevacizumab in pediatrics

Bevacizumab (Avastin) has rapidly gained status as a broadly active agent for malignancies of several different histologies in adults. This activity has spawned a range of uses in pediatrics for both oncologic and non-oncologic indications. Early analyses indicate that pediatric cancers exhibit a spectrum of responses to bevacizumab that suggest its activity may be more limited than in adult oncology. Most exciting, is that for low-grade tumors that threaten vision and hearing, there is not only evidence for objective tumor response but for recovery of lost function as well. In addition to oncological indications, there is a range of uses for non-oncologic disease for which bevacizumab has clear activity. Finally, a number of mechanisms have been identified as contributing to bevacizumab resistance in cancer. Elucidating these mechanisms will guide the development of future clinical trials of bevacizumab in pediatric oncology.

MRI multiparametric hemodynamic characterization of the normal brain

Characterization of the brain's vascular system is of major clinical importance in the assessment of patients with cerebrovascular disease. The aim of this study was to characterize brain hemodynamics using multiparametric methods and to obtain reference values from the healthy brain. A multimodal magnetic resonance imaging (MRI) study was performed in twenty healthy subjects, including dynamic susceptibility contrast imaging and blood oxygen level dependence (BOLD) during hypercapnia and carbogen challenges. Brain tissues were defined using unsupervised cluster analysis based on these three methods, and several hemodynamic parameters were calculated for gray matter (GM), white matter (WM), blood vessels and dura (BVD); the three main vascular territories within the GM; and arteries and veins defined within the BVD cluster. The carbogen challenge produced a BOLD signal twice as high as the hypercapnia challenge, in all brain tissues. The three brain tissues differed significantly from each other in their hemodynamic characteristics, supporting the graded vascularity of the tissues, with BVD>GM>WM. Within the GM cluster, a significant delay of ∼1.2 s of the bolus arrival time was detected within the posterior cerebral artery territory relative to the middle and anterior cerebral artery territories. No differences were detected between right and left middle cerebral artery territories for all hemodynamic parameters. Within the BVD cluster, a significant delay of ∼1.9 s of the bolus arrival time was detected within the veins relative to the arteries. This parameter enabled to differentiate between the various blood vessels, including arteries, veins and choroid plexus. This study provides reference values for several hemodynamic parameters, obtained from healthy brains, and may be clinically important in the assessment of patients with various vascular pathologies.