The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Image-based Phenotyping.

Background Radiomic features may quantify characteristics present in medical imaging. However, the lack of standardized definitions and validated reference values have hampered clinical use. Purpose To standardize a set of 174 radiomic features. Materials and Methods Radiomic features were assessed in three phases. In phase I, 487 features were derived from the basic set of 174 features. Twenty-five research teams with unique radiomics software implementations computed feature values directly from a digital phantom, without any additional image processing. In phase II, 15 teams computed values for 1347 derived features using a CT image of a patient with lung cancer and predefined image processing configurations. In both phases, consensus among the teams on the validity of tentative reference values was measured through the frequency of the modal value and classified as follows: less than three matches, weak; three to five matches, moderate; six to nine matches, strong; 10 or more matches, very strong. In the final phase (phase III), a public data set of multimodality images (CT, fluorine 18 fluorodeoxyglucose PET, and T1-weighted MRI) from 51 patients with soft-tissue sarcoma was used to prospectively assess reproducibility of standardized features. Results Consensus on reference values was initially weak for 232 of 302 features (76.8%) at phase I and 703 of 1075 features (65.4%) at phase II. At the final iteration, weak consensus remained for only two of 487 features (0.4%) at phase I and 19 of 1347 features (1.4%) at phase II. Strong or better consensus was achieved for 463 of 487 features (95.1%) at phase I and 1220 of 1347 features (90.6%) at phase II. Overall, 169 of 174 features were standardized in the first two phases. In the final validation phase (phase III), most of the 169 standardized features could be excellently reproduced (166 with CT; 164 with PET; and 164 with MRI). Conclusion A set of 169 radiomics features was standardized, which enabled verification and calibration of different radiomics software. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Kuhl and Truhn in this issue.

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.

CT radiomics and PET radiomics: ready for clinical implementation?

INTRODUCTION: Today, rapid technical and clinical developments result in an increasing number of treatment options for oncological diseases. Thus, decision support systems are needed to offer the right treatment to the right patient. Imaging biomarkers hold great promise in patient-individual treatment guidance. Routinely performed for diagnosis and staging, imaging datasets are expected to hold more information than used in the clinical practice. Radiomics describes the extraction of a large number of meaningful quantitative features from medical images, such as computed tomography (CT) and positron emission tomography (PET). Due to the non-invasive nature and ability to capture 3D image-based heterogeneity, radiomic features are potential surrogate markers of the cancer phenotype. Several radiomic studies are published per day, owing to encouraging results of many radiomics-based patient outcome models. Despite this comparably large number of studies, radiomics is mainly studied in proof of principle concept. Hence, a translation of radiomics from a hot topic research field into an essential clinical decision-making tool is lacking, but of high clinical interest. EVIDENCE ACQUISITION: Herein, we present a literature review addressing the clinical evidence of CT and PET radiomics. An extensive literature review was conducted in PubMed, including papers on robustness and clinical applications. EVIDENCE SYNTHESIS: We summarize image-modality related influences on the robustness of radiomic features and provide an overview of clinical evidence reported in the literature. Today, more evidence has been provided for CT imaging, however, PET imaging offers the promise of direct imaging of biological processes and functions. We provide a summary of future research directions, which needs to be addressed in order to successfully introduce radiomics into clinical medicine. In comparison to CT, more focus should be directed towards harmonization of PET acquisition and reconstruction protocols, which is important for transferable modelling. CONCLUSIONS: Both CT and PET radiomics are promising pre-treatment and intra-treatment biomarkers for outcome prediction. Most studies are performed in retrospective setting, however their validation in prospective data collections is ongoing.

Influence of inter-observer delineation variability on radiomics stability in different tumor sites.

BACKGROUND: Radiomics is a promising methodology for quantitative analysis and description of radiological images using advanced mathematics and statistics. Tumor delineation, which is still often done manually, is an essential step in radiomics, however, inter-observer variability is a well-known uncertainty in radiation oncology. This study investigated the impact of inter-observer variability (IOV) in manual tumor delineation on the reliability of radiomic features (RF). METHODS: Three different tumor types (head and neck squamous cell carcinoma (HNSCC), malignant pleural mesothelioma (MPM) and non-small cell lung cancer (NSCLC)) were included. For each site, eleven individual tumors were contoured on CT scans by three experienced radiation oncologists. Dice coefficients (DC) were calculated for quantification of delineation variability. RF were calculated with an in-house developed software implementation, which comprises 1404 features: shape (n = 18), histogram (n = 17), texture (n = 137) and wavelet (n = 1232). The IOV of RF was studied using the intraclass correlation coefficient (ICC). An ICC >0.8 indicates a good reproducibility. For the stable RF, an average linkage hierarchical clustering was performed to identify classes of uncorrelated features. RESULTS: Median DC was high for NSCLC (0.86, range 0.57-0.90) and HNSCC (0.72, 0.21-0.89), whereas it was low for MPM (0.26, 0-0.9) indicating substantial IOV. Stability rate of RF correlated with DC and depended on tumor site, showing a high stability in NSCLC (90% of total parameters), acceptable stability in HNSCC (59% of total parameters) and low stability in MPM (36% of total parameters). Shape features showed the weakest stability across all tumor types. Hierarchical clustering revealed 14 groups of correlated and stable features for NSCLC and 6 groups for both HNSCC and MPM. CONCLUSION: Inter-observer delineation variability has a relevant influence on radiomics analysis and is strongly influenced by tumor type. This leads to a reduced number of suitable imaging features.

Comparison of PET and CT radiomics for prediction of local tumor control in head and neck squamous cell carcinoma.

PURPOSE: An association between radiomic features extracted from CT and local tumor control in the head and neck squamous cell carcinoma (HNSCC) has been shown. This study investigated the value of pretreatment functional imaging (18F-FDG PET) radiomics for modeling of local tumor control. MATERIAL AND METHODS: Data from HNSCC patients (n = 121) treated with definitive radiochemotherapy were used for model training. In total, 569 radiomic features were extracted from both contrast-enhanced CT and 18F-FDG PET images in the primary tumor region. CT, PET and combined PET/CT radiomic models to assess local tumor control were trained separately. Five feature selection and three classification methods were implemented. The performance of the models was quantified using concordance index (CI) in 5-fold cross validation in the training cohort. The best models, per image modality, were compared and verified in the independent validation cohort (n = 51). The difference in CI was investigated using bootstrapping. Additionally, the observed and radiomics-based estimated probabilities of local tumor control were compared between two risk groups. RESULTS: The feature selection using principal component analysis and the classification based on the multivariabale Cox regression with backward selection of the variables resulted in the best models for all image modalities (CICT = 0.72, CIPET = 0.74, CIPET/CT = 0.77). Tumors more homogenous in CT density (decreased GLSZMsize_zone_entropy) and with a focused region of high FDG uptake (higher GLSZMSZLGE) indicated better prognosis. No significant difference in the performance of the models in the validation cohort was observed (CICT = 0.73, CIPET = 0.71, CIPET/CT = 0.73). However, the CT radiomics-based model overestimated the probability of tumor control in the poor prognostic group (predicted = 68%, observed = 56%). CONCLUSIONS: Both CT and PET radiomics showed equally good discriminative power for local tumor control modeling in HNSCC. However, CT-based predictions overestimated the local control rate in the poor prognostic validation cohort, and thus, we recommend to base the local control modeling on the 18F-FDG PET.

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.

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.

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.

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.

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.].

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.

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.

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.

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. .

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.

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.

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.

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.