Development and multicenter validation of a multiparametric imaging model to predict treatment response in rectal cancer.

OBJECTIVES: To develop and validate a multiparametric model to predict neoadjuvant treatment response in rectal cancer at baseline using a heterogeneous multicenter MRI dataset. METHODS: Baseline staging MRIs (T2W (T2-weighted)-MRI, diffusion-weighted imaging (DWI) / apparent diffusion coefficient (ADC)) of 509 patients (9 centres) treated with neoadjuvant chemoradiotherapy (CRT) were collected. Response was defined as (1) complete versus incomplete response, or (2) good (Mandard tumor regression grade (TRG) 1-2) versus poor response (TRG3-5). Prediction models were developed using combinations of the following variable groups: (1) Non-imaging: age/sex/tumor-location/tumor-morphology/CRT-surgery interval (2) Basic staging: cT-stage/cN-stage/mesorectal fascia involvement, derived from (2a) original staging reports, or (2b) expert re-evaluation (3) Advanced staging: variables from 2b combined with cTN-substaging/invasion depth/extramural vascular invasion/tumor length (4) Quantitative imaging: tumour volume + first-order histogram features (from T2W-MRI and DWI/ADC) Models were developed with data from 6 centers (n = 412) using logistic regression with the Least Absolute Shrinkage and Selector Operator (LASSO) feature selection, internally validated using repeated (n = 100) random hold-out validation, and externally validated using data from 3 centers (n = 97). RESULTS: After external validation, the best model (including non-imaging and advanced staging variables) achieved an area under the curve of 0.60 (95%CI=0.48-0.72) to predict complete response and 0.65 (95%CI=0.53-0.76) to predict a good response. Quantitative variables did not improve model performance. Basic staging variables consistently achieved lower performance compared to advanced staging variables. CONCLUSIONS: Overall model performance was moderate. Best results were obtained using advanced staging variables, highlighting the importance of good-quality staging according to current guidelines. Quantitative imaging features had no added value (in this heterogeneous dataset). CLINICAL RELEVANCE STATEMENT: Predicting tumour response at baseline could aid in tailoring neoadjuvant therapies for rectal cancer. This study shows that image-based prediction models are promising, though are negatively affected by variations in staging quality and MRI acquisition, urging the need for harmonization. KEY POINTS: This multicenter study combining clinical information and features derived from MRI rendered disappointing performance to predict response to neoadjuvant treatment in rectal cancer. Best results were obtained with the combination of clinical baseline information and state-of-the-art image-based staging variables, highlighting the importance of good quality staging according to current guidelines and staging templates. No added value was found for quantitative imaging features in this multicenter retrospective study. This is likely related to acquisition variations, which is a major problem for feature reproducibility and thus model generalizability.

Sources of variation in multicenter rectal MRI data and their effect on radiomics feature reproducibility.

OBJECTIVES: To investigate sources of variation in a multicenter rectal cancer MRI dataset focusing on hardware and image acquisition, segmentation methodology, and radiomics feature extraction software. METHODS: T2W and DWI/ADC MRIs from 649 rectal cancer patients were retrospectively acquired in 9 centers. Fifty-two imaging features (14 first-order/6 shape/32 higher-order) were extracted from each scan using whole-volume (expert/non-expert) and single-slice segmentations using two different software packages (PyRadiomics/CapTk). Influence of hardware, acquisition, and patient-intrinsic factors (age/gender/cTN-stage) on ADC was assessed using linear regression. Feature reproducibility was assessed between segmentation methods and software packages using the intraclass correlation coefficient. RESULTS: Image features differed significantly (p < 0.001) between centers with more substantial variations in ADC compared to T2W-MRI. In total, 64.3% of the variation in mean ADC was explained by differences in hardware and acquisition, compared to 0.4% by patient-intrinsic factors. Feature reproducibility between expert and non-expert segmentations was good to excellent (median ICC 0.89-0.90). Reproducibility for single-slice versus whole-volume segmentations was substantially poorer (median ICC 0.40-0.58). Between software packages, reproducibility was good to excellent (median ICC 0.99) for most features (first-order/shape/GLCM/GLRLM) but poor for higher-order (GLSZM/NGTDM) features (median ICC 0.00-0.41). CONCLUSIONS: Significant variations are present in multicenter MRI data, particularly related to differences in hardware and acquisition, which will likely negatively influence subsequent analysis if not corrected for. Segmentation variations had a minor impact when using whole volume segmentations. Between software packages, higher-order features were less reproducible and caution is warranted when implementing these in prediction models. KEY POINTS: • Features derived from T2W-MRI and in particular ADC differ significantly between centers when performing multicenter data analysis. • Variations in ADC are mainly (> 60%) caused by hardware and image acquisition differences and less so (< 1%) by patient- or tumor-intrinsic variations. • Features derived using different image segmentations (expert/non-expert) were reproducible, provided that whole-volume segmentations were used. When using different feature extraction software packages with similar settings, higher-order features were less reproducible.

Studying local tumour heterogeneity on MRI and FDG-PET/CT to predict response to neoadjuvant chemoradiotherapy in rectal cancer.

OBJECTIVE: To investigate whether quantifying local tumour heterogeneity has added benefit compared to global tumour features to predict response to chemoradiotherapy using pre-treatment multiparametric PET and MRI data. METHODS: Sixty-one locally advanced rectal cancer patients treated with chemoradiotherapy and staged at baseline with MRI and FDG-PET/CT were retrospectively analyzed. Whole-tumour volumes were segmented on the MRI and PET/CT scans from which global tumour features (T2Wvolume/T2Wentropy/ADCmean/SUVmean/TLG/CTmean-HU) and local texture features (histogram features derived from local entropy/mean/standard deviation maps) were calculated. These respective feature sets were combined with clinical baseline parameters (e.g. age/gender/TN-stage) to build multivariable prediction models to predict a good (Mandard TRG1-2) versus poor (Mandard TRG3-5) response to chemoradiotherapy. Leave-one-out cross-validation (LOOCV) with bootstrapping was performed to estimate performance in an 'independent' dataset. RESULTS: When using only imaging features, local texture features showed an AUC = 0.81 versus AUC = 0.74 for global tumour features. After internal cross-validation (LOOCV), AUC to predict a good response was the highest for the combination of clinical baseline variables + global tumour features (AUC = 0.83), compared to AUC = 0.79 for baseline + local texture and AUC = 0.76 for all combined (baseline + global + local texture). CONCLUSION: In imaging-based prediction models, local texture analysis has potential added value compared to global tumour features to predict response. However, when combined with clinical baseline parameters such as cTN-stage, the added value of local texture analysis appears to be limited. The overall performance to predict response when combining baseline variables with quantitative imaging parameters is promising and warrants further research. KEY POINTS: • Quantification of local tumour texture on pre-therapy FDG-PET/CT and MRI has potential added value compared to global tumour features to predict response to chemoradiotherapy in rectal cancer. • However, when combined with clinical baseline parameters such as cTN-stage, the added value of local texture over global tumour features is limited. • Predictive performance of our optimal model-combining clinical baseline variables with global quantitative tumour features-was encouraging (AUC 0.83), warranting further research in this direction on a larger scale.

Value of combined multiparametric MRI and FDG-PET/CT to identify well-responding rectal cancer patients before the start of neoadjuvant chemoradiation.

OBJECTIVES: To explore the value of multiparametric MRI combined with FDG-PET/CT to identify well-responding rectal cancer patients before the start of neoadjuvant chemoradiation. METHODS: Sixty-one locally advanced rectal cancer patients who underwent a baseline FDG-PET/CT and MRI (T2W + DWI) and received long-course neoadjuvant chemoradiotherapy were retrospectively analysed. Tumours were delineated on MRI and PET/CT from which the following quantitative parameters were calculated: T2W volume and entropy, ADC mean and entropy, CT density (mean-HU), SUV maximum and mean, metabolic tumour volume (MTV42%) and total lesion glycolysis (TLG). These features, together with sex, age, mrTN-stage ("baseline parameters") and the CRT-surgery interval were analysed using multivariable stepwise logistic regression. Outcome was a good (TRG 1-2) versus poor histopathological response. Performance (AUC) to predict response was compared for different combinations of baseline ± quantitative imaging parameters and performance in an 'independent' dataset was estimated using bootstrapped leave-one-out cross-validation (LOOCV). RESULTS: The optimal multivariable prediction model consisted of a combination of baseline + quantitative imaging parameters and included mrT-stage (OR 0.004, p < 0.001), T2W-signal entropy (OR 7.81, p = 0.0079) and T2W volume (OR 1.028, p = 0.0389) as the selected predictors. AUC in the study dataset was 0.88 and 0.83 after LOOCV. No PET/CT features were selected as predictors. CONCLUSIONS: A multivariable model incorporating mrT-stage and quantitative parameters from baseline MRI can aid in identifying well-responding patients before the start of treatment. Addition of FDG-PET/CT is not beneficial. KEY POINTS: • A multivariable model incorporating the mrT-stage and quantitative features derived from baseline MRI can aid in identifying well-responding patients before the start of neoadjuvant chemoradiotherapy. • mrT-stage was the strongest predictor in the model and was complemented by the tumour volume and signal entropy calculated from T2W-MRI. • Adding quantitative features derived from pre-treatment PET/CT or DWI did not contribute to the model's predictive performance.

A multicentric randomized controlled phase III trial of adaptive and 18F-FDG-PET-guided dose-redistribution in locally advanced head and neck squamous cell carcinoma (ARTFORCE).

BACKGROUND AND PURPOSE: This multicenter randomized phase III trial evaluated whether locoregional control of patients with LAHNSCC could be improved by fluorodeoxyglucose-positron emission tomography (FDG-PET)-guided dose-escalation while minimizing the risk of increasing toxicity using a dose-redistribution and scheduled adaptation strategy. MATERIALS AND METHODS: Patients with T3-4-N0-3-M0 LAHNSCC were randomly assigned (1:1) to either receive a dose distribution ranging from 64-84 Gy/35 fractions with adaptation at the 10thfraction (rRT) or conventional 70 Gy/35 fractions (cRT). Both arms received concurrent three-cycle 100 mg/m2cisplatin. Primary endpoints were 2-year locoregional control (LRC) and toxicity. Primary analysis was based on the intention-to-treat principle. RESULTS: Due to slow accrual, the study was prematurely closed (at 84 %) after randomizing 221 eligible patients between 2012 and 2019 to receive rRT (N = 109) or cRT (N = 112). The 2-year LRC estimate difference of 81 % (95 %CI 74-89 %) vs. 74 % (66-83 %) in the rRT and cRT arm, respectively, was not found statistically significant (HR 0.75, 95 %CI 0.43-1.31,P=.31). Toxicity prevalence and incidence rates were similar between trial arms, with exception for a significant increased grade ≥ 3 pharyngolaryngeal stenoses incidence rate in the rRT arm (0 versus 4 %,P=.05). In post-hoc subgroup analyses, rRT improved LRC for patients with N0-1 disease (HR 0.21, 95 %CI 0.05-0.93) and oropharyngeal cancer (0.31, 0.10-0.95), regardless of HPV. CONCLUSION: Adaptive and dose redistributed radiotherapy enabled dose-escalation with similar toxicity rates compared to conventional radiotherapy. While FDG-PET-guided dose-escalation did overall not lead to significant tumor control or survival improvements, post-hoc results showed improved locoregional control for patients with N0-1 disease or oropharyngeal cancer treated with rRT.

Mitigating differential baseline shifts in locally advanced lung cancer patients using an average anatomy model.

BACKGROUND AND PURPOSE: Differential baseline shifts between primary tumor and involved lymph nodes in locally advanced lung cancer patients compromise the accuracy of radiotherapy. The purpose of this study was to evaluate the performance of an average anatomy model (AAM) derived from repeat imaging and deformable registration to reduce these geometrical uncertainties. METHODS AND MATERIALS: An in-house implementation of a B-Spline deformable image registration (DIR) algorithm was first validated using three different validation approaches: (a) a circle method to test the consistency of the DIR, (b) fiducial marker target registration error, and (c) the recovery of a known deformation vector field (DVF). Subsequently, AAM was generated by first averaging five DVFs resulting from cone beam CT (CBCT) to planning CT (pCT) DIR and second by applying the inverse of the average DVF to the pCT. The proposed method was evaluated on 15 locally advanced lung cancer patients receiving daily motion compensated CBCT and a repeat CT (rCT) for adaptive radiotherapy. Reduction of systematic baseline shifts of the primary tumor were quantified for the fractions used to build the AAM as well as over the whole treatment and compared to the performance of the rCT. RESULTS: The deformable registration accuracy was ≤ 2 mm vector length for the first two validation methods and about 3 mm for the third method. The systematic baseline shifts over the five fractions prior to the rCT used to build the AAM reduced from 5.9 mm vector length relative to the pCT to 2.3 and 4.2 mm relative to the AAM and rCT, respectively. The overall systematic errors in the left-right, cranio-caudal, and anterior-posterior directions were [3.4, 3.8, 3.3] mm, [2.3, 2.9, 2.6] mm, and [2.3, 3.1, 2.7] mm for the pCT, AAM, and rCT, respectively. CONCLUSIONS: The AAM mitigates systematic errors occurring during treatment due to differential baseline shifts between the primary tumor and involved lymph nodes similar to (or even better than) rCT. The superior performance of the AAM in terms of the systematic error derived from the initial fractions indicates that further analysis of the optimum intervention time is required. This model has the potential to be used as an efficient and accurate alternative for rCT in adaptive radiotherapy of locally advanced lung cancer patients, obviating the need for rescanning and recontouring.

Analysis of GTV reduction during radiotherapy for oropharyngeal cancer: Implications for adaptive radiotherapy.

BACKGROUND AND PURPOSE: Adaptive field size reduction based on gross tumor volume (GTV) shrinkage imposes risk on coverage. Fiducial markers were used as surrogate for behavior of tissue surrounding the GTV edge to assess this risk by evaluating if GTVs during treatment are dissolving or actually shrinking. MATERIALS AND METHODS: Eight patients with oropharyngeal tumors treated with chemo-radiation were included. Before treatment, fiducial markers (0.035×0.2cm2, n=40) were implanted at the edge of the primary tumor. All patients underwent planning-CT, daily cone beam CT (CBCT) and MRIs (pre-treatment, weeks 3 and 6). Marker displacement on CBCT was compared to local GTV surface displacement on MRIs. Additionally, marker displacement relative to the GTV surfaces during treatment was measured. RESULTS: GTV surface displacement derived from MRI was larger than derived from fiducial markers (average difference: 0.1cm in week 3). During treatment, the distance between markers and GTV surface on MRI in week 3 increased in 33%>0.3cm and in 10%>0.5cm. The MRI-GTV shrank faster than the surrounding tissue represented by the markers, i.e. adapting to GTV shrinkage may cause under-dosage of microscopic disease. CONCLUSIONS: We showed that adapting to primary tumor GTV shrinkage on MRI mid-treatment is potentially not safe since at least part of the GTV is likely to be dissolving. Adjustment to clear anatomical boundaries, however, may be done safely.

Effects of anatomical changes on pencil beam scanning proton plans in locally advanced NSCLC patients.

BACKGROUND AND PURPOSE: Daily anatomical variations can cause considerable differences between delivered and planned dose. This study simulates and evaluates these effects in spot-scanning proton therapy for lung cancer patients. MATERIALS AND METHODS: Robust intensity modulated treatment plans were designed on the mid-position CT scan for sixteen locally advanced lung cancer patients. To estimate dosimetric uncertainty, deformable registration was performed on their daily CBCTs to generate 4DCT equivalent scans for each fraction and to map recomputed dose to a common frame. RESULTS: Without adaptive planning, eight patients had an undercoverage of the targets of more than 2GyE (maximum of 14.1GyE) on the recalculated treatment dose from the daily anatomy variations including respiration. In organs at risk, a maximum increase of 4.7GyE in the D1 was found in the mediastinal structures. The effect of respiratory motion alone is smaller: 1.4GyE undercoverage for targets and less than 1GyE for organs at risk. CONCLUSIONS: Daily anatomical variations over the course of treatment can cause considerable dose differences in the robust planned dose distribution. An advanced planning strategy including knowledge of anatomical uncertainties would be recommended to improve plan robustness against interfractional variations. For large anatomical changes, adaptive therapy is mandatory.

Geometric changes of parotid glands caused by hydration during chemoradiotherapy.

BACKGROUND: Plan adaptation during the course of (chemo)radiotherapy of H&N cancer requires repeat CT scanning to capture anatomy changes such as parotid gland shrinkage. Hydration, applied to prevent nephrotoxicity from cisplatin, could temporarily alter the hydrogen balance and hence the captured anatomy. The aim of this study was to determine geometric changes of parotid glands as function of hydration during chemoradiotherapy compared to a control group treated with radiotherapy only. METHODS: This study included an experimental group (n = 19) receiving chemoradiotherapy, and a control group (n = 19) receiving radiotherapy only. Chemoradiotherapy patients received cisplatin with 9 l of saline solution during hydration in the first, fourth and seventh week. The delineations of the parotid glands on the planning CT scan were automatically propagated to Cone Beam CT scans using deformable image registration. Relative volume and position of the parotid glands were determined at the second chemotherapy cycle (week four) and at fraction 35. RESULTS: When saline solution was administrated, the volume temporarily increased on the first day (7.2 %, p < 0.001), second day (10.8 %, p < 0.001) and third day (7.0 %, p = 0.016). The gland positions shifted lateral, the distance between glands increased on the first day with 1.5 mm (p < 0.001), on the second day 2.2 mm (p < 0.001). At fraction 35, with both groups the mean shrinkage was 24 % ± 11 % (1SD) and the mean medial distance between the parotid glands decreased by 0.47 cm ± 0.27 cm. CONCLUSIONS: Hydration significantly modulates parotid gland geometry. Unless, in the context of adaptive RT, a repeat CT scan is timed during a chemotherapy cycle, these effects are of minor clinical relevance.