Identification and validation of radiomic features from computed tomography for preoperative classification of neuroblastic tumors in children
BMC Pediatrics volume 23, Article number: 262 (2023)
To identify radiomic features that can predict the pathological type of neuroblastic tumor in children.
Data on neuroblastic tumors in 104 children were retrospectively analyzed. There were 14 cases of ganglioneuroma, 24 cases of ganglioneuroblastoma, and 65 cases of neuroblastoma. Stratified sampling was used to randomly allocate the cases into the training and validation sets in a ratio of 3:1. The maximum relevance–minimum redundancy algorithm was used to identify the top 10 of two clinical features and 851 radiomic features in portal venous–phase contrast-enhanced computed tomography images. Least absolute shrinkage and selection operator regression was used to classify tumors in two binary steps: first as ganglioneuroma compared to the other two types, then as ganglioneuroblastoma compared to neuroblastoma.
Based on 10 clinical-radiomic features, the classifier identified ganglioneuroma compared to the other two tumor types in the validation dataset with sensitivity of 100.0%, specificity of 81.8%, and an area under the receiver operating characteristic curve (AUC) of 0.875. The classifier identified ganglioneuroblastoma versus neuroblastoma with a sensitivity of 83.3%, a specificity of 87.5%, and an AUC of 0.854. The overall accuracy of the classifier across all three types of tumors was 80.8%.
Radiomic features can help predict the pathological type of neuroblastic tumors in children.
Neuroblastic tumors, which originate from sympathetic natural crest cells in the early stage of neural development , are the most frequent type of extracranial solid tumors in children . Although such tumors can occur in any part of the sympathetic nervous system, most often they include the abdominal spinal sympathetic ganglia (in 60% of cases) and adrenal gland (in 30% of cases) .
Three types of neuroblastic tumors have been described, which differ significantly in the course of their disease and in optimal treatment : ganglioneuroma, ganglioneuroblastoma, and neuroblastoma. Neuroblastoma, which is highly malignant, contains undifferentiated neuroblasts . This tumor type accounts for 10% of malignant tumors and 15% of tumor-related mortality in children. Ganglioneuroblastoma, which is considerably less malignant than neuroblastoma, contains neuroblasts, glial fibers, proliferative nerve sheath cells, and ganglion cells in different degrees of differentiation. Ganglioneuroma is benign.
To optimize treatment, it is necessary to accurately determine the type of neuroblastic tumors in children . At present, tumors are assigned to one of the three types based on pathological examination of the biopsies, which is invasive and carries a risk of complications . Moreover, such typing depends on the experience of the clinician and the exact location of the biopsy in the heterogeneous tumor tissue. A more objective and noninvasive approach could avoid these shortcomings.
With the constant development and cross-integration of medical imaging, computer science, informatics, and other disciplines, traditional image diagnostics is undergoing a new round of change. Radiomics realizes high-throughput feature extraction, analysis, and quantification of image data using automated algorithms, thus providing new interpretations of the potential features of the images. Radiomics can extract microscopic details of a large number of tumor lesions that are difficult to identify with human eyes and may quantify their internal subtle structures to obtain a set of image markers related to the disease. It has been proven that many radiomic features are useful for predicting the stage of an abdominal tumor, evaluation of effectiveness, and prognosis . Computed tomography (CT) is widely used to detect and diagnose neuroblastic tumors in the clinical environment, and certain radiomic features from images obtained using CT or other methods are useful in the diagnosis and typing of various diseases [9,10,11,12]. CT is noninvasive and effectively records the heterogeneity of the tumor [13, 14], so it is important to determine whether we could identify CT-based radiomic features for typing neuroblastic tumors in children. Based on this idea, such features have proven useful in predicting the proliferation of MYCN in neuroblastoma and ganglioneuroblastoma .
Therefore, here we have defined and validated a radiomic classifier based on contrast-enhanced CT to distinguish the three types of neuroblastic tumors.
Methods and materials
Data were analyzed retrospectively for all patients who had been treated for pathology-confirmed neuroblastic tumors of any type at the Children’s Hospital of Soochow University (Suzhou, China) between January 2015 and December 2021. This study was approved by the institutional review board of the Children’s Hospital of Soochow University. Written informed consent was provided by the parents or legal guardians of the children.
CT and image segmentation
Contrast-enhanced CT was conducted using a 64-slice system (GE Optima CT660, GE Healthcare, Optima 660; GE Medical System, Milwaukee, WI, USA). The scanning parameters were as follows: tube voltage, 70–120 kVp; tube current, 10–1041 mA; rotation time, 0.35–4.53 s; pixel spacing, 0.32–0.73 mm; slice thickness, 5 mm; and slice interval, 5 mm.
Patients who could not cooperate were sedated with 10% chloral hydrate (0.5 ml/kg) for a total of no more than 10ml. Subjects whose CT image motion artifacts were too large and could influence clinical diagnosis were again given an enhanced CT scan at the scheduled time.
Using the 3D Slicer software (https://www.slicer.org/), a radiologist with 5 years of experience manually described all the tumors in the CT images of the portal vein phase (Fig. 1a). Segmentation was supervised and confirmed by a radiologist with 10 years of experience.
Image preprocessing and feature extraction
Portal venous–phase contrast-enhanced CT images, which were resampled to a voxel size of 0.5 mm ⋅ 0.5 mm ⋅ 5 mm to reduce inter-subject variation in the reconstruction parameters, were re-binned using a bin width of 20 to minimize bias due to sparsely populated matrices [16,17,18].
The following four types of radiomic features were extracted from the images using the Pyradiomics package (https://pyradiomics.readthedocs.io/en/latest/index.html) in Python 3.6: shape features, first-order features, texture features, and wavelet features (Fig. 1b) [9, 19]. Shape features, such as sphericity and compactness, describe the three-dimensional shape and size of the tumor. First-order features, such as kurtosis and median, are first-order statistics that quantify the distribution of voxel intensities. Texture features, which characterize the spatial distribution of voxel intensities in a neighborhood, are calculated based on the gray level co-occurrence matrix (“glcm”), gray level dependence matrix (“gldm”), gray level run length matrix (“glrlm”), gray level size zone matrix (“glszm”), and neighboring gray tone difference matrix (“ngtdm”). Wavelet features are the first-order and texture features that emerge after the wavelet filter is applied to the images. Finally, the extracted features were normalized using a robust method that scaled the measured parameters and excluded outliers based on the observed means and standard deviations .
Stratified sampling was used to randomly distribute the data into training and validation sets in a 3:1 ratio while achieving a similar distribution of the three tumor types between the two sets.
Next, the potentially diagnostic radiomic features were selected from all of the features in the training set (Fig. 1c). In the selection process, the radiomic features were compared among the three tumor types using one-way analysis of variance (ANOVA) for normally distributed variables, or the Kruskal–Wallis test for variables with a skewed distribution. The radiomics features associated with P < 0.05 and two clinical features (age and gender) were ranked using the maximum relevance–minimum redundancy algorithm [13, 14, 21], and the top 10 features were selected to create a classifier.
Derivation of the classifier
The final classifier was defined using least absolute shrinkage and selection operator (LASSO) regression (Fig. 1c) [14, 21, 22]. LASSO regularization of linear models of the input variables (radiomic features) and the response variable (tumor type) generated least-squares regression coefficients, which were then used to build the LASSO classifier. The model was optimized based on the minimum five-fold cross-validated mean squared error (MSE).
The final model for classifying the tumors into one of the three types proceeded through two binary classifications: first, ganglioneuroma compared to the other two types of tumors, followed by ganglioneuroblastoma compared to neuroblastoma. To compare the classification performance of conventional clinical features and radiomics features, two models were developed, one based solely on clinical data (clinical model) and the other on the same clinical data together with the top radiomic features mentioned in Sect. 2.4 (clinical-radiomics model).
Assessment of the classifier’s performance
The performance of the model in classifying the tumors into one type compared to the other two types was evaluated in terms of sensitivity, specificity, and the area under the receiver operating characteristic curve (AUC). To assess the overall accuracy, the number of correctly classified subjects with any tumor type was divided by the total number of subjects. To evaluate the accuracy for each type of tumor, the number of correctly classified subjects with that tumor type was divided by the total number of subjects with that type of tumor.
The analysis included 103 subjects, of whom 65 had neuroblastoma, 24 had ganglioneuroblastoma, and 14 had ganglioneuroma based on pathological data (Table 1).
For each subject, data were extracted for 851 shape, first-order, texture, and wavelet features and two clinical features (age and gender), and the following top 10 features emerged from the maximum relevance–minimum redundancy algorithm: “age,” “original_glcm_Correlation,” “wavelet-HLH_glszm_SmallAreaEmphasis,” “wavelet-LHL_glcm_Imc1,” “wavelet-LLH_glszm_LowGrayLevelZoneEmphasis,” “wavelet-LLL_glcm_Imc2,” “original_glcm_Imc1,” “original_glcm_Imc2,” “wavelet-LLL_glcm_Correlation,” and “wavelet-HLH_glszm_SizeZoneNonUniformityNormalized” (Fig. 2). In this way, three of the top 10 features were texture features, while six were wavelets.
By combining radiomic features and clinical data, the classifier showed a balanced accuracy of 80.8% in the validation set (Table 2), compared to only 46.2% for the classifier that included only clinical data (Table 3). In the training set, the radiomic-clinical classifier distinguished ganglioneuroma from the other two tumor types with a sensitivity of 90.0%, a specificity of 94.0%, and an AUC of 0.969 (Table 2; Fig. 3). The corresponding values in the validation set were 100.0%, 81.8%, and 0.875. In the training set, the radiomic-clinical classifier distinguished ganglioneuroblastoma from neuroblastoma with a sensitivity of 83.3%, a specificity of 91.8%, and an AUC of 0.931. The corresponding values in the validation set were 83.3%, 87.5%, and 0.854.
The overall accuracy of the radiomic-clinical classifier was 85.7% in the training set and 80.8% in the validation set (Fig. 4).
Although the three types of neuroblastic tumors occur in similar locations and cause similar symptoms in children, they differ substantially in prognosis and optimal treatment . Currently, the tumors are typed on the basis of pathological examination of biopsies, but this is invasive and carries the risk of complications . Here, we provide evidence that noninvasive tumor assessment based on clinical data and preoperative CT can predict the type of tumor and therefore help guide treatment.
The overall accuracy of our radiomic-clinical classifier was 80.8% in the validation set, suggesting the need for further optimization. Nevertheless, the model was reasonably accurate at identifying ganglioneuroma, showing sensitivity of 100.0% and specificity of 81.8%. Model performance can be improved by using large samples from multiple centers and by exploring machine learning algorithms other than LASSO. Taking into account our small, single-center sample, the radiomic-clinical model should be tested for generalizability and further optimized on the basis of larger samples from multiple centers.
We selected the radiomic features for our classifier using the maximum relevance–minimum redundancy algorithm, which has increasingly been used in radiomic studies since its first description in 2005 . Unlike earlier methods for feature selection, this algorithm takes into account the correlation between features, thereby reducing the risk of false positives and improving the generalizability of the model [24, 25].
Our results emphasize the usefulness of texture features in differentiating the types of neuroblastic tumors, most likely reflecting that such features can capture the heterogeneity of the tumor structure . Among the three tumor types, 10 features were distinguished, the overall accuracy of which was 85.7% in the training set and 80.8% in the validation set, and the performance was even better in the two binary classifications within the model. Thus, the ability to distinguish between neuroblastoma, ganglioneuroma, and ganglioneuroblastoma depends on the identification of differences in the structure of tissues.
Six of the 10 features extracted in the present study belong to spatial gray level co-occurrence matrix (GLCM) features. Previous studies have shown that GLCM features are helpful in the pretreatment prediction of pathological complete response to no special type(NST) in breast cancer . GLSZM texture features have been useful in differentiating between two different tumors in studies using CT imaging omics to distinguish between pelvic rhabdomyosarcoma and yolk cystoma in children . In the present study, GLCM and GLSMZ were useful in predicting the pathological type of neuroblastic tumors in children, suggesting that the texture characteristics of different tumor types differ. At the same time, most of the valuable features in the present study were based on wavelet and LoG. Both the wavelet and the LoG are high-order statistical methods for placing filter grids on the image to extract repetitive or nonrepetitive patterns, both of which help reveal more valuable information that is not visible in the lesion .
The present study determines the feasibility of using contrast-enhanced CT to diagnose neuroblastic tumors before surgery. Further validation and optimization studies may help develop an accurate, noninvasive tool that can guide treatment.
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Whittle SB, Smith V, Doherty E, Zhao S, McCarty S, Zage PE. Overview and recent advances in the treatment of neuroblastoma. Expert Rev Anticancer Ther. 2017;17:369–86.
Pearson ADJ, P.R.N.I.P R, Plowman PN, Pieters R, editors. Paediatric oncology. London, England: Arnold; 2004. pp. 386–414.
Louis CU, Shohet JM. Neuroblastoma: molecular pathogenesis and therapy. Annu Rev Med. 2015;66:49–63.
Brisse HJ, McCarville MB, Granata C, Krug KB, Wootton-Gorges SL, Kanegawa K, et al. Guidelines for imaging and staging of neuroblastic tumors: consensus report from the International Neuroblastoma Risk Group Project. Radiology. 2011;261:243–57.
Pinto NR, Applebaum MA, Volchenboum SL, Matthay KK, London WB, Ambros PF, et al. Advances in risk classification and treatment strategies for Neuroblastoma. J Clin Oncol. 2015;33:3008–17.
Simon T, Hero B, Schulte JH, Deubzer H, Hundsdoerfer P, von Schweinitz D, et al. 2017 GPOH Guidelines for diagnosis and treatment of patients with neuroblastic tumors. Klin Padiatr. 2017;229:147–67.
Hassan SF, Mathur S, Magliaro TJ, Larimer EL, Ferrell LB, Vasudevan SA, et al. Needle core vs open biopsy for diagnosis of intermediate- and high-risk neuroblastoma in children. J Pediatr Surg. 2012;47:1261–6.
Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278(2):563–77.
Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. 2014;5:4006.
Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Reviews Clin Oncol. 2017;14:749–62.
Liu C, Jiao D, Liu Z. Artificial Intelligence (AI)-aided Disease Prediction. BIO Integr. 2020;1:130–6.
Wang J, Wu D, Sun M, Peng Z, Lin Y, Lin H et al. Deep Segmentation Feature-Based Radiomics Improves Recurrence Prediction of Hepatocellular Carcinoma. BME Frontiers. 2022; 2022:1–12.
Jiang W, Song Y, Sun Z, Qiu J, Shi L. Dosimetric factors and Radiomics features within different regions of interest in planning CT images for improving the prediction of Radiation Pneumonitis. Int J Radiat Oncol Biol Phys. 2021;110:1161–70.
Xia W, Hu B, Li H, Geng C, Wu Q, Yang L, et al. Multiparametric-MRI-Based Radiomics Model for differentiating primary Central Nervous System Lymphoma from Glioblastoma: Development and Cross-Vendor Validation. J Magn Reson Imaging. 2021;53:242–50.
Di Giannatale A, Di Paolo PL, Curione D, Lenkowicz J, Napolitano A, Secinaro A, et al. Radiogenomics prediction for MYCN amplification in neuroblastoma: a hypothesis generating study. Pediatr Blood Cancer. 2021;68:e29110.
Loi S, Mori M, Benedetti G, Partelli S, Broggi S, Cattaneo GM, et al. Robustness of CT radiomic features against image discretization and interpolation in characterizing pancreatic neuroendocrine neoplasms. Phys Med. 2020;76:125–33.
Fave X, Mackin D, Yang J, Zhang J, Fried D, Balter P, et al. Can radiomics features be reproducibly measured from CBCT images for patients with non-small cell lung cancer? Med Phys. 2015;42:6784–97.
Ferreira Junior JR, Koenigkam-Santos M, Cipriano FEG, Fabro AT, Azevedo-Marques PM. Radiomics-based features for pattern recognition of lung cancer histopathology and metastases. Comput Methods Programs Biomed. 2018;159:23–30.
Zwanenburg A, Vallières M, Abdalah MA, Aerts HJWL, Andrearczyk V, Apte A, et al. The image Biomarker Standardization Initiative: standardized quantitative Radiomics for High-Throughput Image-based phenotyping. Radiology. 2020;295:328–38.
Cao XH, Stojkovic I, Obradovic Z. A robust data scaling algorithm to improve classification accuracies in biomedical data. BMC Bioinformatics. 2016;17:359.
Meng X, Xia W, Xie P, Zhang R, Li W, Wang M. Preoperative radiomic signature based on multiparametric magnetic resonance imaging for noninvasive evaluation of biological characteristics in rectal cancer. Eur Radiol. 2019;29(6):3200–9.
Jian J, Li Y, Pickhardt PJ, Xia W, He Z, Zhang R, et al. MR image-based radiomics to differentiate type Ι and type ΙΙ epithelial ovarian cancers. Eur Radiol. 2021;31:403–10.
Peng H, Long F, Ding C. Feature selection based on mutual information: criteria of max-dependency, max-relevance, and min-redundancy. IEEE Trans Pattern Anal Mach Intell. 2005;27:1226–38.
Destyanto TYR, Lin RF. Detecting computer activities using eye-movement features. J Ambient Intell Human Comput. 2020. https://doi.org/10.1007/s12652-020-02683-8.
Zhou X, Gao X, Wang J, Yu H, Wang Z, Chi Z. Eye tracking data guided feature selection for image classification. Pattern Recognit. 2017;63:56–70.
Li Z, Mao Y, Li H, Yu G, Wan H, Li B. Differentiating brain metastases from different pathological types of lung cancers using texture analysis of T1 postcontrast MR. Magn Reson Med. 2016;76:1410–9.
Tsai HY, Tsai TY, Wu CH, Chung WS, Wang JC, Hsu JS, et al. Integration of clinical and CT-Based Radiomic features for pretreatment prediction of pathologic complete response to Neoadjuvant systemic therapy in breast Cancer. Cancers (Basel). 2022;14(24):6261.
Chen X, Huang Y, He L, Zhang T, Zhang L, Ding H. CT-Based Radiomics to differentiate pelvic Rhabdomyosarcoma from Yolk Sac Tumors in Children. Front Oncol. 2020;10:584272.
This study was supported by the Key Research and Development Program of Shandong Province (2021SFGC0104), Key Research and Development Program of Jiangsu Province (BE2021663), Key Research and Development Program (Social Development) of Jiangsu Province (BE2020659), Introductory Project in the Suzhou “Clinical Medicine Expert Team” Program (SZYJTD201706), Suzhou Science and Technology Development (Minsheng Technology) Project (SS2019011), and Suzhou Clinical Medicine Center (Szlcyxzx202104).
Ethics approval and consent to participate
The authors ensured that questions related to the accuracy or integrity of any part of this work were appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (2013 revision), and the study protocol was approved by the Ethics Committee of the Children’s Hospital of Soochow University. Written informed consent was provided by the parents or legal guardians of the children.
Consent for publication
Conflict of interest
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhao, L., Shi, L., Huang, Sg. et al. Identification and validation of radiomic features from computed tomography for preoperative classification of neuroblastic tumors in children. BMC Pediatr 23, 262 (2023). https://doi.org/10.1186/s12887-023-04057-3