Machine-learning with region-level radiomic and dosimetric features for predicting radiotherapy-induced rectal toxicities in prostate cancer patients.

BACKGROUND AND PURPOSE: This study aims to build machine learning models to predict radiation-induced rectal toxicities for three clinical endpoints and explore whether the inclusion of radiomic features calculated on radiotherapy planning computerised tomography (CT) scans combined with dosimetric features can enhance the prediction performance. MATERIALS AND METHODS: 183 patients recruited to the VoxTox study (UK-CRN-ID-13716) were included. Toxicity scores were prospectively collected after 2 years with grade ≥ 1 proctitis, haemorrhage (CTCAEv4.03); and gastrointestinal (GI) toxicity (RTOG) recorded as the endpoints of interest. The rectal wall on each slice was divided into 4 regions according to the centroid, and all slices were divided into 4 sections to calculate region-level radiomic and dosimetric features. The patients were split into a training set (75%, N = 137) and a test set (25%, N = 46). Highly correlated features were removed using four feature selection methods. Individual radiomic or dosimetric or combined (radiomic + dosimetric) features were subsequently classified using three machine learning classifiers to explore their association with these radiation-induced rectal toxicities. RESULTS: The test set area under the curve (AUC) values were 0.549, 0.741 and 0.669 for proctitis, haemorrhage and GI toxicity prediction using radiomic combined with dosimetric features. The AUC value reached 0.747 for the ensembled radiomic-dosimetric model for haemorrhage. CONCLUSIONS: Our preliminary results show that region-level pre-treatment planning CT radiomic features have the potential to predict radiation-induced rectal toxicities for prostate cancer. Moreover, when combined with region-level dosimetric features and using ensemble learning, the model prediction performance slightly improved.

Sub-regional analysis of the parotid glands: model development for predicting late xerostomia with radiomics features in head and neck cancer patients.

BACKGROUND: The irradiation of sub-regions of the parotid has been linked to xerostomia development in patients with head and neck cancer (HNC). In this study, we compared the xerostomia classification performance of radiomics features calculated on clinically relevant and de novo sub-regions of the parotid glands of HNC patients. MATERIAL AND METHODS: All patients (N = 117) were treated with TomoTherapy in 30-35 fractions of 2-2.167 Gy per fraction with daily mega-voltage-CT (MVCT) acquisition for image-guidance purposes. Radiomics features (N = 123) were extracted from daily MVCTs for the whole parotid gland and nine sub-regions. The changes in feature values after each complete week of treatment were considered as predictors of xerostomia (CTCAEv4.03, grade ≥ 2) at 6 and 12 months. Combinations of predictors were generated following the removal of statistically redundant information and stepwise selection. The classification performance of the logistic regression models was evaluated on train and test sets of patients using the Area Under the Curve (AUC) associated with the different sub-regions at each week of treatment and benchmarked with the performance of models solely using dose and toxicity at baseline. RESULTS: In this study, radiomics-based models predicted xerostomia better than standard clinical predictors. Models combining dose to the parotid and xerostomia scores at baseline yielded an AUCtest of 0.63 and 0.61 for xerostomia prediction at 6 and 12 months after radiotherapy while models based on radiomics features extracted from the whole parotid yielded a maximum AUCtest of 0.67 and 0.75, respectively. Overall, across sub-regions, maximum AUCtest was 0.76 and 0.80 for xerostomia prediction at 6 and 12 months. Within the first two weeks of treatment, the cranial part of the parotid systematically yielded the highest AUCtest. CONCLUSION: Our results indicate that variations of radiomics features calculated on sub-regions of the parotid glands can lead to earlier and improved prediction of xerostomia in HNC patients.

Assessing the generalisability of radiomics features previously identified as predictive of radiation-induced sticky saliva and xerostomia.

Background and purpose: While core to the scientific approach, reproducibility of experimental results is challenging in radiomics studies. A recent publication identified radiomics features that are predictive of late irradiation-induced toxicity in head and neck cancer (HNC) patients. In this study, we assessed the generalisability of these findings. Materials and Methods: The procedure described in the publication in question was applied to a cohort of 109 HNC patients treated with 50-70 Gy in 20-35 fractions using helical radiotherapy although there were inherent differences between the two patient populations and methodologies. On each slice of the planning CT with delineated parotid and submandibular glands, the imaging features that were previously identified as predictive of moderate-to-severe xerostomia and sticky saliva 12 months post radiotherapy (Xer12m and SS12m) were calculated. Specifically, Short Run Emphasis (SRE) and maximum CT intensity (maxHU) were evaluated for improvement in prediction of Xer12m and SS12m respectively, compared to models solely using baseline toxicity and mean dose to the salivary glands. Results: None of the associations previously identified as statistically significant and involving radiomics features in univariate or multivariate models could be reproduced on our cohort. Conclusion: The discrepancies observed between the results of the two studies delineate limits to the generalisability of the previously reported findings. This may be explained by the differences in the approaches, in particular the imaging characteristics and subsequent methodological implementation. This highlights the importance of external validation, high quality reporting guidelines and standardisation protocols to ensure generalisability, replication and ultimately clinical implementation.

Predicting radiotherapy-induced xerostomia in head and neck cancer patients using day-to-day kinetics of radiomics features.

Background and purpose: The images acquired during radiotherapy for image-guidance purposes could be used to monitor patient-specific response to irradiation and improve treatment personalisation. We investigated whether the kinetics of radiomics features from daily mega-voltage CT image-guidance scans (MVCT) improve prediction of moderate-to-severe xerostomia compared to dose/volume parameters in radiotherapy of head-and-neck cancer (HNC). Materials and Methods: All included HNC patients (N = 117) received 30 or more fractions of radiotherapy with daily MVCTs. Radiomics features were calculated on the contra-lateral parotid glands of daily MVCTs. Their variations over time after each complete week of treatment were used to predict moderate-to-severe xerostomia (CTCAEv4.03 grade ≥ 2) at 6, 12 and 24 months post-radiotherapy. After dimensionality reduction, backward/forward selection was used to generate combinations of predictors.Three types of logistic regression model were generated for each follow-up time: 1) a pre-treatment reference model using dose/volume parameters, 2) a combination of dose/volume and radiomics-based predictors, and 3) radiomics-based predictors. The models were internally validated by cross-validation and bootstrapping and their performance evaluated using Area Under the Curve (AUC) on separate training and testing sets. Results: Moderate-to-severe xerostomia was reported by 46 %, 33 % and 26 % of the patients at 6, 12 and 24 months respectively. The selected models using radiomics-based features extracted at or before mid-treatment outperformed the dose-based models with an AUCtrain/AUCtest of 0.70/0.69, 0.76/0.74, 0.86/0.86 at 6, 12 and 24 months, respectively. Conclusion: Our results suggest that radiomics features calculated on MVCTs from the first half of the radiotherapy course improve prediction of moderate-to-severe xerostomia in HNC patients compared to a dose-based pre-treatment model.

Estimating the percentage of patients who might benefit from proton beam therapy instead of X-ray radiotherapy.

OBJECTIVES: High-energy Proton Beam Therapy (PBT) commenced in England in 2018 and NHS England commissions PBT for 1.5% of patients receiving radical radiotherapy. We sought expert opinion on the level of provision. METHODS: Invitations were sent to 41 colleagues working in PBT, most at one UK centre, to contribute by completing a spreadsheet. 39 responded: 23 (59%) completed the spreadsheet; 16 (41%) declined, arguing that clinical outcome data are lacking, but joined six additional site-specialist oncologists for two consensus meetings. The spreadsheet was pre-populated with incidence data from Cancer Research UK and radiotherapy use data from the National Cancer Registration and Analysis Service. 'Mechanisms of Benefit' of reduced growth impairment, reduced toxicity, dose escalation and reduced second cancer risk were examined. RESULTS: The most reliable figure for percentage of radical radiotherapy patients likely to benefit from PBT was that agreed by 95% of the 23 respondents at 4.3%, slightly larger than current provision. The median was 15% (range 4-92%) and consensus median 13%. The biggest estimated potential benefit was from reducing toxicity, median benefit to 15% (range 4-92%), followed by dose escalation median 3% (range 0 to 47%); consensus values were 12 and 3%. Reduced growth impairment and reduced second cancer risk were calculated to benefit 0.5% and 0.1%. CONCLUSIONS: The most secure estimate of percentage benefit was 4.3% but insufficient clinical outcome data exist for confident estimates. The study supports the NHS approach of using the evidence base and developing it through randomised trials, non-randomised studies and outcomes tracking. ADVANCES IN KNOWLEDGE: Less is known about the percentage of patients who may benefit from PBT than is generally acknowledged. Expert opinion varies widely. Insufficient clinical outcome data exist to provide robust estimates. Considerable further work is needed to address this, including international collaboration; much is already underway but will take time to provide mature data.

50 years of radiotherapy research: Evolution, trends and lessons for the future.

Rapid and relentless technological advances in an ever-more globalized world have shaped the field of radiation oncology in which we practise today. These developments have drastically modified the habitus1 of health professionals and researchers at an individual and organisational level. In this article we present an analysis of trends in radiation oncology research over the last half a century. To do so, the data from >350,000 scientific publications pertaining to a yearly search of the PubMed database with the keywords cancer radiotherapy was analysed. This analysis revealed that, over the years, radiotherapy research output has declined relative to alternative cancer therapies, representing 64% in 1970 it decreased to 31% in 2019. Also, the pace of research has significantly accelerated with, in the last 15 years, a doubling in the number of articles published by the 10% most productive researchers. Researchers are also facing stronger competition today with a proportion of first authors that will never get to publish as a last author increasing steadily from 58% in 1970 to 84% in 2000. Additionally, radiotherapy research output is extremely unequally distributed in the world, with Africa and South America contributing to ∼3% of radiotherapy articles in 2019 while representing 23% of the world's population. This disparity, reflecting economic situations and radiotherapy capabilities, has a knock-on effect for the provision of routine clinical treatment. Since research activity is inherent to delivery of high quality clinical care, this contributes to the global inequity of radiotherapy services. Learning from these trends is crucial for the future not only of radiation oncology research but also for effective and equitable cancer care.

Associations between voxel-level accumulated dose and rectal toxicity in prostate radiotherapy.

BACKGROUND AND PURPOSE: Associations between dose and rectal toxicity in prostate radiotherapy are generally poorly understood. Evaluating spatial dose distributions to the rectal wall (RW) may lead to improvements in dose-toxicity modelling by incorporating geometric information, masked by dose-volume histograms. Furthermore, predictive power may be strengthened by incorporating the effects of interfraction motion into delivered dose calculations.Here we interrogate 3D dose distributions for patients with and without toxicity to identify rectal subregions at risk (SRR), and compare the discriminatory ability of planned and delivered dose. MATERIAL AND METHODS: Daily delivered dose to the rectum was calculated using image guidance scans, and accumulated at the voxel level using biomechanical finite element modelling. SRRs were statistically determined for rectal bleeding, proctitis, faecal incontinence and stool frequency from a training set (n = 139), and tested on a validation set (n = 47). RESULTS: SRR patterns differed per endpoint. Analysing dose to SRRs improved discriminative ability with respect to the full RW for three of four endpoints. Training set AUC and OR analysis produced stronger toxicity associations from accumulated dose than planned dose. For rectal bleeding in particular, accumulated dose to the SRR (AUC 0.76) improved upon dose-toxicity associations derived from planned dose to the RW (AUC 0.63). However, validation results could not be considered significant. CONCLUSIONS: Voxel-level analysis of dose to the RW revealed SRRs associated with rectal toxicity, suggesting non-homogeneous intra-organ radiosensitivity. Incorporating spatial features of accumulated delivered dose improved dose-toxicity associations. This may be an important tool for adaptive radiotherapy in the future.

Anatomical change during radiotherapy for head and neck cancer, and its effect on delivered dose to the spinal cord.

BACKGROUND AND PURPOSE: The impact of weight loss and anatomical change during head and neck (H&N) radiotherapy on spinal cord dosimetry is poorly understood, limiting evidence-based adaptive management strategies. MATERIALS AND METHODS: 133 H&N patients treated with daily mega-voltage CT image-guidance (MVCT-IG) on TomoTherapy, were selected. Elastix software was used to deform planning scan SC contours to MVCT-IG scans, and accumulate dose. Planned (DP) and delivered (DA) spinal cord D2% (SCD2%) were compared. Univariate relationships between neck irradiation strategy (unilateral vs bilateral), T-stage, N-stage, weight loss, and changes in lateral separation (LND) and CT slice surface area (SSA) at C1 and the superior thyroid notch (TN), and ΔSCD2% [(DA - DP) D2%] were examined. RESULTS: The mean value for (DA - DP) D2% was -0.07 Gy (95%CI -0.28 to 0.14, range -5.7 Gy to 3.8 Gy), and the mean absolute difference between DP and DA (independent of difference direction) was 0.9 Gy (95%CI 0.76-1.04 Gy). Neck treatment strategy (p = 0.39) and T-stage (p = 0.56) did not affect ΔSCD2%. Borderline significance (p = 0.09) was seen for higher N-stage (N2-3) and higher ΔSCD2%. Mean reductions in anatomical metrics were substantial: weight loss 6.8 kg; C1LND 12.9 mm; C1SSA 12.1 cm2; TNLND 5.3 mm; TNSSA 11.2 cm2, but no relationship between weight loss or anatomical change and ΔSCD2% was observed (all r2 < 0.1). CONCLUSIONS: Differences between delivered and planned spinal cord D2% are small in patients treated with daily IG. Even patients experiencing substantial weight loss or anatomical change during treatment do not require adaptive replanning for spinal cord safety.

Identifying risk factors for L'Hermitte's sign after IMRT for head and neck cancer.

BACKGROUND: L'Hermitte's sign (LS) after chemoradiotherapy for head and neck cancer appears related to higher spinal cord doses. IMRT plans limit spinal cord dose, but the incidence of LS remains high. METHODS: One hundred seventeen patients treated with TomoTherapy™ between 2008 and 2015 prospectively completed a side-effect questionnaire (VoxTox Trial Registration: UK CRN ID 13716). Baseline patient and treatment data were collected. Radiotherapy plans were analysed; mean and maximum spinal cord dose and volumes receiving 10, 20, 30 and 40 Gy were recorded. Dose variation across the cord was examined. These data were included in a logistic regression model. RESULTS: Forty two patients (35.9%) reported LS symptoms. Concurrent weekly cisplatin did not increase LS risk (p = 0.70, OR = 1.23 {95% CI 0.51-2.34}). Of 13 diabetic participants (9 taking metformin), only 1 developed LS (p = 0.025, OR = 0.13 {95% CI 0.051-3.27}). A refined binary logistic regression model showed that patients receiving unilateral radiation (p = 0.019, OR = 2.06 {95% CI 0.15-0.84}) were more likely to develop LS. Higher V40Gy (p = 0.047, OR = 1.06 {95% CI 1.00-1.12}), and younger age (mean age 56.6 vs 59.7, p = 0.060, OR = 0.96 {95% CI 0.92-1.00}) were associated with elevated risk of LS, with borderline significance. CONCLUSIONS: In this cohort, concomitant cisplatin did not increase risk, and LS incidence was lower in diabetic patients. Patient age and dose gradients across the spinal cord may be important factors.

The future of image-guided radiotherapy-is image everything?

MR-based image-guided (IG) radiotherapy via all-in-one MR treatment units (MR-linacs) is one of the hottest topics in contemporary radiotherapy research. From ingenious engineering solutions to complex physical problems, researchers have developed machines with the promise of superior image quality, and all the advantages this may confer. Benefits include better tumour visualisation, online adaptation and the potential for image biomarker-based personalised RT. However, it is important to remember that the technical challenges are real. In many instances, they are skillfully managed rather than abolished, a point illustrated by the wide variety of MR-linac designs. The proposed benefits also deserve careful inspection. Better visibility of the primary tumour on an IG scan cannot be bad, but does not automatically equate to better IG, which often depends on a more generalised match to daily anatomy. MR-linac will undoubtedly be a rich milieu to search for IMBs, but these will need to be carefully validated, and similar work with CT-based biomarkers using existing, cheaper, and more widely available hardware is currently ongoing. Online adaptation is an attractive concept, but practicalities are complex, and more work is required to understand which patients will benefit from plan adaptation, and when. Finally, the issue of cost cannot be overlooked, nor can the research community's responsibilities to global healthcare inequalities. MR-linac is an exciting and ingenious technology, which merits both investment and research. It may not, however, have the future to itself.

Applying physical science techniques and CERN technology to an unsolved problem in radiation treatment for cancer: the multidisciplinary 'VoxTox' research programme.

The VoxTox research programme has applied expertise from the physical sciences to the problem of radiotherapy toxicity, bringing together expertise from engineering, mathematics, high energy physics (including the Large Hadron Collider), medical physics and radiation oncology. In our initial cohort of 109 men treated with curative radiotherapy for prostate cancer, daily image guidance computed tomography (CT) scans have been used to calculate delivered dose to the rectum, as distinct from planned dose, using an automated approach. Clinical toxicity data have been collected, allowing us to address the hypothesis that delivered dose provides a better predictor of toxicity than planned dose.

Long-term tumour control in sacral chordoma following high-dose palliative image-guided intensity-modulated radiotherapy (IG-IMRT).

Chordomas can present a challenge to the radiation oncologist and surgeon owing to the proximity to neurological structures. We describe a case of long-term tumour control in a chordoma of the lumbar spine following high-dose palliative radiotherapy. Image-guided intensity-modulated radiotherapy with photons provided a good solution to deliver 65 Gy to the tumour in a technically challenging case, and local control has been sustained over a period of years.

Fast imaging employing steady-state acquisition (FIESTA) MRI to investigate cerebrospinal fluid (CSF) within dural reflections of posterior fossa cranial nerves.

OBJECTIVE: There is no consensus approach to covering skull base meningeal reflections-and cerebrospinal fluid (CSF) therein-of the posterior fossa cranial nerves (CNs VII-XII) when planning radiotherapy (RT) for medulloblastoma and ependymoma. We sought to determine whether MRI and specifically fast imaging employing steady-state acquisition (FIESTA) sequences can answer this anatomical question and guide RT planning. METHODS: 96 posterior fossa FIESTA sequences were reviewed. Following exclusions, measurements were made on the following scans for each foramen respectively (left, right); internal acoustic meatus (IAM) (86, 84), jugular foramen (JF) (83, 85) and hypoglossal canal (HC) (42, 45). A protocol describes measurement procedure. Two observers measured distances for five cases and agreement was assessed. One observer measured all the remaining cases. RESULTS: IAM and JF measurement interobserver variability was compared. Mean measurement difference between observers was -0.275 mm (standard deviation 0.557). IAM and JF measurements were normally distributed. Mean IAM distance was 12.2 mm [95% confidence interval (CI) 8.8-15.6]; JF was 7.3 mm (95% CI 4.0-10.6). The HC was difficult to visualize on many images and data followed a bimodal distribution. CONCLUSION: Dural reflections of posterior fossa CNs are well demonstrated by FIESTA MRI. Measuring CSF extension into these structures is feasible and robust; mean CSF extension into IAM and JF was measured. We plan further work to assess coverage of these structures with photon and proton RT plans. Advances in knowledge: We have described CSF extension beyond the internal table of the skull into the IAM, JF and HC. Oncologists planning RT for patients with medulloblastoma and ependymoma may use these data to guide contouring.