Longitudinal Correlations Between Intravoxel Incoherent Motion (IVIM) and Dynamic Contrast-Enhanced (DCE) MRI During Radiotherapy in Prostate Cancer Patients

Multiple diffusion models for predicting pathologic response of esophageal squamous cell carcinoma to neoadjuvant chemotherapy

BACKGROUND

To assess the feasibility and diagnostic performance of the fractional order calculus (FROC), continuous-time random-walk (CTRW), diffusion kurtosis imaging (DKI), intravoxel incoherent motion (IVIM), mono-exponential (MEM) and stretched exponential models (SEM) for predicting response to neoadjuvant chemotherapy (NACT) in patients with esophageal squamous cell carcinoma (ESCC).

MATERIALS AND METHODS

This study prospectively included consecutive ESCC patients with baseline and follow up MR imaging and pathologically confirmed cT1-4aN + M0 or T3-4aN0M0 and underwent radical resection after neoadjuvant chemotherapy (NACT) between July 2019 and January 2023. Patients were divided into pCR (TRG 0) and non-pCR (TRG1 + 2 + 3) groups according to tumor regression grading (TRG). The Pre-, Post- and Delta-treatment models were built. 18 predictive models were generated according to different feature categories, based on six models by five-fold cross-validation. Areas under the curve (AUCs) of the models were compared by using DeLong method.

RESULTS

Overall, 90 patients (71 men, 19 women; mean age, 64 years ± 6 [SD]) received NACT and underwent baseline and Post-NACT esophageal MRI, with 29 patients in the pCR group and 61 patients in the non-pCR group. Among 18 predictive models, The Pre-, Post-, and Delta-CTRW model showed good predictive efficacy (AUC = 0.722, 0.833 and 0.790). Additionally, the Post-FROC model (AUC = 0.907) also exhibited good diagnostic performance.

CONCLUSIONS

Our study indicates that the CTRW model, along with the Post-FROC model, holds significant promise for the future of NACT efficacy prediction in ESCC patients.

Longitudinal Image Data for Outcome Modeling

In oncology, medical imaging is crucial for diagnosis, treatment planning and therapy execution. Treatment responses can be complex and varied and are known to involve factors of treatment, patient characteristics and tumor microenvironment. Longitudinal image analysis is able to track temporal changes, aiding in disease monitoring, treatment evaluation, and outcome prediction. This allows for the enhancement of personalized medicine. However, analyzing longitudinal 2D and 3D images presents unique challenges, including image registration, reliable segmentation, dealing with variable imaging intervals, and sparse data. This review presents an overview of techniques and methodologies in longitudinal image analysis, with a primary focus on outcome modeling in radiation oncology.

Clinical Workflow of Cone Beam Computer Tomography-Based Daily Online Adaptive Radiotherapy with Offline Magnetic Resonance Guidance: The Modular Adaptive Radiotherapy System (MARS)

PURPOSE

The Ethos (Varian Medical Systems) radiotherapy device combines semi-automated anatomy detection and plan generation for cone beam computer tomography (CBCT)-based daily online adaptive radiotherapy (oART). However, CBCT offers less soft tissue contrast than magnetic resonance imaging (MRI). This work aims to present the clinical workflow of CBCT-based oART with shuttle-based offline MR guidance.

METHODS

From February to November 2023, 31 patients underwent radiotherapy on the Ethos (Varian, Palo Alto, CA, USA) system with machine learning (ML)-supported daily oART. Moreover, patients received weekly MRI in treatment position, which was utilized for daily plan adaptation, via a shuttle-based system. Initial and adapted treatment plans were generated using the Ethos treatment planning system. Patient clinical data, fractional session times (MRI + shuttle transport + positioning, adaptation, QA, RT delivery) and plan selection were assessed for all fractions in all patients.

RESULTS

In total, 737 oART fractions were applied and 118 MRIs for offline MR guidance were acquired. Primary sites of tumors were prostate (n = 16), lung (n = 7), cervix (n = 5), bladder (n = 1) and endometrium (n = 2). The treatment was completed in all patients. The median MRI acquisition time including shuttle transport and positioning to initiation of the Ethos adaptive session was 53.6 min (IQR 46.5-63.4). The median total treatment time without MRI was 30.7 min (IQR 24.7-39.2). Separately, median adaptation, plan QA and RT times were 24.3 min (IQR 18.6-32.2), 0.4 min (IQR 0.3-1,0) and 5.3 min (IQR 4.5-6.7), respectively. The adapted plan was chosen over the scheduled plan in 97.7% of cases.

CONCLUSION

This study describes the first workflow to date of a CBCT-based oART combined with a shuttle-based offline approach for MR guidance. The oART duration times reported resemble the range shown by previous publications for first clinical experiences with the Ethos system.

Towards tumour hypoxia imaging: Incorporating relative oxygen extraction fraction mapping of prostate with multi-parametric quantitative MRI on a 1.5T MR-linac

Hypoxia plays a central role in tumour radioresistance. Reliable tumour hypoxia imaging would allow the monitoring of tumour response and a more personalized adaptation of radiotherapy planning. Here, we showed a proof of concept of the feasibility and repeatability of relative oxygen extraction fraction (rOEF) mapping of prostate using multi-parametric quantitative MRI (qMRI) achieved for the first time on a 1.5T MR-linac. T2, T2* relaxation times maps, and intra-voxel incoherent motion (IVIM) parametric maps mapping were computed on a 29 years old healthy volunteer. R2' and rOEF maps were calculated based on a multi-parametric model. Long-term repeatability and repeatability coefficient (RC) were determined for each parameter according to QIBA recommendations. Mean values for the entire healthy prostate were 0.99 ± 0.14 × 10-3  mm/s2 , 81 ± 2.1 × 10-3  mm/s2 , 21.6 ± 3.6%, 92.7 ± 19.7 ms and 62.4 ± 17.3 ms for Dslow , Dfast , f, T2 and T2*, respectively. R2' and rOEF in the prostate were 6.1 ± 3.4 s-1 and 18.2 ± 10.1% respectively. The RC of rOEF was 4.43%. Long-term repeatability of quantitative parameters based on a test-retest ranged from 2 to 18%. qMRI parameters are measurable and repeatable on 1.5T MR LINAC. From T2, T2* and IVIM parameters maps, we were able to obtain a rOEF mapping of the prostate. These results are the first step to a non-invasive imaging of tumour hypoxia during radiotherapy leading to a biological image-guided adaptive radiotherapy.

Magnetic Resonance Perfusion Imaging of Prostate

Magnetic resonance (MR) perfusion imaging, both with and without exogenous contrast agents, has the potential to assess tissue perfusion and vascularity in prostate cancer. Dynamic contrast-enhanced (DCE) MRI is an important element of the clinical non-invasive multiparametric MRI, which can be used to differentiate benign from malignant lesions, to stage tumors, and to monitor response to therapy. The arterial spin labeled (ASL) and intravoxel incoherent motion (IVIM) diffusion-weighted MRI have the advantage of quantitative perfusion measurements without the concerns of gadolinium-based contrast agent safety and retention issues. The adoption of these non-contrast techniques in clinical practice needs more research and clinical evaluation.

Reproducibility of diffusion-weighted magnetic resonance imaging in head and neck cancer assessed on a 1.5 T MR-Linac and comparison to parallel measurements on a 3 T diagnostic scanner

BACKGROUND AND PURPOSE

Before quantitative imaging biomarkers (QIBs) acquired with magnetic resonance imaging (MRI) can be used for interventional trials in radiotherapy (RT), technical validation of these QIBs is necessary. The aim of this study was to assess the reproducibility of apparent diffusion coefficient (ADC) values, derived from diffusion-weighted (DW) MRI, in head and neck cancer using a 1.5 T MR-Linac (MRL) by comparison to a 3 T diagnostic scanner (DS).

MATERIAL AND METHODS

DW-MRIs were acquired on MRL and DS for 15 head and neck cancer patients before RT and in week 2 and rigidly registered to the planning computed tomography. Mean ADC values were calculated for submandibular (SG) and parotid (PG) glands as well as target volumes (TV, gross tumor volume and lymph nodes), which were delineated based on computed tomography. Mean absolute ADC differences as well as within-subject coefficient of variation (wCV) and intraclass correlation coefficients (ICCs) were calculated for all volumes of interest.

RESULTS

A total of 23 datasets were analyzed. Mean ADC difference (DS-MRL) for SG, PG and TV resulted in 142, 254 and 93·10-6 mm2/s. wCVs/ICCs, comparing MRL and DS, were determined as 13.7 %/0.26, 24.4 %/0.23 and 16.1 %/0.73 for SG, PG and TV, respectively.

CONCLUSION

ADC values, measured on the 1.5 T MRL, showed reasonable reproducibility with an ADC underestimation in contrast to the DS. This ADC shift must be validated in further experiments and considered for future translation of QIB candidates from DS to MRL for response adaptive RT.

MRI-Guided Adaptive Radiation Therapy

Magnetic resonance imaging-guided radiation therapy (MRIgRT) has improved soft tissue contrast over computed tomography (CT) based image-guided RT. Superior visualization of the target and surrounding radiosensitive structures has the potential to improve oncological outcomes partly due to safer dose-escalation and adaptive planning. In this review, we highlight the workflow of adaptive MRIgRT planning, which includes simulation imaging, daily MRI, identifying isocenter shifts, contouring, plan optimization, quality control, and delivery. Increased utilization of MRIgRT will depend on addressing technical limitations of this technology, while addressing treatment efficacy, cost-effectiveness, and workflow training.

Quantitative MRI on MR-Linacs: Towards Biological Image-Guided Adaptive Radiotherapy

Recognizing the potential of quantitative imaging biomarkers (QIBs) in radiotherapy, many studies have investigated the prognostic value of quantitative MRI (qMRI). With the introduction of MRI-guided radiotherapy systems, the practical challenges of repeated imaging have been substantially reduced. Since patients are treated inside an MRI scanner, acquisition of qMRI can be done during each fraction with limited or no prolongation of the fraction duration. In this review paper, we identify the steps that need been taken to move from MR as an imaging technique to a useful biomarker for MRI-guided radiotherapy (MRgRT).

Longitudinal monitoring of Apparent Diffusion Coefficient (ADC) in patients with prostate cancer undergoing MR-guided radiotherapy on an MR-Linac at 1.5 T: a prospective feasibility study

BACKGROUND

Hybrid MRI linear accelerators (MR-Linac) might enable individualized online adaptation of radiotherapy using quantitative MRI sequences as diffusion-weighted imaging (DWI). The purpose of this study was to investigate the dynamics of lesion apparent diffusion coefficient (ADC) in patients with prostate cancer undergoing MR-guided radiation therapy (MRgRT) on a 1.5T MR-Linac. The ADC values at a diagnostic 3T MRI scanner were used as the reference standard.

PATIENTS AND AND METHODS

In this prospective single-center study, patients with biopsy-confirmed prostate cancer who underwent both an MRI exam at a 3T scanner (MRI_3T) and an exam at a 1.5T MR-Linac (MRL) at baseline and during radiotherapy were included. Lesion ADC values were measured by a radiologist and a radiation oncologist on the slice with the largest lesion. ADC values were compared before vs. during radiotherapy (during the second week) on both systems via paired t-tests. Furthermore, Pearson correlation coefficient and inter-reader agreement were computed.

RESULTS

A total of nine male patients aged 67 ± 6 years [range 60 - 67 years] were included. In seven patients, the cancerous lesion was in the peripheral zone, and in two patients the lesion was in the transition zone. Inter-reader reliability regarding lesion ADC measurement was excellent with an intraclass correlation coefficient of (ICC) > 0.90 both at baseline and during radiotherapy. Thus, the results of the first reader will be reported. In both systems, there was a statistically significant elevation of lesion ADC during radiotherapy (mean MRL-ADC at baseline was 0.97 ± 0.18 × 10^-3 mm²/s vs. mean MRL-ADC during radiotherapy 1.38 ± 0.3 × 10^-3 mm²/s, yielding a mean lesion ADC elevation of 0.41 ± 0.20 × 10^-3 mm²/s, p < 0.001). Mean MRI_3T-ADC at baseline was 0.78 ± 0.165 × 10^-3 mm²/s vs. mean MRI_3T-ADC during radiotherapy 0.99 ± 0.175 × 10^-3 mm²/s, yielding a mean lesion ADC elevation of 0.21 ± 0.96 × 10^-3 mm²/s p < 0.001). The absolute ADC values from MRL were consistently significantly higher than those from MRI_3T at baseline and during radiotherapy (p < = 0.001). However, there was a strong positive correlation between MRL-ADC and MRI_3T-ADC at baseline (r = 0.798, p = 0.01) and during radiotherapy (r = 0.863, p = 0.003).

CONCLUSIONS

Lesion ADC as measured on MRL increased significantly during radiotherapy and ADC measurements of lesions on both systems showed similar dynamics. This indicates that lesion ADC as measured on the MRL may be used as a biomarker for evaluation of treatment response. In contrast, absolute ADC values as calculated by the algorithm of the manufacturer of the MRL showed systematic deviations from values obtained on a diagnostic 3T MRI system. These preliminary findings are promising but need large-scale validation. Once validated, lesion ADC on MRL might be used for real-time assessment of tumor response in patients with prostate cancer undergoing MR-guided radiation therapy.

Three-dimensional amide proton transfer-weighted and intravoxel incoherent motion imaging for predicting bone metastasis in patients with prostate cancer: A pilot study

PURPOSE

To explore the value of 3-dimensional amide proton transfer-weighted (APTw) and intravoxel incoherent motion (IVIM) imaging in predicting bone metastasis (BM) of prostate cancer (PCa) in addition to routine diffusion-weighted imaging (DWI).

METHODS

The clinical and imaging data of 39 PCa patients who were pathologically confirmed in our hospital from March 2019 to February 2022 were retrospectively analyzed, and they were divided into BM-negative (27 patients) and BM-positive (12 patients) groups. MR examination included APTw, DWI and IVIM imaging. The IVIM data was fitted by single-exponential IVIM model (IVIM_mono) and double-exponential IVIM model (IVIM_bi), respectively. The APTw, ADC, IVIM_mono (D_mono, D*_mono, and f_mono), and IVIM_bi (D_bi, D*_bi, and f_bi) parameters were independently measured by two radiologists. The synthetic minority oversampling technique (SMOTE) was conducted to balance the minority group. Mann-Whitney U test or Student's t-test was used to compare above values between the BM-negative and BM-positive groups. The diagnostic performance was evaluated with receiver operating characteristic (ROC) analysis of each parameter and their combination. The Delong test was used for ROC curve comparison.The relationship between APTw and IVIM was explored through Spearman's rank correlation analysis.

RESULTS

The APTw and D*_mono values were higher, and the ADC, f_mono, and f_bi values were lower in the BM-positive group than in the BM-negative group (all P < 0.05). Among the individual parameters, the AUC of f_mono was the highest (AUC = 0.865), and AUC (f_mono) was significantly higher than AUC (f_bi), AUC (D*_mono), and AUC (ADC) (all P < 0.05). The AUC (IVIM_mono) was higher than the AUC (IVIM_bi) (P = 0.0068). The combination of APTw and IVIM_mono further improved diagnostic capability, and the AUC of APTw+IVIM_mono was significantly higher than those of APTw and DWI (all P < 0.05). No correlation was found between IVIM-derived parameters and APTw value.

CONCLUSION

Both 3D APTw and IVIM imaging could predict BM of PCa. IVIM showed better performance than APTw and DWI, and the single-exponential IVIM model was superior to the double-exponential IVIM model. The combination of APTw and IVIM could further improve diagnostic performance.