Radiomic detection of abnormal brain regions in tuberous sclerosis complex

Radiomics refers to the extraction of quantitative information from medical images and is most commonly utilized in oncology to provide ancillary information for solid tumor diagnosis, prognosis, and treatment response. The traditional radiomic pipeline involves segmentation of volumes of interest with comparison to normal brain. In other neurologic disorders, such as epilepsy, lesion delineation may be difficult or impossible due to poor anatomic definition, small size, and multifocal or diffuse distribution. Tuberous sclerosis complex (TSC) is a rare genetic disease in which brain magnetic resonance imaging (MRI) demonstrates multifocal abnormalities with variable imaging and epileptogenic features.

The purpose of this study was to develop a radiomic workflow for identification of abnormal brain regions in TSC, using a whole-brain atlas-based approach with generation of heatmaps based on signal deviation from normal controls.

This was a retrospective pilot study utilizing high-resolution whole-brain 3D FLAIR MRI datasets from retrospective enrollment of tuberous sclerosis complex (TSC) patients and normal controls. Subjects underwent MRI including high-resolution 3D FLAIR sequences. Preprocessing included skull stripping, coregistration, and intensity normalization. Using the Brainnetome and Harvard-Oxford atlases, brain regions were parcellated into 318 discrete regions. Expert neuroradiologists spatially labeled all tubers in TSC patients using ITK-SNAP. The pyradiomics toolbox was used to extract 88 radiomic features based on IBSI guidelines, comparing tuber-affected and non-tuber-affected parenchyma in TSC patients, as well as normal brain tissue in control patients. For model training and validation, regions with tubers from 20 TSC patients and 30 normal control subjects were randomly divided into two training sets (80%) and two validation sets (20%). Additional model testing was performed on a separate group of 20 healthy controls. LASSO (least absolute shrinkage and selection operator) was used to perform variable selection and regularization to identify regions containing tubers. Relevant radiomic features selected by LASSO were combined to produce a radiomic score ω, defined as the sum of squared differences from average control group values. Region-specific ω scores were converted to heat maps and spatially coregistered with brain MRI to reflect overall radiomic deviation from normal.

The proposed radiomic workflow allows for quantification of deviation from normal in 318 regions of the brain with the use of a summative radiomic score ω. This score can be used to generate spatially registered heatmaps to identify brain regions with radiomic abnormalities. The pilot study of TSC showed radiomic scores ω that were statistically different in regions containing tubers from regions without tubers/normal brain (p < 0.0001). Our model exhibits an AUC of 0.81 (95% confidence interval: 0.78-0.84) on the testing set, and the best threshold obtained on the training set, when applied to the testing set, allows us to identify regions with tubers with a specificity of 0.91 and a sensitivity of 0.60.

We describe a whole-brain atlas-based radiomic approach to identify abnormal brain regions in TSC patients. This approach may be helpful for identifying specific regions of interest based on relatively greater signal deviation, particularly in clinical scenarios with numerous or poorly defined anatomic lesions.