Identification of Multiclass Pediatric Low-Grade Neuroepithelial Tumor Molecular Subtype with ADC MR Imaging and Machine Learning

BACKGROUND AND PURPOSE: Molecular biomarker identification increasingly influences the treatment planning of pediatric low-grade neuroepithelial tumors (PLGNTs). We aimed to develop and validate a radiomics-based ADC signature predictive of the molecular status of PLGNTs. MATERIALS AND METHODS: In this retrospective bi-institutional study, we searched the PACS for baseline brain MRIs from children with PLGNTs. Semiautomated tumor segmentation on ADC maps was performed using the semiautomated level tracing effect tool with 3D Slicer. Clinical variables, including age, sex, and tumor location, were collected from chart review. The molecular status of tumors was derived from biopsy. Multiclass random forests were used to predict the molecular status and fine-tuned using a grid search on the validation sets. Models were evaluated using independent and unseen test sets based on the combined data, and the area under the receiver operating characteristic curve (AUC) was calculated for the prediction of 3 classes: KIAA1549-BRAF fusion, BRAF V600E mutation, and non-BRAF cohorts. Experiments were repeated 100 times using different random data splits and model initializations to ensure reproducible results. RESULTS: Two hundred ninety-nine children from the first institution and 23 children from the second institution were included (53.6% male; mean, age 8.01 years; 51.8% supratentorial; 52.2% with KIAA1549-BRAF fusion). For the 3-class prediction using radiomics features only, the average test AUC was 0.74 (95% CI, 0.73–0.75), and using clinical features only, the average test AUC was 0.67 (95% CI, 0.66–0.68). The combination of both radiomics and clinical features improved the AUC to 0.77 (95% CI, 0.75–0.77). The diagnostic performance of the per-class test AUC was higher in identifying KIAA1549-BRAF fusion tumors among the other subgroups (AUC = 0.81 for the combined radiomics and clinical features versus 0.75 and 0.74 for BRAF V600E mutation and non-BRAF, respectively). CONCLUSIONS: ADC values of tumor segmentations have differentiative signals that can be used for training machine learning classifiers for molecular biomarker identification of PLGNTs. ADC-based pretherapeutic differentiation of the BRAF status of PLGNTs has the potential to avoid invasive tumor biopsy and enable earlier initiation of targeted therapy.