Diagnostic Value of Fractal Analysis for the Differentiation of Brain Tumors Using 3-Tesla Magnetic Resonance Susceptibility-Weighted Imaging

Fractal dimension and lacunarity measures of glioma subcomponents are discriminative of the grade of gliomas and IDH status

Since the overall glioma mass and its subcomponents-enhancing region (malignant part of the tumor), non-enhancing (less aggressive tumor cells), necrotic core (dead cells), and edema (water deposition)-are complex and irregular structures, non-Euclidean geometric measures such as fractal dimension (FD or "fractality") and lacunarity are needed to quantify their structural complexity. Fractality measures the extent of structural irregularity, while lacunarity measures the spatial distribution or gaps. The complex geometric patterns of the glioma subcomponents may be closely associated with the grade and molecular landscape. Therefore, we measured FD and lacunarity in the glioma subcomponents and developed machine learning models to discriminate between tumor grades and isocitrate dehydrogenase (IDH) gene status. 3D fractal dimension (FD3D) and lacunarity (Lac3D) were measured for the enhancing, non-enhancing plus necrotic core, and edema-subcomponents using preoperative structural-MRI obtained from the The Cancer Genome Atlas (TCGA) and University of California San Francisco Preoperative Diffuse Glioma MRI (UCSF-PDGM) glioma cohorts. The FD3D and Lac3D measures of the tumor-subcomponents were then compared across glioma grades (HGGs: high-grade gliomas vs. LGGs: low-grade gliomas) and IDH status (mutant vs. wild type). Using these measures, machine learning platforms discriminative of glioma grade and IDH status were developed. Kaplan-Meier survival analysis was used to assess the prognostic significance of FD3D and Lac3D measurements. HGG exhibited significantly higher fractality and lower lacunarity in the enhancing subcomponent, along with lower fractality in the non-enhancing subcomponent compared to LGG. This suggests that a highly irregular and complex geometry in the enhancing-subcomponent is a characteristic feature of HGGs. A comparison of FD3D and Lac3D between IDH-wild type and IDH-mutant gliomas revealed that mutant gliomas had ~2.5-fold lower FD3D in the enhancing-subcomponent and higher FD3D with lower Lac3D in the non-enhancing subcomponent, indicating a less complex and smooth enhancing subcomponent, and a more continuous non-enhancing subcomponent as features of IDH-mutant gliomas. Supervised ML models using FD3D from both the enhancing and non-enhancing subcomponents together demonstrated high-sensitivity in discriminating glioma grades (~97.9%) and IDH status (~94.4%). A combined fractal estimation of the enhancing and non-enhancing subcomponents using MR images could serve as a non-invasive, precise, and quantitative measure for discriminating glioma grade and IDH status. The combination of 2-hydroxyglutarate-magnetic resonance spectroscopy (2HG-MRS) with FD3D and Lac3D quantification may be established as a robust imaging signature for glioma subtyping.

Computational Fractal-Based Analysis of MR Susceptibility-Weighted Imaging (SWI) in Neuro-Oncology and Neurotraumatology

Susceptibility-weighted imaging (SWI) is a magnetic resonance imaging (MRI) technique able to depict the magnetic susceptibility produced by different substances, such as deoxyhemoglobin, calcium, and iron. The main application of SWI in clinical neuroimaging is detecting microbleedings and venous vasculature. Quantitative analyses of SWI have been developed over the last few years, aimed to offer new parameters, which could be used as neuroimaging biomarkers. Each technique has shown pros and cons, but no gold standard exists yet. The fractal dimension (FD) has been investigated as a novel potential objective parameter for monitoring intratumoral space-filling properties of SWI patterns. We showed that SWI patterns found in different tumors or different glioma grades can be represented by a gradient in the fractal dimension, thereby enabling each tumor to be assigned a specific SWI fingerprint. Such results were especially relevant in the differentiation of low-grade versus high-grade gliomas, as well as from high-grade gliomas versus lymphomas.Therefore, FD has been suggested as a potential image biomarker to analyze intrinsic neoplastic architecture in order to improve the differential diagnosis within clinical neuroimaging, determine appropriate therapy, and improve outcome in patients.These promising preliminary findings could be extended into the field of neurotraumatology, by means of the application of computational fractal-based analysis for the qualitative and quantitative imaging of microbleedings in traumatic brain injury patients. In consideration of some evidences showing that SWI signals are correlated with trauma clinical severity, FD might offer some objective prognostic biomarkers.In conclusion, fractal-based morphometrics of SWI could be further investigated to be used in a complementary way with other techniques, in order to form a holistic understanding of the temporal evolution of brain tumors and follow-up response to treatment, with several further applications in other fields, such as neurotraumatology and cerebrovascular neurosurgery as well.

Fractal Dimension Analysis in Neurological Disorders: An Overview

Fractal analysis has emerged as a powerful tool for characterizing irregular and complex patterns found in the nervous system. This characterization is typically applied by estimating the fractal dimension (FD), a scalar index that describes the topological complexity of the irregular components of the nervous system, both at the macroscopic and microscopic levels, that may be viewed as geometric fractals. Moreover, temporal properties of neurophysiological signals can also be interpreted as dynamic fractals. Given its sensitivity for detecting changes in brain morphology, FD has been explored as a clinically relevant marker of brain damage in several neuropsychiatric conditions as well as in normal and pathological cerebral aging. In this sense, evidence is accumulating for decreases in FD in Alzheimer's disease, frontotemporal dementia, Parkinson's disease, multiple sclerosis, and many other neurological disorders. In addition, it is becoming increasingly clear that fractal analysis in the field of clinical neurology opens the possibility of detecting structural alterations in the early stages of the disease, which highlights FD as a potential diagnostic and prognostic tool in clinical practice.

Fractals, Pattern Recognition, Memetics, and AI: A Personal Journal in the Computational Neurosurgery

In this chapter, the personal journey of the author in many countries, including Italy, Germany, Austria, the United Kingdom, Switzerland, the United States, Canada, and Australia, is summarized, aimed to merge different translational fields (such as neurosurgery and the clinical neurosciences in general, biomedical engineering, mathematics, computer science, and cognitive sciences) and lay the foundations of a new field defined computational neurosurgery, with fractals, pattern recognition, memetics, and artificial intelligence as the common key words of the journey.

Computational Fractal-Based Neurosurgery

Fractal geometry is a branch of mathematics used to characterize and quantify the geometrical complexity of natural objects, with many applications in different fields, including physics, astronomy, geology, meteorology, finances, social sciences, and computer graphics. In the biomedical sciences, the use of fractal parameters has allowed the introduction of novel morphometric parameters, which have been shown to be useful to characterize any biomedical images as well as any time series within different domains of applications. Specifically, in the neurosciences and neurosurgery, the use of the fractal dimension and other computationally inferred fractal parameters has offered robust morphometric quantitators to characterize the brain in its wholeness, from neurons to the cortical structure and connections, and introduced new prognostic, diagnostic, and eventually therapeutic markers of many diseases of neurosurgical interest, including brain tumors and cerebrovascular malformations, as summarized in this chapter.

Fractal dimension: analyzing its potential as a neuroimaging biomarker for brain tumor diagnosis using machine learning

Purpose: The main purpose of this study was to comprehensively investigate the potential of fractal dimension (FD) measures in discriminating brain gliomas into low-grade glioma (LGG) and high-grade glioma (HGG) by examining tumor constituents and non-tumorous gray matter (GM) and white matter (WM) regions. Methods: Retrospective magnetic resonance imaging (MRI) data of 42 glioma patients (LGG, n = 27 and HGG, n = 15) were used in this study. Using MRI, we calculated different FD measures based on the general structure, boundary, and skeleton aspects of the tumorous and non-tumorous brain GM and WM regions. Texture features, namely, angular second moment, contrast, inverse difference moment, correlation, and entropy, were also measured in the tumorous and non-tumorous regions. The efficacy of FD features was assessed by comparing them with texture features. Statistical inference and machine learning approaches were used on the aforementioned measures to distinguish LGG and HGG patients. Results: FD measures from tumorous and non-tumorous regions were able to distinguish LGG and HGG patients. Among the 15 different FD measures, the general structure FD values of enhanced tumor regions yielded high accuracy (93%), sensitivity (97%), specificity (98%), and area under the receiver operating characteristic curve (AUC) score (98%). Non-tumorous GM skeleton FD values also yielded good accuracy (83.3%), sensitivity (100%), specificity (60%), and AUC score (80%) in classifying the tumor grades. These measures were also found to be significantly (p < 0.05) different between LGG and HGG patients. On the other hand, among the 25 texture features, enhanced tumor region features, namely, contrast, correlation, and entropy, revealed significant differences between LGG and HGG. In machine learning, the enhanced tumor region texture features yielded high accuracy, sensitivity, specificity, and AUC score. Conclusion: A comparison between texture and FD features revealed that FD analysis on different aspects of the tumorous and non-tumorous components not only distinguished LGG and HGG patients with high statistical significance and classification accuracy but also provided better insights into glioma grade classification. Therefore, FD features can serve as potential neuroimaging biomarkers for glioma.

Foundations of Multiparametric Brain Tumour Imaging Characterisation Using Machine Learning

The heterogeneity of brain tumours at the molecular, metabolic and structural levels poses significant challenge for accurate tissue characterisation. Artificial intelligence and radiomics have emerged as valuable tools to analyse quantitative features extracted from medical images which capture the complex microenvironment of brain tumours. In particular, a number of computational tools including machine learning algorithms have been proposed for image preprocessing, tumour segmentation, feature extraction, classification, and prognostic stratifications as well. In this chapter, we explore the fundamentals of multiparametric brain tumour characterisation, as an understanding of the strengths, limitations and applications of these tools allows clinicians to better develop and evaluate models with improved diagnostic and prognostic value in brain tumour patients.

Shape matters: morphological metrics of glioblastoma imaging abnormalities as biomarkers of prognosis

Lacunarity, a quantitative morphological measure of how shapes fill space, and fractal dimension, a morphological measure of the complexity of pixel arrangement, have shown relationships with outcome across a variety of cancers. However, the application of these metrics to glioblastoma (GBM), a very aggressive primary brain tumor, has not been fully explored. In this project, we computed lacunarity and fractal dimension values for GBM-induced abnormalities on clinically standard magnetic resonance imaging (MRI). In our patient cohort (n = 402), we connect these morphological metrics calculated on pretreatment MRI with the survival of patients with GBM. We calculated lacunarity and fractal dimension on necrotic regions (n = 390), all abnormalities present on T1Gd MRI (n = 402), and abnormalities present on T2/FLAIR MRI (n = 257). We also explored the relationship between these metrics and age at diagnosis, as well as abnormality volume. We found statistically significant relationships to outcome for all three imaging regions that we tested, with the shape of T2/FLAIR abnormalities that are typically associated with edema showing the strongest relationship with overall survival. This link between morphological and survival metrics could be driven by underlying biological phenomena, tumor location or microenvironmental factors that should be further explored.

Spherical coordinates transformation pre-processing in Deep Convolution Neural Networks for brain tumor segmentation in MRI

Magnetic Resonance Imaging (MRI) is used in everyday clinical practice to assess brain tumors. Deep Convolutional Neural Networks (DCNN) have recently shown very promising results in brain tumor segmentation tasks; however, DCNN models fail the task when applied to volumes that are different from the training dataset. One of the reasons is due to the lack of data standardization to adjust for different models and MR machines. In this work, a 3D spherical coordinates transform during the pre-processing phase has been hypothesized to improve DCNN models' accuracy and to allow more generalizable results even when the model is trained on small and heterogeneous datasets and translated into different domains. Indeed, the spherical coordinate system avoids several standardization issues since it works independently of resolution and imaging settings. The model trained on spherical transform pre-processed inputs resulted in superior performance over the Cartesian-input trained model on predicting gliomas' segmentation on Tumor Core and Enhancing Tumor classes, achieving a further improvement in accuracy by merging the two models together. The proposed model is not resolution-dependent, thus improving segmentation accuracy and theoretically solving some transfer learning problems related to the domain shifting, at least in terms of image resolution in the datasets.

The diagnostic value of postcontrast susceptibility-weighted imaging in the assessment of intracranial brain neoplasm at 3T

BACKGROUND

Susceptibility-weighted imaging (SWI) is occasionally performed with intravenous gadolinium (Gd). It was reported that SWI can be performed after Gd injection without information loss or signal change.

PURPOSE

To investigate the diagnostic value of contrast-enhanced SWI (CE-SWI) in the assessment of intracranial brain neoplasm.

MATERIAL AND METHODS

After obtaining the approval of the local ethics committee, 35 brain neoplasm patients (24 with metastasis and 11 with glioblastoma multiforme [GBM]) were enrolled. In order to investigate the value of using CE-SWI, two neuroradiologists performed an evaluation of the frequency of the intratumoral susceptibility signals (ITSS) in SWI and CE-SWI with visual assessment using 5-grade scales. We evaluated the visibility of the tumor margins and the internal architecture of tumors on T1-weighted imaging (T1WI), contrast-enhanced T1 (CE-T1), SWI, and CE-SWI.

RESULTS

The contrast-enhanced scans (CE-T1 and CE-SWI) showed statistically significant higher scores compared to non-enhanced scans (T1WI and SWI) for the analysis of the tumor margin in GBM and metastasis (P < 0.05, Wilcoxon signed rank test). Statistically significant higher scores are noted in GBMs compared to metastases in the visibility of the internal architecture of tumors on CE-SWI and the visibility of the tumor margin on CE-T1 (P < 0.05, Mann-Whitney test).

CONCLUSION

Based on our results, SWI can be performed after gadolinium injection without information loss or signal change. CE-SWI is useful in evaluating intracranial neoplasm due to its ability to simultaneously demonstrate both ITSS that are not visible with conventional magnetic resonance sequences and contrast enhancement.

Radiomic Based Machine Learning Performance for a Three Class Problem in Neuro-Oncology: Time to Test the Waters?

Prior radiomics studies have focused on two-class brain tumor classification, which limits generalizability. The performance of radiomics in differentiating the three most common malignant brain tumors (glioblastoma (GBM), primary central nervous system lymphoma (PCNSL), and metastatic disease) is assessed; factors affecting the model performance and usefulness of a single sequence versus multiparametric MRI (MP-MRI) remain largely unaddressed. This retrospective study included 253 patients (120 metastatic (lung and brain), 40 PCNSL, and 93 GBM). Radiomic features were extracted for whole a tumor mask (enhancing plus necrotic) and an edema mask (first pipeline), as well as for separate enhancing and necrotic and edema masks (second pipeline). Model performance was evaluated using MP-MRI, individual sequences, and the T1 contrast enhanced (T1-CE) sequence without the edema mask across 45 model/feature selection combinations. The second pipeline showed significantly high performance across all combinations (Brier score: 0.311-0.325). GBRM fit using the full feature set from the T1-CE sequence was the best model. The majority of the top models were built using a full feature set and inbuilt feature selection. No significant difference was seen between the top-performing models for MP-MRI (AUC 0.910) and T1-CE sequence with (AUC 0.908) and without edema masks (AUC 0.894). T1-CE is the single best sequence with comparable performance to that of multiparametric MRI (MP-MRI). Model performance varies based on tumor subregion and the combination of model/feature selection methods.

Brain volumetric and fractal analysis of synthetic MRI: A comparative study with conventional 3D T1-weighted images

PURPOSE

The estimation of brain volumetric measurements based on Synthetic MRI (SyMRI) is easy and fast, however, the consistency of brain volumetric and morphologic measurements based on SyMRI and 3D T1WI should be further addressed. The current study evaluated the impact of spatial resolution on brain volumetric and morphologic measurements using SyMRI, and test whether the brain measurements derived from SyMRI were consistent with those resulted from 3D T1WI.

METHOD

Brain volumetric and fractal analysis were applied to thirty healthy subjects, each underwent four SyMRI acquisitions with different spatial resolutions (1 × 1 × 2 mm, 1 × 1x3mm, 1 × 1 × 4 mm, 2 × 2 × 2 mm) and a 3D T1WI (1 × 1 × 1 mm isotropic). The consistency of the SyMRI measurements was tested using one-way non-parametric Kruskal-Wallis test and post hoc Dwass-Steel-Critchlow-Fligner test. The association between SyMRI and 3D T1WI derived measurements was evaluated using linear regression models.

RESULTS

Our results demonstrated that both in- and through-plane resolutions show an impact on brain volumetric measurements, while brain parenchymal volume showed high consistency across the SyMRI acquisitions, and high association with the measurements from 3D T1WI. In addition, SyMRI with 1 × 1 × 4 mm resolution showed the strongest association with 3D T1WI compared to other SyMRI acquisitions in both volumetric and fractal analyses. Moreover, substantial differences were found in fractal dimension of both gray and white matter between the SyMRI and 3D T1WI tissue segmentations.

CONCLUSIONS

Our results suggested that the measurements from SyMRI with relatively higher in-plane and lower through-plane resolution (1 × 1 × 4 mm) are much closer to 3D T1WI.

Radiomics-Based Differentiation between Glioblastoma, CNS Lymphoma, and Brain Metastases: Comparing Performance across MRI Sequences and Machine Learning Models

Prior radiomics studies have focused on two-class brain tumor classification, which limits generalizability. The performance of radiomics in differentiating the three most common malignant brain tumors (glioblastoma (GBM), primary central nervous system lymphoma (PCNSL), and metastatic disease) is assessed; factors affecting the model performance and usefulness of a single sequence versus multiparametric MRI (MP-MRI) remain largely unaddressed. This retrospective study included 253 patients (120 metastatic (lung and brain), 40 PCNSL, and 93 GBM). Radiomic features were extracted for whole a tumor mask (enhancing plus necrotic) and an edema mask (first pipeline), as well as for separate enhancing and necrotic and edema masks (second pipeline). Model performance was evaluated using MP-MRI, individual sequences, and the T1 contrast enhanced (T1-CE) sequence without the edema mask across 45 model/feature selection combinations. The second pipeline showed significantly high performance across all combinations (Brier score: 0.311–0.325). GBRM fit using the full feature set from the T1-CE sequence was the best model. The majority of the top models were built using a full feature set and inbuilt feature selection. No significant difference was seen between the top-performing models for MP-MRI (AUC 0.910) and T1-CE sequence with (AUC 0.908) and without edema masks (AUC 0.894). T1-CE is the single best sequence with comparable performance to that of multiparametric MRI (MP-MRI). Model performance varies based on tumor subregion and the combination of model/feature selection methods.

Machine Learning for the Prediction of Molecular Markers in Glioma on Magnetic Resonance Imaging: A Systematic Review and Meta-Analysis

BACKGROUND

Molecular characterization of glioma has implications for prognosis, treatment planning, and prediction of treatment response. Current histopathology is limited by intratumoral heterogeneity and variability in detection methods. Advances in computational techniques have led to interest in mining quantitative imaging features to noninvasively detect genetic mutations.

OBJECTIVE

To evaluate the diagnostic accuracy of machine learning (ML) models in molecular subtyping gliomas on preoperative magnetic resonance imaging (MRI).

METHODS

A systematic search was performed following PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analysis) guidelines to identify studies up to April 1, 2020. Methodological quality of studies was assessed using the Quality Assessment for Diagnostic Accuracy Studies (QUADAS)-2. Diagnostic performance estimates were obtained using a bivariate model and heterogeneity was explored using metaregression.

RESULTS

Forty-four original articles were included. The pooled sensitivity and specificity for predicting isocitrate dehydrogenase (IDH) mutation in training datasets were 0.88 (95% CI 0.83-0.91) and 0.86 (95% CI 0.79-0.91), respectively, and 0.83 to 0.85 in validation sets. Use of data augmentation and MRI sequence type were weakly associated with heterogeneity. Both O6-methylguanine-DNA methyltransferase (MGMT) gene promoter methylation and 1p/19q codeletion could be predicted with a pooled sensitivity and specificity between 0.76 and 0.83 in training datasets.

CONCLUSION

ML application to preoperative MRI demonstrated promising results for predicting IDH mutation, MGMT methylation, and 1p/19q codeletion in glioma. Optimized ML models could lead to a noninvasive, objective tool that captures molecular information important for clinical decision making. Future studies should use multicenter data, external validation and investigate clinical feasibility of ML models.

Intratumoral Susceptibility Signals Reflect Biomarker Status in Gliomas

Susceptibility-weighted imaging (SWI) can be a useful tool to depict vascular structures in brain tumors as well as micro-bleedings, which represent tumor invasion to blood vessels and could also be representative of tumoral angiogenesis. In this study, we investigated the relationship between SWI features and glioma grades, and the expression of key molecular markers isocitrate dehydrogenase 1 (IDH1), O-6-methylguanine-DNA methyltransferase (MGMT), and 1p19q. The gliomas were graded according to the intratumoral susceptibility signals (ITSS). We used the Mann-Whitney test to analyze the relationship between ITSS grades and the pathological level and status of these markers. Additionally, the area under the curve (AUC) was used to determine the predictive value of glioma SWI characteristics for the molecular marker status. In these cases, the ITSS grades of low-grade gliomas (LGG) were significantly lower than those of high-grade gliomas (HGG). Similarly, the ITSS grades of gliomas with IDH1 mutations and MGMT methylation were significantly lower than those of gliomas with Wild-type IDH1 and unmethylated MGMT. However, ITSS grades showed no relationship with 1p19q deletion status, while they did show significant predictive ability for glioma grade, IDH1 mutation, and MGMT methylation. These findings indicate an association between some molecular markers and cerebral microbleeds in gliomas, providing a new avenue for non-invasive prediction of molecular genetics in gliomas and an important basis for preoperative personalized surgical treatment based on molecular pathology.

Clinical potential of meningioma genomic insights: a practical review for neurosurgeons

Meningiomas are among the most common intracranial pathological conditions, accounting for 36% of intracranial lesions treated by neurosurgeons. Although the majority of these lesions are benign, the classical categorization of tumors by histological type or World Health Organization (WHO) grade has not fully captured the potential for meningioma progression and recurrence. Many targeted treatments have failed to generate a long-lasting effect on these tumors. Recently, several seminal studies evaluating the genomics of intracranial meningiomas have rapidly changed the understanding of the disease. The importance of NF2 (neurofibromin 2), TRAF7 (tumor necrosis factor [TNF] receptor-associated factor 7), KLF4 (Kruppel-like factor 4), AKT1, SMO (smoothened), PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha), and POLR2 (RNA polymerase II subunit A) demonstrates that there are at least 6 distinct mutational classes of meningiomas. In addition, 6 methylation classes of meningioma have been appreciated, enabling improved prediction of prognosis compared with traditional WHO grades. Genomic studies have shed light on the nature of recurrent meningioma, distinct intracranial locations and mutational patterns, and a potential embryonic cancer stem cell-like origin. However, despite these exciting findings, the clinical relevance of these findings remains elusive. The authors review the key findings from recent genomic studies in meningiomas, specifically focusing on how these findings relate to clinical insights for the practicing neurosurgeon.

A healthy dose of chaos: Using fractal frameworks for engineering higher-fidelity biomedical systems

Optimal levels of chaos and fractality are distinctly associated with physiological health and function in natural systems. Chaos is a type of nonlinear dynamics that tends to exhibit seemingly random structures, whereas fractality is a measure of the extent of organization underlying such structures. Growing bodies of work are demonstrating both the importance of chaotic dynamics for proper function of natural systems, as well as the suitability of fractal mathematics for characterizing these systems. Here, we review how measures of fractality that quantify the dose of chaos may reflect the state of health across various biological systems, including: brain, skeletal muscle, eyes and vision, lungs, kidneys, tumours, cell regulation, skin and wound repair, bone, vasculature, and the heart. We compare how reports of either too little or too much chaos and fractal complexity can be damaging to normal biological function, and suggest that aiming for the healthy dose of chaos may be an effective strategy for various biomedical applications. We also discuss rising examples of the implementation of fractal theory in designing novel materials, biomedical devices, diagnostics, and clinical therapies. Finally, we explain important mathematical concepts of fractals and chaos, such as fractal dimension, criticality, bifurcation, and iteration, and how they are related to biology. Overall, we promote the effectiveness of fractals in characterizing natural systems, and suggest moving towards using fractal frameworks as a basis for the research and development of better tools for the future of biomedical engineering.

Computational quantitative MR image features - a potential useful tool in differentiating glioblastoma from solitary brain metastasis

PURPOSE

Glioblastomas (GBM) and metastases are the most frequent malignant brain tumors in the adult population. Their presentation on conventional MRI is quite similar, but treatment strategy and prognosis are substantially different. Even with advanced MR techniques, in some cases diagnostic uncertainty remains. The main objective of this study was to determine whether fractal, texture, or both MR image analyses could aid in differentiating glioblastoma from solitary brain metastasis.

METHOD

In a retrospective study of 55 patients (30 glioblastomas and 25 solitary metastases) who underwent T2W/SWI/CET1 MRI, quantitative parameters of fractal and texture analysis were estimated, using box-counting and gray level co-occurrence matrix (GLCM) methods.

RESULTS

All five GLCM parameters obtained from T2W images showed significant difference between glioblastomas and solitary metastases, as well as on CET1 images except correlation (S_COR), contrary to SWI images which showed different values of two parameters (angular second moment-S_ASM and contrast-S_CON). Only three fractal features (binary box dimension-D_bin, normalized box dimension-D_norm and lacunarity-λ) measured on T2W and D_norm measured on CET1 images significantly differed GBMs from solitary metastases. The highest sensitivity and specificity were obtained from inverse difference moment (S_IDM) on T2W and S_IDM on CET1 images, respectively. Combination of several GLCM parameters yielded better results. The processing of T2W images provided the most significantly different parameters between the groups, followed by CET1 and SWI images.

CONCLUSIONS

Computational-aided quantitative image analysis may potentially improve diagnostic accuracy. According to our results texture features are more significant than fractal-based features in differentiation glioblastoma from solitary metastasis.

Multi-parametric effect in predicting tumor histological grade by using susceptibility weighted magnetic resonance imaging in tongue squamous cell carcinoma

BACKGROUND

Susceptibility weighted imaging (SWI) is helpful for depicting hemorrhage, calcification, and increased vascularity in some neoplasms, which may reflect tumor grade. In this study, we aimed to apply SWI in patients with oral tongue squamous cell carcinomas (OTSCCs) and relate multi-parametric effect to tumor histological grade prediction.

METHODS

Preoperative MR examinations were performed on a 1 .5T MRI scanner with T1-, T2- and contrast-enhanced (CE) T1-weighted imaging. In addition to routine head and neck MRI sequences, SWI was performed. Tumor thickness and volume were measured. Intratumoral susceptibility signal intensities (ITSSs), ITSS score and ITSS ratio on SWI were evaluated and recorded. Subjects were sub-grouped into low- and high-grade according to the histological findings post operation. Parameters such as tumor thickness, tumor volume and three ITSS related parameters were compared between low- and high-grade groups. ROC analysis was performed on above parameters to access the capability in predicting tumor histological grade. Different multi-parametric models were run to access multi-parametric combination effect.

RESULTS

Thirty patients with OTSCC were finally included in the study. Twenty of them were categorized as low-grade SCC and the other ten subjects were high-grade SCC according to the pathologic findings. No significant difference was seen for tumor thickness or tumor volume between two sub-groups. ITSSs were seen in 23/30 patients. Significant difference of ITSS scores between low- and high-grade OTSCCs was observed, with mean value of 0.95 ± 0.83 and 1.70 ± 0.95, respectively. Univariate ROC analysis demonstrated ITSSs, ITSS score and ITSS ratio were valuable parameters for predicting tumor histological grade and ITSSs was superior to the other two parameters, with an area under ROC curve of 0.790. Multi-parametric model using combination of ITSSs and tumor thickness would greatly improve the predictive capability in comparison with a univariate approach, yielding the area under ROC curve of 0.84(0.69,0.99). On contrast-enhanced SWI (CE-SWI), ITSSs were shown more clearly delineated in comparison with non-contrast enhanced SWI.

CONCLUSIONS

In conclusion, SWI was superior in depiction of internal characteristics of OTSCCs, which would potentially provide more diagnostic information. Multi-parametric model using combination of ITSSs and tumor thickness would be valuable in predicting tumor histological grade.

Relationship between necrotic patterns in glioblastoma and patient survival: fractal dimension and lacunarity analyses using magnetic resonance imaging

Necrosis is a hallmark feature of glioblastoma (GBM). This study investigated the prognostic role of necrotic patterns in GBM using fractal dimension (FD) and lacunarity analyses of magnetic resonance imaging (MRI) data and evaluated the role of lacunarity in the biological processes leading to necrosis. We retrospectively reviewed clinical and MRI data of 95 patients with GBM. FD and lacunarity of the necrosis on MRI were calculated by fractal analysis and subjected to survival analysis. We also performed gene ontology analysis in 32 patients with available RNA-seq data. Univariate analysis revealed that FD < 1.56 and lacunarity > 0.46 significantly correlated with poor progression-free survival (p = 0.006 and p = 0.012, respectively) and overall survival (p = 0.008 and p = 0.005, respectively). Multivariate analysis revealed that both parameters were independent factors for unfavorable progression-free survival (p = 0.001 and p = 0.015, respectively) and overall survival (p = 0.002 and p = 0.007, respectively). Gene ontology analysis revealed that genes positively correlated with lacunarity were involved in the suppression of apoptosis and necrosis-associated biological processes. We demonstrate that the fractal parameters of necrosis in GBM can predict patient survival and are associated with the biological processes of tumor necrosis.