Innovative Education Method for a More Effective, Faster, and Valued Training in Radiation Therapy Treatment Planning.

Purpose: Because of the automation of radiation therapy, competencies of radiation technologists (RTTs) change, and training methods are challenged. This study aims to develop, and pilot test an innovative training method based on lean management principles. Methods and Materials: A new training method was developed for lung cancer treatment planning (TP). The novelty is summarized by including a stable environment and an increased focus on the how and why of key decision making. Trainees have to motivate their decisions during TP process, and to argue their choices with peers. Six students and 6 RTTs completed this training for lung cancer TP. Effects of the training were measured by (1) quality of TP, using doses in organs at risk and target volumes, (2) perceived experiences (survey), measured at baseline (T0); after peer session (T1); and 6 months later (T2). Finally, training throughput time was measured. Results: At T0, RTTs showed a larger intragroup interquartile range (IIR) (2.63Gy vs 1.51Gy), but lower mean doses to heart and esophagus than students (6.79Gy vs 8.49Gy; 20.87Gy vs 24.62Gy). At T1, quality of TPs was similar between RTTs and students (IIR: 1.39Gy vs 1.33Gy) and no significant differences in mean dose to heart and esophagus (4.48Gy vs 4.69Gy; 17.75Gy vs 18.47Gy). At T2, students still performed equal to RTTs (IIR: 1.07Gy vs 1.45Gy) and achieved lower maximum dose to esophagus (44.75Gy vs 46.45Gy). The training method and peer sessions were experienced positive: at baseline (T0): 8 score on a scale 1-10, directly after the peer sessions; (T1): 8 by the students and 7 by the RTTs, after 9 months; (T2): 9 by the students and 7 by the RTTs. Training throughput time decreased from 12 to 3 months. Conclusions: This training method based on lean management principles was successfully applied to training of RTTs for lung cancer TP. Training throughput time was reduced dramatically and TP quality sustained after 6 months. This method can potentially improve training efficiency in diverse situations with complex decision-making.

Evaluation of a cone-beam computed tomography system calibrated for accurate radiotherapy dose calculation.

Background and purpose: Dose calculation on cone-beam computed tomography (CBCT) images has been less accurate than on computed tomography (CT) images due to lower image quality and discrepancies in CT numbers for CBCT. As increasing interest arises in offline and online re-planning, dose calculation accuracy was evaluated for a novel CBCT imager integrated into a ring gantry treatment machine. Materials and methods: The new CBCT system allowed fast image acquisition (5.9 s) by using new hardware, including a large-size flat panel detector, and incorporated image-processing algorithms with iterative reconstruction techniques, leading to accurate CT numbers allowing dose calculation. In this study, CBCT- and CT-based dose calculations were compared based on three anthropomorphic phantoms, after CBCT-to-mass-density calibration was performed. Six plans were created on the CT scans covering various target locations and complexities, followed by CBCT to CT registrations, copying of contours, and re-calculation of the plans on the CBCT scans. Dose-volume histogram metrics for target volumes and organs-at-risk (OARs) were evaluated, and global gamma analyses were performed. Results: Target coverage differences were consistently below 1.2 %, demonstrating the agreement between CT and re-calculated CBCT dose distributions. Differences in Dmean for OARs were below 0.5 Gy for all plans, except for three OARs, which were below 0.8 Gy (<1.1 %). All plans had a 3 %/1mm gamma pass rate > 97 %. Conclusions: This study demonstrated comparable results between dose calculations performed on CBCT and CT acquisitions. The new CBCT system with enhanced image quality and CT number accuracy opens possibilities for off-line and on-line re-planning.

A Framework for Assessing the Effect of Cardiac and Respiratory Motion for Stereotactic Arrhythmia Radioablation Using a Digital Phantom With a 17-Segment Model: A STOPSTORM.eu Consortium Study.

PURPOSE: The optimal motion management strategy for patients receiving stereotactic arrhythmia radioablation (STAR) for the treatment of ventricular tachycardia (VT) is not fully known. We developed a framework using a digital phantom to simulate cardiorespiratory motion in combination with different motion management strategies to gain insight into the effect of cardiorespiratory motion on STAR. METHODS AND MATERIALS: The 4-dimensional (4D) extended cardiac-torso (XCAT) phantom was expanded with the 17-segment left ventricular (LV) model, which allowed placement of STAR targets in standardized ventricular regions. Cardiac- and respiratory-binned 4D computed tomography (CT) scans were simulated for free-breathing, reduced free-breathing, respiratory-gating, and breath-hold scenarios. Respiratory motion of the heart was set to population-averaged values of patients with VT: 6, 2, and 1 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction was adjusted by reducing LV ejection fraction to 35%. Target displacement was evaluated for all segments using envelopes encompassing the cardiorespiratory motion. Envelopes incorporating only the diastole plus respiratory motion were created to simulate the scenario where cardiac motion is not fully captured on 4D respiratory CT scans used for radiation therapy planning. RESULTS: The average volume of the 17 segments was 6 cm3 (1-9 cm3). Cardiac contraction-relaxation resulted in maximum segment (centroid) motion of 4, 6, and 3.5 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction-relaxation resulted in a motion envelope increase of 49% (24%-79%) compared with individual segment volumes, whereas envelopes increased by 126% (79%-167%) if respiratory motion also was considered. Envelopes incorporating only the diastole and respiration motion covered on average 68% to 75% of the motion envelope. CONCLUSIONS: The developed LV-segmental XCAT framework showed that free-wall regions display the most cardiorespiratory displacement. Our framework supports the optimization of STAR by evaluating the effect of (cardio)respiratory motion and motion management strategies for patients with VT.

Clinical benefit of range uncertainty reduction in proton treatment planning based on dual-energy CT for neuro-oncological patients.

OBJECTIVE: Several studies have shown that dual-energy CT (DECT) can lead to improved accuracy for proton range estimation. This study investigated the clinical benefit of reduced range uncertainty, enabled by DECT, in robust optimisation for neuro-oncological patients. METHODS: DECT scans for 27 neuro-oncological patients were included. Commercial software was applied to create stopping-power ratio (SPR) maps based on the DECT scan. Two plans were robustly optimised on the SPR map, keeping the beam and plan settings identical to the clinical plan. One plan was robustly optimised and evaluated with a range uncertainty of 3% (as used clinically; denoted 3%-plan); the second plan applied a range uncertainty of 2% (2%-plan). Both plans were clinical acceptable and optimal. The dose-volume histogram parameters were compared between the two plans. Two experienced neuro-radiation oncologists determined the relevant dose difference for each organ-at-risk (OAR). Moreover, the OAR toxicity levels were assessed. RESULTS: For 24 patients, a dose reduction >0.5/1 Gy (relevant dose difference depending on the OAR) was seen in one or more OARs for the 2%-plan; e.g. for brainstem D0.03cc in 10 patients, and hippocampus D40% in 6 patients. Furthermore, 12 patients had a reduction in toxicity level for one or two OARs, showing a clear benefit for the patient. CONCLUSION: Robust optimisation with reduced range uncertainty allows for reduction of OAR toxicity, providing a rationale for clinical implementation. Based on these results, we have clinically introduced DECT-based proton treatment planning for neuro-oncological patients, accompanied with a reduced range uncertainty of 2%. ADVANCES IN KNOWLEDGE: This study shows the clinical benefit of range uncertainty reduction from 3% to 2% in robustly optimised proton plans. A dose reduction to one or more OARs was seen for 89% of the patients, and 44% of the patients had an expected toxicity level decrease.

Transferability of radiomic signatures from experimental to human interstitial lung disease.

Background: Interstitial lung disease (ILD) defines a group of parenchymal lung disorders, characterized by fibrosis as their common final pathophysiological stage. To improve diagnosis and treatment of ILD, there is a need for repetitive non-invasive characterization of lung tissue by quantitative parameters. In this study, we investigated whether CT image patterns found in mice with bleomycin induced lung fibrosis can be translated as prognostic factors to human patients diagnosed with ILD. Methods: Bleomycin was used to induce lung fibrosis in mice (n_control = 36, n_experimental = 55). The patient cohort consisted of 98 systemic sclerosis (SSc) patients (n_ILD = 65). Radiomic features (n_histogram = 17, n_texture = 137) were extracted from microCT (mice) and HRCT (patients) images. Predictive performance of the models was evaluated with the area under the receiver-operating characteristic curve (AUC). First, predictive performance of individual features was examined and compared between murine and patient data sets. Second, multivariate models predicting ILD were trained on murine data and tested on patient data. Additionally, the models were reoptimized on patient data to reduce the influence of the domain shift on the performance scores. Results: Predictive power of individual features in terms of AUC was highly correlated between mice and patients (r = 0.86). A model based only on mean image intensity in the lung scored AUC = 0.921 ± 0.048 in mice and AUC = 0.774 (CI95% 0.677-0.859) in patients. The best radiomic model based on three radiomic features scored AUC = 0.994 ± 0.013 in mice and validated with AUC = 0.832 (CI95% 0.745-0.907) in patients. However, reoptimization of the model weights in the patient cohort allowed to increase the model's performance to AUC = 0.912 ± 0.058. Conclusion: Radiomic signatures of experimental ILD derived from microCT scans translated to HRCT of humans with SSc-ILD. We showed that the experimental model of BLM-induced ILD is a promising system to test radiomic models for later application and validation in human cohorts.

Corrigendum to "Robustness of radiomic features in magnetic resonance imaging for patients with glioblastoma: Multi-center study" [Phys. Imaging Radiat. Oncol. 22 (2022) 131-136].

[This corrects the article DOI: 10.1016/j.phro.2022.05.006.].

Robustness of radiomic features in magnetic resonance imaging for patients with glioblastoma: Multi-center study.

Background and purpose: Radiomics offers great potential in improving diagnosis and treatment for patients with glioblastoma multiforme. However, in order to implement radiomics in clinical routine, the features used for prognostic modelling need to be stable. This comprises significant challenge in multi-center studies. The aim of this study was to evaluate the impact of different image normalization methods on MRI features robustness in multi-center study. Methods: Radiomics stability was checked on magnetic resonance images of eleven patients. The images were acquired in two different hospitals using contrast-enhanced T1 sequences. The images were normalized using one of five investigated approaches including grey-level discretization, histogram matching and z-score. Then, radiomic features were extracted and features stability was evaluated using intra-class correlation coefficients. In the second part of the study, improvement in the prognostic performance of features was tested on 60 patients derived from publicly available dataset. Results: Depending on the normalization scheme, the percentage of stable features varied from 3.4% to 8%. The histogram matching based on the tumor region showed the highest amount of the stable features (113/1404); while normalization using fixed bin size resulted in 48 stable features. The histogram matching also led to better prognostic value (median c-index increase of 0.065) comparing to non-normalized images. Conclusions: MRI normalization plays an important role in radiomics. Appropriate normalization helps to select robust features, which can be used for prognostic modelling in multicenter studies. In our study, histogram matching based on tumor region improved both stability of radiomic features and their prognostic value.

Improved Survival Prediction by Combining Radiological Imaging and S-100B Levels Into a Multivariate Model in Metastatic Melanoma Patients Treated With Immune Checkpoint Inhibition.

Purpose: We explored imaging and blood bio-markers for survival prediction in a cohort of patients with metastatic melanoma treated with immune checkpoint inhibition. Materials and Methods: 94 consecutive metastatic melanoma patients treated with immune checkpoint inhibition were included into this study. PET/CT imaging was available at baseline (Tp0), 3 months (Tp1) and 6 months (Tp2) after start of immunotherapy. Radiological response at Tp2 was evaluated using iRECIST. Total tumor burden (TB) at each time-point was measured and relative change of TB compared to baseline was calculated. LDH, CRP and S-100B were also analyzed. Cox proportional hazards model and logistic regression were used for survival analysis. Results: iRECIST at Tp2 was significantly associated with overall survival (OS) with C-index=0.68. TB at baseline was not associated with OS, whereas TB at Tp1 and Tp2 provided similar predictive power with C-index of 0.67 and 0.71, respectively. Appearance of new metastatic lesions during follow-up was an independent prognostic factor (C-index=0.73). Elevated LDH and S-100B ratios at Tp2 were significantly associated with worse OS: C-index=0.73 for LDH and 0.73 for S-100B. Correlation of LDH with TB was weak (r=0.34). A multivariate model including TB change, S-100B, and appearance of new lesions showed the best predictive performance with C-index=0.83. Conclusion: Our analysis shows only a weak correlation between LDH and TB. Additionally, baseline TB was not a prognostic factor in our cohort. A multivariate model combining early blood and imaging biomarkers achieved the best predictive power with regard to survival, outperforming iRECIST.

A 2.5D convolutional neural network for HPV prediction in advanced oropharyngeal cancer.

BACKGROUND: Infection with human papilloma virus (HPV) is one of the most relevant prognostic factors in advanced oropharyngeal cancer (OPC) treatment. In this study we aimed to assess the diagnostic accuracy of a deep learning-based method for HPV status prediction in computed tomography (CT) images of advanced OPC. METHOD: An internal dataset and three public collections were employed (internal: n = 151, HNC1: n = 451; HNC2: n = 80; HNC3: n = 110). Internal and HNC1 datasets were used for training, whereas HNC2 and HNC3 collections were used as external test cohorts. All CT scans were resampled to a 2 mm3 resolution and a sub-volume of 72x72x72 pixels was cropped on each scan, centered around the tumor. Then, a 2.5D input of size 72x72x3 pixels was assembled by selecting the 2D slice containing the largest tumor area along the axial, sagittal and coronal planes, respectively. The convolutional neural network employed consisted of the first 5 modules of the Xception model and a small classification network. Ten-fold cross-validation was applied to evaluate training performance. At test time, soft majority voting was used to predict HPV status. RESULTS: A final training mean [range] area under the curve (AUC) of 0.84 [0.76-0.89], accuracy of 0.76 [0.64-0.83] and F1-score of 0.74 [0.62-0.83] were achieved. AUC/accuracy/F1-score values of 0.83/0.75/0.69 and 0.88/0.79/0.68 were achieved on the HNC2 and HNC3 test sets, respectively. CONCLUSION: Deep learning was successfully applied and validated in two external cohorts to predict HPV status in CT images of advanced OPC, proving its potential as a support tool in cancer precision medicine.

MR-Guided Adaptive Radiotherapy for Head and Neck Cancer: Prospective Evaluation of Migration and Anatomical Changes of the Major Salivary Glands.

The aim of this study was to quantify anatomical changes of parotids and submandibular glands and evaluate potential dosimetric advantages during weekly adaptive MR-guided radiotherapy (MRgRT) for the definitive treatment of head and neck cancer (HNC). The data and plans of 12 patients treated with bilateral intensity-modulated radiotherapy for HNC using MR-linac, with weekly offline adaptations, were prospectively evaluated. The positional and volumetric changes of the salivary glands were analyzed by manual segmentation in weekly MRI images and the dosimetric impact of these anatomical changes on the adapted treatment plans was assessed. The mean volume change in parotid and submandibular gland volume was -31.9% (p < 0.0001) and -29.7% (p < 0.0001) after five weeks, respectively. The volume change was significantly correlated with the cumulative dose for the respective gland at the time of volume measurement. Inter-parotid distance changed by -5.4% (6.5 mm) on average after five weeks (p = 0.0005). The distance became significantly smaller only in the left-right direction. The inter-submandibular gland distance changed by 0.7 mm (p = 0.38). This study demonstrated significant changes in salivary gland volumes and position following daily MR guidance and weekly plan adaptation. Ongoing clinical trials will provide data on the clinical impact of these changes and novel MR-based adaptation strategies.

Quantification of the spatial distribution of primary tumors in the lung to develop new prognostic biomarkers for locally advanced NSCLC.

The anatomical location and extent of primary lung tumors have shown prognostic value for overall survival (OS). However, its manual assessment is prone to interobserver variability. This study aims to use data driven identification of image characteristics for OS in locally advanced non-small cell lung cancer (NSCLC) patients. Five stage IIIA/IIIB NSCLC patient cohorts were retrospectively collected. Patients were treated either with radiochemotherapy (RCT): RCT1* (n = 107), RCT2 (n = 95), RCT3 (n = 37) or with surgery combined with radiotherapy or chemotherapy: S1* (n = 135), S2 (n = 55). Based on a deformable image registration (MIM Vista, 6.9.2.), an in-house developed software transferred each primary tumor to the CT scan of a reference patient while maintaining the original tumor shape. A frequency-weighted cumulative status map was created for both exploratory cohorts (indicated with an asterisk), where the spatial extent of the tumor was uni-labeled with 2 years OS. For the exploratory cohorts, a permutation test with random assignment of patient status was performed to identify regions with statistically significant worse OS, referred to as decreased survival areas (DSA). The minimal Euclidean distance between primary tumor to DSA was extracted from the independent cohorts (negative distance in case of overlap). To account for the tumor volume, the distance was scaled with the radius of the volume-equivalent sphere. For the S1 cohort, DSA were located at the right main bronchus whereas for the RCT1 cohort they further extended in cranio-caudal direction. In the independent cohorts, the model based on distance to DSA achieved performance: AUCRCT2 [95% CI] = 0.67 [0.55-0.78] and AUCRCT3 = 0.59 [0.39-0.79] for RCT patients, but showed bad performance for surgery cohort (AUCS2 = 0.52 [0.30-0.74]). Shorter distance to DSA was associated with worse outcome (p = 0.0074). In conclusion, this explanatory analysis quantifies the value of primary tumor location for OS prediction based on cumulative status maps. Shorter distance of primary tumor to a high-risk region was associated with worse prognosis in the RCT cohort.

Preselection of robust radiomic features does not improve outcome modelling in non-small cell lung cancer based on clinical routine FDG-PET imaging.

BACKGROUND: Radiomics is a promising tool for identifying imaging-based biomarkers. Radiomics-based models are often trained on single-institution datasets; however, multi-centre imaging datasets are preferred for external generalizability owing to the influence of inter-institutional scanning differences and acquisition settings. The study aim was to determine the value of preselection of robust radiomic features in routine clinical positron emission tomography (PET) images to predict clinical outcomes in locally advanced non-small cell lung cancer (NSCLC). METHODS: A total of 1404 primary tumour radiomic features were extracted from pre-treatment [18F]fluorodeoxyglucose (FDG)-PET scans of stage IIIA/N2 or IIIB NSCLC patients using a training cohort (n = 79; prospective Swiss multi-centre randomized phase III trial SAKK 16/00; 16 centres) and an internal validation cohort (n = 31; single centre). Robustness studies investigating delineation variation, attenuation correction and motion were performed (intraclass correlation coefficient threshold > 0.9). Two 12-/24-month event-free survival (EFS) and overall survival (OS) logistic regression models were trained using standardized imaging: (1) with robust features alone and (2) with all available features. Models were then validated using fivefold cross-validation, and validation on a separate single-centre dataset. Model performance was assessed using area under the receiver operating characteristic curve (AUC). RESULTS: Robustness studies identified 179 stable features (13%), with 25% stable features for 3D versus 4D acquisition, 31% for attenuation correction and 78% for delineation. Univariable analysis found no significant robust features predicting 12-/24-month EFS and 12-month OS (p value > 0.076). Prognostic models without robust preselection performed well for 12-month EFS in training (AUC = 0.73) and validation (AUC = 0.74). Patient stratification into two risk groups based on 12-month EFS was significant for training (p value = 0.02) and validation cohorts (p value = 0.03). CONCLUSIONS: A PET-based radiomics model using a standardized, multi-centre dataset to predict EFS in locally advanced NSCLC was successfully established and validated with good performance. Prediction models with robust feature preselection were unsuccessful, indicating the need for a standardized imaging protocol.

Targeting Treatment Resistance in Head and Neck Squamous Cell Carcinoma - Proof of Concept for CT Radiomics-Based Identification of Resistant Sub-Volumes.

PURPOSE: Radiomics has already been proposed as a prognostic biomarker in head and neck cancer (HNSCC). However, its predictive power in radiotherapy has not yet been studied. Here, we investigated a local radiomics approach to distinguish between tumor sub-volumes with different levels of radiosensitivity as a possible target for radiation dose intensification. MATERIALS AND METHODS: Of 40 patients (n=28 training and n=12 validation) with biopsy confirmed locally recurrent HNSCC, pretreatment contrast-enhanced CT images were registered with follow-up PET/CT imaging allowing identification of controlled (GTVcontrol) vs non-controlled (GTVrec) tumor sub-volumes on pretreatment imaging. A bi-regional model was built using radiomic features extracted from pretreatment CT in the GTVrec and GTVcontrol to differentiate between those regions. Additionally, concept of local radiomics was implemented to perform detection task. The original tumor volume was divided into sub-volumes with no prior information on the location of recurrence. Radiomic features from those sub-volumes were then used to detect recurrent sub-volumes using multivariable logistic regression. RESULTS: Radiomic features extracted from non-controlled regions differed significantly from those in controlled regions (training AUC = 0.79 CI 95% 0.66 - 0.91 and validation AUC = 0.88 CI 95% 0.72 - 1.00). Local radiomics analysis allowed efficient detection of non-controlled sub-volumes both in the training AUC = 0.66 (CI 95% 0.56 - 0.75) and validation cohort 0.70 (CI 95% 0.53 - 0.86), however performance of this model was inferior to bi-regional model. Both models indicated that sub-volumes characterized by higher heterogeneity were linked to tumor recurrence. CONCLUSION: Local radiomics is able to detect sub-volumes with decreased radiosensitivity, associated with location of tumor recurrence in HNSCC in the pre-treatment CT imaging. This proof of concept study, indicates that local CT radiomics can be used as predictive biomarker in radiotherapy and potential target for dose intensification.

Radiomic Analysis to Predict Outcome in Recurrent Glioblastoma Based on Multi-Center MR Imaging From the Prospective DIRECTOR Trial.

BACKGROUND: Based on promising results from radiomic approaches to predict O6-methylguanine DNA methyltransferase promoter methylation status (MGMT status) and clinical outcome in patients with newly diagnosed glioblastoma, the current study aimed to evaluate radiomics in recurrent glioblastoma patients. METHODS: Pre-treatment MR-imaging data of 69 patients enrolled into the DIRECTOR trial in recurrent glioblastoma served as a training cohort, and 49 independent patients formed an external validation cohort. Contrast-enhancing tumor and peritumoral volumes were segmented on MR images. 180 radiomic features were extracted after application of two MR intensity normalization techniques: fixed number of bins and linear rescaling. Radiomic feature selection was performed via principal component analysis, and multivariable models were trained to predict MGMT status, progression-free survival from first salvage therapy, referred to herein as PFS2, and overall survival (OS). The prognostic power of models was quantified with concordance index (CI) for survival data and area under receiver operating characteristic curve (AUC) for the MGMT status. RESULTS: We established and validated a radiomic model to predict MGMT status using linear intensity interpolation and considering features extracted from gadolinium-enhanced T1-weighted MRI (training AUC = 0.670, validation AUC = 0.673). Additionally, models predicting PFS2 and OS were found for the training cohort but were not confirmed in our validation cohort. CONCLUSIONS: A radiomic model for prediction of MGMT promoter methylation status from tumor texture features in patients with recurrent glioblastoma was successfully established, providing a non-invasive approach to anticipate patient's response to chemotherapy if biopsy cannot be performed. The radiomic approach to predict PFS2 and OS failed.

Comparison of beam segment versus full plan re-optimization in daily magnetic resonance imaging-guided online-adaptive radiotherapy.

The optimal approach for magnetic resonance imaging-guided online adaptive radiotherapy is currently unknown and needs to consider patient on-couch time constraints. The aim of this study was to compare two different plan optimization approaches. The comparison was performed in 238 clinically applied online-adapted treatment plans from 55 patients, in which the approach of re-optimization was selected based on the physician's choice. For 33 patients where both optimization approaches were used at least once, the median treatment planning dose metrics of both target and organ at risk differed less than 1%. Therefore, we concluded that beam segment weight optimization was chosen adequately for most patients without compromising plan quality.

Impact of CT convolution kernel on robustness of radiomic features for different lung diseases and tissue types.

OBJECTIVES: In this study, we aimed to assess the impact of different CT reconstruction kernels on the stability of radiomic features and the transferability between different diseases and tissue types. Three lung diseases were evaluated, i.e. non-small cell lung cancer (NSCLC), malignant pleural mesothelioma (MPM) and interstitial lung disease related to systemic sclerosis (SSc-ILD) as well as four different tissue types, i.e. primary tumor, largest involved lymph node ipsilateral and contralateral lung. METHODS: Pre-treatment non-contrast enhanced CT scans from 23 NSCLC, 10 MPM and 12 SSc-ILD patients were collected retrospectively. For each patient, CT scans were reconstructed using smooth and sharp kernel in filtered back projection. The regions of interest (ROIs) were contoured on the smooth kernel-based CT and transferred to the sharp kernel-based CT. The voxels were resized to the largest voxel dimension of each cohort. In total, 1386 features were analyzed. Feature stability was assessed using the intraclass correlation coefficient. Features above the stability threshold >0.9 were considered stable. RESULTS: We observed a strong impact of the reconstruction method on stability of the features (at maximum 26% of the 1386 features were stable). Intensity features were the most stable followed by texture and wavelet features. The wavelet features showed a positive correlation between percentage of stable features and size of the ROI (R2 = 0.79, p = 0.005). Lymph node radiomics showed poorest stability (<10%) and lung radiomics the largest stability (26%). Robustness analysis done on the contralateral lung could to a large extent be transferred to the ipsilateral lung, and the overlap of stable lung features between different lung diseases was more than 50%. However, results of robustness studies cannot be transferred between tissue types, which was investigated in NSCLC and MPM patients; the overlap of stable features for lymph node and lung, as well as for primary tumor and lymph node was very small in both disease types. CONCLUSION: The robustness of radiomic features is strongly affected by different reconstruction kernels. The effect is largely influenced by the tissue type and less by the disease type. ADVANCES IN KNOWLEDGE: The study presents to our knowledge the most complete analysis on the impact of convolution kernel on the robustness of CT-based radiomics for four relevant tissue types in three different lung diseases. .

Radiomics Feature Activation Maps as a New Tool for Signature Interpretability.

INTRODUCTION: In the field of personalized medicine, radiomics has shown its potential to support treatment decisions. However, the limited feature interpretability hampers its introduction into the clinics. Here, we propose a new methodology to create radiomics feature activation maps, which allows to identify the spatial-anatomical locations responsible for signature activation based on local radiomics. The feasibility of this technique will be studied for histological subtype differentiation (adenocarcinoma versus squamous cell carcinoma) in non-small cell lung cancer (NSCLC) using computed tomography (CT) radiomics. MATERIALS AND METHODS: Pre-treatment CT scans were collected from a multi-centric Swiss trial (training, n=73, IIIA/N2 NSCLC, SAKK 16/00) and an independent cohort (validation, n=32, IIIA/N2/IIIB NSCLC). Based on the gross tumor volume (GTV), four peritumoral region of interests (ROI) were defined: lung_exterior (expansion into the lung), iso_exterior (expansion into lung and soft tissue), gradient (GTV border region), GTV+Rim (GTV and iso_exterior). For each ROI, 154 radiomic features were extracted using an in-house developed software implementation (Z-Rad, Python v2.7.14). Features robust against delineation variability served as an input for a multivariate logistic regression analysis. Model performance was quantified using the area under the receiver operating characteristic curve (AUC) and verified using five-fold cross validation and internal validation. Local radiomic features were extracted from the GTV+Rim ROI using non-overlapping 3x3x3 voxel patches previously marked as GTV or rim. A binary activation map was created for each patient using the median global feature value from the training. The ratios of activated/non-activated patches of GTV and rim regions were compared between histological subtypes (Wilcoxon test). RESULTS: Iso_exterior, gradient, GTV+Rim showed good performances for histological subtype prediction (AUCtraining=0.68-0.72 and AUCvalidation=0.73-0.74) whereas GTV and lung_exterior models failed validation. GTV+Rim model feature activation maps showed that local texture feature distribution differed significantly between histological subtypes in the rim (p=0.0481) but not in the GTV (p=0.461). CONCLUSION: In this exploratory study, radiomics-based prediction of NSCLC histological subtypes was predominantly based on the peritumoral region indicating that radiomics activation maps can be useful for tracing back the spatial location of regions responsible for signature activation.

Treatment plan quality during online adaptive re-planning.

BACKGROUND: Online adaptive radiotherapy is intended to prevent plan degradation caused by inter-fractional tumor volume and shape changes, but time limitations make online re-planning challenging. The aim of this study was to compare the quality of online-adapted plans to their respective reference treatment plans. METHODS: Fifty-two patients treated on a ViewRay MRIdian Linac were included in this retrospective study. In total 238 online-adapted plans were analyzed, which were optimized with either changing of the segment weights (n = 85) or full re-optimization (n = 153). Five different treatment sites were evaluated: prostate, abdomen, liver, lung and pelvis. Dosimetric parameters of gross tumor volume (GTV), planning target volume (PTV), 2 cm ring around the PTV and organs at risk (OARs) were considered. The Wilcoxon signed-rank test was used to assess differences between online-adapted and reference treatment plans, p < 0.05 was considered significant. RESULTS: The average duration of the online adaptation, consisting of contour editing, plan optimization and quality assurance (QA), was 24 ± 6 min. The GTV was slightly larger (average ± SD: 1.9% ± 9.0%) in the adapted plans than in the reference plans (p < 0.001). GTV-D95% exhibited no significant changes when considering all plans, but GTV-D2% increased by 0.40% ± 1.5% on average (p < 0.001). There was a very small yet significant decrease in GTV-coverage for the abdomen plans. The ring Dmean increased on average by 1.0% ± 3.6% considering all plans (p < 0.001). There was a significant reduction of the dose to the rectum of 4.7% ± 16% on average (p < 0.001) for prostate plans. CONCLUSIONS: Dosimetric quality of online-adapted plans was comparable to reference treatment plans and OAR dose was either comparable or decreased, depending on treatment site. However, dose spillage was slightly increased.

Radiomic biomarkers for head and neck squamous cell carcinoma.

Tumor heterogeneity is a well-known prognostic factor in head and neck squamous cell carcinoma (HNSCC). A major limitation of tissue- and blood-derived tumor markers is the lack of spatial resolution to image tumor heterogeneity. Tissue markers derived from tumor biopsies usually represent only a small tumor subregion at a single timepoint and are therefore often not representative of the tumors' biology or the biological alterations during and after treatment. Similarly, liquid biopsies give an overall picture of the tumors' secreted factors but completely lack any spatial resolution. Radiomics has the potential to give complete three-dimensional information about the tumor. We conducted a comprehensive literature search to assess the correlation of radiomics to tumor biology and treatment outcome in HNSCC and to assess current limitations of the radiomic biomarkers. In total, 25 studies that explored the ability of radiomics to predict tumor biology and phenotype in HNSCC and 28 studies that explored radiomics to predict post-treatment events were identified. Out of these 53 studies, only three failed to show a significant correlation. The major technical challenges are currently artifacts due to metal implants, non-standardized contrast injection, and delineation uncertainties. All studies to date were retrospective and none of the above-mentioned radiomics signatures have been validated in an independent cohort using an independent software implementation, which shows that transferability due to the numerous technical challenges is currently a major limitation. However, radiomics is a very young field and these studies hopefully pave the way for clinical implementation of radiomics for HNSCC in the future.

Consolidation cetuximab after concurrent triplet radiochemotherapy+cetuximab in patients with advanced head and neck cancer: A randomized phase II study.

BACKGROUND AND PURPOSE: Preclinical data suggest that cetuximab should be continued after end of concurrent radiotherapy+cetuximab due to its efficacy against residual tumor cells in the irradiated tumor bed. Based on this concept the phase II add-on cetuximab (AOC) study was designed. MATERIALS AND METHODS: Altogether 63 patients with advanced head and neck cancer were treated with radiochemotherapy (70 Gy, cisplatin 40 mg/m2 weekly) in combination with concurrent cetuximab (loading dose 400 mg/m2, then 250 mg/m2 weekly). Thereafter patients were randomized to cetuximab consolidation (500 mg/m2 biweekly × 6) or no further treatment. The primary endpoint was the 2-year locoregional control (LRC) rate. As translational research endpoints serum markers were analyzed before and during treatment and CT-based quantitative image analysis (radiomics) was performed. RESULTS: Median follow-up was 24 months. The 2-year LRC rates were 67.9% and 67.7% in the treatment arms with and without consolidation cetuximab, respectively. Higher than median levels of three serum markers were negatively associated with the 2-year LRC rate in the overall patient cohort: Osteopontin, IL8 and FasL2 (p ≤ 0.05). A radiomics model consisting of two radiomics features could be built showing that higher entropy and higher complexity of tumor Hounsfield unit distribution indicates worse LRC (concordance index 0.66). No correlation was found between biological and imaging markers. CONCLUSIONS: There was no evidence that consolidation cetuximab would improve the 2-year LRC rate. Prognostic biological and imaging markers could be identified for the overall patient cohort. Studies with larger patient numbers are needed to correlate biological and imaging markers.