A common [18F]-FDG PET radiomic signature to predict survival in patients with HPV-induced cancers.

Locally advanced cervical cancer (LACC) and anal and oropharyngeal squamous cell carcinoma (ASCC and OPSCC) are mostly caused by oncogenic human papillomaviruses (HPV). In this paper, we developed machine learning (ML) models based on clinical, biological, and radiomic features extracted from pre-treatment fluorine-18-fluorodeoxyglucose positron emission tomography ([18F]-FDG PET) images to predict the survival of patients with HPV-induced cancers. For this purpose, cohorts from five institutions were used: two cohorts of patients treated for LACC including 104 patients from Gustave Roussy Campus Cancer (Center 1) and 90 patients from Leeds Teaching Hospitals NHS Trust (Center 2), two datasets of patients treated for ASCC composed of 66 patients from Institut du Cancer de Montpellier (Center 3) and 67 patients from Oslo University Hospital (Center 4), and one dataset of 45 OPSCC patients from the University Hospital of Zurich (Center 5). Radiomic features were extracted from baseline [18F]-FDG PET images. The ComBat technique was applied to mitigate intra-scanner variability. A modified consensus nested cross-validation for feature selection and hyperparameter tuning was applied on four ML models to predict progression-free survival (PFS) and overall survival (OS) using harmonized imaging features and/or clinical and biological variables as inputs. Each model was trained and optimized on Center 1 and Center 3 cohorts and tested on Center 2, Center 4, and Center 5 cohorts. The radiomic-based CoxNet model achieved C-index values of 0.75 and 0.78 for PFS and 0.76, 0.74, and 0.75 for OS on the test sets. Radiomic feature-based models had superior performance compared to the bioclinical ones, and combining radiomic and bioclinical variables did not improve the performances. Metabolic tumor volume (MTV)-based models obtained lower C-index values for a majority of the tested configurations but quite equivalent performance in terms of time-dependent AUCs (td-AUC). The results demonstrate the possibility of identifying common PET-based image signatures for predicting the response of patients with induced HPV pathology, validated on multi-center multiconstructor data.

Reinventing radiation therapy with machine learning and imaging bio-markers (radiomics): State-of-the-art, challenges and perspectives.

Radiation therapy is a pivotal cancer treatment that has significantly progressed over the last decade due to numerous technological breakthroughs. Imaging is now playing a critical role on deployment of the clinical workflow, both for treatment planning and treatment delivery. Machine-learning analysis of predefined features extracted from medical images, i.e. radiomics, has emerged as a promising clinical tool for a wide range of clinical problems addressing drug development, clinical diagnosis, treatment selection and implementation as well as prognosis. Radiomics denotes a paradigm shift redefining medical images as a quantitative asset for data-driven precision medicine. The adoption of machine-learning in a clinical setting and in particular of radiomics features requires the selection of robust, representative and clinically interpretable biomarkers that are properly evaluated on a representative clinical data set. To be clinically relevant, radiomics must not only improve patients' management with great accuracy but also be reproducible and generalizable. Hence, this review explores the existing literature and exposes its potential technical caveats, such as the lack of quality control, standardization, sufficient sample size, type of data collection, and external validation. Based upon the analysis of 165 original research studies based on PET, CT-scan, and MRI, this review provides an overview of new concepts, and hypotheses generating findings that should be validated. In particular, it describes evolving research trends to enhance several clinical tasks such as prognostication, treatment planning, response assessment, prediction of recurrence/relapse, and prediction of toxicity. Perspectives regarding the implementation of an AI-based radiotherapy workflow are presented.

Elucidation of roles for vitamin B12 in regulation of folate, ubiquinone, and methionine metabolism.

Only a small fraction of vitamin B12-requiring organisms are able to synthesize B12 de novo, making it a common commodity in microbial communities. Initially recognized as an enzyme cofactor of a few enzymes, recent studies have revealed additional B12-binding enzymes and regulatory roles for B12 Here we report the development and use of a B12-based chemical probe to identify B12-binding proteins in a nonphototrophic B12-producing bacterium. Two unexpected discoveries resulted from this study. First, we identified a light-sensing B12-binding transcriptional regulator and demonstrated that it controls folate and ubiquinone biosynthesis. Second, our probe captured proteins involved in folate, methionine, and ubiquinone metabolism, suggesting that it may play a role as an allosteric effector of these processes. These metabolic processes produce precursors for synthesis of DNA, RNA, and protein. Thereby, B12 likely modulates growth, and by limiting its availability to auxotrophs, B12-producing organisms may facilitate coordination of community metabolism.

Automatic gross tumor volume segmentation with failure detection for safe implementation in locally advanced cervical cancer.

Background and Purpose: Automatic segmentation methods have greatly changed the RadioTherapy (RT) workflow, but still need to be extended to target volumes. In this paper, Deep Learning (DL) models were compared for Gross Tumor Volume (GTV) segmentation in locally advanced cervical cancer, and a novel investigation into failure detection was introduced by utilizing radiomic features. Methods and materials: We trained eight DL models (UNet, VNet, SegResNet, SegResNetVAE) for 2D and 3D segmentation. Ensembling individually trained models during cross-validation generated the final segmentation. To detect failures, binary classifiers were trained using radiomic features extracted from segmented GTVs as inputs, aiming to classify contours based on whether their Dice Similarity Coefficient ( DSC ) < T and DSC ⩾ T . Two distinct cohorts of T2-Weighted (T2W) pre-RT MR images captured in 2D sequences were used: one retrospective cohort consisting of 115 LACC patients from 30 scanners, and the other prospective cohort, comprising 51 patients from 7 scanners, used for testing. Results: Segmentation by 2D-SegResNet achieved the best DSC, Surface DSC ( SDSC 3 mm ), and 95th Hausdorff Distance (95HD): DSC = 0.72 ± 0.16, SDSC 3 mm =0.66 ± 0.17, and 95HD = 14.6 ± 9.0 mm without missing segmentation ( M =0) on the test cohort. Failure detection could generate precision ( P = 0.88 ), recall ( R = 0.75 ), F1-score ( F = 0.81 ), and accuracy ( A = 0.86 ) using Logistic Regression (LR) classifier on the test cohort with a threshold T = 0.67 on DSC values. Conclusions: Our study revealed that segmentation accuracy varies slightly among different DL methods, with 2D networks outperforming 3D networks in 2D MRI sequences. Doctors found the time-saving aspect advantageous. The proposed failure detection could guide doctors in sensitive cases.

Deep-Learning for Rapid Estimation of the Out-of-Field Dose in External Beam Photon Radiation Therapy - A Proof of Concept.

PURPOSE: The dose deposited outside of the treatment field during external photon beam radiation therapy treatment, also known as out-of-field dose, is the subject of extensive study as it may be associated with a higher risk of developing a second cancer and could have deleterious effects on the immune system that compromise the efficiency of combined radio-immunotherapy treatments. Out-of-field dose estimation tools developed today in research, including Monte Carlo simulations and analytical methods, are not suited to the requirements of clinical implementation because of their lack of versatility and their cumbersome application. We propose a proof of concept based on deep learning for out-of-field dose map estimation that addresses these limitations. METHODS AND MATERIALS: For this purpose, a 3D U-Net, considering as inputs the in-field dose, as computed by the treatment planning system, and the patient's anatomy, was trained to predict out-of-field dose maps. The cohort used for learning and performance evaluation included 3151 pediatric patients from the FCCSS database, treated in 5 clinical centers, whose whole-body dose maps were previously estimated with an empirical analytical method. The test set, composed of 433 patients, was split into 5 subdata sets, each containing patients treated with devices unseen during the training phase. Root mean square deviation evaluated only on nonzero voxels located in the out-of-field areas was computed as performance metric. RESULTS: Root mean square deviations of 0.28 and 0.41 cGy/Gy were obtained for the training and validation data sets, respectively. Values of 0.27, 0.26, 0.28, 0.30, and 0.45 cGy/Gy were achieved for the 6 MV linear accelerator, 16 MV linear accelerator, Alcyon cobalt irradiator, Mobiletron cobalt irradiator, and betatron device test sets, respectively. CONCLUSIONS: This proof-of-concept approach using a convolutional neural network has demonstrated unprecedented generalizability for this task, although it remains limited, and brings us closer to an implementation compatible with clinical routine.

AutoComBat: a generic method for harmonizing MRI-based radiomic features.

The use of multicentric data is becoming essential for developing generalizable radiomic signatures. In particular, Magnetic Resonance Imaging (MRI) data used in brain oncology are often heterogeneous in terms of scanners and acquisitions, which significantly impact quantitative radiomic features. Various methods have been proposed to decrease dependency, including methods acting directly on MR images, i.e., based on the application of several preprocessing steps before feature extraction or the ComBat method, which harmonizes radiomic features themselves. The ComBat method used for radiomics may be misleading and presents some limitations, such as the need to know the labels associated with the "batch effect". In addition, a statistically representative sample is required and the applicability of a signature whose batch label is not present in the train set is not possible. This work aimed to compare a priori and a posteriori radiomic harmonization methods and propose a code adaptation to be machine learning compatible. Furthermore, we have developed AutoComBat, which aims to automatically determine the batch labels, using either MRI metadata or quality metrics as inputs of the proposed constrained clustering. A heterogeneous dataset consisting of high and low-grade gliomas coming from eight different centers was considered. The different methods were compared based on their ability to decrease relative standard deviation of radiomic features extracted from white matter and on their performance on a classification task using different machine learning models. ComBat and AutoComBat using image-derived quality metrics as inputs for batch assignment and preprocessing methods presented promising results on white matter harmonization, but with no clear consensus for all MR images. Preprocessing showed the best results on the T1w-gd images for the grading task. For T2w-flair, AutoComBat, using either metadata plus quality metrics or metadata alone as inputs, performs better than the conventional ComBat, highlighting its potential for data harmonization. Our results are MRI weighting, feature class and task dependent and require further investigations on other datasets.

Imaging approaches and radiomics: toward a new era of ultraprecision radioimmunotherapy?

Strong rationale and a growing number of preclinical and clinical studies support combining radiotherapy and immunotherapy to improve patient outcomes. However, several critical questions remain, such as the identification of patients who will benefit from immunotherapy and the identification of the best modalities of treatment to optimize patient response. Imaging biomarkers and radiomics have recently emerged as promising tools for the non-invasive assessment of the whole disease of the patient, allowing comprehensive analysis of the tumor microenvironment, the spatial heterogeneity of the disease and its temporal changes. This review presents the potential applications of medical imaging and the challenges to address, in order to help clinicians choose the optimal modalities of both radiotherapy and immunotherapy, to predict patient's outcomes and to assess response to these promising combinations.

COMBING: Clustering in Oncology for Mathematical and Biological Identification of Novel Gene Signatures.

Precision medicine is a paradigm shift in healthcare relying heavily on genomics data. However, the complexity of biological interactions, the large number of genes as well as the lack of comparisons on the analysis of data, remain a tremendous bottleneck regarding clinical adoption. In this paper, we introduce a novel, automatic and unsupervised framework to discover low-dimensional gene biomarkers. Our method is based on the LP-Stability algorithm, a high dimensional center-based unsupervised clustering algorithm. It offers modularity as concerns metric functions and scalability, while being able to automatically determine the best number of clusters. Our evaluation includes both mathematical and biological criteria to define a quantitative metric. The recovered signature is applied to a variety of biological tasks, including screening of biological pathways and functions, and characterization relevance on tumor types and subtypes. Quantitative comparisons among different distance metrics, commonly used clustering methods and a referential gene signature used in the literature, confirm state of the art performance of our approach. In particular, our signature, based on 27 genes, reports at least 30 times better mathematical significance (average Dunn's Index) and 25% better biological significance (average Enrichment in Protein-Protein Interaction) than those produced by other referential clustering methods. Finally, our signature reports promising results on distinguishing immune inflammatory and immune desert tumors, while reporting a high balanced accuracy of 92% on tumor types classification and averaged balanced accuracy of 68% on tumor subtypes classification, which represents, respectively 7% and 9% higher performance compared to the referential signature.

Radiomics to evaluate interlesion heterogeneity and to predict lesion response and patient outcomes using a validated signature of CD8 cells in advanced melanoma patients treated with anti-PD1 immunotherapy.

PURPOSE: While there is still a significant need to identify potential biomarkers that can predict which patients are most likely to respond to immunotherapy treatments, radiomic approaches have shown promising results. The objectives of this study were to evaluate whether a previously validated radiomics signature of CD8 T-cells could predict progressions at a lesion level and whether the spatial heterogeneity of this radiomics score could be used at a patient level to assess the clinical response and survival of melanoma patients. METHODS: Clinical data from patients with advanced melanoma treated in our center with immunotherapy were retrieved. Radiomic features were extracted and the CD8 radiomics signature was applied. A progressive lesion was defined by an increase in lesion size of 20% or more. Dispersion metrics of the radiomics signature were estimated to evaluate the impact of interlesion heterogeneity on patient's response. Fine-tuned cut-offs for predicting overall survival were evaluated after splitting data into training and test sets. RESULTS: A total of 136 patients were included in this study, with 1120 segmented lesions at baseline, and 1052 lesions at first evaluation. A low CD8 radiomics score at baseline was associated with a significantly higher risk of lesion progression (AUC=0.55, p=0.0091), especially for lesions larger than >1 mL (AUC=0.59 overall, p=0.0035, with AUC=0.75, p=0.002 for subcutaneous lesions, AUC=0.68, p=0.01, for liver lesions and AUC=0.62, p=0.03 for nodes). The least infiltrated lesion according to the radiomics score of CD8 T-cells was positively associated with overall survival (training set HR=0.31, p=0.00062, test set HR=0.28, p=0.016), which remained significant in a multivariate analysis including clinical and biological variables. CONCLUSIONS: These results confirm the predictive value at a lesion level of the biologically inspired CD8 radiomics score in melanoma patients treated with anti-PD1-based immunotherapy and may be interesting to assess the disease spatial heterogeneity to evaluate the patient prognosis with potential clinical implication such as tumor selection for focal ablative therapies.

Development of a Machine Learning Classifier Based on Radiomic Features Extracted From Post-Contrast 3D T1-Weighted MR Images to Distinguish Glioblastoma From Solitary Brain Metastasis.

OBJECTIVES: To differentiate Glioblastomas (GBM) and Brain Metastases (BM) using a radiomic features-based Machine Learning (ML) classifier trained from post-contrast three-dimensional T1-weighted (post-contrast 3DT1) MR imaging, and compare its performance in medical diagnosis versus human experts, on a testing cohort. METHODS: We enrolled 143 patients (71 GBM and 72 BM) in a retrospective bicentric study from January 2010 to May 2019 to train the classifier. Post-contrast 3DT1 MR images were performed on a 3-Tesla MR unit and 100 radiomic features were extracted. Selection and optimization of the Machine Learning (ML) classifier was performed using a nested cross-validation. Sensitivity, specificity, balanced accuracy, and area under the receiver operating characteristic curve (AUC) were calculated as performance metrics. The model final performance was cross-validated, then evaluated on a test set of 37 patients, and compared to human blind reading using a McNemar's test. RESULTS: The ML classifier had a mean [95% confidence interval] sensitivity of 85% [77; 94], a specificity of 87% [78; 97], a balanced accuracy of 86% [80; 92], and an AUC of 92% [87; 97] with cross-validation. Sensitivity, specificity, balanced accuracy and AUC were equal to 75, 86, 80 and 85% on the test set. Sphericity 3D radiomic index highlighted the highest coefficient in the logistic regression model. There were no statistical significant differences observed between the performance of the classifier and the experts' blinded examination. CONCLUSIONS: The proposed diagnostic support system based on radiomic features extracted from post-contrast 3DT1 MR images helps in differentiating solitary BM from GBM with high diagnosis performance and generalizability.

Deep Learning-Based Concurrent Brain Registration and Tumor Segmentation.

Image registration and segmentation are the two most studied problems in medical image analysis. Deep learning algorithms have recently gained a lot of attention due to their success and state-of-the-art results in variety of problems and communities. In this paper, we propose a novel, efficient, and multi-task algorithm that addresses the problems of image registration and brain tumor segmentation jointly. Our method exploits the dependencies between these tasks through a natural coupling of their interdependencies during inference. In particular, the similarity constraints are relaxed within the tumor regions using an efficient and relatively simple formulation. We evaluated the performance of our formulation both quantitatively and qualitatively for registration and segmentation problems on two publicly available datasets (BraTS 2018 and OASIS 3), reporting competitive results with other recent state-of-the-art methods. Moreover, our proposed framework reports significant amelioration (p < 0.005) for the registration performance inside the tumor locations, providing a generic method that does not need any predefined conditions (e.g., absence of abnormalities) about the volumes to be registered. Our implementation is publicly available online at https://github.com/TheoEst/joint_registration_tumor_segmentation.

Standardization of brain MR images across machines and protocols: bridging the gap for MRI-based radiomics.

Radiomics relies on the extraction of a wide variety of quantitative image-based features to provide decision support. Magnetic resonance imaging (MRI) contributes to the personalization of patient care but suffers from being highly dependent on acquisition and reconstruction parameters. Today, there are no guidelines regarding the optimal pre-processing of MR images in the context of radiomics, which is crucial for the generalization of published image-based signatures. This study aims to assess the impact of three different intensity normalization methods (Nyul, WhiteStripe, Z-Score) typically used in MRI together with two methods for intensity discretization (fixed bin size and fixed bin number). The impact of these methods was evaluated on first- and second-order radiomics features extracted from brain MRI, establishing a unified methodology for future radiomics studies. Two independent MRI datasets were used. The first one (DATASET1) included 20 institutional patients with WHO grade II and III gliomas who underwent post-contrast 3D axial T1-weighted (T1w-gd) and axial T2-weighted fluid attenuation inversion recovery (T2w-flair) sequences on two different MR devices (1.5 T and 3.0 T) with a 1-month delay. Jensen-Shannon divergence was used to compare pairs of intensity histograms before and after normalization. The stability of first-order and second-order features across the two acquisitions was analysed using the concordance correlation coefficient and the intra-class correlation coefficient. The second dataset (DATASET2) was extracted from the public TCIA database and included 108 patients with WHO grade II and III gliomas and 135 patients with WHO grade IV glioblastomas. The impact of normalization and discretization methods was evaluated based on a tumour grade classification task (balanced accuracy measurement) using five well-established machine learning algorithms. Intensity normalization highly improved the robustness of first-order features and the performances of subsequent classification models. For the T1w-gd sequence, the mean balanced accuracy for tumour grade classification was increased from 0.67 (95% CI 0.61-0.73) to 0.82 (95% CI 0.79-0.84, P = .006), 0.79 (95% CI 0.76-0.82, P = .021) and 0.82 (95% CI 0.80-0.85, P = .005), respectively, using the Nyul, WhiteStripe and Z-Score normalization methods compared to no normalization. The relative discretization makes unnecessary the use of intensity normalization for the second-order radiomics features. Even if the bin number for the discretization had a small impact on classification performances, a good compromise was obtained using the 32 bins considering both T1w-gd and T2w-flair sequences. No significant improvements in classification performances were observed using feature selection. A standardized pre-processing pipeline is proposed for the use of radiomics in MRI of brain tumours. For models based on first- and second-order features, we recommend normalizing images with the Z-Score method and adopting an absolute discretization approach. For second-order feature-based signatures, relative discretization can be used without prior normalization. In both cases, 32 bins for discretization are recommended. This study may pave the way for the multicentric development and validation of MR-based radiomics biomarkers.

Dosimetry-Driven Quality Measure of Brain Pseudo Computed Tomography Generated From Deep Learning for MRI-Only Radiation Therapy Treatment Planning.

PURPOSE: This study aims to evaluate the impact of key parameters on the pseudo computed tomography (pCT) quality generated from magnetic resonance imaging (MRI) with a 3-dimensional (3D) convolutional neural network. METHODS AND MATERIALS: Four hundred two brain tumor cases were retrieved, yielding associations between 182 computed tomography (CT) and T1-weighted MRI (T1) scans, 180 CT and contrast-enhanced T1-weighted MRI (T1-Gd) scans, and 40 CT, T1, and T1-Gd scans. A 3D CNN was used to map T1 or T1-Gd onto CT scans and evaluate the importance of different components. First, the training set size's influence on testing set accuracy was assessed. Moreover, we evaluated the MRI sequence impact, using T1-only and T1-Gd-only cohorts. We then investigated 4 MRI standardization approaches (histogram-based, zero-mean/unit-variance, white stripe, and no standardization) based on training, validation, and testing cohorts composed of 242, 81, and 79 patients cases, respectively, as well as a bias field correction influence. Finally, 2 networks, namely HighResNet and 3D UNet, were compared to evaluate the architecture's impact on the pCT quality. The mean absolute error, gamma indices, and dose-volume histograms were used as evaluation metrics. RESULTS: Generating models using all the available cases for training led to higher pCT quality. The T1 and T1-Gd models had a maximum difference in gamma index means of 0.07 percentage point. The mean absolute error obtained with white stripe was 78 ± 22 Hounsfield units, which slightly outperformed histogram-based, zero-mean/unit-variance, and no standardization (P < .0001). Regarding the network architectures, 3%/3 mm gamma indices of 99.83% ± 0.19% and 99.74% ± 0.24% were obtained for HighResNet and 3D UNet, respectively. CONCLUSIONS: Our best pCTs were generated using more than 200 samples in the training data set. Training with T1 only and T1-Gd only did not significantly affect performance. Regardless of the preprocessing applied, the dosimetry quality remained equivalent and relevant for potential use in clinical practice.

Radiomics to predict outcomes and abscopal response of patients with cancer treated with immunotherapy combined with radiotherapy using a validated signature of CD8 cells.

BACKGROUND: Combining radiotherapy (RT) with immuno-oncology (IO) therapy (IORT) may enhance IO-induced antitumor response. Quantitative imaging biomarkers can be used to provide prognosis, predict tumor response in a non-invasive fashion and improve patient selection for IORT. A biologically inspired CD8 T-cells-associated radiomics signature has been developed on previous cohorts. We evaluated here whether this CD8 radiomic signature is associated with lesion response, whether it may help to assess disease spatial heterogeneity for predicting outcomes of patients treated with IORT. We also evaluated differences between irradiated and non-irradiated lesions. METHODS: Clinical data from patients with advanced solid tumors in six independent clinical studies of IORT were investigated. Immunotherapy consisted of 4 different drugs (antiprogrammed death-ligand 1 or anticytotoxic T-lymphocyte-associated protein 4 in monotherapy). Most patients received stereotactic RT to one lesion. Irradiated and non-irradiated lesions were delineated from baseline and the first evaluation CT scans. Radiomic features were extracted from contrast-enhanced CT images and the CD8 radiomics signature was applied. A responding lesion was defined by a decrease in lesion size of at least 30%. Dispersion metrices of the radiomics signature were estimated to evaluate the impact of tumor heterogeneity in patient's response. RESULTS: A total of 94 patients involving multiple lesions (100 irradiated and 189 non-irradiated lesions) were considered for a statistical interpretation. Lesions with high CD8 radiomics score at baseline were associated with significantly higher tumor response (area under the receiving operating characteristic curve (AUC)=0.63, p=0.0020). Entropy of the radiomics scores distribution on all lesions was shown to be associated with progression-free survival (HR=1.67, p=0.040), out-of-field abscopal response (AUC=0.70, p=0.014) and overall survival (HR=2.08, p=0.023), which remained significant in a multivariate analysis including clinical and biological variables. CONCLUSIONS: These results enhance the predictive value of the biologically inspired CD8 radiomics score and suggests that tumor heterogeneity should be systematically considered in patients treated with IORT. This CD8 radiomics signature may help select patients who are most likely to benefit from IORT.

The complexity of tumor shape, spiculatedness, correlates with tumor radiomic shape features.

Radiomics extracts high-throughput quantitative data from medical images to contribute to precision medicine. Radiomic shape features have been shown to correlate with patient outcomes. However, how radiomic shape features vary in function of tumor complexity and tumor volume, as well as with method used for meshing and voxel resampling, remains unknown. The aims of this study are to create tumor models with varying degrees of complexity, or spiculatedness, and evaluate their relationship with quantitatively extracted shape features. Twenty-eight tumor models were mathematically created using spherical harmonics with the spiculatedness degree d being increased by increments of 3 (d = 11 to d = 92). Models were 3D printed with identical bases of 5 cm, imaged with a CT scanner with two different slice thicknesses, and semi-automatically delineated. Resampling of the resulting masks on a 1 × 1 × 1 mm3 grid was performed, and the voxel size of each model was then calculated to eliminate volume differences. Four MATLAB-based algorithms (isosurface (M1), isosurface filter (M2), isosurface remeshing (M3), and boundary (M4)) were used to extract nine 3D features (Volume, Surface area, Surface-to-volume, Compactness1, Compactness2, Compactness3, Spherical Disproportion, Sphericity and Fractional Concavity). To quantify the impact of 3D printing, acquisition, segmentation and meshing, features were computed directly from the stereolithography (STL) file format that was used for 3D printing, and compared to those computed. Changes in feature values between 0.6 and 2 mm slice acquisitions were also compared. Spearman's rank-order correlation coefficients were computed to determine the relationship of each shape feature with spiculatedness for each of the four meshing algorithms. Percent changes were calculated between shape features extracted from the original and resampled contoured images to evaluate the influence of spatial resampling. Finally, the percent change in shape features when the volume was changed from 25% to 150% of their original volume was quantified for three distinct tumor models and compared to the percent change observed when modifying the spiculatedness of the model from d = 11 to d = 92. Values extracted using isosurface remeshing method are the closest to the STL reference ones, with mean differences less than 10.8% (Compactness2) for all features. Seven of the eight features had strong significant correlations with tumor model complexity irrespective of the meshing algorithm (r > 0.98, p < 10-4), with fractional concavity having the lowest correlation coefficient (r = 0.83, p < 10-4, M2). Comparisons of features extracted from the 0.6 and 2 mm slice thicknesses showed that mean differences were from 2.1% (Compactness3) to 12.7% (Compactness2) for the isosurface remeshing method. Resampling on a 1 × 1 × 1 mm3 grid resulted in between 1.3% (Compactness3) to 9.5% (Fractional Concavity) mean changes in feature values. Compactness2, Compactness3, Spherical Disproportion, Sphericity and Fractional Concavity were the features least affected by volume changes. Compactness1 had a 90.4% change with volume, which was greater than the change between the least and most spiculated models. This is the first methodological study that directly demonstrates the relationship of tumor spiculatedness with radiomic shape features, that also produced 3D tumor models, which may serve as reference phantoms for future radiomic studies. Surface Area, Surface-to-volume, and Spherical Disproportion had direct relationships with spiculatedness while the three formulas for Compactness, Sphericity and Fractional Concavity had inverse relationships. The features Compactness2, Compactness3, Spherical Disproportion, and Sphericity should be prioritized as these have minimal variations with volume changes, slice thickness and resampling.

Commentary: Toward an improved understanding of recidivism.

Setting out to identify predictors of recidivism in seriously mentally ill (SMI) persons who return to their community from prison, Hall and colleagues offer to the forensic community a substantive contribution to the organization of interventions aimed to promote both the social and the psychiatric recovery of these individuals and to reduce their rate of re-entry into the penal system. While their work in the field is praiseworthy, I would like to share my thoughts and considerations to stimulate a dialogue about a subject that should be of interest to both clinicians and policy makers.