Skip to main content

ORIGINAL RESEARCH article

Front. Med., 12 November 2021
Sec. Precision Medicine
This article is part of the Research Topic Intelligent Analysis of Biomedical Imaging Data for Precision Medicine View all 19 articles

Differentiation of Brain Abscess From Cystic Glioma Using Conventional MRI Based on Deep Transfer Learning Features and Hand-Crafted Radiomics Features

\nLinlin Bo&#x;Linlin Bo1Zijian Zhang&#x;Zijian Zhang2Zekun JiangZekun Jiang1Chao YangChao Yang1Pu HuangPu Huang1Tingyin ChenTingyin Chen3Yifan WangYifan Wang3Gang YuGang Yu1Xiao TanXiao Tan2Quan Cheng,
Quan Cheng4,5*Dengwang Li
Dengwang Li1*Zhixiong Liu,
Zhixiong Liu4,5*
  • 1Shandong Key Laboratory of Medical Physics and Image Processing, Shandong Institute of Industrial Technology for Health Sciences and Precision Medicine, School of Physics and Electronics, Shandong Normal University, Jinan, China
  • 2Department of Oncology, Xiangya Hospital, Central South University, Changsha, China
  • 3Department of Network Information Center, Xiangya Hospital, Centra South University, Changsha, China
  • 4Department of Neurosurgery, Xiangya Hospital, Central South University, Changsha, China
  • 5National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, Changsha, China

Objectives: To develop and validate the model for distinguishing brain abscess from cystic glioma by combining deep transfer learning (DTL) features and hand-crafted radiomics (HCR) features in conventional T1-weighted imaging (T1WI) and T2-weighted imaging (T2WI).

Methods: This single-center retrospective analysis involved 188 patients with pathologically proven brain abscess (102) or cystic glioma (86). One thousand DTL and 105 HCR features were extracted from the T1WI and T2WI of the patients. Three feature selection methods and four classifiers, such as k-nearest neighbors (KNN), random forest classifier (RFC), logistic regression (LR), and support vector machine (SVM), for distinguishing brain abscess from cystic glioma were compared. The best feature combination and classifier were chosen according to the quantitative metrics including area under the curve (AUC), Youden Index, and accuracy.

Results: In most cases, deep learning-based radiomics (DLR) features, i.e., DTL features combined with HCR features, contributed to a higher accuracy than HCR and DTL features alone for distinguishing brain abscesses from cystic gliomas. The AUC values of the model established, based on the DLR features in T2WI, were 0.86 (95% CI: 0.81, 0.91) in the training cohort and 0.85 (95% CI: 0.75, 0.95) in the test cohort, respectively.

Conclusions: The model established with the DLR features can distinguish brain abscess from cystic glioma efficiently, providing a useful, inexpensive, convenient, and non-invasive method for differential diagnosis. This is the first time that conventional MRI radiomics is applied to identify these diseases. Also, the combination of HCR and DTL features can lead to get impressive performance.

Introduction

Brain glioma is the most common intracranial brain tumor that is extremely difficult to treat. Currently, surgical resection is the standard treatment of resectable diseases, followed by postoperative radiotherapy and chemotherapy (1). The majority of gliomas are solid tumors, but some present cystic changes, such as cystic glioma, which has different clinicopathological features from other tumors. Brain abscess is an infectious disease that has high morbidity and mortality (2, 3). Though the treatment and prognosis of these two diseases are different, accurate and timely differential diagnosis is crucial. In many cases, CT and MR images lack specificity for cystic glioma and brain abscess, especially when the medical history and clinical manifestations of the diseases cannot provide a differential diagnosis for timely treatment measures. At present, the two diseases are mainly distinguished by pathological examination, with the caveat of invasive procedure and intra-operator variability. To accurately distinguish the two diseases, previous studies have proposed advanced MR images diagnosis techniques (2, 4), such as susceptibility-weighted imaging and apparent diffusion coefficients (ADC). However, these diagnosis techniques cannot obtain high accuracy, and they rely on the experience of radiologists (5). The use of the most rudimentary imaging modalities of T1-weighted imaging (T1WI) and T2-weighted imaging (T2WI) for a training model with a large sample size contributes to more universality and fewer errors.

As a method of machine learning, radiomics is used for quantitative image feature extraction from tumor regions of interest. It has great potential for oncology practice, including differential diagnosis, prediction of pathological classification, lymph node metastasis, and survival (610). Radiomics has been applied to brain tumor diseases (1115), especially in differentiating brain tumors (1621). For example, Qian et al. investigated the ability of radiomic analysis to distinguish between isolated brain metastases and glioblastoma (16); Dong et al. used the radiomic features derived from the areas of peripheral enhancing edema to differentiate glioblastoma from supratentorial single brain metastasis (17); Zhang et al. investigated the feasibility of contrast-enhanced T1WI radiomics features extracted by machine-learning algorithms to distinguish between low-grade oligodendroglioma and atypical anaplastic oligodendroglioma (18); Chen et al. applied radiomics analysis to distinguish between metastatic brain tumors and glioblastomas based on contrast-enhanced T1WI, and they validated the discriminative performance of this method (19); Artzi et al. used radiomics-based machine learning to differentiate between brain metastasis subtypes and glioblastoma based on conventional postcontrast T1WI (20). However, the radiomics features are mainly the texture, size, volume, shape, and intensity characteristics of the tumor, limiting the potential of this method. Therefore, extracting more complex features and fusing them with radiomic features may improve the prediction and generalization capabilities of the model (2123).

In recent years, deep convolutional neural networks (CNNs) (24) with complex network structures have achieved remarkable results in the field of computer vision, such as tumor grade prediction, patient prognosis, pathology classification, and organ segmentation (25, 26). The successful application of deep learning requires a large number of training cohort sets. Since the available medical data sets have a limited size, a pretrained CNN known as “transfer learning” can be employed to avoid overfitting and replace deep learning in many practical applications (21, 27, 28).

It is not clear whether T1WI and T2WI, as conventional routine images in hospitals, also have diagnostic values for distinguishing brain abscess from cystic glioma. In this study, we hypothesized that conventional T1WI and T2WI would also be valuable in distinguishing between these two diseases. To this end, the DLR features extracted from patients with brain abscesses or cystic gliomas were used to validate the diagnostic capability of T1WI and T2WI.

Materials and Methods

Patients

This study was reviewed and approved by the Institutional Review Board of Xiangya Hospital and informed consent was provided by the patient participating in this study. From January 2017 to October 2020, 216 patients who met requirements and underwent T1WI and T2WI MRI were included in the cohort after an initial case screening. Twenty-eight patients were excluded due to poor MRI quality caused by technical operations or inspection processes. Finally, 188 patients were enrolled in this study, among which, 102 patients were diagnosed with brain abscesses (age [mean ± SD], 47.8 ± 17.6 years; 33 males and 69 females) and 86 were diagnosed with cystic gliomas (age [mean ± SD], 46.2 ± 15.1 years; 27 males and 59 females). The training cohort and test cohort were divided by stratified sampling with a ratio of 7:3, and the distribution of the two diseases was almost the same as that of the overall data set. Then, a nested 5-fold cross-validation was performed on the training cohort.

Figure 1 showed the flowchart of our study, consisting of image preprocessing, feature extraction, feature analysis, and model construction.

FIGURE 1
www.frontiersin.org

Figure 1. Study workflow overview. (A) Imaging processing; (B) Feature extraction; (C) Feature analysis and model construction.

Image Acquisition

All MRI examinations were conducted in the radiology department of Xiangya hospital with a 3.0T MR Scanner. High-quality MR images were obtained under the following configurations: ①axial T1WI: layer thickness = 5 mm, layer spacing = 1.5 mm, matrix = 512 × 416, field of view = 24 × 24 cm. ②axial T2WI: layer thickness = 5 mm, layer spacing = 1.5 mm, matrix = 416 × 512, field of view = 24 × 24 cm. All MR images were retrieved from the picture archiving and communication system for further image feature extraction.

Image Preprocessing and Tumor Segmentation

To convert T1WI images to the space of T2WI images, automatic rigid registration was performed with the ITK-SNAP software (Version 3.8.0, http://www.itksnap.org/) to segment the structures in the 3D medical image. Meanwhile, manual segmentation of the lesions of all subjects was performed on registered T2WI and T1WI images by a neuroradiologist with 10 years of experience. Then, a radiologist (with 10-years experience) segmented 50 cases, consisting of 25 pathologically proven brain abscess cases and 25 cystic glioma cases randomly selected from all samples. In this way, the consistency of the extracted HCR and DTL features of the neuroradiologist and radiologist was evaluated, and the impact of inter-operator variation on model stability and generalizability was reduced. Besides, the intra-class correlation coefficient of each feature of these 50 cases is calculated.

Image preprocessing was performed as follows. Before DTL features were extracted, an image of the largest cross-sectional area and its upper and lower layers were chosen as a three-channel image. Then, a rectangular region of interest around the tumor contour was used to crop the MR image. Next, the size of the tumor patch was resized to 224 × 224 to meet the input size requirement of the pretrained CNN model. Before the HCR features were extracted, B-spline interpolation was adopted to resample all images to the same voxel size of 1 × 1 × 1 mm3. To avoid the influence of different MR image machine scanners on feature extraction, all images were normalized. Moreover, it seems that deep learning is less affected by different MR machine types than radiomics.

Feature Extraction

In line with the Imaging Biomarker Standardization Initiative, two kinds of features were extracted, i.e., DTL feature and HCR feature. As for the extraction of the DTL feature, ResNet-50 (29) and VGG-19 (30) pretrained on the natural image dataset ImageNet (http://www.image-net.org/) were taken as our base models. (Visual Geometry Group) VGG-19 contains 19 hidden layers (16 convolution layers and 3 full connection layers). It uses 3 × 3 convolutional kernels in all layers to deepen the number of layers and avoid excessive parameters. As for ResNet, it integrates residual learning to avoid gradient dispersion and accuracy reduction in deep networks, improving the network efficiency, accuracy, and execution speed. The internal deep learning features in the image are also visualized while the convolutional layer receives the input features and generates the output feature mapping. As for the extraction of the HCR feature, 105 original HCR features were extracted from each of the axial T1WI and T2WI images using PyRadiomics (Version 2.1.0, https://pyradiomics.readthedocs.io/).

Feature Selection

To prevent overfitting, the multistep feature selection method was adopted to select the best features for distinguishing brain abscesses from cystic gliomas. First, all the HCR features were analyzed in order by the Spearman rank correlation test and mutual information method. The Spearman rank correlation test was used to investigate the internal linear correlation between individual features. The higher the absolute value of the correlation coefficient, the stronger the correlation. As for the non-redundant features (with a linear correlation coefficient <0.95), the mutual information method was used to capture arbitrary relationships (both linear and non-linear) between each feature and object variable. Then, according to the imaging modality and feature category, the remaining radiomics features and all DLR features were divided into four feature groups, i.e., T1WI-HCR, T2WI-HCR, T1WI-DTL, and T2WI-DTL group. Finally, the least absolute shrinkage and selection operator (LASSO) and recursive feature elimination (RFE) method based on LR and SVM were adopted to repeatedly create the model and select the best feature subset in each feature group.

Feature Fusion

Feature fusion indicates that two feature groups are put together. The T1WI-HCR and T2WI-HCR feature groups were fused to combined HCR (comb-HCR), and 10 optimal MR image features were selected; similarly, T1WI-DTL and T2WI-DTL were fused to combined DTL (comb-DTL); the T1WI-HCR and T1WI-DTL feature groups were fused to T1WI-DLR, and 10 optimal MR image features were selected; the T2WI-HCR and T2WI-DTL feature groups were fused to T2WI-DLR, and 10 optimal MR image features were selected; the comb-HCR and comb-DTL feature groups were fused to combined DLR (comb-DLR), and 10 optimal MR image features were selected. Refer Figure 2 for details of feature selection and fusion flow chart.

FIGURE 2
www.frontiersin.org

Figure 2. Feature selection and fusion flow chart. The red box is the original extracted feature. The black box is the feature group. Each feature group contains 10 features. A new feature group was formed after feature fusion and screening of two feature groups.

Feature Analysis/Model Construction and Validation

After feature fusion and selection, we used each feature group separately to build machine learning classification models, including LR, RFC, KNN, and SVM, implemented by Python Scikit-learn (https://scikit-learn.org/stable/user_guide.html). The performance of different classifiers was compared. To prevent overfitting, we performed 1,000 iterations of nested 5-fold cross-validation to select the best parameters for the classifier in the training cohorts. The discriminative power of the model was assessed by AUC values, Youden Index, and receiver operating characteristic (ROC) curves. Accuracy, precision, recall, specificity, and F1-score were also used as quantitative metrics. The AUC values for comparative disease identification were carried out using DeLong test.

Clinical ADC Maps vs. Our Model

In this study, two experienced radiologists (with more than 10 years of experience in brain tumor MRI) were assigned to jointly perform image ADC diagnosis, but they were not involved in the quantitative image analysis described above. All clinicopathological information were removed, and the radiologists distinguished brain abscesses from cystic gliomas based on ADC images only. The proportion of ADC maps in our study cohort was counted. Besides, the diagnostic performance of the two radiologists on the same dataset following the current clinical practice (including the use of ADC maps) was compared to that of the established classifier.

Statistical Analysis

The comparison of categorical variables was performed through chi-square tests or Fisher tests, and the comparison between quantitative variables was performed through t-tests or Mann-Whitney U-test. Meanwhile, the Spearman rank correlation test was adopted to evaluate the correlation and executed in Python. A p < 0.05 (two-sided) indicates a significant difference in distinguishing cystic gliomas from brain abscesses. Statistical analysis was performed with IBM SPSS Statistics (version 25; IBM Corporation, Armonk, NY, USA), R (https://cran.r-project.org/), and Python (version 3.6.6, https://www.python.org). A pretrained CNN model was run using Keras with a Tensorflow backend (https://keras.io/applications/#Resnet-50 and https://keras.io/applications/#VGG-19).

Results

Patient Characteristics

A total of 131 and 57 patients were involved in the training and test cohort of this study, respectively. Specifically, the training cohort involved 71 patients with brain abscesses and 60 patients with gliomas, while the test cohort involved 31 patients with brain abscesses and 26 patients with gliomas. The patient characteristics are provided in Supplementary Table 1. The gold standard for distinguishing between brain abscess and cystic glioma was confirmed by pathologists through pathological examination.

Results of the Feature Extraction

To extract the DTL features, the tumor patch images were input to the pretrained CNN to extract 1,000 features from each MR image modality, and the extracted features were outputs from the last fully connected layer of VGG-19 and ResNet-50. The extracted HCR features included First Order Statistics (18 features), Shape-based (3D) (14 features), Gray Level Cooccurrence Matrix (22 features), Gray Level Run Length Matrix (16 features), Gray Level Size Zone Matrix (16 features), Neighboring Gray Tone Difference Matrix (5 features), and Gray Level Dependence Matrix (14 features).

Results of the Feature Selection and Fusion

After feature selection and fusion, only the features with intra-class correlation coefficients >0.95 were retained, indicating that these features are not affected by multiple tumor segmentation operators and present high reproducibility. For single image modality analysis, a feature selection method was used to preserve 10 optimal features in each group. In multimodality analysis, the two groups of each modality were combined, and 10 features were filtered out. Refer Supplementary Table 2 for detailed feature selection results.

DLR vs. DTL vs. HCR

The modeling effects of the combined modality feature with different category were compared. On the test cohort, the feature obtained from comb-DLR presented higher diagnostic accuracy than those from comb-HCR and comb-DTL for distinguishing brain abscess from cystic glioma. The AUC values of the models established with the comb-DLR, comb-HCR, and comb-DTL models features were 0.82, 0.79, and 0.76, respectively; the AUC values of the models established with the T2WI-DLR, T2WI-HCR, and T2WI-DTL features were 0.85, 0.80, and 0.71, respectively; the AUC values of the models established with the T1WI-DLR, T1WI-HCR, and T1WI-DTL features were 0.74, 0.77, and 0.80, respectively. Refer Figure 3, Table 1 for details of model comparison results.

FIGURE 3
www.frontiersin.org

Figure 3. ROC comparison. (A,B), comb-HCR vs. comb-DTL vs. comb-DLR in the training cohort (A) and test cohort (B). (C,D) T1WI-DLR vs. T2WI-DLR vs. comb-DLR in the training cohort (C) and test cohort (D).

TABLE 1
www.frontiersin.org

Table 1. Model representation.

Multimodality vs. Single Modality

Since the AUC value of the model established with the DLR features was statistically higher than that of the models established with the HCR and DTL features in most cases, the DLR features were used in the multimodal experiments. On the test cohort, the AUC value of the model established with the T2WI-DLR features was statistically higher than that of the model established with the T1WI-DLR and com-DLR features. Model comparison results are provided in Figure 3, Table 1.

Construction and Validation of the Final Model

As can be seen from Table 1, the optimal model was obtained by using the T2WI-DLR features combined with an SVM-based RFE feature selection method, and an SVM classifier. The AUC value of the model on the training and test cohort reached 0.86 (95% CI: 0.81, 0.91) and 0.85 (95% CI: 0.75, 0.95) for distinguishing brain abscess from cystic glioma, respectively. Figure 4 presents the performance of the optimal model. It can be seen from the figure that on the training set, the AUC values of the nested 5-fold were 0.86, 0.86, 0.92, 0.88, and 0.80, respectively. Besides, the standard deviation of the mean AUC value was 0.04, indicating that our model has good stability and robustness and reduces overfitting. The optimal cutoff value of the model was determined by Youden Index. On the test cohort, the sensitivity and specificity of the model were 73.1 and 87.1%, respectively, with an optimal critical value of 0.512 and a Jorden index of 0.601. The details of HCR and DTL feature selection and model construction are listed in Tables 2, 3, respectively. The selection of the hyperparameters for each model is listed in Supplementary Table 3. Moreover, the ROC curves for each model were compared, and the Delong test results are detailed in Supplementary Table 4. It can be seen from the table, that T2WI-DLR features (AUC, 0.85) are superior to T2WI-DTL (AUC, 0.71; P = 0.0058; DeLong test).

FIGURE 4
www.frontiersin.org

Figure 4. ROC curves of the optimal model (T2WI-DLR). (A) training cohort; (B) test cohort.

TABLE 2
www.frontiersin.org

Table 2. HCR model construction.

TABLE 3
www.frontiersin.org

Table 3. DTL model construction.

DTL Feature Visualization

As shown in Figure 5, the feature maps output by the last convolutional layer in the VGG-19 and Resnet-50 model are visualized. The feature map of the visually perceptible tumor region captures most of the details in the image. To a certain extent, it confirms the reliability of transfer learning for feature extraction. By visualizing the learning of the features in the image by these two networks, more insight into the working of the networks can be obtained, and the reasons why the disease may be correctly identified by transfer learning models can be understood.

FIGURE 5
www.frontiersin.org

Figure 5. Feature visualization. (A,D) The grayscale T1WI, T2WI, and the corresponding heat map are shown, and the red areas indicating greater weighting, with color bars on the right side of the plot. (B,E), the feature map of the last convolutional layer in the VGG-19 model. (C,F), the feature map of the last convolutional layer in the Resnet-50 model.

Clinical ADC Maps vs. Our Model

The cases containing ADC maps in the cohort were counted. There were 33 cases in total, including 21 cases of brain abscess and 12 cases of cystic glioma. A detailed comparison of the distinguishment performance between the clinical ADC maps diagnosis and our model is shown in Table 4. For each patient, the diagnosis of the radiologist and the model prediction are listed in Supplementary Table 5. It can be found that our model has the same accuracy as the clinical ADC assessment, so it has great potential for distinguishing between the two diseases.

TABLE 4
www.frontiersin.org

Table 4. Clinical ADC maps vs. our model.

Discussions

Accurate identification of brain abscesses and cystic gliomas is essential to planning appropriate treatment, assessing outcome, and future prognosis. However, due to the similarity in the conventional MR images, i.e., ring enhancement, it is difficult to distinguish between the two diseases. In this study, a deep learning-based statistical analysis method based on multistep feature selection and fusion was demonstrated and verified. The experimental results indicate that the method can be used to distinguish between brain abscess and cystic glioma in conventional T1WI and T2WI. The previous literature on disease prediction prognosis and classification differential diagnosis for quantitative image analysis has shown that deep learning contributes to better performance of radiomics analysis (21, 27, 31, 32). Our study demonstrates that by extracting DTL features with VGG-19, a model with excellent feature learning and feature representation abilities can be obtained. Besides, as shown in Figure 5, VGG-19 can better focus on the details of the tumor region than Resnet-50.

According to the feature selection results of the optimal model, two “good” HCR features were selected for statistical significance analysis. It can be seen from the box plots in Figure 6 that there is not much difference in the distribution of the features between brain abscess and cystic glioma. All corresponding p-values of the statistical tests for distinguishing the two diseases are presented in the figure. These results indicate that these two features have a good identification ability in this work, showing the reproducibility and usefulness of feature engineering. Besides, the two “good” HCR features are all texture features, which reflect the homogeneous phenomenon in the image, which once again demonstrates the superiority of texture features in distinguishing brain abscess from cystic glioma. Previous studies have also shown that texture features are highly predictive in many tasks, which is consistent with the results obtained in this study (33, 34).

FIGURE 6
www.frontiersin.org

Figure 6. The box plots of two “good” T2WI-HCR of T2WI-DLR with the results of the statistical test. (A) original_glszm_SmallAreaEmphasis; (B) original_gldm_DependenceNonUniformityNormalized. Brain abscess:0; Cystic glioma:1.

Compared with using T1WI and T2WI alone in DLR, the model based on combined modality does not achieve improved performance. This indicates that the single modality of T2WI is also a good predictor for distinguishing between brain abscesses and gliomas, and this is consistent with the fact that the T2 modality is more commonly used in imaging diagnosis of brain diseases in clinical practice. When the performance of the two models is stable and the results are complementary, model fusion can lead to better performance (35).

Previous studies (5, 3638) have demonstrated that the advanced MRI techniques, such as magnetic resonance spectroscopy, susceptibility-weighted imaging, ADC, and dynamic susceptibility contrast-enhanced, can distinguish brain abscesses from gliomas, but these techniques have some limitations. Refer Supplementary Table 6 for model performance comparison. First, the sample size of these techniques is small, and only a few cases of pyogenic abscess and glioblastoma are included, hindering the direct application of the results to daily clinical practice. Second, the model based on the combination of intralesional susceptibility signal and ADC achieved a good AUC value (38). However, the combination does not lead to obvious improvement of differential diagnosis, because only a small number of patients with abscess/glioblastoma show atypical high/low ADC. Finally, none of the techniques were based on conventional MR images, increasing the image acquisition difficulty and cost. In this study, some measures were taken to overcome these limitations. The relatively larger sample size contributes to a better performance of conventional MR images for distinguishing brain abscesses and gliomas. The comparison of dataset size is listed in Supplementary Table 6. Also, the HCR and DTL features were extracted from conventional MR images. To our knowledge, there is no report on integrating HCR and DTL features for distinguishing brain abscesses from cystic gliomas. Besides, our research is based on some ordinary image data and does not require special training, so it has significant potential. In addition, DTL feature extraction uses a fixed-size bounding box for the tumor region, which not only provides information about intertumoral heterogeneity but also provides tumor microenvironment information to a certain extent.

To promote the development of radiomics as imaging biomarkers, a plethora of studies have used radiomics quality score (RQS) to evaluate and standardize radiomics (8, 39, 40). The RQS of our study was satisfactory at 15 points (41.7% of the ideal quality score), and the detailed result is listed in Supplementary Table 7. The RQSs of the relevant work (1620) were analyzed in our study, but only our study is open to science and data, only one study conducted a multivariable analysis with non-radiomics features (20), and only one study based on multicenter validation (16). Besides, no research has conducted a phantom study, collected images of individuals at additional time points, discussed biological correlates, conducted a prospective study, or reported the cost-effectiveness of the clinical application.

The limitations of our work are as follows. First, due to the difficulty of obtaining external validation data, the patients in our study were single-center. The effects of clinical ADC maps diagnosis were compared with our proposed model, but this sequence was not added to our model due to the insufficient sample size of ADC maps. Multicenter validation, multi-MRI sequences, and prospective studies will be involved in our future work. Meanwhile, additional features such as proteomics, transcriptomics, pathomics, and genomic features were not considered in our study. Multi-omics joint analysis that integrates complex structural systems with multiple layers, levels, and functions may enhance the performance to identify brain abscess and cystic glioma and overcome the limitations of a single theoretical model. Besides, this study only used the image of the largest cross-section area with the upper and lower layers as the input to the VGG-19 or ResNet-50 model. The use of the 3D volume of the tumor/region of interest should be investigated in future research. Finally, the application of our study to identify other brain tumors and the enhancement of the algorithm will be explored.

Conclusions

This paper first reports a model combining DTL features and HCR features from conventional MRI for distinguishing brain abscesses from cystic glioma. The study results provide an effective, inexpensive, convenient, and non-invasive method for differential diagnosis.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, the main code are available from https://github.com/better0123/radiomics-analysis.git. Further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board (IRB) of Xiangya Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

LB, ZZ, and ZJ: study design. ZL, QC, ZZ, XT, YW, and TC: data collection. LB, PH, ZZ, and CY: data analysis. DL, GY, and ZL: supervision. LB, QC, ZZ, and ZJ: manuscript writing. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Nature Science Foundation of China (NO. 82073893, NO. 81873635, NO. 81703622 and NO. 81472693); the China Postdoctoral Science Foundation (NO. 2018M63302); the Natural Science Foundation of Hunan Province (NO. 2018JJ3838 and NO. 2018SK2101); the Hunan Provincial Health and Health Committee Foundation of China (C2019186); Xiangya Hospital Central South University postdoctoral foundation; the National Natural Science Foundation of China (61971271); the Taishan Scholars Project of Shandong Province (Tsqn20161023) and the Primary Research; Development Plan of Shandong Province (No. 2018GGX101018 and No. 2019QYTPY020).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank the Network and Information Center and Department of Radiology of Xiangya Hospital that have given support in data extraction and analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.748144/full#supplementary-material

Abbreviations

DTL, deep transfer learning; HCR, hand-crafted radiomics; T1WI, T1-weighted image; T2WI, T2-weighted image; KNN, k-nearest neighbors; RFC, random forest classifier; LR, logistic regression; SVM, support vector machine; AUC, Area Under the Curve; DLR, deep learning-based radiomics; ADC, apparent diffusion coefficients; CNN, convolutional neural networks; LASSO, least absolute shrinkage and selection operator; RFE, recursive feature elimination; comb-HCR, combined HCR; comb-DTL, combined DTL; comb-DLR, combined DLR; ROC, receiver operating characteristic; RQS, radiomics quality score.

References

1. Ganau L, Paris M, Ligarotti GK, Ganau M. Management of gliomas: overview of the latest technological advancements and related behavioral drawbacks. Behav Neurol. (2015) 2015:862634. doi: 10.1155/2015/862634

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Afridi M, Asadollahi S, Nath K, Chawla SJJNPN. Differentiation of brain abscesses from necrotic high-grade gliomas using advanced MR imaging techniques: a mini review. J Neuroimaging Psychiatry Neurol. (2019) 4:17–22. doi: 10.17756/jnpn.2019-030

CrossRef Full Text | Google Scholar

3. Alvis Miranda H, Castellar-Leones SM, Elzain MA, Moscote-Salazar LR. Brain abscess: current management. J Neurosci Rural Pract. (2013) 4(Suppl 1):S67–81. doi: 10.4103/0976-3147.116472

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Toh CH, Wei KC, Ng SH, Wan YL, Lin CP, Castillo M. Differentiation of brain abscesses from necrotic glioblastomas and cystic metastatic brain tumors with diffusion tensor imaging. Am J Neuroradiol. (2011) 32:1646. doi: 10.3174/ajnr.A2581

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Toh CH, Wei K-C, Chang C-N, Ng S-H, Wong H-F, Lin C.-P. Differentiation of brain abscesses from glioblastomas and metastatic brain tumors: comparisons of diagnostic performance of dynamic susceptibility contrast-enhanced perfusion MR imaging before and after mathematic contrast leakage correction. PLoS ONE. (2014) 9:e109172. doi: 10.1371/journal.pone.0109172

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Mayerhoefer ME, Materka A, Langs G, Häggström I, Szczypiński P, Gibbs P, et al. Introduction to radiomics. J Nucl Med. (2020) 61:488–95. doi: 10.2967/jnumed.118.222893

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Liu Z, Wang S, Di Dong JW, Fang C, Zhou X, Sun K, et al. The applications of radiomics in precision diagnosis and treatment of oncology: opportunities and challenges. Theranostics. (2019) 9:1303. doi: 10.7150/thno.30309

PubMed Abstract | CrossRef Full Text | Google Scholar

8. 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 Rev Clin Oncol. (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Giraud P, Giraud P, Gasnier A, El Ayachy R, Kreps S, Foy J-P, et al. Radiomics and machine learning for radiotherapy in head and neck cancers. Front Oncol. (2019) 9:174. doi: 10.3389/fonc.2019.00174

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Song J, Yin Y, Wang H, Chang Z, Liu Z, Cui L. A review of original articles published in the emerging field of radiomics. Eur J Radiol. (2020) 127:108991. doi: 10.1016/j.ejrad.2020.108991

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Xie H, Ma S, Wang X, Zhang XJ. Noncontrast computer tomography–based radiomics model for predicting intracerebral hemorrhage expansion: preliminary findings and comparison with conventional radiological model. Eur Radiol. (2020) 30:87–98. doi: 10.1007/s00330-019-06378-3

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Qiu Q, Duan J, Deng H, Han Z, Gu J, Yue NJ, et al. Development and validation of a radiomics nomogram model for predicting postoperative recurrence in patients with esophageal squamous cell cancer who achieved pCR after neoadjuvant chemoradiotherapy followed by surgery. Front Oncol. (2020) 10:1398. doi: 10.3389/fonc.2020.01398

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Dong F, Li Q, Xu D, Xiu W, Zeng Q, Zhu X, et al. Differentiation between pilocytic astrocytoma and glioblastoma: a decision tree model using contrast-enhanced magnetic resonance imaging-derived quantitative radiomic features. Eur Radiol. (2019) 29:3968–75. doi: 10.1007/s00330-018-5706-6

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Bae S, An C, Ahn SS, Kim H, Han K, Kim SW, et al. Robust performance of deep learning for distinguishing glioblastoma from single brain metastasis using radiomic features: model development and validation. Sci Rep. (2020) 10:12110. doi: 10.1038/s41598-020-68980-6

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zhang Z, Yang J, Ho A, Jiang W, Logan J, Wang X, et al. A predictive model for distinguishing radiation necrosis from tumour progression after gamma knife radiosurgery based on radiomic features from MR images. Eur Radiol. (2018) 28:2255–63. doi: 10.1007/s00330-017-5154-8

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Qian Z, Li Y, Wang Y, Li L, Li R, Wang K, et al. Differentiation of glioblastoma from solitary brain metastases using radiomic machine-learning classifiers. Cancer Lett. (2019) 451:128–35. doi: 10.1016/j.canlet.2019.02.054

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dong F, Li Q, Jiang B, Zhu X, Zeng Q, Huang P, et al. Differentiation of supratentorial single brain metastasis and glioblastoma by using peri-enhancing oedema region–derived radiomic features and multiple classifiers. Eur Radiol. (2020) 30:3015–22. doi: 10.1007/s00330-019-06460-w

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zhang Y, Chen C, Cheng Y, Teng Y, Guo W, Xu H, et al. Ability of radiomics in differentiation of anaplastic oligodendroglioma from atypical low-grade oligodendroglioma using machine-learning approach. Front Oncol. (2019) 9:1371. doi: 10.3389/fonc.2019.01371

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chen C, Ou X, Wang J, Guo W, Ma XJF. Radiomics-based machine learning in differentiation between glioblastoma and metastatic brain tumors. Front Oncol. (2019) 9:806. doi: 10.3389/fonc.2019.00806

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Artzi M, Bressler I, Ben Bashat D. Differentiation between glioblastoma, brain metastasis and subtypes using radiomics analysis. J Magn Reson Imaging. (2019) 50:519–28. doi: 10.1002/jmri.26643

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Han W, Qin L, Bay C, Chen X, Yu KH, Miskin N, et al. Deep transfer learning and radiomics feature prediction of survival of patients with high-grade gliomas. Am J Neuroradiol. (2020) 41:40. doi: 10.3174/ajnr.A6365

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Lao J, Chen Y, Li Z-C, Li Q, Zhang J, Liu J, et al. A deep learning-based radiomics model for prediction of survival in glioblastoma multiforme. Sci Rep. (2017) 7:10353. doi: 10.1038/s41598-017-10649-8

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wu X, Li Y, Chen X, Huang Y, He L, Zhao K, et al. Deep learning features improve the performance of a radiomics signature for predicting kras status in patients with colorectal cancer. Acad Radiol. (2020) 27:e254–62. doi: 10.1016/j.acra.2019.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Litjens G, Kooi T, Bejnordi BE, Setio AAA, Ciompi F, Ghafoorian M, et al. A survey on deep learning in medical image analysis. Med Image Anal. (2017) 42:60–88. doi: 10.1016/j.media.2017.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhang W, Yang G, Huang H, Yang W, Xu X, Liu Y, et al. ME-Net: multi-encoder net framework for brain tumor segmentation. Int J Imaging Syst Technol. (2021) 1–15. doi: 10.1002/ima.22571

CrossRef Full Text | Google Scholar

26. Jin Y, Yang G, Fang Y, Li R, Xu X, Liu Y, et al. 3D PBV-Net: an automated prostate MRI data segmentation method. Comput Biol Med. (2021) 128:104160. doi: 10.1016/j.compbiomed.2020.104160

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Xue L-Y, Jiang Z-Y, Fu T-T, Wang Q-M, Zhu Y-L, Dai M, et al. Transfer learning radiomics based on multimodal ultrasound imaging for staging liver fibrosis. Eur Radiol. (2020) 30:2973–83. doi: 10.1007/s00330-019-06595-w

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hu Y, Xie C, Yang H, Ho JWK, Wen J, Han L, et al. Computed tomography-based deep-learning prediction of neoadjuvant chemoradiotherapy treatment response in esophageal squamous cell carcinoma. Radiother Oncol. (2021) 154:6–13. doi: 10.1016/j.radonc.2020.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

29. He K, Zhang X, Ren S, Sun J. Deep residual learning for image recognition[C]. In: Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition. (2016) 770–8.

Google Scholar

30. Simonyan K, Zisserman A. Very deep convolutional networks for large-scale image recognition. arXiv [Preprint] arXiv:1409.1556. (2014).

Google Scholar

31. Yu J, Deng Y, Liu T, Zhou J, Jia X, Xiao T, et al. Lymph node metastasis prediction of papillary thyroid carcinoma based on transfer learning radiomics. Nat Commun. (2020) 11:4807. doi: 10.1038/s41467-020-18497-3

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhang Y, Lobo-Mueller EM, Karanicolas P, Gallinger S, Haider MA, Khalvati F. Improving prognostic performance in resectable pancreatic ductal adenocarcinoma using radiomics and deep learning features fusion in CT images. Sci Rep. (2021) 11:1378. doi: 10.1038/s41598-021-80998-y

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Tian Q, Yan LF, Zhang X, Zhang X, Hu YC, Han Y, et al. Radiomics strategy for glioma grading using texture features from multiparametric MRI. J Magn Reson Imaging. (2018) 48:1518–28. doi: 10.1002/jmri.26010

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhang X, Xu X, Tian Q, Li B, Wu Y, Yang Z, et al. Radiomics assessment of bladder cancer grade using texture features from diffusion-weighted imaging. J Magn Reson Imaging. (2017) 46:1281–8. doi: 10.1002/jmri.25669

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wu X, Dong D, Zhang L, Fang M, Zhu Y, He B, et al. Exploring the predictive value of additional peritumoral regions based on deep learning and radiomics: a multicenter study. Med Phys. (2021) 48:2374–85. doi: 10.1002/mp.14767

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Fu J-H, Chuang T-C, Chung H-W, Chang H-C, Lin H-S, Hsu S-S, et al. Discriminating pyogenic brain abscesses, necrotic glioblastomas, and necrotic metastatic brain tumors by means of susceptibility-weighted imaging. Eur Radiol. (2015) 25:1413–20. doi: 10.1007/s00330-014-3518-x

PubMed Abstract | CrossRef Full Text | Google Scholar

37. De Simone M, Brogna B, Sessa G, Oliva G, Guida B, Magliulo M. Valuable contribution of magnetic resonance spectroscopy in differentiation of brain abscess from glioma. Infect Dis. (2017) 49:871–3. doi: 10.1080/23744235.2017.1331464

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Lai P-H, Chung H-W, Chang H-C, Fu J-H, Wang P-C, Hsu S-H, et al. Susceptibility-weighted imaging provides complementary value to diffusion-weighted imaging in the differentiation between pyogenic brain abscesses, necrotic glioblastomas, and necrotic metastatic brain tumors. Eur J Radiol. (2019) 117:56–61. doi: 10.1016/j.ejrad.2019.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Park JE, Kim D, Kim HS, Park SY, Kim JY, Cho SJ, et al. Quality of science and reporting of radiomics in oncologic studies: room for improvement according to radiomics quality score and TRIPOD statement. Eur Radiol. (2020) 30:523–36. doi: 10.1007/s00330-019-06360-z

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Park JE, Kim HS, Kim D, Park SY, Kim JY, Cho SJ, et al. A systematic review reporting quality of radiomics research in neuro-oncology: toward clinical utility and quality improvement using high-dimensional imaging features. BMC Cancer. (2020) 20:29. doi: 10.1186/s12885-019-6504-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: brain abscess, deep transfer learning, radiomics, convolutional neural network, cystic glioma

Citation: Bo L, Zhang Z, Jiang Z, Yang C, Huang P, Chen T, Wang Y, Yu G, Tan X, Cheng Q, Li D and Liu Z (2021) Differentiation of Brain Abscess From Cystic Glioma Using Conventional MRI Based on Deep Transfer Learning Features and Hand-Crafted Radiomics Features. Front. Med. 8:748144. doi: 10.3389/fmed.2021.748144

Received: 27 July 2021; Accepted: 06 October 2021;
Published: 12 November 2021.

Edited by:

Yong Xu, Harbin Institute of Technology, Shenzhen, China

Reviewed by:

Qi Zhu, Nanjing University of Aeronautics and Astronautics, China
Chunshan Li, Harbin Institute of Technology, China

Copyright © 2021 Bo, Zhang, Jiang, Yang, Huang, Chen, Wang, Yu, Tan, Cheng, Li and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhixiong Liu, emhpeGlvbmdsaXUmI3gwMDA0MDtjc3UuZWR1LmNu; Dengwang Li, ZGVuZ3dhbmcmI3gwMDA0MDtzZG51LmVkdS5jbg==; Quan Cheng, Y2hlbmdxdWFuJiN4MDAwNDA7Y3N1LmVkdS5jbg==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.