The Image Biomarker Standardization Initiative: Standardized Convolutional Filters for Reproducible Radiomics and Enhanced Clinical Insights.

Filters are commonly used to enhance specific structures and patterns in images, such as vessels or peritumoral regions, to enable clinical insights beyond the visible image using radiomics. However, their lack of standardization restricts reproducibility and clinical translation of radiomics decision support tools. In this special report, teams of researchers who developed radiomics software participated in a three-phase study (September 2020 to December 2022) to establish a standardized set of filters. The first two phases focused on finding reference filtered images and reference feature values for commonly used convolutional filters: mean, Laplacian of Gaussian, Laws and Gabor kernels, separable and nonseparable wavelets (including decomposed forms), and Riesz transformations. In the first phase, 15 teams used digital phantoms to establish 33 reference filtered images of 36 filter configurations. In phase 2, 11 teams used a chest CT image to derive reference values for 323 of 396 features computed from filtered images using 22 filter and image processing configurations. Reference filtered images and feature values for Riesz transformations were not established. Reproducibility of standardized convolutional filters was validated on a public data set of multimodal imaging (CT, fluorodeoxyglucose PET, and T1-weighted MRI) in 51 patients with soft-tissue sarcoma. At validation, reproducibility of 486 features computed from filtered images using nine configurations × three imaging modalities was assessed using the lower bounds of 95% CIs of intraclass correlation coefficients. Out of 486 features, 458 were found to be reproducible across nine teams with lower bounds of 95% CIs of intraclass correlation coefficients greater than 0.75. In conclusion, eight filter types were standardized with reference filtered images and reference feature values for verifying and calibrating radiomics software packages. A web-based tool is available for compliance checking.

Investigating the Robustness Neighborhood Gray Tone Difference Matrix and Gray Level Co-occurrence Matrix Radiomic Features on Clinical Computed Tomography Systems Using Anthropomorphic Phantoms: Evidence From a Multivendor Study.

OBJECTIVE: The aim of this study was to determine if optimized imaging protocols across multiple computed tomography (CT) vendors could result in reproducible radiomic features calculated from an anthropomorphic phantom. METHODS: Materials with varying degrees of heterogeneity were placed throughout the lungs of the phantom. Twenty scans of the phantom were acquired on 3 CT manufacturers with chest CT protocols that had optimized protocol parameters. Scans were reconstructed using vendor-specific standards and lung kernels. The concordance correlation coefficient (CCC) was used to calculate reproducibility between features. For features with high CCC values, Bland-Altman analysis was also used to quantify agreement. RESULTS: The mean Hounsfield unit (HU) was 32.93 HU (141.7 to -26.5 HU) for the rubber insert and 347.2 HU (-320.9 to -347.7 HU) for the wood insert. Low CCC values of less than 0.9 were calculated for all features across all scans. CONCLUSIONS: Radiomic features that are derived from the spatial distribution of voxel intensities should be particularly scrutinized for reproducibility in a multivendor environment.

Breast cancer molecular subtype classifier that incorporates MRI features.

PURPOSE: To use features extracted from magnetic resonance (MR) images and a machine-learning method to assist in differentiating breast cancer molecular subtypes. MATERIALS AND METHODS: This retrospective Health Insurance Portability and Accountability Act (HIPAA)-compliant study received Institutional Review Board (IRB) approval. We identified 178 breast cancer patients between 2006-2011 with: 1) ERPR + (n = 95, 53.4%), ERPR-/HER2 + (n = 35, 19.6%), or triple negative (TN, n = 48, 27.0%) invasive ductal carcinoma (IDC), and 2) preoperative breast MRI at 1.5T or 3.0T. Shape, texture, and histogram-based features were extracted from each tumor contoured on pre- and three postcontrast MR images using in-house software. Clinical and pathologic features were also collected. Machine-learning-based (support vector machines) models were used to identify significant imaging features and to build models that predict IDC subtype. Leave-one-out cross-validation (LOOCV) was used to avoid model overfitting. Statistical significance was determined using the Kruskal-Wallis test. RESULTS: Each support vector machine fit in the LOOCV process generated a model with varying features. Eleven out of the top 20 ranked features were significantly different between IDC subtypes with P < 0.05. When the top nine pathologic and imaging features were incorporated, the predictive model distinguished IDC subtypes with an overall accuracy on LOOCV of 83.4%. The combined pathologic and imaging model's accuracy for each subtype was 89.2% (ERPR+), 63.6% (ERPR-/HER2+), and 82.5% (TN). When only the top nine imaging features were incorporated, the predictive model distinguished IDC subtypes with an overall accuracy on LOOCV of 71.2%. The combined pathologic and imaging model's accuracy for each subtype was 69.9% (ERPR+), 62.9% (ERPR-/HER2+), and 81.0% (TN). CONCLUSION: We developed a machine-learning-based predictive model using features extracted from MRI that can distinguish IDC subtypes with significant predictive power. J. Magn. Reson. Imaging 2016;44:122-129.

Predictors of acute toxicities during definitive chemoradiation using intensity-modulated radiotherapy for anal squamous cell carcinoma.

PURPOSE: To identify clinical and dosimetric factors associated with acute hematologic and gastrointestinal (GI) toxicities during definitive therapy using intensity-modulated radiotherapy (IMRT) for anal squamous cell carcinoma (ASCC). MATERIALS AND METHODS: We retrospectively analyzed 108 ASCC patients treated with IMRT. Clinical information included age, gender, stage, concurrent chemotherapy, mitomycin (MMC) chemotherapy and weekly hematologic and GI toxicity during IMRT. From contours of the bony pelvis and bowel, dose-volume parameters were extracted. Logistic regression models were used to test associations between toxicities and clinical or dosimetric predictors. RESULTS: The median age was 59 years, 81 patients were women and 84 patients received concurrent MMC and 5-fluorouracil (5FU). On multivariate analysis (MVA), the model most predictive of Grade 2 + anemia included the maximum bony pelvis dose (Dmax), female gender, and T stage [p = 0.035, cross validation area under the curve (cvAUC) = 0.66]. The strongest model of Grade 2 + leukopenia included V10 (percentage of pelvic bone volume receiving ≥ 10 Gy) and number of MMC cycles (p = 0.276, cvAUC = 0.57). The model including MMC cycle number and T stage correlated best with Grade 2 + neutropenia (p = 0.306, cvAUC = 0.57). The model predictive of combined Grade 2 + hematologic toxicity (HT) included V10 and T stage (p = 0.016, cvAUC = 0.66). A model including VA45 (absolute bowel volume receiving ≥ 45 Gy) and MOH5 (mean dose to hottest 5% of bowel volume) best predicted diarrhea (p = 0.517, cvAUC = 0.56). CONCLUSION: Dosimetric constraints to the pelvic bones should be integrated into IMRT planning to reduce toxicity, potentially reducing treatment interruptions and improving disease outcomes in ASCC. Specifically, our results indicate that Dmax should be confined to ≤ 57 Gy to minimize anemia and that V10 should be restricted to ≤ 87% to reduce incidence of all HT.

Breast cancer subtype intertumor heterogeneity: MRI-based features predict results of a genomic assay.

PURPOSE: To investigate the association between a validated, gene-expression-based, aggressiveness assay, Oncotype Dx RS, and morphological and texture-based image features extracted from magnetic resonance imaging (MRI). MATERIALS AND METHODS: This retrospective study received Internal Review Board approval and need for informed consent was waived. Between 2006-2012, we identified breast cancer patients with: 1) ER+, PR+, and HER2- invasive ductal carcinoma (IDC); 2) preoperative breast MRI; and 3) Oncotype Dx RS test results. Extracted features included morphological, histogram, and gray-scale correlation matrix (GLCM)-based texture features computed from tumors contoured on pre- and three postcontrast MR images. Linear regression analysis was performed to investigate the association between Oncotype Dx RS and different clinical, pathologic, and imaging features. P < 0.05 was considered statistically significant. RESULTS: Ninety-five patients with IDC were included with a median Oncotype Dx RS of 16 (range: 0-45). Using stepwise multiple linear regression modeling, two MR-derived image features, kurtosis in the first and third postcontrast images and histologic nuclear grade, were found to be significantly correlated with the Oncotype Dx RS with P = 0.0056, 0.0005, and 0.0105, respectively. The overall model resulted in statistically significant correlation with Oncotype Dx RS with an R-squared value of 0.23 (adjusted R-squared = 0.20; P = 0.0002) and a Spearman's rank correlation coefficient of 0.49 (P < 0.0001). CONCLUSION: A model for IDC using imaging and pathology information correlates with Oncotype Dx RS scores, suggesting that image-based features could also predict the likelihood of recurrence and magnitude of chemotherapy benefit.

Clinical and Dosimetric Risk Factors Associated With Radiation-Induced Lung Toxicities After Multiple Courses of Lung Stereotactic Body Radiation Therapy.

Purpose: Data are limited on radiation-induced lung toxicities (RILT) after multiple courses of lung stereotactic body radiation therapy (SBRT). We herein analyze a large cohort of patients to explore the clinical and dosimetric risk factors associated with RILT in such settings. Methods and Materials: A single institutional database of patients treated with multiple courses of lung SBRT between January 2014 and December 2019 was analyzed. Grade 2 or higher (G2+) RILT after the last course of SBRT was the primary endpoint. Composite plans were generated with advanced algorithms including deformable registration and equivalent dose adjustment. Logistic regression analyses were performed to examine correlations between patient or treatment factors including dosimetry and G2+ RILT. Risk stratification of patients and lung constraints based on acceptable normal tissue complication probability were calculated based on risk factors identified. Results: Among 110 eligible patients (56 female and 54 male), there were 64 synchronous (58.2%; defined as 2 courses of SBRT delivered within 30 days) and 46 metachronous (41.8%) courses of SBRT. The composite median lung V20, lung V5, and mean lung dose were 9.9% (interquartile range [IQR], 7.3%-12.4%), 32.2% (IQR, 25.5%-40.1%), and 7.0 Gy (IQR, 5.5 Gy-8.6 Gy), respectively. With a median follow-up of 21.1 months, 30 patients (27.3%) experienced G2+ RILT. Five patients (4.5%) developed G3 RILT, and 1 patient (0.9%) developed G4 RILT, and no patients developed G5 RILT. On multivariable regression analysis, female sex (odds ratio [OR], 4.35; 95% CI, 1.49%-14.3%; P = .01), synchronous SBRT (OR, 8.78; 95% CI, 2.27%-47.8%; P = .004), prior G2+ RILT (OR, 29.8; 95% CI, 2.93%-437%; P = .007) and higher composite lung V20 (OR, 1.18; 95% CI, 1.02%-1.38%; P = .030) were associated with significantly higher likelihood of G2+ RILT. Conclusions: Our data suggest an acceptable incidence of G2+ RILT after multiple courses of lung SBRT. Female sex, synchronous SBRT, prior G2+ RILT, and higher composite lung V20 may be risk factors for G2+ RILT.

Interpretable Machine Learning for Choosing Radiation Dose-volume Constraints on Cardio-pulmonary Substructures Associated with Overall Survival in NRG Oncology RTOG 0617.

PURPOSE: Our objective was to use interpretable machine learning for choosing dose-volume constraints on cardiopulmonary substructures (CPSs) associated with overall survival (OS) in radiation therapy for locally advanced non-small cell lung cancer. METHODS AND MATERIALS: A total of 428 patients with non-small cell lung cancer were randomly divided into training/validation/test subsets (n = 230/149/49) in Radiation Therapy Oncology Group 0617. Manual or automated contouring was performed to segment CPSs, including heart, atria, ventricles, aorta, left/right ventricle/atrium (LV+RV+LA+RA), inferior/superior vena cava, pulmonary artery, and pericardium. Peri (pericardium-heart), rest (heart-[LV+RV+LA+RA]), clinical target volume (CTV), and lungs-CTV contours were also obtained. Dose-volume histogram features were extracted, including minimum/mean dose to the hottest x% volume (Dx%[Gy]/MOHx%[Gy]), minimum/mean/maximum dose, percent volume receiving at least xGy (VxGy[%]), and overlapping volume of each CPS with planning target volume (PTV_Voverlap[%]). Clinical parameters were collected from the National Clinical Trials Network/Community oncology research program data archive. Feature selection was performed using a series of multiblock sparse partial least squares regression, stability selection supervised principal component analysis, and Boruta. Explainable boosting machine (EBM) was trained using a conditional survival distribution-based approach for imputing censored data, treating survival analysis as a regression problem. Harrell's C-index was used to evaluate OS discrimination performance of EBM, Cox proportional hazards (CPH), random survival forest, extreme gradient boosting survival embeddings, and CPH deep neural network (DeepSurv) models in the test set. Dose-volume constraints were selected using the binary change point detection algorithm in Shapley additive explanations-based partial dependence functions. RESULTS: Selected features included LA_V60Gy(%), pericardium_D30%(Gy), lungs-CTV_PTV_Voverlap(%), RA_V55Gy(%), and received_cons_chemo. All models ranked LA_V60Gy(%) as the most important feature. EBM achieved the best performance for predicting OS, followed by extreme gradient boosting survival embeddings, random survival forest, DeepSurv, and CPH (C-index = 0.653, 0.646, 0.642, 0.638, and 0.632). EBM global explanations suggested that LA_V60Gy(%) < 25.6, lungs-CTV_PTV_Voverlap(%) < 1.1, pericardium_D30%(Gy) < 18.9, RA_V55Gy(%) < 19.5, and received_cons_chemo = 'Yes' for improved OS. CONCLUSIONS: EBM can be used to discriminate OS while also guiding dose-volume constraint selection for optimal management of cardiac toxicity in lung cancer radiation therapy.

ROE (Radiotherapy Outcomes Estimator): An open-source tool for optimizing radiotherapy prescriptions.

BACKGROUND AND OBJECTIVES: Radiotherapy prescriptions currently derive from population-wide guidelines established through large clinical trials. We provide an open-source software tool for patient-specific prescription determination using personalized dose-response curves. METHODS: We developed ROE, a plugin to the Computational Environment for Radiotherapy Research to visualize predicted tumor control and normal tissue complication simultaneously, as a function of prescription dose. ROE can be used natively with MATLAB and is additionally made accessible in GNU Octave and Python, eliminating the need for commercial licenses. It provides a curated library of published and validated predictive models and incorporates clinical restrictions on normal tissue outcomes. ROE additionally provides batch-mode tools to evaluate and select among different fractionation schemes and analyze radiotherapy outcomes across patient cohorts. CONCLUSION: ROE is an open-source, GPL-copyrighted tool for interactive exploration of the dose-response relationship to aid in radiotherapy planning. We demonstrate its potential clinical relevance in (1) improving patient awareness by quantifying the risks and benefits of a given treatment protocol (2) assessing the potential for dose escalation across patient cohorts and (3) estimating accrual rates of new protocols.

CT and cone-beam CT of ablative radiation therapy for pancreatic cancer with expert organ-at-risk contours.

We describe a dataset from patients who received ablative radiation therapy for locally advanced pancreatic cancer (LAPC), consisting of computed tomography (CT) and cone-beam CT (CBCT) images with physician-drawn organ-at-risk (OAR) contours. The image datasets (one CT for treatment planning and two CBCT scans at the time of treatment per patient) were collected from 40 patients. All scans were acquired with the patient in the treatment position and in a deep inspiration breath-hold state. Six radiation oncologists delineated the gastrointestinal OARs consisting of small bowel, stomach and duodenum, such that the same physician delineated all image sets belonging to the same patient. Two trained medical physicists further edited the contours to ensure adherence to delineation guidelines. The image and contour files are available in DICOM format and are publicly available from The Cancer Imaging Archive ( https://doi.org/10.7937/TCIA.ESHQ-4D90 , Version 2). The dataset can serve as a criterion standard for evaluating the accuracy and reliability of deformable image registration and auto-segmentation algorithms, as well as a training set for deep-learning-based methods.

Prospectively-validated deep learning model for segmenting swallowing and chewing structures in CT.

Objective.Delineating swallowing and chewing structures aids in radiotherapy (RT) treatment planning to limit dysphagia, trismus, and speech dysfunction. We aim to develop an accurate and efficient method to automate this process.Approach.CT scans of 242 head and neck (H&N) cancer patients acquired from 2004 to 2009 at our institution were used to develop auto-segmentation models for the masseters, medial pterygoids, larynx, and pharyngeal constrictor muscle using DeepLabV3+. A cascaded framework was used, wherein models were trained sequentially to spatially constrain each structure group based on prior segmentations. Additionally, an ensemble of models, combining contextual information from axial, coronal, and sagittal views was used to improve segmentation accuracy. Prospective evaluation was conducted by measuring the amount of manual editing required in 91 H&N CT scans acquired February-May 2021.Main results. Medians and inter-quartile ranges of Dice similarity coefficients (DSC) computed on the retrospective testing set (N = 24) were 0.87 (0.85-0.89) for the masseters, 0.80 (0.79-0.81) for the medial pterygoids, 0.81 (0.79-0.84) for the larynx, and 0.69 (0.67-0.71) for the constrictor. Auto-segmentations, when compared to two sets of manual segmentations in 10 randomly selected scans, showed better agreement (DSC) with each observer than inter-observer DSC. Prospective analysis showed most manual modifications needed for clinical use were minor, suggesting auto-contouring could increase clinical efficiency. Trained segmentation models are available for research use upon request viahttps://github.com/cerr/CERR/wiki/Auto-Segmentation-models.Significance.We developed deep learning-based auto-segmentation models for swallowing and chewing structures in CT and demonstrated its potential for use in treatment planning to limit complications post-RT. To the best of our knowledge, this is the only prospectively-validated deep learning-based model for segmenting chewing and swallowing structures in CT. Segmentation models have been made open-source to facilitate reproducibility and multi-institutional research.

An Incognizant Glitch-pinnacle Polyliner Dissociation: A Case Report.

Introduction: Liner dissociation of the pinnacle total hip arthroplasty system is a rare but documented complication. Although a few reports are published internationally, to the best of our knowledge no cases have been documented from India so far. Case Report: A 31-year-old male presented with failed femoral head fracture fixation for which total hip replacement was done. Postoperatively at 18 months, he was diagnosed with pinnacle liner dissociation and liner exchange was performed. Conclusion: This report aims to raise awareness about the incidence of pinnacle liner dissociation.

Reproducibility of radiomic features using network analysis and its application in Wasserstein k-means clustering.

Purpose: The goal of this study is to develop innovative methods for identifying radiomic features that are reproducible over varying image acquisition settings. Approach: We propose a regularized partial correlation network to identify reliable and reproducible radiomic features. This approach was tested on two radiomic feature sets generated using two different reconstruction methods on computed tomography (CT) scans from a cohort of 47 lung cancer patients. The largest common network component between the two networks was tested on phantom data consisting of five cancer samples. To further investigate whether radiomic features found can identify phenotypes, we propose a k -means clustering algorithm coupled with the optimal mass transport theory. This approach following the regularized partial correlation network analysis was tested on CT scans from 77 head and neck squamous cell carcinoma (HNSCC) patients in the Cancer Imaging Archive (TCIA) and validated using an independent dataset. Results: A set of common radiomic features was found in relatively large network components between the resultant two partial correlation networks resulting from a cohort of lung cancer patients. The reliability and reproducibility of those radiomic features were further validated on phantom data using the Wasserstein distance. Further analysis using the network-based Wasserstein k -means algorithm on the TCIA HNSCC data showed that the resulting clusters separate tumor subsites as well as HPV status, and this was validated on an independent dataset. Conclusion: We showed that a network-based analysis enables identifying reproducible radiomic features and use of the selected set of features can enhance clustering results.

Registering Study Analysis Plans (SAPs) Before Dissecting Your Data-Updating and Standardizing Outcome Modeling.

Public preregistration of study analysis plans (SAPs) is widely recognized for clinical trials, but adopted to a much lesser extent in observational studies. Registration of SAPs prior to analysis is encouraged to not only increase transparency and exactness but also to avoid positive finding bias and better standardize outcome modeling. Efforts to generally standardize outcome modeling, which can be based on clinical trial and/or observational data, have recently spurred. We suggest a three-step SAP concept in which investigators are encouraged to (1) Design the SAP and circulate it among the co-investigators, (2) Log the SAP with a public repository, which recognizes the SAP with a digital object identifier (DOI), and (3) Cite (using the DOI), briefly summarize and motivate any deviations from the SAP in the associated manuscript. More specifically, the SAP should include the scope (brief data and study description, co-investigators, hypotheses, primary outcome measure, study title), in addition to step-by-step details of the analysis (handling of missing data, resampling, defined significance level, statistical function, validation, and variables and parameterization).

A novel kernel Wasserstein distance on Gaussian measures: An application of identifying dental artifacts in head and neck computed tomography.

The Wasserstein distance is a powerful metric based on the theory of optimal mass transport. It gives a natural measure of the distance between two distributions with a wide range of applications. In contrast to a number of the common divergences on distributions such as Kullback-Leibler or Jensen-Shannon, it is (weakly) continuous, and thus ideal for analyzing corrupted and noisy data. Until recently, however, no kernel methods for dealing with nonlinear data have been proposed via the Wasserstein distance. In this work, we develop a novel method to compute the L2-Wasserstein distance in reproducing kernel Hilbert spaces (RKHS) called kernel L2-Wasserstein distance, which is implemented using the kernel trick. The latter is a general method in machine learning employed to handle data in a nonlinear manner. We evaluate the proposed approach in identifying computed tomography (CT) slices with dental artifacts in head and neck cancer, performing unsupervised hierarchical clustering on the resulting Wasserstein distance matrix that is computed on imaging texture features extracted from each CT slice. We further compare the performance of kernel Wasserstein distance with alternatives including kernel Kullback-Leibler divergence we previously developed. Our experiments show that the kernel approach outperforms classical non-kernel approaches in identifying CT slices with artifacts.

Radiomic analysis identifies tumor subtypes associated with distinct molecular and microenvironmental factors in head and neck squamous cell carcinoma.

PURPOSE: To identify whether radiomic features from pre-treatment computed tomography (CT) scans can predict molecular differences between head and neck squamous cell carcinoma (HNSCC) using The Cancer Imaging Archive (TCIA) and The Cancer Genome Atlas (TCGA). METHODS: 77 patients from the TCIA with HNSCC had imaging suitable for analysis. Radiomic features were extracted and unsupervised consensus clustering was performed to identify subtypes. Genomic data was extracted from the matched patients in the TCGA database. We explored relationships between radiomic features and molecular profiles of tumors, including the tumor immune microenvironment. A machine learning method was used to build a model predictive of CD8 + T-cells. An independent cohort of 83 HNSCC patients was used to validate the radiomic clusters. RESULTS: We initially extracted 104 two-dimensional radiomic features, and after feature stability tests and removal of volume dependent features, reduced this to 67 features for subsequent analysis. Consensus clustering based on these features resulted in two distinct clusters. The radiomic clusters differed by primary tumor subsite (p = 0.0096), HPV status (p = 0.0127), methylation-based clustering results (p = 0.0025), and tumor immune microenvironment. A random forest model using radiomic features predicted CD8 + T-cells independent of HPV status with R2 = 0.30 (p < 0.0001) on cross validation. Consensus clustering on the validation cohort resulted in two distinct clusters that differ in tumor subsite (p = 1.3 × 10-7) and HPV status (p = 4.0 × 10-7). CONCLUSION: Radiomic analysis can identify biologic features of tumors such as HPV status and T-cell infiltration and may be able to provide other information in the near future to help with patient stratification.

Predictive Modeling of Thoracic Radiotherapy Toxicity and the Potential Role of Serum Alpha-2-Macroglobulin.

Background: To investigate the impact of alpha-2-macroglobulin (A2M), a suspected intrinsic radioprotectant, on radiation pneumonitis and esophagitis using multifactorial predictive models. Materials and Methods: Baseline A2M levels were obtained for 258 patients prior to thoracic radiotherapy (RT). Dose-volume characteristics were extracted from treatment plans. Spearman's correlation (Rs) test was used to correlate clinical and dosimetric variables with toxicities. Toxicity prediction models were built using least absolute shrinkage and selection operator (LASSO) logistic regression on 1,000 bootstrapped datasets. Results: Grade ≥2 esophagitis and pneumonitis developed in 61 (23.6%) and 36 (14.0%) patients, respectively. The median A2M level was 191 mg/dL (range: 94-511). Never/former/current smoker status was 47 (18.2%)/179 (69.4%)/32 (12.4%). We found a significant negative univariate correlation between baseline A2M levels and esophagitis (Rs = -0.18/p = 0.003) and between A2M and smoking status (Rs = 0.13/p = 0.04). Further significant parameters for grade ≥2 esophagitis included age (Rs = -0.32/p < 0.0001), chemotherapy use (Rs = 0.56/p < 0.0001), dose per fraction (Rs = -0.57/p < 0.0001), total dose (Rs = 0.35/p < 0.0001), and several other dosimetric variables with Rs > 0.5 (p < 0.0001). The only significant non-dosimetric parameter for grade ≥2 pneumonitis was sex (Rs = -0.32/p = 0.037) with higher risk for women. For pneumonitis D15 (lung) (Rs = 0.19/p = 0.006) and D45 (heart) (Rs = 0.16/p = 0.016) had the highest correlation. LASSO models applied on the validation data were statistically significant and resulted in areas under the receiver operating characteristic curve of 0.84 (esophagitis) and 0.78 (pneumonitis). Multivariate predictive models did not require A2M to reach maximum predictive power. Conclusion: This is the first study showing a likely association of higher baseline A2M values with lower risk of radiation esophagitis and with smoking status. However, the baseline A2M level was not a significant risk factor for radiation pneumonitis.

Library of deep-learning image segmentation and outcomes model-implementations.

An open-source library of implementations for deep-learning-based image segmentation and outcomes models based on radiotherapy and radiomics is presented. As oncology treatment planning becomes increasingly driven by automation, such a library of model implementations is crucial to (i) validate existing models on datasets collected at different institutions, (ii) automate segmentation, (iii) create ensembles for improving performance and (iv) incorporate validated models in the clinical workflow. Inclusion of deep-learning-based image segmentation and outcomes models in the same library provides a fully automated and reproduceable pipeline to estimate prognosis. The library was developed with the Computational Environment for Radiological Research (CERR) software platform. Centralizing model implementations in CERR builds upon its rich set of radiotherapy and radiomics tools and caters to the world-wide user base. CERR provides well-validated feature extraction pipelines for radiotherapy dosimetry and radiomics with fine control over the calculation settings, allowing users to select appropriate parameters used in model derivation. Models for automatic image segmentation are distributed via containers, allowing them to be deployed with a variety of scientific computing architectures. The library includes implementations of popular DVH-based models outlined in the Quantitative Analysis of Normal Tissue Effects in the Clinic effort and recently published literature. Radiomics models include features from the Image Biomarker Standardization Initiative and application-specific features found to be relevant across multiple sites and image modalities. The library is distributed as a module within CERR at https://www.github.com/cerr/CERR under the GNU-GPL copyleft with additional restrictions on clinical and commercial use and provision to dual license in future.

MRI features predictive of negative surgical margins in patients with HER2 overexpressing breast cancer undergoing breast conservation.

Here we develop a tool to predict resectability of HER2+ breast cancer at breast conservation surgery (BCS) utilizing features identified on preoperative breast MRI. We identified patients with HER2+ breast cancer who obtained pre-operative breast MRI and underwent BCS between 2002-2013. From the contoured tumor on pre-operative MRI, shape, histogram, and co-occurrence and size zone matrix texture features were extracted. In univariate analysis, Spearman's correlation coefficient (Rs) was used to assess the correlation between each image feature and an endpoint (surgical re-excision). For multivariate modeling, we employed a support vector machine (SVM) method in a manner of leave-one-out cross-validation (LOOCV). Of 109 patients with HER2+breast cancer who underwent BCS, 39% underwent surgical re-excision. 62% had residual cancer at re-excision. In univariate analysis, solidity (Rs = -0.32, p = 0.009) and extent (Rs = -0.29, p = 0.019) were significantly associated with re-excision. Skewness in post-contrast 1, 2, and 3 (Rs = 0.25, p = 0.045; Rs = 0.30, p = 0.015; Rs = 0.28, p = 0.026) and kurtosis in post-contrast 1 (Rs = 0.26, p = 0.035) were also statistically significant. LOOCV-based SVM test achieved 74.4% specificity and 71.4% sensitivity when 21 features were used. Thus, tumor texture, histogram and morphological MRI features may assist surgical planning, encouraging wide margins or mastectomy in patients who may otherwise go on to re-excision.

Image-guided radiotherapy reduces the risk of under-dosing high-risk prostate cancer extra-capsular disease and improves biochemical control.

BACKGROUND: To determine if reduced dose delivery uncertainty is associated with daily image-guidance (IG) and Prostate Specific Antigen Relapse Free Survival (PRFS) in intensity-modulated radiotherapy (IMRT) of high-risk prostate cancer (PCa). METHODS: Planning data for consecutive PCa patients treated with IMRT (n = 67) and IG-IMRT (n = 35) was retrieved. Using computer simulations of setup errors, we estimated the patient-specific uncertainty in accumulated treatment dose distributions for the prostate and for posterolateral aspects of the gland that are at highest risk for extra-capsular disease. Multivariate Cox regression for PRFS considering Gleason score, T-stage, pre-treatment PSA, number of elevated clinical risk factors (T2c+, GS7+ and PSA10+), nomogram-predicted risk of extra-capsular disease (ECD), and dose metrics was performed. RESULTS: For IMRT vs. IG-IMRT, plan dosimetry values were similar, but simulations revealed uncertainty in delivered dose external to the prostate was significantly different, due to positioning uncertainties. A patient-specific interaction term of the risk of ECD and risk of low dose to the ECD (p = 0.005), and the number of elevated clinical risk factors (p = 0.008), correlate with reduced PRFS. CONCLUSIONS: Improvements in PSA outcomes for high-risk PCa using IG-IMRT vs. IMRT without IG may be due to improved dosimetry for ECD.

Technical Note: Extension of CERR for computational radiomics: A comprehensive MATLAB platform for reproducible radiomics research.

PURPOSE: Radiomics is a growing field of image quantitation, but it lacks stable and high-quality software systems. We extended the capabilities of the Computational Environment for Radiological Research (CERR) to create a comprehensive, open-source, MATLAB-based software platform with an emphasis on reproducibility, speed, and clinical integration of radiomics research. METHOD: The radiomics tools in CERR were designed specifically to quantitate medical images in combination with CERR's core functionalities of radiological data import, transformation, management, image segmentation, and visualization. CERR allows for batch calculation and visualization of radiomics features, and provides a user-friendly data structure for radiomics metadata. All radiomics computations are vectorized for speed. Additionally, a test suite is provided for reconstruction and comparison with radiomics features computed using other software platforms such as the Insight Toolkit (ITK) and PyRadiomics. CERR was evaluated according to the standards defined by the Image Biomarker Standardization Initiative. CERR's radiomics feature calculation was integrated with the clinically used MIM software using its MATLAB® application programming interface. RESULTS: The CERR provides a comprehensive computational platform for radiomics analysis. Matrix formulations for the compute-intensive Haralick texture resulted in speeds that are superior to the implementation in ITK 4.12. For an image discretized into 32 bins, CERR achieved a speedup of 3.5 times over ITK. The CERR test suite enabled the successful identification of programming errors as well as genuine differences in radiomics definitions and calculations across the software packages tested. CONCLUSION: The CERR's radiomics capabilities are comprehensive, open-source, and fast, making it an attractive platform for developing and exploring radiomics signatures across institutions. The ability to both choose from a wide variety of radiomics implementations and to integrate with a clinical workflow makes CERR useful for retrospective as well as prospective research analyses.