A radiomics model from joint FDG-PET and MRI texture features for the prediction of lung metastases in soft-tissue sarcomas of the extremities

This study aims at developing a joint FDG-PET and MRI texture-based model for the early evaluation of lung metastasis risk in soft-tissue sarcomas (STSs). We investigate if the creation of new composite textures from the combination of FDG-PET and MR imaging information could better identify aggressive tumours. Towards this goal, a cohort of 51 patients with histologically proven STSs of the extremities was retrospectively evaluated. All patients had pre-treatment FDG-PET and MRI scans comprised of T1-weighted and T2-weighted fat-suppression sequences (T2FS). Nine non-texture features (SUV metrics and shape features) and forty-one texture features were extracted from the tumour region of separate (FDG-PET, T1 and T2FS) and fused (FDG-PET/T1 and FDG-PET/T2FS) scans. Volume fusion of the FDG-PET and MRI scans was implemented using the wavelet transform. The influence of six different extraction parameters on the predictive value of textures was investigated. The incorporation of features into multivariable models was performed using logistic regression. The multivariable modeling strategy involved imbalance-adjusted bootstrap resampling in the following four steps leading to final prediction model construction: (1) feature set reduction; (2) feature selection; (3) prediction performance estimation; and (4) computation of model coefficients. Univariate analysis showed that the isotropic voxel size at which texture features were extracted had the most impact on predictive value. In multivariable analysis, texture features extracted from fused scans significantly outperformed those from separate scans in terms of lung metastases prediction estimates. The best performance was obtained using a combination of four texture features extracted from FDG-PET/T1 and FDG-PET/T2FS scans. This model reached an area under the receiver-operating characteristic curve of 0.984 ± 0.002, a sensitivity of 0.955 ± 0.006, and a specificity of 0.926 ± 0.004 in bootstrapping evaluations. Ultimately, lung metastasis risk assessment at diagnosis of STSs could improve patient outcomes by allowing better treatment adaptation.

[1]  E. van den Berg,et al.  Metastasis in Soft Tissue Sarcomas: Prognostic Criteria and Treatment Perspectives , 2004, Cancer and Metastasis Reviews.

[2]  M. Hatt,et al.  Intratumor Heterogeneity Characterized by Textural Features on Baseline 18F-FDG PET Images Predicts Response to Concomitant Radiochemotherapy in Esophageal Cancer , 2011, The Journal of Nuclear Medicine.

[3]  Robert M. Haralick,et al.  Textural Features for Image Classification , 1973, IEEE Trans. Syst. Man Cybern..

[4]  M. Hatt,et al.  18F-FDG PET Uptake Characterization Through Texture Analysis: Investigating the Complementary Nature of Heterogeneity and Functional Tumor Volume in a Multi–Cancer Site Patient Cohort , 2015, The Journal of Nuclear Medicine.

[5]  Vicky Goh,et al.  Are Pretreatment 18F-FDG PET Tumor Textural Features in Non–Small Cell Lung Cancer Associated with Response and Survival After Chemoradiotherapy? , 2013, The Journal of Nuclear Medicine.

[6]  Belur V. Dasarathy,et al.  Image characterizations based on joint gray level-run length distributions , 1991, Pattern Recognit. Lett..

[7]  F. O’Sullivan,et al.  Sarcoma tumor FDG uptake measured by PET and patient outcome: a retrospective analysis , 2002, European Journal of Nuclear Medicine and Molecular Imaging.

[8]  Robert King,et al.  Textural features corresponding to textural properties , 1989, IEEE Trans. Syst. Man Cybern..

[9]  Lubomir M. Hadjiiski,et al.  Classifier performance prediction for computer-aided diagnosis using a limited dataset. , 2008, Medical physics.

[10]  Joel Max,et al.  Quantizing for minimum distortion , 1960, IRE Trans. Inf. Theory.

[11]  Mary M. Galloway,et al.  Texture analysis using gray level run lengths , 1974 .

[12]  I El Naqa,et al.  Dose response explorer: an integrated open-source tool for exploring and modelling radiotherapy dose–volume outcome relationships , 2006, Physics in medicine and biology.

[13]  Joseph O Deasy,et al.  CERR: a computational environment for radiotherapy research. , 2003, Medical physics.

[14]  A. Jemal,et al.  Cancer statistics, 2014 , 2014, CA: a cancer journal for clinicians.

[15]  Issam El-Naqa,et al.  Exploring feature-based approaches in PET images for predicting cancer treatment outcomes , 2009, Pattern Recognit..

[16]  M. Brennan Soft tissue sarcoma: advances in understanding and management. , 2005, The surgeon : journal of the Royal Colleges of Surgeons of Edinburgh and Ireland.

[17]  M. Burt,et al.  Pulmonary metastases from soft tissue sarcoma: analysis of patterns of diseases and postmetastasis survival. , 1999, Annals of surgery.

[18]  I. Buvat,et al.  A review of partial volume correction techniques for emission tomography and their applications in neurology, cardiology and oncology , 2012, Physics in medicine and biology.

[19]  D H Leung,et al.  Multifactorial analysis of the survival of patients with distant metastasis arising from primary extremity sarcoma , 1999, Cancer.

[20]  James F. Greenleaf,et al.  Use of gray value distribution of run lengths for texture analysis , 1990, Pattern Recognit. Lett..

[21]  John G. Griffiths,et al.  Least squares ellipsoid specific fitting , 2004, Geometric Modeling and Processing, 2004. Proceedings.

[22]  S. P. Lloyd,et al.  Least squares quantization in PCM , 1982, IEEE Trans. Inf. Theory.

[23]  Andre Dekker,et al.  Radiomics: the process and the challenges. , 2012, Magnetic resonance imaging.

[24]  Patrick Granton,et al.  Radiomics: extracting more information from medical images using advanced feature analysis. , 2012, European journal of cancer.

[25]  C. Freeman,et al.  Metabolic activity measured on PET/CT correlates with clinical outcomes in patients with limb and girdle sarcomas , 2011, Journal of surgical oncology.

[26]  P. Lambin,et al.  Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach , 2014, Nature Communications.

[27]  K. Aldape,et al.  Identification of noninvasive imaging surrogates for brain tumor gene-expression modules , 2008, Proceedings of the National Academy of Sciences.

[28]  G. Collewet,et al.  Influence of MRI acquisition protocols and image intensity normalization methods on texture classification. , 2004, Magnetic resonance imaging.

[29]  P. Lambin,et al.  Predicting outcomes in radiation oncology—multifactorial decision support systems , 2013, Nature Reviews Clinical Oncology.

[30]  Ronald Boellaard,et al.  Evaluation of a cumulative SUV-volume histogram method for parameterizing heterogeneous intratumoural FDG uptake in non-small cell lung cancer PET studies , 2011, European Journal of Nuclear Medicine and Molecular Imaging.

[31]  Michael Mitzenmacher,et al.  Detecting Novel Associations in Large Data Sets , 2011, Science.

[32]  M. D. Ernst Permutation Methods: A Basis for Exact Inference , 2004 .

[33]  D. Longo,et al.  Tumor heterogeneity and personalized medicine. , 2012, The New England journal of medicine.

[34]  Howard Y. Chang,et al.  Decoding global gene expression programs in liver cancer by noninvasive imaging , 2007, Nature Biotechnology.

[35]  J. Bradley,et al.  Combined PET/CT image characteristics for radiotherapy tumor response in lung cancer. , 2012, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[36]  Bernard Fertil,et al.  Texture indexes and gray level size zone matrix. Application to cell nuclei classification , 2009 .