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

Validation of single reference variable flip angle (SR-VFA) dynamic T1 mapping with T2 * correction using a novel rotating phantom

PURPOSE

To validate single reference variable flip angle (SR-VFA) dynamic T1 mapping with and without T2 * correction against inversion recovery (IR) T1 measurements.

METHODS

A custom cylindrical phantom with three concentric compartments was filled with variably doped agar to produce a smooth spatial gradient of the T1 relaxation rate as a function of angle across each compartment. IR T1 , VFA T1 , and B1 + measurements were made on the phantom before rotation, and multi-echo stack-of-radial dynamic images were acquired during rotation via an MRI-compatible motor. B1 + -corrected SR-VFA and SR-VFA-T2 * T1 maps were computed from the sliding window reconstructed images and compared against rotationally registered IR and VFA T1 maps to determine the percentage error.

RESULTS

Both VFA and SR-VFA-T2 * T1 maps fell within 10% of IR T1 measurements for a low rotational speed, with a mean accuracy of 2.3% ± 2.6% and 2.8% ± 2.6%, respectively. Increasing rotational speed was found to decrease the accuracy due to increasing temporal smoothing over ranges where the T1 change had a nonconstant slope. SR-VFA T1 mapping was found to have similar accuracy as the SR-VFA-T2 * and VFA methods at low TEs (˜<2 ms), whereas accuracy degraded strongly with later TEs. T2 * correction of the SR-VFA T1 maps was found to consistently improve accuracy and precision, especially at later TEs.

CONCLUSION

SR-VFA-T2 * dynamic T1 mapping was found to be accurate against reference IR T1 measurements within 10% in an agar phantom. Further validation is needed in mixed fat-water phantoms and in vivo.

Correlation of MRI quantitative perfusion parameters with EGFR, VEGF and EGFR gene mutations in non-small cell cancer

To explore the relationship between quantitative perfusion histogram parameters of dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) with the expression of tumor tissue epidermal growth factor receptor (EGFR), vascular endothelial growth factor (VEGF) and EGFR gene mutations in non-small cell lung cancer (NSCLC). A total of 44 consecutive patients with known NSCLC were recruited from March 2018 to August 2021. Histogram parameters (mean, uniformity, skewness, energy, kurtosis, entropy, percentile) of each (Ktrans, Kep, Ve, Vp, Fp) were obtained by Omni Kinetics software. Immunohistochemistry staining was used in the detection of the expression of VEGF and EGFR protein, and the mutation of EGFR gene was detected by PCR. Corresponding statistical test was performed to compare the parameters and protein expression between squamous cell carcinoma (SCC) and adenocarcinoma (AC), as well as EGFR mutations and wild-type. Correlation analysis was used to evaluate the correlation between parameters with the expression of VEGF and EGFR protein. Fp (skewness, kurtosis, energy) were statistically significant between SCC and AC, and the area under the ROC curve were 0.733, 0.700 and 0.675, respectively. The expression of VEGF in AC was higher than in SCC. Fp (skewness, kurtosis, energy) were negatively correlated with VEGF (r =  - 0.527, - 0.428, - 0.342); Ktrans (Q50) was positively correlated with VEGF (r = 0.32); Kep (energy), Ktrans (skewness, kurtosis) were positively correlated with EGFR (r = 0.622, r = 0.375, 0.358), some histogram parameters of Kep, Ktrans (uniformity, entropy) and Ve (kurtosis) were negatively correlated with EGFR (r =  - 0.312 to - 0.644). Some perfusion histogram parameters were statistically significant between EGFR mutations and wild-type, they were higher in wild-type than mutated (P < 0.05). Quantitative perfusion histogram parameters of DCE-MRI have a certain value in the differential diagnosis of NSCLC, which have the potential to non-invasively evaluate the expression of cell signaling pathway-related protein.

Machine learning identifies multi-parametric functional PET/MR imaging cluster to predict radiation resistance in preclinical head and neck cancer models

PURPOSE

Tumor hypoxia and other microenvironmental factors are key determinants of treatment resistance. Hypoxia positron emission tomography (PET) and functional magnetic resonance imaging (MRI) are established prognostic imaging modalities to identify radiation resistance in head-and-neck cancer (HNC). The aim of this preclinical study was to develop a multi-parametric imaging parameter specifically for focal radiotherapy (RT) dose escalation using HNC xenografts of different radiation sensitivities.

METHODS

A total of eight human HNC xenograft models were implanted into 68 immunodeficient mice. Combined PET/MRI using dynamic [18F]-fluoromisonidazole (FMISO) hypoxia PET, diffusion-weighted (DW), and dynamic contrast-enhanced MRI was carried out before and after fractionated RT (10 × 2 Gy). Imaging data were analyzed on voxel-basis using principal component (PC) analysis for dynamic data and apparent diffusion coefficients (ADCs) for DW-MRI. A data- and hypothesis-driven machine learning model was trained to identify clusters of high-risk subvolumes (HRSs) from multi-dimensional (1-5D) pre-clinical imaging data before and after RT. The stratification potential of each 1D to 5D model with respect to radiation sensitivity was evaluated using Cohen's d-score and compared to classical features such as mean/peak/maximum standardized uptake values (SUVmean/peak/max) and tumor-to-muscle-ratios (TMRpeak/max) as well as minimum/valley/maximum/mean ADC.

RESULTS

Complete 5D imaging data were available for 42 animals. The final preclinical model for HRS identification at baseline yielding the highest stratification potential was defined in 3D imaging space based on ADC and two FMISO PCs ([Formula: see text]). In 1D imaging space, only clusters of ADC revealed significant stratification potential ([Formula: see text]). Among all classical features, only ADCvalley showed significant correlation to radiation resistance ([Formula: see text]). After 2 weeks of RT, FMISO_c1 showed significant correlation to radiation resistance ([Formula: see text]).

CONCLUSION

A quantitative imaging metric was described in a preclinical study indicating that radiation-resistant subvolumes in HNC may be detected by clusters of ADC and FMISO using combined PET/MRI which are potential targets for future functional image-guided RT dose-painting approaches and require clinical validation.

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.

Estimating quantitative physiological and morphological tissue parameters of murine tumor models using hyperspectral imaging and optical profilometry

Understanding tumors and their microenvironment are essential for successful and accurate disease diagnosis. Tissue physiology and morphology are altered in tumors compared to healthy tissues, and there is a need to monitor tumors and their surrounding tissues, including blood vessels, non-invasively. This preliminary study utilizes a multimodal optical imaging system combining hyperspectral imaging (HSI) and three-dimensional (3D) optical profilometry (OP) to capture hyperspectral images and surface shapes of subcutaneously grown murine tumor models. Hyperspectral images are corrected with 3D OP data and analyzed using the inverse-adding doubling (IAD) method to extract tissue properties such as melanin volume fraction and oxygenation. Blood vessels are segmented using the B-COSFIRE algorithm from oxygenation maps. From 3D OP data, tumor volumes are calculated and compared to manual measurements using a vernier caliper. Results show that tumors can be distinguished from healthy tissue based on most extracted tissue parameters ( p < 0.05 ). Furthermore, blood oxygenation is 50% higher within the blood vessels than in the surrounding tissue, and tumor volumes calculated using 3D OP agree within 26% with manual measurements using a vernier caliper. Results suggest that combining HSI and OP could provide relevant quantitative information about tumors and improve the disease diagnosis.

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.

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.

High Doses of Photons and Carbon Ions Comparably Increase Vascular Permeability in R3327-HI Prostate Tumors: A Dynamic Contrast-Enhanced MRI Study

Carbon- (12C-) ion radiotherapy exhibits enhanced biological effectiveness compared to photon radiotherapy, however, the contribution of its interaction with the vasculature remains debatable. The effect of high-dose 12C-ion and photon irradiation on vascular permeability in moderately differentiated rat prostate tumors was compared using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Syngeneic R3327-HI rat prostate tumors were irradiated with a single dose of either 18 or 37 Gy 12C ions, or 37 or 75 Gy 6-MV photons (sub-curative and curative dose levels, respectively). DCE-MRI was performed one day prior to and 3, 7, 14 and 21 days postirradiation. Voxel-based tumor concentration-time curves were clustered based on their curve shape and treatment response was assessed as the longitudinal changes in the relative abundance per cluster. Radiation-induced vascular damage and increased permeability occurred at day 7 postirradiation for all treatment groups except for the 75 Gy photon-irradiated group, where the onset of vascular damage was delayed until day 14. No differences between irradiation modalities were found. Therefore, early vascular damage cannot explain the higher effectiveness of 12C ions relative to photons in terms of local tumor control for this moderately differentiated prostate tumor and the applied single high doses.

Exploring MRI based radiomics analysis of intratumoral spatial heterogeneity in locally advanced nasopharyngeal carcinoma treated with intensity modulated radiotherapy

Background

We hypothesized that spatial heterogeneity exists between recurrent and non-recurrent regions within a tumor. The aim of this study was to determine if there is a difference between radiomics features derived from recurrent versus non recurrent regions within the tumor based on pre-treatment MRI.

Methods

A total of 14 T4NxM0 NPC patients with histologically proven “in field” recurrence in the post nasal space following curative intent IMRT were included in this study. Pretreatment MRI were co-registered with MRI at the time of recurrence for the delineation of gross tumor volume at diagnosis(GTV) and at recurrence(GTVr). A total of 7 histogram features and 40 texture features were computed from the recurrent(GTVr) and non-recurrent region(GTV-GTVr). Paired t-tests and Wilcoxon signed-rank tests were carried out on the 47 quantified radiomics features.

Results

A total of 7 features were significantly different between recurrent and non-recurrent regions. Other than the variance from intensity-based histogram, the remaining six significant features were either from the gray-level size zone matrix (GLSZM) or the neighbourhood gray-tone difference matrix (NGTDM).

Conclusions

The radiomic features extracted from pre-treatment MRI can potentially reflect the difference between recurrent and non-recurrent regions within a tumor and has a potential role in pre-treatment identification of intra-tumoral radio-resistance for selective dose escalation.

Multiparametric Analysis of Longitudinal Quantitative MRI data to Identify Distinct Tumor Habitats in Preclinical Models of Breast Cancer

This study identifies physiological tumor habitats from quantitative magnetic resonance imaging (MRI) data and evaluates their alterations in response to therapy. Two models of breast cancer (BT-474 and MDA-MB-231) were imaged longitudinally with diffusion-weighted MRI and dynamic contrast-enhanced MRI to quantify tumor cellularity and vascularity, respectively, during treatment with trastuzumab or albumin-bound paclitaxel. Tumors were stained for anti-CD31, anti-Ki-67, and H&E. Imaging and histology data were clustered to identify tumor habitats and percent tumor volume (MRI) or area (histology) of each habitat was quantified. Histological habitats were correlated with MRI habitats. Clustering of both the MRI and histology data yielded three clusters: high-vascularity high-cellularity (HV-HC), low-vascularity high-cellularity (LV-HC), and low-vascularity low-cellularity (LV-LC). At day 4, BT-474 tumors treated with trastuzumab showed a decrease in LV-HC (p = 0.03) and increase in HV-HC (p = 0.03) percent tumor volume compared to control. MDA-MB-231 tumors treated with low-dose albumin-bound paclitaxel showed a longitudinal decrease in LV-HC percent tumor volume at day 3 (p = 0.01). Positive correlations were found between histological and imaging-derived habitats: HV-HC (BT-474: p = 0.03), LV-HC (MDA-MB-231: p = 0.04), LV-LC (BT-474: p = 0.04; MDA-MB-231: p < 0.01). Physiologically distinct tumor habitats associated with therapeutic response were identified with MRI and histology data in preclinical models of breast cancer.

An Automated Segmentation Pipeline for Intratumoural Regions in Animal Xenografts Using Machine Learning and Saturation Transfer MRI

Saturation transfer MRI can be useful in the characterization of different tumour types. It is sensitive to tumour metabolism, microstructure, and microenvironment. This study aimed to use saturation transfer to differentiate between intratumoural regions, demarcate tumour boundaries, and reduce data acquisition times by identifying the imaging scheme with the most impact on segmentation accuracy. Saturation transfer-weighted images were acquired over a wide range of saturation amplitudes and frequency offsets along with T₁ and T₂ maps for 34 tumour xenografts in mice. Independent component analysis and Gaussian mixture modelling were used to segment the images and identify intratumoural regions. Comparison between the segmented regions and histopathology indicated five distinct clusters: three corresponding to intratumoural regions (active tumour, necrosis/apoptosis, and blood/edema) and two extratumoural (muscle and a mix of muscle and connective tissue). The fraction of tumour voxels segmented as necrosis/apoptosis quantitatively matched those calculated from TUNEL histopathological assays. An optimal protocol was identified providing reasonable qualitative agreement between MRI and histopathology and consisting of T₁ and T₂ maps and 22 magnetization transfer (MT)-weighted images. A three-image subset was identified that resulted in a greater than 90% match in positive and negative predictive value of tumour voxels compared to those found using the entire 24-image dataset. The proposed algorithm can potentially be used to develop a robust intratumoural segmentation method.

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.

Multiparametric MRI and Coregistered Histology Identify Tumor Habitats in Breast Cancer Mouse Models

It is well-recognized that solid tumors are genomically, anatomically, and physiologically heterogeneous. In general, more heterogeneous tumors have poorer outcomes, likely due to the increased probability of harboring therapy-resistant cells and regions. It is hypothesized that the genomic and physiologic heterogeneity are related, because physiologically distinct regions will exert variable selection pressures leading to the outgrowth of clones with variable genomic/proteomic profiles. To investigate this, methods must be in place to interrogate and define, at the microscopic scale, the cytotypes that exist within physiologically distinct subregions ("habitats") that are present at mesoscopic scales. MRI provides a noninvasive approach to interrogate physiologically distinct local environments, due to the biophysical principles that govern MRI signal generation. Here, we interrogate different physiologic parameters, such as perfusion, cell density, and edema, using multiparametric MRI (mpMRI). Signals from six different acquisition schema were combined voxel-by-voxel into four clusters identified using a Gaussian mixture model. These were compared with histologic and IHC characterizations of sections that were coregistered using MRI-guided 3D printed tumor molds. Specifically, we identified a specific set of MRI parameters to classify viable-normoxic, viable-hypoxic, nonviable-hypoxic, and nonviable-normoxic tissue types within orthotopic 4T1 and MDA-MB-231 breast tumors. This is the first coregistered study to show that mpMRI can be used to define physiologically distinct tumor habitats within breast tumor models. SIGNIFICANCE: This study demonstrates that noninvasive imaging metrics can be used to distinguish subregions within heterogeneous tumors with histopathologic correlation.

Data-driven identification of tumor subregions based on intravoxel incoherent motion reveals association with proliferative activity

PURPOSE

Intravoxel incoherent motion (IVIM) analysis gives information on tissue diffusion and perfusion and may thus have a potential for e.g. tumor tissue characterization. This work aims to study if clustering based on IVIM parameter maps can identify tumor subregions, and to assess the relevance of obtained subregions by histological analysis.

METHODS

Fourteen mice with human neuroendocrine tumors were examined with diffusion-weighted imaging to obtain IVIM parameter maps. Gaussian mixture models with IVIM maps from all tumors as input were used to partition voxels into k clusters, where k = 2 was chosen for further analysis based on goodness of fit. Clustering was performed with and without the perfusion-related IVIM parameter D * , and with and without including spatial information. The validity of the clustering was assessed by comparison with corresponding histologically stained tumor sections. A Ki-67-based index quantifying the degree of tumor proliferation was considered appropriate for the comparison based on the obtained cluster characteristics.

RESULTS

The clustering resulted in one class with low diffusion and high perfusion and another with slightly higher diffusion and low perfusion. Strong agreement was found between tumor subregions identified by clustering and subregions identified by histological analysis, both regarding size and spatial agreement. Neither D * nor spatial information had substantial effects on the clustering results.

CONCLUSIONS

The results of this study show that IVIM parameter maps can be used to identify tumor subregions using a data-driven framework based on Gaussian mixture models. In the studied tumor model, the obtained subregions showed agreement with proliferative activity.