Dynamic and simultaneous MR measurement of R1 and R2* changes during respiratory challenges for the assessment of blood and tissue oxygenation

Navigator-based motion compensation for liver BOLD measurement with five-echo SAGE EPI and breath-hold task

PURPOSE

The purpose of this work is to apply multi-echo spin- and gradient-echo (SAGE) echo-planar imaging (EPI) combined with a navigator-based (NAV) prospective motion compensation method for a quantitative liver blood oxygen level dependent (BOLD) measurement with a breath-hold (BH) task.

METHODS

A five-echo SAGE sequence was developed to quantitatively measure T2 and T2* to depict function with sufficient signal-to-noise ratio, spatial resolution and sensitivity to BOLD changes induced by the BH task. To account for respiratory motion, a navigator was employed in the form of a single gradient-echo projection readout, located at the diaphragm along the inferior-superior direction. Prior to each transverse imaging slice of the spin-echo EPI-based readouts, navigator acquisition and fat suppression were incorporated. Motion data was obtained from the navigator and transmitted back to the sequence, allowing real-time adjustments to slice positioning. Six healthy volunteers and three patients with liver carcinoma were included in this study. Quantitative T2 and T2* were calculated at each time point of the BH task. Parameters of t value from first-level analysis using a general linear model and hepatovascular reactivity (HVR) of Echo1, T2 and T2* were calculated.

RESULTS

The motion caused by respiratory activity was successfully compensated using the navigator signal. The average changes of T2 and T2* during breath-hold were about 1% and 0.7%, respectively. With the help of NAV prospective motion compensation whole liver t values could be obtained without motion artifacts. The quantified liver T2 (34.7 ± 0.7 ms) and T2* (29 ± 1.2 ms) values agreed with values from literature. In healthy volunteers, the distribution of statistical t value and HVR was homogeneous throughout the whole liver. In patients with liver carcinoma, the distribution of t value and HVR was inhomogeneous due to metastases or therapy.

CONCLUSIONS

This study demonstrates the feasibility of using a NAV prospective motion compensation technique in conjunction with five-echo SAGE EPI for the quantitative measurement of liver BOLD with a BH task.

Effect of inhaled oxygen level on dynamic glucose-enhanced MRI in mouse brain

PURPOSE

To investigate the effect of inhaled oxygen level on dynamic glucose enhanced (DGE) MRI in mouse brain tissue and CSF at 3 T.

METHODS

DGE data of brain tissue and CSF from mice under normoxia or hyperoxia were acquired in independent and interleaved experiments using on-resonance variable delay multi-pulse (onVDMP) MRI. A bolus of 0.15 mL filtered 50% D-glucose was injected through the tail vein over 1 min during DGE acquisition. MRS was acquired before and after DGE experiments to confirm the presence of D-glucose.

RESULTS

A significantly higher DGE effect under normoxia than under hyperoxia was observed in brain tissue (p = 0.0001 and p = 0.0002 for independent and interleaved experiments, respectively), but not in CSF (p > 0.3). This difference is attributed to the increased baseline MR tissue signal under hyperoxia induced by a shortened T1 and an increased BOLD effect. When switching from hyperoxia to normoxia without glucose injection, a signal change of ˜3.0% was found in brain tissue and a signal change of ˜1.5% was found in CSF.

CONCLUSIONS

DGE signal was significantly lower under hyperoxia than that under normoxia in brain tissue, but not in CSF. The reason is that DGE effect size of brain tissue is affected by the baseline signal, which could be influenced by T1 change and BOLD effect. Therefore, DGE experiments in which the oxygenation level is changed from baseline need to be interpreted carefully.

The effects of normobaric and hyperbaric oxygenation on MRI signal intensities in T1 -weighted, T2 -weighted and FLAIR images in human brain

BACKGROUND

Dissolved oxygen has known paramagnetic effects in magnetic resonance imaging (MRI). The aim of this study was to compare the effects of normobaric oxygenation (NBO) and hyperbaric oxygenation (HBO) on human brain MRI signal intensities.

PATIENTS AND METHODS

Baseline brain MRI was performed in 17 healthy subjects (mean age 27.8 ± 3.2). MRI was repeated after exposure to the NBO and HBO at different time points (0 min, 25 min, 50 min). Signal intensities in T 1-weighted, T 2-weighted images and fluid attenuated inversion recovery (FLAIR) signal intensities of several intracranial structures were compared between NBO and HBO.

RESULTS

Increased T 1-weighted signal intensities were observed in white and deep grey brain matter, cerebrospinal fluid (CSF), venous blood and vitreous body after exposure to NBO as well as to HBO compared to baseline (Dunnett's test, p < 0.05) without significant differences between both protocols. There was also no significant difference in T 2-weighted signal intensities between NBO and HBO. FLAIR signal intensities were increased only in the vitreous body after NBO and HBO and FLAIR signal of caudate nucleus was decreased after NBO (Dunnett's test, p < 0.05). The statistically significant differences in FLAIR signal intensities were found between NBO and HBO (paired t-test, p < 0.05) in most observed brain structures (paired t-test, p < 0.05).

CONCLUSIONS

Our results show that NBO and HBO alters signal intensities T 1-weighted and FLAIR images of human brain. The differences between NBO and HBO are most pronounced in FLAIR imaging.

Functional Imaging of Hypoxia: PET and MRI

Molecular and functional imaging have critical roles in cancer care. Existing evidence suggests that noninvasive detection of hypoxia within a particular type of cancer can provide new information regarding the relationship between hypoxia, cancer aggressiveness and altered therapeutic responses. Following the identification of hypoxia inducible factor (HIF), significant progress in understanding the regulation of hypoxia-induced genes has been made. These advances have provided the ability to therapeutically target HIF and tumor-associated hypoxia. Therefore, by utilizing the molecular basis of hypoxia, hypoxia-based theranostic strategies are in the process of being developed which will further personalize care for cancer patients. The aim of this review is to provide an overview of the significance of tumor hypoxia and its relevance in cancer management as well as to lay out the role of imaging in detecting hypoxia within the context of cancer.

Advances in PET and MRI imaging of tumor hypoxia

Tumor hypoxia is a complex and evolving phenomenon both in time and space. Molecular imaging allows to approach these variations, but the tracers used have their own limitations. PET imaging has the disadvantage of low resolution and must take into account molecular biodistribution, but has the advantage of high targeting accuracy. The relationship between the signal in MRI imaging and oxygen is complex but hopefully it would lead to the detection of truly oxygen-depleted tissue. Different ways of imaging hypoxia are discussed in this review, with nuclear medicine tracers such as [¹⁸F]-FMISO, [¹⁸F]-FAZA, or [⁶⁴Cu]-ATSM but also with MRI techniques such as perfusion imaging, diffusion MRI or oxygen-enhanced MRI. Hypoxia is a pejorative factor regarding aggressiveness, tumor dissemination and resistance to treatments. Therefore, having accurate tools is particularly important.

Oxygen-enhanced MRI and radiotherapy in patients with oropharyngeal squamous cell carcinoma

Background and purpose

This study aimed to assess the role of T1 mapping and oxygen-enhanced MRI in patients undergoing radical dose radiotherapy for HPV positive oropharyngeal cancer, which has not yet been examined in an OE-MRI study.

Materials and methods

Variable Flip Angle T1 maps were acquired on a 3T MRI scanner while patients (n = 12) breathed air and/or 100 % oxygen, before and after fraction 10 of the planned 30 fractions of chemoradiotherapy ('visit 1' and 'visit 2', respectively). The analysis aimed to assess to what extent (1) native R1 relates to patient outcome; (2) OE-MRI response relates to patient outcome; (3) changes in mean R1 before and after radiotherapy related to clinical outcome in patients with oropharyngeal squamous cell carcinoma.

Results

Due to the radiotherapy being largely successful, the sample sizes of non-responder groups were small, and therefore it was not possible to properly assess the predictive nature of OE-MRI. The tumour R1 increased in some patients while decreasing in others, in a pattern that was overall consistent with the underlying OE-MRI theory and previously reported tumour OE-MRI responses. In addition, we discuss some practical challenges faced when integrating this technique into a clinical trial, with the aim that sharing this is helpful to researchers planning to use OE-MRI in future clinical studies.

Conclusion

Altogether, these results suggest that further clinical OE-MRI studies to assess hypoxia and radiotherapy response are worth pursuing, and that there is important work to be done to improve the robustness of the OE-MRI technique in human applications in order for it to be useful as a widespread clinical technique.

The Role of Imaging Biomarkers to Guide Pharmacological Interventions Targeting Tumor Hypoxia

Hypoxia is a common feature of solid tumors that contributes to angiogenesis, invasiveness, metastasis, altered metabolism and genomic instability. As hypoxia is a major actor in tumor progression and resistance to radiotherapy, chemotherapy and immunotherapy, multiple approaches have emerged to target tumor hypoxia. It includes among others pharmacological interventions designed to alleviate tumor hypoxia at the time of radiation therapy, prodrugs that are selectively activated in hypoxic cells or inhibitors of molecular targets involved in hypoxic cell survival (i.e., hypoxia inducible factors HIFs, PI3K/AKT/mTOR pathway, unfolded protein response). While numerous strategies were successful in pre-clinical models, their translation in the clinical practice has been disappointing so far. This therapeutic failure often results from the absence of appropriate stratification of patients that could benefit from targeted interventions. Companion diagnostics may help at different levels of the research and development, and in matching a patient to a specific intervention targeting hypoxia. In this review, we discuss the relative merits of the existing hypoxia biomarkers, their current status and the challenges for their future validation as companion diagnostics adapted to the nature of the intervention.

Target Definition in MR-Guided Adaptive Radiotherapy for Head and Neck Cancer

In recent years, MRI-guided radiotherapy (MRgRT) has taken an increasingly important position in image-guided radiotherapy (IGRT). Magnetic resonance imaging (MRI) offers superior soft tissue contrast in anatomical imaging compared to computed tomography (CT), but also provides functional and dynamic information with selected sequences. Due to these benefits, in current clinical practice, MRI is already used for target delineation and response assessment in patients with head and neck squamous cell carcinoma (HNSCC). Because of the close proximity of target areas and radiosensitive organs at risk (OARs) during HNSCC treatment, MRgRT could provide a more accurate treatment in which OARs receive less radiation dose. With the introduction of several new radiotherapy techniques (i.e., adaptive MRgRT, proton therapy, adaptive cone beam computed tomography (CBCT) RT, (daily) adaptive radiotherapy ensures radiation dose is accurately delivered to the target areas. With the integration of a daily adaptive workflow, interfraction changes have become visible, which allows regular and fast adaptation of target areas. In proton therapy, adaptation is even more important in order to obtain high quality dosimetry, due to its susceptibility for density differences in relation to the range uncertainty of the protons. The question is which adaptations during radiotherapy treatment are oncology safe and at the same time provide better sparing of OARs. For an optimal use of all these new tools there is an urgent need for an update of the target definitions in case of adaptive treatment for HNSCC. This review will provide current state of evidence regarding adaptive target definition using MR during radiotherapy for HNSCC. Additionally, future perspectives for adaptive MR-guided radiotherapy will be discussed.

A General Model to Calculate the Spin-Lattice Relaxation Rate (R1) of Blood, Accounting for Hematocrit, Oxygen Saturation, Oxygen Partial Pressure, and Magnetic Field Strength Under Hyperoxic Conditions

BACKGROUND

Under normal physiological conditions, the spin-lattice relaxation rate (R1) in blood is influenced by many factors, including hematocrit, field strength, and the paramagnetic effects of deoxyhemoglobin and dissolved oxygen. In addition, techniques such as oxygen-enhanced magnetic resonance imaging (MRI) require high fractions of inspired oxygen to induce hyperoxia, which complicates the R1 signal further. A quantitative model relating total blood oxygen content to R1 could help explain these effects.

PURPOSE

To propose and assess a general model to estimate the R1 of blood, accounting for hematocrit, SO2 , PO2 , and B0 under both normal physiological and hyperoxic conditions.

STUDY TYPE

Mathematical modeling.

POPULATION

One hundred and twenty-six published values of R1 from phantoms and animal models.

FIELD STRENGTH/SEQUENCE

5-8.45 T.

ASSESSMENT

We propose a two-compartment nonlinear model to calculate R1 as a function of hematocrit, PO2 , and B0. The Akaike Information Criterion (AIC) was used to select the best-performing model with the fewest parameters. A previous model of R1 as a function of hematocrit, SO2 , and B0 has been proposed by Hales et al, and our work builds upon this work to make the model applicable under hyperoxic conditions (SO2  > 0.99). Models were assessed using the AIC, mean squared error (MSE), coefficient of determination (R2 ), and Bland-Altman analysis. The effect of volume fraction constants W RBC and W plasma was assessed by the SD of resulting R1. The range of the model was determined by the maximum and minimum B0, hematocrit, SO2 , and PO2 of the literature data points.

STATISTICAL TESTS

Bland-Altman, AIC, MSE, coefficient of determination (R2 ), SD.

RESULTS

The model estimates agreed well with the literature values of R1 of blood (R2  = 0.93, MSE = 0.0013 s-2 ), and its performance was consistent across the range of parameters: B0 = 1.5-8.45 T, SO2  = 0.40-1, PO2  = 30-700 mmHg.

DATA CONCLUSION

Using the results from this model, we have quantified and explained the contradictory decrease in R1 reported in oxygen-enhanced MRI and oxygen-delivery experiments.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

Using Variable Flip Angle (VFA) and Modified Look-Locker Inversion Recovery (MOLLI) T1 mapping in clinical OE-MRI

BACKGROUND AND PURPOSE

The imaging technique known as Oxygen-Enhanced MRI is under development as a noninvasive technique for imaging hypoxia in tumours and pulmonary diseases. While promising results have been shown in preclinical experiments, clinical studies have mentioned experiencing difficulties with patient motion, image registration, and the limitations of single-slice images compared to 3D volumes. As clinical studies begin to assess feasibility of using OE-MRI in patients, it is important for researchers to communicate about the practical challenges experienced when using OE-MRI on patients to help the technique advance.

MATERIALS AND METHODS

We report on our experience with using two types of T1 mapping (MOLLI and VFA) for a recently completed OE-MRI clinical study on oropharyngeal squamous cell carcinoma.

RESULTS

We report: (1) the artefacts and practical difficulties encountered in this study; (2) the difference in estimated T1 from each method used - the VFA T1 estimation was higher than the MOLLI estimation by 27% on average; (3) the standard deviation within the tumour ROIs - there was no significant difference in the standard deviation seen within the tumour ROIs from the VFA versus MOLLI; and (4) the OE-MRI response collected from either method. Lastly, we collated the MRI acquisition details from over 45 relevant manuscripts as a convenient reference for researchers planning future studies.

CONCLUSION

We have reported our practical experience from an OE-MRI clinical study, with the aim that sharing this is helpful to researchers planning future studies. In this study, VFA was a more useful technique for using OE-MRI in tumours than MOLLI T1 mapping.

Oxygen-sensitive MRI assessment of tumor response to hypoxic gas breathing challenge

Oxygen-sensitive MRI has been extensively used to investigate tumor oxygenation based on the response (R₂ * and/or R₁ ) to a gas breathing challenge. Most studies have reported response to hyperoxic gas indicating potential biomarkers of hypoxia. Few studies have examined hypoxic gas breathing and we have now evaluated acute dynamic changes in rat breast tumors. Rats bearing syngeneic subcutaneous (n = 15) or orthotopic (n = 7) 13762NF breast tumors were exposed to a 16% O₂ gas breathing challenge and monitored using blood oxygen level dependent (BOLD) R₂ * and tissue oxygen level dependent (TOLD) T₁ -weighted measurements at 4.7 T. As a control, we used a traditional hyperoxic gas breathing challenge with 100% O₂ on a subset of the subcutaneous tumor bearing rats (n = 6). Tumor subregions identified as responsive on the basis of R₂ * dynamics coincided with the viable tumor area as judged by subsequent H&E staining. As expected, R₂ * decreased and T₁ -weighted signal increased in response to 100% O₂ breathing challenge. Meanwhile, 16% O₂ breathing elicited an increase in R₂ *, but divergent response (increase or decrease) in T₁ -weighted signal. The T₁ -weighted signal increase may signify a dominating BOLD effect triggered by 16% O₂ in the relatively more hypoxic tumors, whereby the influence of increased paramagnetic deoxyhemoglobin outweighs decreased pO₂ . The results emphasize the importance of combined BOLD and TOLD measurements for the correct interpretation of tumor oxygenation properties.

Imaging tumour hypoxia with oxygen-enhanced MRI and BOLD MRI

Hypoxia is known to be a poor prognostic indicator for nearly all solid tumours and also is predictive of treatment failure for radiotherapy, chemotherapy, surgery and targeted therapies. Imaging has potential to identify, spatially map and quantify tumour hypoxia prior to therapy, as well as track changes in hypoxia on treatment. At present no hypoxia imaging methods are available for routine clinical use. Research has largely focused on positron emission tomography (PET)-based techniques, but there is gathering evidence that MRI techniques may provide a practical and more readily translational alternative. In this review we focus on the potential for imaging hypoxia by measuring changes in longitudinal relaxation [R1; termed oxygen-enhanced MRI or tumour oxygenation level dependent (TOLD) MRI] and effective transverse relaxation [R2*; termed blood oxygenation level dependent (BOLD) MRI], induced by inhalation of either 100% oxygen or the radiosensitising hyperoxic gas carbogen. We explain the scientific principles behind oxygen-enhanced MRI and BOLD and discuss significant studies and their limitations. All imaging biomarkers require rigorous validation in order to translate into clinical use and the steps required to further develop oxygen-enhanced MRI and BOLD MRI into decision-making tools are discussed.

The effect of flow on blood oxygen level dependent (R* 2 ) MRI of orthotopic lung tumors

PURPOSE

Blood oxygen level dependent (BOLD) MRI based on R 2 * measurements can provide insights into tumor vascular oxygenation. However, measurements are susceptible to blood flow, which may vary accompanying a hyperoxic gas challenge. We investigated flow sensitivity by comparing R 2 * measurements with and without flow suppression (fs) in 2 orthotopic lung xenograft tumor models.

METHODS

H460 (n = 20) and A549 (n = 20) human lung tumor xenografts were induced by surgical implantation of cancer cells in the right lung of nude rats. MRI was performed at 4.7T after tumors reached 5 to 8 mm in diameter. A multiecho gradient echo MRI sequence was acquired with and without spatial saturation bands on each side of the imaging plane to evaluate the effect of flow on R 2 * . fs and non-fs R 2 * MRI measurements were interleaved during an oxygen breathing challenge (from air to 100% O₂ ). T 2 * -weighted signal intensity changes (ΔSI(%)) and R 2 * measurements were obtained for regions of interest and on a voxel-by-voxel basis and discrepancies quantified with Bland-Altman analysis.

RESULTS

Flow suppression affected ΔSI(%) and R 2 * measurements in each tumor model. Average discrepancy and limits of agreement from Bland-Altman analyses revealed greater flow-related bias in A549 than H460.

CONCLUSION

The effect of flow on R 2 * , and hence BOLD, was tumor model dependent with measurements being more sensitive in well-perfused A549 tumors.

Brain Morphometry and Longitudinal Relaxation Time of Spontaneously Hypertensive Rats (SHRs) in Early and Intermediate Stages of Hypertension Investigated by 3D VFA-SPGR MRI

Cerebral small vessel disease(s) (SVD) results from pathological changes of the small blood vessels in the brain and is common in older people. The diagnostic features by which SVD manifests in brain includes white matter hyperintensities, lacunes, dilated perivascular spaces, microbleeds, and atrophy. In the present study, we use in vivo magnetic resonance imaging (MRI) to characterize brain morphometry and longitudinal relaxation time (T1) of spontaneously hypertensive rats (SHRs) to study the contribution of chronic hypertension to SVD relevant pathology. Male SHR and Wistar-Kyoto (WKY) rats underwent 3D variable flip angle spoiled gradient echo brain MRI at 9.4 T at early (seven weeks old) and established (19 weeks old) stages of hypertension. The derived proton density weighted and T1 images were utilized for morphometry and to characterize T1 properties in gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF). Custom tissue probability maps were constructed for accurate computerized whole brain tissue segmentations and voxel-wise analyses. Characteristic morphological differences between the two strains included enlarged ventricles, smaller corpus callosum (CC) volumes and general 'thinning' of CC in SHR compared to WKY rats at both age groups. While we did not observe parenchymal T1 differences, the T1 of CSF was elevated in SHR compared to controls. Collectively these findings indicate that SHRs develop WM atrophy which is a clinically robust MRI biomarker associated with WM degeneration.

Mapping Hypoxia in Renal Carcinoma with Oxygen-enhanced MRI: Comparison with Intrinsic Susceptibility MRI and Pathology

Purpose To cross-validate T1-weighted oxygen-enhanced (OE) MRI measurements of tumor hypoxia with intrinsic susceptibility MRI measurements and to demonstrate the feasibility of translation of the technique for patients. Materials and Methods Preclinical studies in nine 786-0-R renal cell carcinoma (RCC) xenografts and prospective clinical studies in eight patients with RCC were performed. Longitudinal relaxation rate changes (∆R1) after 100% oxygen inhalation were quantified, reflecting the paramagnetic effect on tissue protons because of the presence of molecular oxygen. Native transverse relaxation rate (R2*) and oxygen-induced R2* change (∆R2*) were measured, reflecting presence of deoxygenated hemoglobin molecules. Median and voxel-wise values of ∆R1 were compared with values of R2* and ∆R2*. Tumor regions with dynamic contrast agent-enhanced MRI perfusion, refractory to signal change at OE MRI (referred to as perfused Oxy-R), were distinguished from perfused oxygen-enhancing (perfused Oxy-E) and nonperfused regions. R2* and ∆R2* values in each tumor subregion were compared by using one-way analysis of variance. Results Tumor-wise and voxel-wise ∆R1 and ∆R2* comparisons did not show correlative relationships. In xenografts, parcellation analysis revealed that perfused Oxy-R regions had faster native R2* (102.4 sec-1 vs 81.7 sec-1) and greater negative ∆R2* (-22.9 sec-1 vs -5.4 sec-1), compared with perfused Oxy-E and nonperfused subregions (all P < .001), respectively. Similar findings were present in human tumors (P < .001). Further, perfused Oxy-R helped identify tumor hypoxia, measured at pathologic analysis, in both xenografts (P = .002) and human tumors (P = .003). Conclusion Intrinsic susceptibility biomarkers provide cross validation of the OE MRI biomarker perfused Oxy-R. Consistent relationship to pathologic analyses was found in xenografts and human tumors, demonstrating biomarker translation. Published under a CC BY 4.0 license. Online supplemental material is available for this article.

Imaging biomarkers in oncology: Basics and application to MRI

Cancer remains a global killer alongside cardiovascular disease. A better understanding of cancer biology has transformed its management with an increasing emphasis on a personalized approach, so-called "precision cancer medicine." Imaging has a key role to play in the management of cancer patients. Imaging biomarkers that objectively inform on tumor biology, the tumor environment, and tumor changes in response to an intervention complement genomic and molecular diagnostics. In this review we describe the key principles for imaging biomarker development and discuss the current status with respect to magnetic resonance imaging (MRI).

LEVEL OF EVIDENCE

5 TECHNICAL EFFICACY: Stage 5 J. Magn. Reson. Imaging 2018;48:13-26.

Data-driven mapping of hypoxia-related tumor heterogeneity using DCE-MRI and OE-MRI

PURPOSE

Previous work has shown that combining dynamic contrast-enhanced (DCE)-MRI and oxygen-enhanced (OE)-MRI binary enhancement maps can identify tumor hypoxia. The current work proposes a novel, data-driven method for mapping tissue oxygenation and perfusion heterogeneity, based on clustering DCE/OE-MRI data.

METHODS

DCE-MRI and OE-MRI were performed on nine U87 (glioblastoma) and seven Calu6 (non-small cell lung cancer) murine xenograft tumors. Area under the curve and principal component analysis features were calculated and clustered separately using Gaussian mixture modelling. Evaluation metrics were calculated to determine the optimum feature set and cluster number. Outputs were quantitatively compared with a previous non data-driven approach.

RESULTS

The optimum method located six robustly identifiable clusters in the data, yielding tumor region maps with spatially contiguous regions in a rim-core structure, suggesting a biological basis. Mean within-cluster enhancement curves showed physiologically distinct, intuitive kinetics of enhancement. Regions of DCE/OE-MRI enhancement mismatch were located, and voxel categorization agreed well with the previous non data-driven approach (Cohen's kappa = 0.61, proportional agreement = 0.75).

CONCLUSION

The proposed method locates similar regions to the previous published method of binarization of DCE/OE-MRI enhancement, but renders a finer segmentation of intra-tumoral oxygenation and perfusion. This could aid in understanding the tumor microenvironment and its heterogeneity. Magn Reson Med 79:2236-2245, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Tumor physiological changes during hypofractionated stereotactic body radiation therapy assessed using multi-parametric magnetic resonance imaging

Radiation therapy is a primary treatment for non-resectable lung cancer and hypoxia is thought to influence tumor response. Hypoxia is expected to be particularly relevant to the evolving new radiation treatment scheme of hypofractionated stereotactic body radiation therapy (SBRT). As such, we sought to develop non-invasive tools to assess tumor pathophysiology and response to irradiation. We applied blood oxygen level dependent (BOLD) and tissue oxygen level dependent (TOLD) MRI, together with dynamic contrast enhanced (DCE) MRI to explore the longitudinal effects of SBRT on tumor oxygenation and vascular perfusion using A549 human lung cancer xenografts in a subcutaneous rat model. Intra-tumor heterogeneity was seen on multi-parametric maps, especially in BOLD, T2* and DCE. At baseline, most tumors showed a positive BOLD signal response (%ΔSI) and increased T2* in response to oxygen breathing challenge, indicating increased vascular oxygenation. Control tumors showed similar response 24 hours and 1 week later. Twenty-four hours after a single dose of 12 Gy, the irradiated tumors showed a significantly decreased T2* (-2.9±4.2 ms) and further decrease was observed (-4.0±6.0 ms) after 1 week, suggesting impaired vascular oxygenation. DCE revealed tumor heterogeneity, but showed minimal changes following irradiation. Rats were cured of the primary tumors by 3x12 Gy, providing long term survival, though with ultimate metastatic recurrence.

Lesion magnetic susceptibility response to hyperoxic challenge: A biomarker for malignant brain tumor microenvironment?

BACKGROUND AND PURPOSE

Quantitative susceptibility mapping has been previously used to differentiate lesions in patients with brain tumors. The aim of this work was to characterize the response of magnetic susceptibility differences in malignant brain tumors and surrounding edema to hyperoxic and hypercapnic respiratory challenges.

METHODS

Images of malignant brain tumor patients (2 glioblastoma multiforme, 2 anaplastic astrocytoma, 1 brain metastasis) with clinical MRI exams (contrast-enhanced T1w) were acquired at 3T. 3D multi-gradient-echo data sets were acquired while the patients inhaled medical-air (21% O2), oxygen (100% O2), and carbogen (95% O2, 5% CO2). Susceptibility maps were generated from real and imaginary data. Regions of interest were analyzed with respect to respiration-gas-induced susceptibility changes.

RESULTS

Contrast-enhancing tumor regions with high baseline magnetic susceptibility exhibited a marked susceptibility reduction under hyperoxic challenges, with a stronger effect (-0.040 to -0.100ppm) under hypercapnia compared to hyperoxia (-0.010 to -0.067ppm). In contrast, regions attributed to necrotic tissue and to edema showed smaller changes of opposite sign, i.e. paramagnetic shift. There was a correlation between malignant tumor tissue magnetic susceptibility at baseline under normoxia and the corresponding susceptibility reduction under hypercapnia and - to a lesser degree - under hyperoxia.

CONCLUSION

In this small cohort of analysis, quantification of susceptibility changes in response to respiratory challenges allowed a complementary, functional differentiation of tumorous sub-regions. Those changes, together with the correlations observed between baseline susceptibility under normoxia and susceptibility reduction with challenges, could prove helpful for a non-invasive characterization of local tumor microenvironment.

Combined endogenous MR biomarkers to predict basal tumor oxygenation and response to hyperoxic challenge

Hypoxia is a common feature of solid tumors, which translates into increased angiogenesis, malignant phenotype cell selection, change in gene expression and greater resistance to radiotherapy and chemotherapy. Therefore, there is a need for markers of hypoxia to stratify patients, in order to personalize treatment to improve therapeutic outcome. However, no modality has yet been validated for the screening of hypoxia in routine clinical practice. Magnetic resonance imaging (MRI) R₁ and R₂ * relaxation parameters are sensitive to tissue oxygenation: R₁ is sensitive to dissolved oxygen and R₂ * is sensitive to intravascular deoxyhemoglobin content. Two rat tumor models with distinct levels of hypoxia, 9L-glioma and rhabdomyosarcoma, were imaged for R₁ and R₂ * under air and carbogen (95% O₂ and 5% CO₂ ) breathing conditions. It was observed that the basal tumor oxygenation level had an impact on the amplitude of response to carbogen in the vascular compartment (R₂ *), but not in the tissue compartment (R₁ ). In addition, the change in tissue oxygenation estimated by ΔR₁ correlated with the change in vascular oxygenation estimated by ΔR₂ *, which is consistent with an increase in oxygen supply generating an elevated tumor pO₂ . At the intra-tumoral level, we identified four types of voxel to which a hypoxic feature was attributed (mild hypoxia, severe hypoxia, normoxia and vascular steal), depending on the carbogen-induced change in R₁ and R₂ * values for each voxel. The results showed that 9L-gliomas present more normoxic fractions, whereas rhabdomyosarcomas present more hypoxic fractions, which is in accordance with a previous study using ¹⁸ F-fluoroazomycin arabinoside (¹⁸ F-FAZA) and electron paramagnetic resonance (EPR) oximetry. The response of the combined endogenous MRI contrasts to carbogen challenge could be a useful tool to predict different tumor hypoxic fractions.

Incorporating Oxygen-Enhanced MRI into Multi-Parametric Assessment of Human Prostate Cancer

Hypoxia is associated with prostate tumor aggressiveness, local recurrence, and biochemical failure. Magnetic resonance imaging (MRI) offers insight into tumor pathophysiology and recent reports have related transverse relaxation rate (R₂*) and longitudinal relaxation rate (R₁) measurements to tumor hypoxia. We have investigated the inclusion of oxygen-enhanced MRI for multi-parametric evaluation of tumor malignancy. Multi-parametric MRI sequences at 3 Tesla were evaluated in 10 patients to investigate hypoxia in prostate cancer prior to radical prostatectomy. Blood oxygen level dependent (BOLD), tissue oxygen level dependent (TOLD), dynamic contrast enhanced (DCE), and diffusion weighted imaging MRI were intercorrelated and compared with the Gleason score. The apparent diffusion coefficient (ADC) was significantly lower in tumor than normal prostate. Baseline R₂* (BOLD-contrast) was significantly higher in tumor than normal prostate. Upon the oxygen breathing challenge, R₂* decreased significantly in the tumor tissue, suggesting improved vascular oxygenation, however changes in R₁ were minimal. R₂* of contralateral normal prostate decreased in most cases upon oxygen challenge, although the differences were not significant. Moderate correlation was found between ADC and Gleason score. ADC and R₂* were correlated and trends were found between Gleason score and R₂*, as well as maximum-intensity-projection and area-under-the-curve calculated from DCE. Tumor ADC and R₂* have been associated with tumor hypoxia, and thus the correlations are of particular interest. A multi-parametric approach including oxygen-enhanced MRI is feasible and promises further insights into the pathophysiological information of tumor microenvironment.

Comparison between model-predicted tumor oxygenation dynamics and vascular-/flow-related Doppler indices

PURPOSE

Mathematical modeling is a powerful and flexible method to investigate complex phenomena. It discloses the possibility of reproducing expensive as well as invasive experiments in a safe environment with limited costs. This makes it suitable to mimic tumor evolution and response to radiotherapy although the reliability of the results remains an issue. Complexity reduction is therefore a critical aspect in order to be able to compare model outcomes to clinical data. Among the factors affecting treatment efficacy, tumor oxygenation is known to play a key role in radiotherapy response. In this work, we aim at relating the oxygenation dynamics, predicted by a macroscale model trained on tumor volumetric data of uterine cervical cancer patients, to vascularization and blood flux indices assessed on Ultrasound Doppler images.

METHODS

We propose a macroscale model of tumor evolution based on three dynamics, namely active portion, necrotic portion, and oxygenation. The model parameters were assessed on the volume size of seven cervical cancer patients administered with 28 fractions of intensity modulated radiation therapy (IMRT) (1.8 Gy/fraction). For each patient, five Doppler ultrasound tests were acquired before, during, and after the treatment. The lesion was manually contoured by an expert physician using 4D View^® (General Electric Company - Fairfield, Connecticut, United States), which automatically provided the overall tumor volume size along with three vascularization and/or blood flow indices. Volume data only were fed to the model for training purpose, while the predicted oxygenation was compared a posteriori to the measured Doppler indices.

RESULTS

The model was able to fit the tumor volume evolution within 8% error (range: 3-8%). A strong correlation between the intrapatient longitudinal indices from Doppler measurements and oxygen predicted by the model (about 90% or above) was found in three cases. Two patients showed an average correlation value (50-70%) and the remaining two presented poor correlations. The latter patients were the ones featuring the smallest tumor reduction throughout the treatment, typical of hypoxic conditions. Moreover, the average oxygenation value predicted by the model was close to the average vascularization-flow index (average difference: 7%).

CONCLUSIONS

The results suggest that the modeled relation between tumor evolution and oxygen dynamics was reasonable enough to provide realistic oxygenation curves in five cases (correlation greater than 50%) out of seven. In case of nonresponsive tumors, the model failed in predicting the oxygenation trend while succeeded in reproducing the average oxygenation value according to the mean vascularization-flow index. Despite the need for deeper investigations, the outcomes of the present work support the hypothesis that a simple macroscale model of tumor response to radiotherapy is able to predict the tumor oxygenation. The possibility of an objective and quantitative validation on imaging data discloses the possibility to translate them as decision support tools in clinical practice and to move a step forward in the treatment personalization.

Cancer heterogeneity and imaging

There is interest in identifying and quantifying tumor heterogeneity at the genomic, tissue pathology and clinical imaging scales, as this may help better understand tumor biology and may yield useful biomarkers for guiding therapy-based decision making. This review focuses on the role and value of using x-ray, CT, MRI and PET based imaging methods that identify, measure and map tumor heterogeneity. In particular we highlight the potential value of these techniques and the key challenges required to validate and qualify these biomarkers for clinical use.

Quantitative T1 mapping under precisely controlled graded hyperoxia at 7T

Increasing the concentration of oxygen dissolved in water is known to increase the recovery rate (R1 = 1/T1) of longitudinal magnetization (T1 relaxation). Direct T1 changes in response to precise hyperoxic gas challenges have not yet been quantified and the actual effect of increasing arterial oxygen concentration on the T1 of brain parenchyma remains unclear. The aim of this work was to use quantitative T1 mapping to measure tissue T1 changes in response to precisely targeted hyperoxic respiratory challenges ranging from baseline end-tidal oxygen (PetO₂) to approximately 500 mmHg. We did not observe measureable T1 changes in either gray matter or white matter parenchymal tissue. The T1 of peripheral cerebrospinal fluid located within the sulci, however, was reduced as a function of PetO₂. No significant T1 changes were observed in the ventricular cerebrospinal fluid under hyperoxia. Our results indicate that care should be taken to distinguish actual T1 changes from those which may be related to partial volume effects with cerebrospinal fluid, or regions with increased fluid content such as edema when examining hyperoxia-induced changes in T1 using methods based on T1-weighted imaging.

Developing oxygen-enhanced magnetic resonance imaging as a prognostic biomarker of radiation response

Oxygen-Enhanced Magnetic Resonance Imaging (OE-MRI) techniques were evaluated as potential non-invasive predictive biomarkers of radiation response. Semi quantitative blood-oxygen level dependent (BOLD) and tissue oxygen level dependent (TOLD) contrast, and quantitative responses of relaxation rates (ΔR1 and ΔR2*) to an oxygen breathing challenge during hypofractionated radiotherapy were applied. OE-MRI was performed on subcutaneous Dunning R3327-AT1 rat prostate tumors (n=25) at 4.7 T prior to each irradiation (2F × 15 Gy) to the gross tumor volume. Response to radiation, while inhaling air or oxygen, was assessed by tumor growth delay measured up to four times the initial irradiated tumor volume (VQT). Radiation-induced hypoxia changes were confirmed using a double hypoxia marker assay. Inhaling oxygen during hypofractionated radiotherapy significantly improved radiation response. A correlation was observed between the difference in the 2nd and 1st ΔR1 (ΔΔR1) and VQT for air breathing rats. The TOLD response before the 2nd fraction showed a moderate correlation with VQT for oxygen breathing rats. The correlations indicate useful prognostic factors to predict tumor response to hypofractionation and could readily be applied for patient stratification and personalized radiotherapy treatment planning.

Vascular Magnetic Resonance Imaging in Brain Tumors During Antiangiogenic Therapy--Are We There Yet?

Abnormal tumor vasculature is a potent mediator of treatment resistance because it results in heterogeneous perfusion, hypoxia, increased interstitial fluid pressure, and incomplete penetration of cytotoxic chemotherapies. Targeting this abnormal tumor vasculature is a promising therapeutic strategy, but results with antiangiogenic drugs in brain cancer have been mixed. Vasculature's response to treatment is a dynamic physiological process that can change rapidly throughout treatment, so it requires noninvasive techniques to serially monitor these changes in order to improve outcome. We review the role of vascular magnetic resonance imaging to measure tumor response to treatment and highlight opportunities and future avenues for expanding these promising techniques.

Tumor radio-sensitivity assessment by means of volume data and magnetic resonance indices measured on prostate tumor bearing rats

PURPOSE

Radiation therapy is one of the most common treatments in the fight against prostate cancer, since it is used to control the tumor (early stages), to slow its progression, and even to control pain (metastasis). Although many factors (e.g., tumor oxygenation) are known to influence treatment efficacy, radiotherapy doses and fractionation schedules are often prescribed according to the principle "one-fits-all," with little personalization. Therefore, the authors aim at predicting the outcome of radiation therapy a priori starting from morphologic and functional information to move a step forward in the treatment customization.

METHODS

The authors propose a two-step protocol to predict the effects of radiation therapy on individual basis. First, one macroscopic mathematical model of tumor evolution was trained on tumor volume progression, measured by caliper, of eighteen Dunning R3327-AT1 bearing rats. Nine rats inhaled 100% O2 during irradiation (oxy), while the others were allowed to breathe air. Second, a supervised learning of the weight and biases of two feedforward neural networks was performed to predict the radio-sensitivity (target) from the initial volume and oxygenation-related information (inputs) for each rat group (air and oxygen breathing). To this purpose, four MRI-based indices related to blood and tissue oxygenation were computed, namely, the variation of signal intensity ΔSI in interleaved blood oxygen level dependent and tissue oxygen level dependent (IBT) sequences as well as changes in longitudinal ΔR1 and transverse ΔR2(*) relaxation rates.

RESULTS

An inverse correlation of the radio-sensitivity parameter, assessed by the model, was found with respect the ΔR2(*) (-0.65) for the oxy group. A further subdivision according to positive and negative values of ΔR2(*) showed a larger average radio-sensitivity for the oxy rats with ΔR2(*)<0 and a significant difference in the two distributions (p < 0.05). Finally, a leave-one-out procedure yielded a radio-sensitivity error lower than 20% in both neural networks.

CONCLUSIONS

While preliminary, these specific results suggest that subjects affected by the same pathology can benefit differently from the same irradiation modalities and support the usefulness of IBT in discriminating between different responses.

Oxygen-Enhanced MRI Accurately Identifies, Quantifies, and Maps Tumor Hypoxia in Preclinical Cancer Models

There is a clinical need for noninvasive biomarkers of tumor hypoxia for prognostic and predictive studies, radiotherapy planning, and therapy monitoring. Oxygen-enhanced MRI (OE-MRI) is an emerging imaging technique for quantifying the spatial distribution and extent of tumor oxygen delivery in vivo. In OE-MRI, the longitudinal relaxation rate of protons (ΔR1) changes in proportion to the concentration of molecular oxygen dissolved in plasma or interstitial tissue fluid. Therefore, well-oxygenated tissues show positive ΔR1. We hypothesized that the fraction of tumor tissue refractory to oxygen challenge (lack of positive ΔR1, termed "Oxy-R fraction") would be a robust biomarker of hypoxia in models with varying vascular and hypoxic features. Here, we demonstrate that OE-MRI signals are accurate, precise, and sensitive to changes in tumor pO2 in highly vascular 786-0 renal cancer xenografts. Furthermore, we show that Oxy-R fraction can quantify the hypoxic fraction in multiple models with differing hypoxic and vascular phenotypes, when used in combination with measurements of tumor perfusion. Finally, Oxy-R fraction can detect dynamic changes in hypoxia induced by the vasomodulator agent hydralazine. In contrast, more conventional biomarkers of hypoxia (derived from blood oxygenation-level dependent MRI and dynamic contrast-enhanced MRI) did not relate to tumor hypoxia consistently. Our results show that the Oxy-R fraction accurately quantifies tumor hypoxia noninvasively and is immediately translatable to the clinic.

Feasibility and reproducibility of BOLD and TOLD measurements in the liver with oxygen and carbogen gas challenge in healthy volunteers and patients with hepatocellular carcinoma

PURPOSE

To quantify baseline relaxation rates R2* and R1 in the abdomen, their changes after respiratory challenges, and their reproducibility in healthy volunteers and patients with hepatocellular carcinoma (HCC) at 1.5T and 3.0T.

MATERIALS AND METHODS

R2* measurements were acquired in the liver in 8 volunteers and 27 patients with 34 HCCs using multiecho T2* at baseline and after respiratory challenges with 100% oxygen (O2 ) and carbogen (CB = 95%O2 /5%CO2 ). R1 was measured at 1.5T in one volunteer and 21 patients with 23 HCCs. Test-retest coefficient of variation (CV) was assessed in 10 subjects. Intra- and interobserver variability of R2* and R1 measurements was assessed in 12 and 10 patients, respectively. Parameters for HCC, liver, and muscle were compared between baseline and after gas challenges.

RESULTS

We observed that R2* and R1 imaging of HCCs with O2 and CB is feasible and reproducible (test-retest CV R2*<15%/R1 <5%; intra- and interobserver intraclass correlation coefficient R2*>0.88/R1 >0.7 and CV R2*<7%/R1 <3% at 1.5T). R2* measurements were observed to be less reproducible at 3.0T (CV<35%). There was a statistically significant decrease in R2* values in HCC before and after O2 (P = 0.02) and increase in R1 after O2 (P = 0.004). CB had no significant effect (P R2* = 0.47/R1  = 0.278).

CONCLUSION

R2* measurements in HCC and liver parenchyma are more reproducible at 1.5T than at 3.0T, and with O2 than with CB challenge. We observed a decrease in R2* and an increase in R1 of HCCs from baseline in response to O2 challenge, as expected with increased tissue and blood oxygenation.

Dynamic oxygen challenge evaluated by NMR T1 and T2*--insights into tumor oxygenation

There is intense interest in developing non-invasive prognostic biomarkers of tumor response to therapy, particularly with regard to hypoxia. It has been suggested that oxygen sensitive MRI, notably blood oxygen level-dependent (BOLD) and tissue oxygen level-dependent (TOLD) contrast, may provide relevant measurements. This study examined the feasibility of interleaved T2*- and T1-weighted oxygen sensitive MRI, as well as R2* and R1 maps, of rat tumors to assess the relative sensitivity to changes in oxygenation. Investigations used cohorts of Dunning prostate R3327-AT1 and R3327-HI tumors, which are reported to exhibit distinct size-dependent levels of hypoxia and response to hyperoxic gas breathing. Proton MRI R1 and R2* maps were obtained for tumors of anesthetized rats (isoflurane/air) at 4.7 T. Then, interleaved gradient echo T2*- and T1-weighted images were acquired during air breathing and a 10 min challenge with carbogen (95% O2 -5% CO2). Signals were stable during air breathing, and each type of tumor showed a distinct signal response to carbogen. T2* (BOLD) response preceded T1 (TOLD) responses, as expected. Smaller HI tumors (reported to be well oxygenated) showed the largest BOLD and TOLD responses. Larger AT1 tumors (reported to be hypoxic and resist modulation by gas breathing) showed the smallest response. There was a strong correlation between BOLD and TOLD signal responses, but ΔR2* and ΔR1 were only correlated for the HI tumors. The magnitude of BOLD and TOLD signal responses to carbogen breathing reflected expected hypoxic fractions and oxygen dynamics, suggesting potential value of this test as a prognostic biomarker of tumor hypoxia.

Quantitative MR Imaging of Brain Tissue and Brain Pathologies

Measurement of basic quantitative magnetic resonance (MR) parameters (e.g., relaxation times T1, T2*, T2 or respective rates R (1/T)) corrected for radiofrequency (RF) coil bias yields different conventional and new tissue contrasts as well as volumes for tissue segmentation. This approach also provides quantitative measures of microstructural and functional tissue changes. We herein demonstrate some prospects of quantitative MR imaging in neurological diagnostics and science.

O2 -sensitive MRI distinguishes brain tumor versus radiation necrosis in murine models

PURPOSE

The goal of this study was to quantify the relationship between the (1) H longitudinal relaxation rate constant, R1 , and oxygen (O2 ) concentration (relaxivity, r1 ) in tissue and to quantify O2 -driven changes in R1 (ΔR1 ) during a breathing gas challenge in normal brain, radiation-induced lesions, and tumor lesions.

METHODS

R1 data were collected in control-state mice (n = 4) during three different breathing gas (and thus tissue O2 ) conditions. In parallel experiments, pO2 was measured in the thalamus of control-state mice (n = 4) under the same breathing gas conditions using an O2 -sensitive microprobe. The relaxivity of tissue O2 was calculated using the R1 and pO2 data. R1 data were collected in control-state (n = 4) mice, a glioma model (n = 7), and a radiation necrosis model (n = 6) during two breathing gas (thus tissue O2 ) conditions. R1 and ΔR1 were calculated for each cohort.

RESULTS

O2 r1 in the brain was 9 × 10(-4)  ± 3 × 10(-4) mm Hg(-1) · s(-1) at 4.7T. R1 and ΔR1 measurements distinguished radiation necrosis from tumor (P< 0.03 and P< 0.01, respectively).

CONCLUSION

The relaxivity of O2 in the brain is determined. R1 and ΔR1 measurements differentiate tumor lesions from radiation necrosis lesions in the mouse models. These pathologies are difficult to distinguish by traditional imaging techniques; O2 -driven changes in R1 holds promise in this regard. Magn Reson Med 75:2442-2447, 2016. © 2015 Wiley Periodicals, Inc.

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.

Noninvasive tumor hypoxia measurement using magnetic resonance imaging in murine U87 glioma xenografts and in patients with glioblastoma

PURPOSE

There is a clinical need for noninvasive, nonionizing imaging biomarkers of tumor hypoxia and oxygenation. We evaluated the relationship of T1 -weighted oxygen-enhanced magnetic resonance imaging (OE-MRI) measurements to histopathology measurements of tumor hypoxia in a murine glioma xenograft and demonstrated technique translation in human glioblastoma multiforme.

METHODS

Preclinical evaluation was performed in a subcutaneous murine human glioma xenograft (U87MG). Animals underwent OE-MRI followed by dynamic contrast-enhanced MRI (DCE-MRI) and histological measurement including reduced pimonidazole adducts and CD31 staining. Area under the curve (AUC) was measured for the R1 curve for OE-MRI and the gadolinium concentration curve for DCE-MRI. Clinical evaluation in five patients used analogous imaging protocols and analyses.

RESULTS

Changes in AUC of OE-MRI (AUCOE ) signal were regionally heterogeneous across all U87MG tumors. Tumor regions with negative AUCOE typically had low DCE-MRI perfusion, had positive correlation with hypoxic area (P = 0.029), and had negative correlation with vessel density (P = 0.004). DCE-MRI measurements did not relate to either hypoxia or vessel density in U87MG tumors. Clinical data confirmed comparable signal changes in patients with glioblastoma.

CONCLUSION

These data support further investigation of T1 -weighted OE-MRI to identify regional tumor hypoxia. The quantification of AUCOE has translational potential as a clinical biomarker of hypoxia.

MRI under hyperbaric air and oxygen: effects on local magnetic field and relaxation times

PURPOSE

Hyperbaric oxygen therapy has shown efficacies in the treatment of a number of diseases. The goal of this study was to develop a rodent hyperbaric chamber for MRI studies and to investigate the effects of hyperbaric air and hyperbaric oxygen on local magnetic field (B0 ) and MRI relaxation parameters in the rat brain.

METHODS

A hyperbaric chamber, constructed to fit inside an animal MRI scanner, was pressurized with air to four atmospheres, while oxygen was delivered locally via nose cone. B0 , T2 , T2 *, and T1 maps in the rat brain were evaluated under normobaric air, hyperbaric air, and hyperbaric oxygen at 7T.

RESULTS

Under hyperbaric oxygen, images exhibited artifacts and temporal instability, attributable to fluctuating oxygen concentration from air and oxygen mixing near the imaging region. Physically shielding the imaging region from fluctuating oxygen concentration resolved the problems. With increasing oxygen at hyperbaric pressure, B0 was shifted downfield with increased inhomogeneity near the ear canals and nose. Brain T2 and T2 * were lengthened, and T1 was shortened.

CONCLUSION

This study establishes the means to perform MRI on rodents under hyperbaric conditions. Hyperbaric air and hyperbaric oxygen have significant effects on B0 and tissue relaxation parameters compared with normobaric air.

Evaluation and immunohistochemical qualification of carbogen-induced ΔR₂ as a noninvasive imaging biomarker of improved tumor oxygenation

PURPOSE

To evaluate and histologically qualify carbogen-induced ΔR2 as a noninvasive magnetic resonance imaging biomarker of improved tumor oxygenation using a double 2-nitroimidazole hypoxia marker approach.

METHODS AND MATERIALS

Multigradient echo images were acquired from mice bearing GH3 prolactinomas, preadministered with the hypoxia marker CCI-103F, to quantify tumor R2 during air breathing. With the mouse remaining positioned within the magnet bore, the gas supply was switched to carbogen (95% O2, 5% CO2), during which a second hypoxia marker, pimonidazole, was administered via an intraperitoneal line, and an additional set of identical multigradient echo images acquired to quantify any changes in tumor R2. Hypoxic fraction was quantified histologically using immunofluorescence detection of CCI-103F and pimonidazole adduct formation from the same whole tumor section. Carbogen-induced changes in tumor pO2 were further validated using the Oxylite fiberoptic probe.

RESULTS

Carbogen challenge significantly reduced mean tumor R2 from 116 ± 13 s(-1) to 97 ± 9 s(-1) (P<.05). This was associated with a significantly lower pimonidazole adduct area (2.3 ± 1%), compared with CCI-103F (6.3 ± 2%) (P<.05). A significant correlation was observed between ΔR2 and Δhypoxic fraction (r=0.55, P<.01). Mean tumor pO2 during carbogen breathing significantly increased from 6.3 ± 2.2 mm Hg to 36.0 ± 7.5 mm Hg (P<.01).

CONCLUSIONS

The combined use of intrinsic susceptibility magnetic resonance imaging with a double hypoxia marker approach corroborates carbogen-induced ΔR2 as a noninvasive imaging biomarker of increased tumor oxygenation.

Correlations of noninvasive BOLD and TOLD MRI with pO2 and relevance to tumor radiation response

PURPOSE

To examine the potential use of blood oxygenation level dependent (BOLD) and tissue oxygenation level dependent (TOLD) contrast MRI to assess tumor oxygenation and predict radiation response.

METHODS

BOLD and TOLD MRI were performed on Dunning R3327-AT1 rat prostate tumors during hyperoxic gas breathing challenge at 4.7 T. Animals were divided into two groups. In Group 1 (n = 9), subsequent (19) F MRI based on spin lattice relaxation of hexafluorobenzene reporter molecule provided quantitative oximetry for comparison. For Group 2 rats (n = 13) growth delay following a single dose of 30 Gy was compared with preirradiation BOLD and TOLD assessments.

RESULTS

Oxygen (100%O2 ) and carbogen (95%O2 /5%CO2 ) challenge elicited similar BOLD, TOLD and pO2 responses. Strong correlations were observed between BOLD or R2* response and quantitative (19) F pO2 measurements. TOLD response showed a general trend with weaker correlation. Irradiation caused a significant tumor growth delay and tumors with larger changes in TOLD and R1 values upon oxygen breathing exhibited significantly increased tumor growth delay.

CONCLUSION

These results provide further insight into the relationships between oxygen sensitive (BOLD/TOLD) MRI and tumor pO2 . Moreover, a larger increase in R1 response to hyperoxic gas challenge coincided with greater tumor growth delay following irradiation.

Simultaneous measurement of tissue oxygen level-dependent (TOLD) and blood oxygenation level-dependent (BOLD) effects in abdominal tissue oxygenation level studies

PURPOSE

To assess oxygenation in abdominal organs with magnetic resonance imaging (MRI), a novel approach is presented to simultaneously measure both T1 - and T2*-maps serially during a single dynamic MRI scan in response to an oxygen challenge.

MATERIALS AND METHODS

The proposed acquisition scheme consists of a multishot multiecho gradient echo planar imaging sequence (ms-GEPI) interleaved with a multishot inversion recovery echo planar imaging (ms-IR-EPI) sequence. Respiratory motion compensation was accomplished with standard belt triggering and by acquiring all image data at the same phase of expiration. This respiratory-triggered, free-breathing, interleaved tissue oxygenation level-dependent (TOLD) and blood oxygenation level-dependent (BOLD) acquisition technique was validated on phantoms and seven healthy volunteers in response to an oxygen challenge.

RESULTS

Measurements of relaxation times both in vitro and in vivo were in good agreement with those obtained using conventional pulse sequences and reported in the literature. The interleaved sequence was able to measure oxygen-induced relaxation time changes in human abdominal organs.

CONCLUSION

The free-breathing respiratory-triggered interleaved T1 and T2* sequence successfully provided relaxation time maps of abdominal organs in a dynamic scan without the need for image registration. The simultaneous monitoring of tissue and blood oxygenation improves time efficiency and should enhance studies comparing dynamic T1 and T2* data within the abdomen.