Research Paper Volume 11, Issue 2 pp 673—696
IDH mutation-specific radiomic signature in lower-grade gliomas
- 1 Beijing Neurosurgical Institute, Capital Medical University, Beijing, China
- 2 Neurological Imaging Center, Beijing Neurosurgical Institute, Capital Medical University, Beijing, China
- 3 Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- 4 Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing, China
- 5 Department of Nuclear Medicine, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- 6 Center of Brain Tumor, Beijing Institute for Brain Disorders, Beijing, China
- 7 China National Clinical Research Center for Neurological Diseases, Beijing, China
- 8 Chinese Glioma Genome Atlas Network (CGGA) and Asian Glioma Genome Atlas Network (AGGA)
- 9 National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, China
Received: November 12, 2018 Accepted: January 6, 2019 Published: January 29, 2019https://doi.org/10.18632/aging.101769
How to Cite
Copyright: Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Unravelling the heterogeneity is the central challenge for glioma precession oncology. In this study, we extracted quantitative image features from T2-weighted MR images and revealed that the isocitrate dehydrogenase (IDH) wild type and mutant lower grade gliomas (LGGs) differed in their expression of 146 radiomic descriptors. The logistic regression model algorithm further reduced these to 86 features. The classification model could discriminate the two types in both the training and validation sets with area under the curve values of 1.0000 and 0.9932, respectively. The transcriptome-radiomic analysis revealed that these features were associated with the immune response, biological adhesion, and several malignant behaviors, all of which are consistent with biological processes that are differentially expressed in IDH wild type and IDH mutant LGGs. Finally, a prognostic signature showed an ability to stratify IDH mutant LGGs into high and low risk groups with distinctive outcomes. By extracting a large number of radiomic features, we identified an IDH mutation-specific radiomic signature with prognostic implications. This radiomic signature may provide a way to non-invasively discriminate lower-grade gliomas as with or without the IDH mutation.
Diffuse gliomas, graded from II to IV according to the World Health Organization (WHO) criteria, are the most common and lethal primary tumors of the central nervous system. Lower grade gliomas (LGGs), designated as astrocytomas, oligodendrogliomas, and mixed oligoastrocytomas of grade II and III gliomas, account for approximately 43.2% of all gliomas diagnosed in adults [1–3]. Although LGGs have a relatively better therapeutic response and longer overall survival (OS) than fully malignant glioblastomas (GBM, WHO grade IV), they eventually transform to higher grade tumors with greater mortality [4, 5].
Isocitrate dehydrogenase (IDH) enzymes are crucial for the tricarboxylic acid cycle, catalyzing the oxidative decarboxylation of isocitrate. Mutations of IDH genes result in production of the oncometabolite 2-hydroxyglutarate (2-HG) instead of α-ketoglutarate . Previous research studies found that IDH mutations are a causative event in gliomagenesis, as well as a diagnostic, classification, and prognostic biomarker for LGG patients [7–9]. Patients harboring these mutations generally have a favorable prognosis, independent of their WHO grade [10, 11]. In light of the crucial role of IDH mutations in glioma management, IDH examination has become a routine diagnostic modality in many neuropathology centers [1, 3]. Currently, immunohistochemistry staining for IDH1R132H using formalin-fixed, paraffin-embedded specimens is the most common approach [12–14]. Alternatively, Sanger sequencing and pyrosequencing analyses are also commonly used, especially for those suspected of harboring mutations other than IDH1R132H [14, 15]. However, the classical determination of IDH status requires surgically removal of tumor tissues. A noninvasive method would be more helpful in the treatment plan and for the prognostic prediction of glioma management.
Previous studies have reported associations between imaging manifestations and IDH mutations. IDH mutant low-grade gliomas occur most frequently in the frontal lobe , especially in the area surrounding the rostral extension of the lateral ventricles . IDH wild type gliomas exhibit more post-contrast enhancement on MR images than their mutant counterparts [17, 18]. Diffusion (the apparent diffusion coefficient and fractional anisotropy) and perfusion (the relative cerebral blood volume and normalized cerebral blood volume) MR imaging can also be used in distinguishing IDH wild type and mutant gliomas [19–21]. Importantly, recent studies showed that the oncometabolite 2-HG can be detected in vivo using magnetic resonance spectroscopy (MRS), providing a better option for IDH testing [22–24]. However, the detection of 2-HG requires a unique MRS sequence device and cannot therefore be feasibly applied in a standard clinical setting . Notably, few of the above approaches are either diagnostic or quantitative.
Radiomics is a quantifying innovation that extracts large numbers of features from radiographic images using automatic data-characterization algorithms [26, 27]. In pioneering work, investigators have applied quantitative radiomics analysis to computer tomography , MR , and positron emission tomography imaging data , deciphering tumor phenotypes of non-small cell lung carcinoma , head and neck cancers , and breast cancers . Gevaert et al. utilized shape, texture, and edge sharpness to divide GBM patients into three clusters with corresponding molecular alterations . These studies highlight the potential of radiomics for quantifying and monitoring tumor-phenotypic characteristics in clinical practice . In the present study, we assessed a total of 431 radiomic features, including first order statistics, shape and size based features, textural features, and wavelet features, from T2-weighted MR images. By comparing radiological and transcriptomic profiles of IDH mutant (IDHMUT) and IDH wild type (IDHWT) LGG patients, diagnostic radiomic features for the IDH mutations were identified and independently validated. Furthermore, transcriptomic differences between the two groups and the biological processes underlying several significant radiomic features were explored. Our results suggested that the radiomic signature can separate the IDHMUT and IDHWT phenotypes of LGG patients and can potentially enable the distinction between molecular subtypes of LGGs and facilitate the design of new treatments.
Demographic and clinical characteristics
A total of 158 patients diagnosed with LGG were enrolled as the training data set. Of these, 118 (74.7%) had IDH mutant (IDHMUT) tumors and 40 (25.3%) had IDH wild type (IDHWT) tumors. No significant differences were observed with respect to age, sex, WHO grade, and tumor location between the two groups. The IDH mutation rate in the validation data set was 75.5% (77 out of 102). The clinical and pathological characteristics of the training and the validation data sets are listed in Table 1.
Table 1. Clinical characteristics of Lower Grade Glioma patients in training and validation set
|Training Set||Validation Set|
|Age median (range), years||37 (18-61)||38 (18-63)||0.403a||44 (22-70)||42 (22-68)||0.754a|
|II||25||89||0.115 b||13||62||0.005 b|
|Left||24||66||0.132 c||10||38||0.565 c|
|aStudent’s t test; bChi-square test; cFisher’s exact test; WT = Wild Type; MUT = Mutation.|
Identification of LGGs with similar radiomic patterns
To assess the radiomic expression patterns, the quantitative radiomic features were extracted from the LGG patients in the training set. An unsupervised hierarchical clustering method with average linkage revealed two clusters of patients with similar radiomic expression patterns (Figure 1). By comparing the clinical parameters of the two clusters, we found that the second cluster was significantly associated with a high frequency of the IDH mutation (P = 0.0020, Fisher’s exact test, Figure 1), which indicates a tight association between IDH mutation status and quantitative radiomic features.
Figure 1. Radiomic patterns of 431 features in LGGs. Each column corresponds to one patient in the training cohort, and each row corresponds to one z-score-normalized radiomic feature. Unsupervised clustering between radiomic features and LGG samples revealed two distinct radiomic patterns. The second cluster showed a higher frequency of the IDH mutation (**, P < 0.01).
Identification of the IDH-mutation specific radiomic signature
Based on previous observations, our goal was to identify a set of radiomic features that would enable the pre diction of the IDH mutation status in LGGs. We first screened the differences in the radiomic features between the IDHWT and IDHMUT LGGs using the SAM algorithm. A total of 146 features, including 5 first order statistics features (energy, entropy, mean, median and root mean square), 1 shape and size based feature (surface to volume ratio), 8 textural features (GLCM: contrast, dissimilarity, energy, entropy, difference entropy, informational measure of correlation 1; GLRLM: short run emphasis, run percentage) and 132 wavelet features (Figure 2A), were found to be expressed differentially. All the 146 screened features were listed in Supplementary Table 1.
Figure 2. Identification and validation of the IDH mutation-specific radiomic signature using the logistic regression. (A) A total of 146 radiomic features were selected using SAM methods. The mean value and the corresponding groups of the differentially expressed features are listed. (B and C) In the training set, the logistic regression-derived radiomic features was able to separate LGGs into two groups with high sensitivity and specificity. The AUCs were 0.86, 0.92, 0.98 and 1.00 for 10, 20, 50 and 86 radiomic features, respectively. (D) Importantly, these 86 features comprised a signature enabling the distinction of LGGs into IDHMUT and IDHWT groups with an AUC of 0.9932. (E) A Western Blot assay confirmed the expression of the mutant IDH1 protein (IDH1R132H, 1:200, DIA-H05, Dianova). (F) The radiogenomic analysis of xenograft gliomas of nude mice. Differential radiomic features between LGGs patients could be used to distinguish the IDH mutation phenotype in the xenograft model as well.
Then, we utilized the logistic regression algorithm to select the IDH mutation-specific signature. A series of ROC curves with relevant AUC was delineated using the given parameters. The results showed that only ten radiomic features were needed to divide 158 LGGs into an IDHWT group and an IDHMUT group with an AUC of 0.86, while the AUC was 0.92 for 20 features and 0.98 for 50 features (Figure 2B and 2C). Notably, using a classification model of 86 radiomic features, the enrolled patients were correctly classified into the IDHWT and the IDHMUT groups (AUC = 1.00, Figure 2B and 2C).
To further validate the classification of this radiomic signature, T2-weighted images from 102 patients and in vivo xenograft glioma models were subjected to feature extraction and logistic regression. Consistently, the results suggested that this signature, which was comprised of 86 radiomic features, could separate the 102 LGG patients into two groups with high sensitivity and specificity (AUC=0.99, Figure 2D). Fifty-four differentially expressed radiomic features were found in xenograft model, and 40 of them (such as Energy, Mean, Median and Root mean square) were shared between patients and xenograft model (Supplementary Table 2). Moreover, by using 40 of the differential features, we could divide the experimental mice into two groups in accordance with their IDH phenotype (Figure 2E and 2F). These in vivo experiments further indicated the robustness of radiomic features in differentiating IDH mutation status. Specifically, the immune associated features were also listed in Supplementary Table 2.
The potential molecular mechanism of IDH-mutation-specific radiomic features
CGGA and TCGA mRNA sequencing data were utilized to identify the biological processes and signaling pathways that differed significantly between the IDHMUT and IDHWT LGGs. Separate gene annotations of the CGGA data and the TCGA data consistently indicated that they primarily involved these biological processes: the immune response, cell adhesion, and vascular development (Supplementary Figure 1 and Supplementary Figure 2). In addition, the TCGA data suggested that the cell cycle phase and cell proliferation were also involved in the radiological manifestation of the IDHWT LGGs (Supplementary Figure 2).
Meanwhile, to attempt to determine the exact biological significance of each radiomic feature, 48 LGGs samples with both transcriptomic and radiomics data were subjected to correlation analysis and gene annotation. The results indicated that the IDH-specific radiomic features, such as the surface to volume ratio (SVR) and entropy, were primarily associated with cell polarity, cell adhesion, cell growth, and immune processes, which were the same biological processes that were found to differ between the IDHMUT and IDHWT LGGs (Supplementary Figures 3–5). Specifically, LGGs with large SVRs were found to be positively associated with a high expression of oncogenes, such as OTX1, HES1, BAG5, and the top 10 genes that were positively associated with SVR were listed in Supplementary Table 3.
Identification of a prognostic-based radiomic signature
Accumulating evidence has revealed that the IDH mutation is a crucial predictor for LGG patient outcomes. The present study also showed that the IDHMUT LGG patients had a longer OS than the IDHWT patients (Figure 3B). To further explore the prognostic ability of the IDH-specific radiomic features, we extracted a compact signature consisting of features with P < 0.05 after a univariate Cox regression analysis (Table 2). The β value of each significant radiomic feature was used for risk evaluation. When they had an elevated risk score, the patients were prone to have a higher mortality and a greater frequency of IDH mutation (Figure 3A).
Figure 3. Identification of a prognostic signature based on differential features between IDHWT and IDHMUT LGGs. (A) The expression pattern of 16 radiomic features along with the elevation of the risk score. The corresponding survival data and IDH status are listed. (B) In 158 LGGs cohort, the IDHMUT patients survived longer than the IDHWT patients (P = 0.0045, HR = 0.4024, 95%CI:0.16–0.70). (C) The risk score divided the LGGs into two groups with distinct outcomes (P = 0.0017, HR = 0.3207, 95%CI:0.19–0.68). (D) IDHMUT LGGs with a low risk score showed a favorable prognosis compared with the IDHWT patients (P = 0.0483, HR = 0.4232, 95%, CI:0.17–0.99). (E) Further, the overall survival time of the IDHMUT patients with a high risk score was not significantly different from that of the IDHWT (P = 0.199).
Table 2. Fourteen prognostic radiomic features identified by Cox regression
|Features||HR||β||95% CI||P Value|
|Mean absolute deviation_HLL||1.011||0.011||1.003||1.019||0.007|
|Surface to Volume Ratio||210.673||5.352||2.872||15475.144||0.015|
|Entropy (group 1)||0.349||-1.052||0.149||0.817||0.015|
|Energy_HLL (group 1 derived)||1.000||0.000||1.000||1.000||0.019|
|Variance_HLL(group 1 derived)||1.000||0.000||1.000||1.000||0.029|
|Root mean square_HLL||1.007||0.007||1.000||1.014||0.037|
|HR = Hazard Ratio; 95% CI = 95% Confidence Interval.|
Using the calculated risk score, we divided the LGG patients into high-risk and low-risk groups. Patients with a high-risk score had a worse prognosis (Figure 3C, HR = 0.3207, P = 0.0017). Intriguingly, IDHMUT patients could be further categorized into IDHMUT-high risk and IDHMUT-low risk groups with significantly different OSs (Figure 3D, HR = 0.4232, P = 0.0483). A multivariate Cox analysis also demonstrated that the radiomic risk score could serve as a prognostic indicator for LGGs (Table 3). Moreover, the OS was not significantly different between the IDHWT–high risk groups and the IDHMUT-high risk groups (Figure 3E, P = 0.199), further emphasizing the prognostic value of the IDH-specific radiomic signature.
Table 3. Univariate and multivariate Cox analysis in Lower Grade Glioma
|Variables||Univariate Cox Regression||Multivariate Cox Regression|
|HR||95% CI||P value||HR||95% CI||P value|
|Female vs. Male||0.722||0.363-1.439||0.355||0.870||0.428-1.767||0.699|
|II vs. III||0.337||0.177-0.643||0.001||0.420||0.213-0.827||0.012|
|MUT vs. WT||0.402||0.209-0.770||0.006||0.538||0.269-1.078||0.081|
|Low vs. High||0.806||0.727-0.894||<0.001||0.875||0.780-0.981||0.022|
A radiogenomic analysis found that a high-risk score was positively associated with genes that included ERCC1, G6PD, SOX9, and EGLN2, which are primarily enriched during the regulation of programmed cell death, cell growth, and metabolic processes. This could partially account for the radiological malignancy of the high-risk group (Figure 4). Representative samples of T2-weighted images with relevant radiomics and clinical features are presented in Figure 5. The first case was a 39-year old male patient with IDH mutant LGG. This patient was classified into the IDHMUT group with a relatively low radiomics risk score. Case 2 was a 46-year-old male with the IDH wildtype LGG who was correctly classified into the IDHWT group with a high-risk score.
Figure 4. Gene annotation of 48 patients with radiomic and transcriptome data. (A) The radiomic features, clinical characteristics, and associated genes are presented. (B) The positively associated genes (blue) that participated in GO in terms (yellow) of apoptosis, cell growth, and metabolic processes.
Figure 5. Case examples of LGG patients with T2-weighted images. Case 1 was a 39-year-old male with an IDH mutant LGG. This patient was classified into the IDHMUT group with a relatively low risk score based on the radiomic features. In contrast, case 2 was a 46-year-old male with an IDH wildtype LGG, who was correctly classified into the IDHWT group with a high risk score.
Additionally, the current study revealed that the radiomic risk score was significantly positively correlated with glioma stem cell markers such as TWIST1 (R = 0.503, P < 0.001) and CD133 (R = 0.346, P = 0.016) (Supplementary Figure 6).
By assessing the comprehensive characteristics of the entire tumor noninvasively, MR imaging is currently an indispensable approach for glioma diagnosis and treatment monitoring. The development of computational methodologies has successfully converted routine MR images to informative descriptors, substituting a quantitative and objective modality for traditionally qualitative and subjective methods. In the present study, we analyzed 431 T2-weighted radiomic features in 158 LGG patients and identified an IDH-specific radiomic signature. An integrated analysis of both radiomic and transcriptomic data indicated that these radiomic features could reflect the tumor immune response, adhesion, and several malignant biological processes, all of which are in accord with behaviors that differentiate between IDHMUT and IDHWT LGGs. Furthermore, these IDH-specific radiomic features could be utilized to establish a prognostic evaluation model. The IDHMUT patients with a low risk score showed a significantly longer OS than the IDHMUT patients with a high-risk score.
Medical imaging holds great promise for monitoring the progression of disease and the therapeutic response because it can noninvasively provide a more comprehensive view of tumors and can be performed repeatedly in routine practice . However, unlike our quantitative radiomics analyses, the conventional evaluation of MR images is subjectively based on the experience of the radiologists and neurosurgeons, leading to a lack of conformity between different clinical centers. In this study, we identified the quantitative radiomics features that were differentially expressed in two IDH phenotypes. The consistency of the findings was validated by an independent cohort of 102 LGGs patients. Intriguingly, some of the differential radiomic features could also be observed in an in vivo glioma model. Several reasons may contribute to this phenomenon. First, IDH1 mutation results in dramatically elevated levels of 2HG, a potential oncometabolite, which could influence the whole metabolic profile. Secondly, IDH1 mutation is sufficient to establish the glioma hypermethylator phenotype, which is a powerful determinant of tumor pathogenicity .
We utilized the integrative analysis of radiomic and transcriptomic data to decipher radiological characteristics that could be associated with biological processes and gene expression. The SVR, also called the surface-area-to-volume ratio, is the amount of surface area per unit volume of an object. For a given volume, the object with the smallest surface area (namely, with the smallest SVR) is a sphere. In contrast, objects with tiny spikes have a very large surface area for a given volume. In keeping with its biological significance, LGGs with large SVRs were found to be positively associated with a high expression of oncogenes, such as OTX1, HES1, BAG5; these genes are involved in biological processes that include the immune process and responses to stimuli. Similarly, a previous study revealed that patients with spherical tumors survived significantly longer than those with irregular tumor surface in glioblastomas, which indicated that tumors with irregular surface could be more malignant than spherical tumors . Another radiomic parameter enriched in the IDHWT LGGs was Entropy, which describes irregularity of pixel values within the tumor area (higher values suggesting increasing irregularity). Entropy is a promising quantitative imaging biomarker for capturing cancer imaging phenotype, and a higher Entropy was found to be associated with higher tumor metabolism, higher tumor grade, worse prognosis, and worse treatment response [37–39]. The current study revealed that Entropy was positively associated with transcription, chromatin organization and other biological functions (Supplementary Figure 3). The current finding is in accordance with previous studies.
The discovery that IDH mutations lead to accumulation of the oncometabolite 2-HG indicates an oncogenic role of the IDH mutation in the genesis of malignant brain tumors . Therefore, the biological functions of the IDH mutation have attracted attention worldwide. Researchers suggested that the IDH mutation is correlated with the regulation of HIF1 , the escape immune system , and aggregation . In the present study, RNA sequencing data from both TCGA and CGGA revealed that IDHMUT and IDHWT LGGs showed differences in immune response, vascular development, adhesion, and even proliferation, which further confirmed the genotypic and phenotypic differences between the two groups. Importantly, these genetic alterations and biological behaviors may be manifested radiologically, supporting a preoperative and noninvasive strategy for IDH prediction.
A recent study showed that IDH mutant high grade gliomas were more amenable to a complete resection of enhancing tumors and had an improved survival with the resection of non-enhanced tissues . Therefore, preoperatively evaluating the IDH status may be beneficial for surgical decision making and for developing selective targeted therapy . An increasing number of studies indicate that patients harboring IDH mutations have a better prognosis. However, some IDHMUT patients had a shorter survival than some IDHWT patients. In the present work, we used radiomic features to promote the prognostic prediction. The results implied that a combination of genetic alterations and radiomic changes could potentially provide a non-invasive methodology for genotype detection. What is more, the worse prognosis of high risk patients could be partially attributable to cell growth, metabolic processes, and programmed cell death, providing a new approach for developing therapeutic targets.
There are several limitations in the present study. We used T2-weighted MR images, a routine type of imaging for clinical glioma management, for the radiomic analysis. Although we normalized the radiological data, weighting the imaging data may not reflect the actual situation in the brain tissue, and the data can and does vary between MR scanners. A quantitative approach, such as T2 mapping, would be more suitable for the analysis of the MR images. Moreover, multi-dimensional data (copy number variance, methylation, proteomic profile, etc.) and a larger cohort would help to further delineate the radiomic and genomic landscape of glioma, and multi-regional samplings would help to obtain more precise results. Additionally, the specific radiomic features of tumors in various pathological types is an interesting topic to be investigated in future studies.
In conclusion, we demonstrated that radiomic features could serve as an alternative approach for IDH phenotype classification in LGG patients. The MR imaging is a routine examination for gliomas, and quantitative radiomic and radiogenomic analyses can potentially provide a noninvasive modality for prognosis prediction and phenotypic monitoring. As the whole process of radiomic analysis takes less than 10 minutes for an individual case and can be practice automatically with out further cost, this newly developed technique is increasingly applied to assisting clinical diagnosis and decision making.
Materials and Methods
Patients and samples
A total of 158 LGG patients who underwent surgical treatment between March 2006 and December 2012 were retrospectively selected as a training set. An additional 102 LGG patients who underwent surgical treatment from January 2016 to June 2016 were prospectively included as a validation set. The inclusion criteria were as follows: (1) pathologically confirmed lower grade glioma; (2) available IDH status; (3) available high-resolution preoperative T2-weighted MR images; (4) available clinical characteristics. The basic clinical characteristics of the patients, including age, gender, WHO grade, and tumor location are summarized in Table 1. The study was approved by the institutional review board of Beijing Tiantan Hospital. The design of the present study is illustrated in Figure 6.
Tumor masking and normalization
The tumor masking was conducted as previously described . Most of the MR images were obtained using a Trio 3.0T scanner (Siemens, Erlangen, Germany), and the remaining clinical structural images were acquired on a Magnetom Verio 3T (Siemens AG, Erlangen, Germany). The T2-weighted image parameters were as follows: repetition time = 5800 ms; echo time = 110 ms; flip angle = 150 degrees; 24 slices; field of view = 240×188 mm2; voxel size = 0.6×0.6×5.0 mm3; matrix = 384×300. Tumors were traced directly using MRIcron (http://www.mccauslandcenter.sc.edu/mricro/mricron). Masks of the brain tumors were drawn on each patient’s T2-weighted images in native space by two board-certified neuroradiologists (K.W. and X.C.), who were blind to the patients’ clinical information. Areas that produced abnormal hyperintense signals on the T2-weighted images were identified as LGG tumor areas. When a greater than 5% discrepancy existed between these two masks, the masks utilized were determined by a senior neuroradiologist (S.L.). The intensities of the voxels in each tumor mask were normalized to the z distribution ([gray value – mean gray value] / SD), on the individual level, to ensure consistency in the distribution of the gray values among the cases in the cohort.
Quantitative radiological features extraction and selection
The extraction of the radiomic features was performed as previously reported . A total of 431 image features were extracted from the tumor masks. The features were categorized into four groups. Group 1 quantitatively described the distribution of voxel intensities in the MR image using 14 descriptors. Group 2 described 8 three-dimensional features based on the shape of the tumor regions. Group 3 described textual features for quantifying intra-tumor heterogeneity, which were calculated from gray level co-occurrence (GLCM, 22 descriptors) and gray level run-length (GLRLM, 11 descriptors) texture matrices. Group 4 calculated the intensity and textural features from wavelet decompositions of the original image. All the algorithms were implemented in MATLAB (2014a).
Transcriptomic comparison between IDHWT and IDHMUT LGG
All available whole genome mRNA sequencing data of LGG patients and clinical information were acquired from the Chinese Glioma Genome Atlas database (http://www.cgga.org.cn) and the Cancer Genome Atlas database (http://cancergenome.nih.gov/). The differentially expressed genes were selected by a significance analysis of the microarray (SAM) algorithm using R programming language (http://cran.r-project.org), and with the criteria of fold change > 20% and false discovery rate (FDR) < 0.05 (Benjamini-Hochberg). Finally, the genes significantly overexpressed in the IDHWT group (CGGA, 1509; TCGA, 1262) or in the IDHMUT group (CGGA, 914; TCGA, 582) were processed using Gene Cluster and Gene Treeview software to construct a heatmap, and the online Database for Annotation, Visualization, and Integrated Discovery (DAVID, http://david.ncifcrf.gov/) program  for the gene ontology (GO).
Xenograft model of glioma and radiomic analysis
Lentiviral vectors carrying IDH1 wildtype or IDH1 R132H mutant cDNA sequences were transduced into the U87MG cells with polybrene (Sigma), as previously described . Stably transduced cells were selected during three days of puromycin (Sigma) treatment. The Western Blot was conducted to verify the expression of the mutant IDH1 enzyme using a cell lysate and the IDH1R132H antibody (1:200, DIA-H05, Dianova). Then, IDH1 wild type or R132H mutant U87MG cells (5 × 105 cells per mouse in 5 µL) were intracranially injected into 5 to 6 weeks old female nude mice (Beijing Vital River Laboratory Animal Technology), as described earlier . After 40 days, the tumors were measured using a 7T MR Image System (Bioclinscan, Bruker). The relevant T2-weighted MR images were subjected to radiomic feature extraction and analysis. After tumor segmentation, a total of 431 radiomic feature were extracted from the T2-weighted MR images from each mouse. The differentially expressed radiomic features were selected using SAM algorithm, with the criteria of false discovery rate (FDR) < 0.05 (Benjamini-Hochberg).
IDH phenotype classification
First, similar to the radiomic analysis in xenograft model, SAM algorithm was conducted on the 431 radiomic features to select the differentially expressed radiomic features between IDHWT and IDHMUT tumors, with the criteria of false discovery rate (FDR) < 0.05 (Benjamini-Hochberg). Next, the selected differentially expressed radiomic features were utilized for IDH phenotype classification. The classification process was conducted using a logistic regression model (Y = expr feature 1 × β feature 1 + expr feature 2 × β feature 2 + … + expr feature n × β feature n + ε) based on the MATLAB (2014a) software. In this model, ‘Y’ represents the IDH status (1 indicates mutation, 0 indicates wildtype), while ‘expr’ represents the expression value of each radiomic feature. β represents the model parameter to be estimated, and ε is the estimated residual. The prediction results were further interpreted using the receiver operating characteristic (ROC) curve. To select the most predictive parameters, the logistic regression algorithm was applied repeatedly. Using the dimensionality reduction principle, the radiomic feature with the highest P value in predicting IDH status was excluded from the model each time until features with the best predictive effect were identified. The signature derived from the training set was subsequently applied to the validation set.
Identification of the prognosis-based signature
The prognostic values for each IDH-specific feature in patients with LGG were calculated using a univariate Cox regression after elimination of the patients with an OS of less than 30 days. The significant radiomic features with P values < 0.05 were selected to develop a prognosis-based signature. The risk score was calculated using the linear combination of the selected features weighted by the regression coefficient derived from the univariate Cox regression analysis (β), referring to previous studies [47, 48]. The risk score for the OS of each individual was calculated as follows: Risk score = expr feature 1 × β feature 1 + expr feature 2 × β feature 2 + … + expr feature n × β feature n. ‘Expr’ represents the expression value of each radiomic feature. The ‘expr’ and β used in methods 2.7 have no relationship with those used in methods 2.6.
We next divided the patients into high-risk and low-risk groups based on the risk score. The cutoff value was determined when the P value was the smallest in log-rank test in the training set, and then the cut-off value was fixed and applied to all Kaplan-Meier curve analyses.
Forty-eight LGG patients involved in the radiomic analysis were subjected to an Agilent Whole Human Genome Array analysis according to the manufacturer’s instructions . The data were acquired using the Agilent G2565BA Microarray Scanner System and Agilent Feature Extraction Software (version 9.1). The probe intensities were normalized using GeneSpring GX 11.0. The IDH mutations of the training data set were assessed by pyrosequencing, while the IDH status for the validation data set was identified by immunohistochemistry (IDH1R132H, DIA-H05, Dianova).
Radiogenomic analysis was further performed. The Pearson correlation coefficients were calculated between the genes and IDH-associated radiomic features. The association was identified to be statistically significant when the absolute value of Pearson correlation coefficient was > 0.4 and the P value was < 0.05. Gene ontology (GO) analysis was conducted to investigate the underlying biological processes of the radiomic features based on the DAVID Bioinformatics Resources (http://david.ncifcrf.gov/). The top 200 positive/negative genes that were significantly associated with each feature were subjected to GO analysis to reveal the underlying biological processes of each feature. Using this method, the underlying biological processes of the risk score and the differentially expressed features in xenograft model were also investigated. In addition, the relationship between the radiomic risk score and glioma stem cell signatures was also assessed with using the Pearson correlation analysis.
The significant differences between the two groups were estimated using a Student’s t-test. Chi-square and Fisher’s exact tests were used to compare the frequencies between the groups. The OS curves were plotted according to the Kaplan–Meier method, with the log-rank test applied for comparison. A Cox regression was used to determine the prognostic value of each radiomic feature for the OS in LGG patients. All the differences were considered statistically significant at the two-sided P < 0.05 level.
Conflicts of Interest
The authors of this manuscript have no conflict of interests to declare.
This work was funded by Beijing Natural Science Foundation (7174295), National Natural Science Foundation of China (81601452), Beijing Postdoctoral Research Foundation (2016ZZ-37), National Basic Research Program of China (No. 2015CB755500), and National Key Research and Development Plan (2016YFC0902500).
- 1. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, Ohgaki H, Wiestler OD, Kleihues P, Ellison DW. The 2016 World Health Organization Classification of Tumors of the Central Nervous System: a summary.Acta neuropathologica.2016;131:803-820. https://doi.org/10.1007/s00401-016-1545-1 [PubMed]
- 2. Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, Scheithauer BW, Kleihues P. The 2007 WHO classification of tumours of the central nervous system.Acta neuropathologica.2007;114:97-109. https://doi.org/10.1007/s00401-007-0243-4 [PubMed]
- 3. Jiang T, Mao Y, Ma W, Mao Q, You Y, Yang X, Jiang C, Kang C, Li X, Chen L, Qiu X, Wang W, Li W, et al. CGCG clinical practice guidelines for the management of adult diffuse gliomas.Cancer letters.2016;375:263-273. https://doi.org/10.1016/j.canlet.2016.01.024 [PubMed]
- 4. Kannan K, Inagaki A, Silber J, Gorovets D, Zhang J, Kastenhuber ER, Heguy A, Petrini JH, Chan TA, Huse JT. Whole-exome sequencing identifies ATRX mutation as a key molecular determinant in lower-grade glioma.Oncotarget.2012;3:1194-1203. https://doi.org/10.18632/oncotarget.689 [PubMed]
- 5. Guan X, Vengoechea J, Zheng S, Sloan AE, Chen Y, Brat DJ, O'Neill BP, de Groot J, Yust-Katz S, Yung WK, Cohen ML, Aldape KD, Rosenfeld S, et al. Molecular subtypes of glioblastoma are relevant to lower grade glioma.PLoS One.2014;9:e91216 https://doi.org/10.1371/journal.pone.0091216 [PubMed]
- 6. Waitkus MS, Diplas BH, Yan H. Isocitrate dehydrogenase mutations in gliomas.Neuro-oncology.2016;18:16-26. https://doi.org/10.1093/neuonc/nov136 [PubMed]
- 7. Yan H, Bigner DD, Velculescu V, Parsons DW. Mutant metabolic enzymes are at the origin of gliomas.Cancer research.2009;69:9157-9159. https://doi.org/10.1158/0008-5472.CAN-09-2650 [PubMed]
- 8. Jiao Y, Killela PJ, Reitman ZJ, Rasheed AB, Heaphy CM, de Wilde RF, Rodriguez FJ, Rosemberg S, Oba-Shinjo SM, Nagahashi Marie SK, Bettegowda C, Agrawal N, Lipp E, et al. Frequent ATRX, CIC, FUBP1 and IDH1 mutations refine the classification of malignant gliomas.Oncotarget.2012;3:709-722. https://doi.org/10.18632/oncotarget.588 [PubMed]
- 9. Sanson M, Marie Y, Paris S, Idbaih A, Laffaire J, Ducray F, El Hallani S, Boisselier B, Mokhtari K, Hoang-Xuan K, Delattre JY. Isocitrate dehydrogenase 1 codon 132 mutation is an important prognostic biomarker in gliomas.J Clin Oncol.2009;27:4150-4154. https://doi.org/10.1200/JCO.2009.21.9832 [PubMed]
- 10. Yan H, Parsons DW, Jin G, McLendon R, Rasheed BA, Yuan W, Kos I, Batinic-Haberle I, Jones S, Riggins GJ, Friedman H, Friedman A, Reardon D, et al. IDH1 and IDH2 mutations in gliomas.N Engl J Med.2009;360:765-773. https://doi.org/10.1056/NEJMoa0808710 [PubMed]
- 11. Zhang CB, Bao ZS, Wang HJ, Yan W, Liu YW, Li MY, Zhang W, Chen L, Jiang T. Correlation of IDH1/2 mutation with clinicopathologic factors and prognosis in anaplastic gliomas: a report of 203 patients from China.J Cancer Res Clin Oncol.2014;140:45-51. https://doi.org/10.1007/s00432-013-1519-9 [PubMed]
- 12. Capper D, Weissert S, Balss J, Habel A, Meyer J, Jager D, Ackermann U, Tessmer C, Korshunov A, Zentgraf H, Hartmann C, von Deimling A. Characterization of R132H mutation-specific IDH1 antibody binding in brain tumors.Brain Pathol.2010;20:245-254. https://doi.org/10.1111/j.1750-3639.2009.00352.x [PubMed]
- 13. Horbinski C, Kofler J, Kelly LM, Murdoch GH, Nikiforova MN. Diagnostic use of IDH1/2 mutation analysis in routine clinical testing of formalin-fixed, paraffin-embedded glioma tissues.J Neuropathol Exp Neurol.2009;68:1319-1325. https://doi.org/10.1097/NEN.0b013e3181c391be [PubMed]
- 14. Reuss DE, Sahm F, Schrimpf D, Wiestler B, Capper D, Koelsche C, Schweizer L, Korshunov A, Jones DT, Hovestadt V, Mittelbronn M, Schittenhelm J, Herold-Mende C, et al. ATRX and IDH1–R132H immunohistochemistry with subsequent copy number analysis and IDH sequencing as a basis for an "integrated" diagnostic approach for adult astrocytoma, oligodendroglioma and glioblastoma.Acta neuropathologica.2015;129:133-146. https://doi.org/10.1007/s00401-014-1370-3 [PubMed]
- 15. Setty P, Hammes J, Rothamel T, Vladimirova V, Kramm CM, Pietsch T, Waha A. A pyrosequencing-based assay for the rapid detection of IDH1 mutations in clinical samples.J Mol Diagn.2010;12:750-756. https://doi.org/10.2353/jmoldx.2010.090237 [PubMed]
- 16. Wang Y, Zhang T, Li S, Fan X, Ma J, Wang L, Jiang T. Anatomical localization of isocitrate dehydrogenase 1 mutation: a voxel-based radiographic study of 146 low-grade gliomas.Eur J Neurol.2015;22:348-354. https://doi.org/10.1111/ene.12578 [PubMed]
- 17. Wang K, Wang Y, Fan X, Wang J, Li G, Ma J, Ma J, Jiang T, Dai J. Radiological features combined with IDH1 status for predicting the survival outcome of glioblastoma patients.Neuro-oncology.2016;18:589-597. https://doi.org/10.1093/neuonc/nov239 [PubMed]
- 18. Wang YY, Wang K, Li SW, Wang JF, Ma J, Jiang T, Dai JP. Patterns of Tumor Contrast Enhancement Predict the Prognosis of Anaplastic Gliomas with IDH1 Mutation.AJNR Am J Neuroradiol.2015;36:2023-2029. https://doi.org/10.3174/ajnr.A4407 [PubMed]
- 19. Kickingereder P, Sahm F, Radbruch A, Wick W, Heiland S, Deimling A, Bendszus M, Wiestler B. IDH mutation status is associated with a distinct hypoxia/angiogenesis transcriptome signature which is non-invasively predictable with rCBV imaging in human glioma.Scientific reports.2015;5:16238 https://doi.org/10.1038/srep16238 [PubMed]
- 20. Lee S, Choi SH, Ryoo I, Yoon TJ, Kim TM, Lee SH, Park CK, Kim JH, Sohn CH, Park SH, Kim IH. Evaluation of the microenvironmental heterogeneity in high-grade gliomas with IDH1/2 gene mutation using histogram analysis of diffusion-weighted imaging and dynamic-susceptibility contrast perfusion imaging.J Neurooncol.2015;121:141-150. https://doi.org/10.1007/s11060-014-1614-z [PubMed]
- 21. Tan WL, Huang WY, Yin B, Xiong J, Wu JS, Geng DY. Can diffusion tensor imaging noninvasively detect IDH1 gene mutations in astrogliomas? A retrospective study of 112 cases.AJNR Am J Neuroradiol.2014;35:920-927. https://doi.org/10.3174/ajnr.A3803 [PubMed]
- 22. Andronesi OC, Kim GS, Gerstner E, Batchelor T, Tzika AA, Fantin VR, Vander Heiden MG, Sorensen AG. Detection of 2-hydroxyglutarate in IDH-mutated glioma patients by in vivo spectral-editing and 2D correlation magnetic resonance spectroscopy.Sci Transl Med.2012;4:116ra114 https://doi.org/10.1126/scitranslmed.3002693 [PubMed]
- 23. Choi C, Ganji SK, DeBerardinis RJ, Hatanpaa KJ, Rakheja D, Kovacs Z, Yang XL, Mashimo T, Raisanen JM, Marin-Valencia I, Pascual JM, Madden CJ, Mickey BE, et al. 2-hydroxyglutarate detection by magnetic resonance spectroscopy in IDH-mutated patients with gliomas.Nat Med.2012;18:624-629. https://doi.org/10.1038/nm.2682 [PubMed]
- 24. Emir UE, Larkin SJ, de Pennington N, Voets N, Plaha P, Stacey R, Al-Qahtani K, McCullagh J, Schofield CJ, Clare S, Jezzard P, Cadoux-Hudson T, Ansorge O. Noninvasive Quantification of 2-Hydroxyglutarate in Human Gliomas with IDH1 and IDH2 Mutations.Cancer research.2016;76:43-49. https://doi.org/10.1158/0008-5472.CAN-15-0934 [PubMed]
- 25. Price SJ. Imaging Markers of Isocitrate Dehydrogenase-1 Mutations in Gliomas.Neurosurgery.2015;62:166-170. https://doi.org/10.1227/NEU.0000000000000812 [PubMed]
- 26. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, Bussink J, Monshouwer R, Haibe-Kains B, Rietveld D, Hoebers F, Rietbergen MM, Leemans CR, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach.Nature communications.2014;5:4006 https://doi.org/10.1038/ncomms5006 [PubMed]
- 27. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, Zegers CM, Gillies R, Boellard R, Dekker A, Aerts HJ. Radiomics: extracting more information from medical images using advanced feature analysis.European journal of cancer.2012;48:441-446. https://doi.org/10.1016/j.ejca.2011.11.036 [PubMed]
- 28. Coroller TP, Grossmann P, Hou Y, Rios Velazquez E, Leijenaar RT, Hermann G, Lambin P, Haibe-Kains B, Mak RH, Aerts HJ. CT-based radiomic signature predicts distant metastasis in lung adenocarcinoma.Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.2015;114:345-350. https://doi.org/10.1016/j.radonc.2015.02.015 [PubMed]
- 29. Itakura H, Achrol AS, Mitchell LA, Loya JJ, Liu T, Westbroek EM, Feroze AH, Rodriguez S, Echegaray S, Azad TD, Yeom KW, Napel S, Rubin DL, et al. Magnetic resonance image features identify glioblastoma phenotypic subtypes with distinct molecular pathway activities.Sci Transl Med.2015;7:303ra138 https://doi.org/10.1126/scitranslmed.aaa7582 [PubMed]
- 30. Leijenaar RT, Carvalho S, Velazquez ER, van Elmpt WJ, Parmar C, Hoekstra OS, Hoekstra CJ, Boellaard R, Dekker AL, Gillies RJ, Aerts HJ, Lambin P. Stability of FDG-PET Radiomics features: an integrated analysis of test-retest and inter-observer variability.Acta Oncol.2013;52:1391-1397. https://doi.org/10.3109/0284186X.2013.812798 [PubMed]
- 31. Parmar C, Leijenaar RT, Grossmann P, Rios Velazquez E, Bussink J, Rietveld D, Rietbergen MM, Haibe-Kains B, Lambin P, Aerts HJ. Radiomic feature clusters and prognostic signatures specific for Lung and Head & Neck cancer.Scientific reports.2015;5:11044 https://doi.org/10.1038/srep11044 [PubMed]
- 32. Grimm LJ, Zhang J, Mazurowski MA. Computational approach to radiogenomics of breast cancer: Luminal A and luminal B molecular subtypes are associated with imaging features on routine breast MRI extracted using computer vision algorithms.Journal of magnetic resonance imaging : JMRI.2015;42:902-907. https://doi.org/10.1002/jmri.24879 [PubMed]
- 33. Parmar C, Grossmann P, Bussink J, Lambin P, Aerts HJ. Machine Learning methods for Quantitative Radiomic Biomarkers.Scientific reports.2015;5:13087 https://doi.org/10.1038/srep13087 [PubMed]
- 34. Dang L, White DW, Gross S, Bennett BD, Bittinger MA, Driggers EM, Fantin VR, Jang HG, Jin S, Keenan MC, Marks KM, Prins RM, Ward PS, et al. Cancer-associated IDH1 mutations produce 2-hydroxyglutarate.Nature.2009;462:739-744. https://doi.org/10.1038/nature08617 [PubMed]
- 35. Turcan S, Rohle D, Goenka A, Walsh LA, Fang F, Yilmaz E, Campos C, Fabius AW, Lu C, Ward PS, Thompson CB, Kaufman A, Guryanova O, et al. IDH1 mutation is sufficient to establish the glioma hypermethylator phenotype.Nature.2012;483:479-483. https://doi.org/10.1038/nature10866 [PubMed]
- 36. Perez-Beteta J, Molina-Garcia D, Ortiz-Alhambra JA, Fernandez-Romero A, Luque B, Arregui E, Calvo M, Borras JM, Melendez B, Rodriguez de Lope A, Moreno de la Presa R, Iglesias Bayo L, Barcia JA, et al. Tumor Surface Regularity at MR Imaging Predicts Survival and Response to Surgery in Patients with Glioblastoma.Radiology.2018;288:218-225. https://doi.org/10.1148/radiol.2018171051 [PubMed]
- 37. Dercle L, Ammari S, Bateson M, Durand PB, Haspinger E, Massard C, Jaudet C, Varga A, Deutsch E, Soria JC, Ferte C. Limits of radiomic-based entropy as a surrogate of tumor heterogeneity: ROI-area, acquisition protocol and tissue site exert substantial influence.Sci Rep.2017;7:7952 https://doi.org/10.1038/s41598-017-08310-5 [PubMed]
- 38. Ganeshan B, Abaleke S, Young RC, Chatwin CR, Miles KA. Texture analysis of non-small cell lung cancer on unenhanced computed tomography: initial evidence for a relationship with tumour glucose metabolism and stage.Cancer Imaging.2010;10:137-143. https://doi.org/10.1102/1470-7330.2010.0021 [PubMed]
- 39. Ganeshan B, Skogen K, Pressney I, Coutroubis D, Miles K. Tumour heterogeneity in oesophageal cancer assessed by CT texture analysis: preliminary evidence of an association with tumour metabolism, stage, and survival.Clin Radiol.2012;67:157-164. https://doi.org/10.1016/j.crad.2011.08.012 [PubMed]
- 40. Zhao S, Lin Y, Xu W, Jiang W, Zha Z, Wang P, Yu W, Li Z, Gong L, Peng Y, Ding J, Lei Q, Guan KL, et al. Glioma-derived mutations in IDH1 dominantly inhibit IDH1 catalytic activity and induce HIF-1alpha.Science.2009;324:261-265. https://doi.org/10.1126/science.1170944 [PubMed]
- 41. Zhang X, Rao A, Sette P, Deibert C, Pomerantz A, Kim WJ, Kohanbash G, Chang Y, Park Y, Engh J, Choi J, Chan T, Okada H, et al. IDH mutant gliomas escape natural killer cell immune surveillance by downregulation of NKG2D ligand expression.Neuro-oncology.2016;18:1402-1412. https://doi.org/10.1093/neuonc/now061 [PubMed]
- 42. Hu H, Wang Z, Liu Y, Zhang C, Li M, Zhang W, Wang K, Cai J, Cheng W, Huang H, Jiang T. Genome-wide transcriptional analyses of Chinese patients reveal cell migration is attenuated in IDH1-mutant glioblastomas.Cancer letters.2015;357:566-574. https://doi.org/10.1016/j.canlet.2014.12.018 [PubMed]
- 43. Beiko J, Suki D, Hess KR, Fox BD, Cheung V, Cabral M, Shonka N, Gilbert MR, Sawaya R, Prabhu SS, Weinberg J, Lang FF, Aldape KD, et al. IDH1 mutant malignant astrocytomas are more amenable to surgical resection and have a survival benefit associated with maximal surgical resection.Neuro-oncology.2014;16:81-91. https://doi.org/10.1093/neuonc/not159 [PubMed]
- 44. Wang Y, Qian T, You G, Peng X, Chen C, You Y, Yao K, Wu C, Ma J, Sha Z, Wang S, Jiang T. Localizing seizure-susceptible brain regions associated with low-grade gliomas using voxel-based lesion-symptom mapping.Neuro-oncology.2015;17:282-288. https://doi.org/10.1093/neuonc/nou130 [PubMed]
- 45. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat Protoc.2009;4:44-57. https://doi.org/10.1038/nprot.2008.211 [PubMed]
- 46. Liu X, Wang X, Du W, Chen L, Wang G, Cui Y, Liu Y, Dou Z, Wang H, Zhang P, Chang L, Yi L, Cai J, et al. Suppressor of fused (Sufu) represses Gli1 transcription and nuclear accumulation, inhibits glioma cell proliferation, invasion and vasculogenic mimicry, improving glioma chemo-sensitivity and prognosis.Oncotarget.2014;5:11681-11694. https://doi.org/10.18632/oncotarget.2585 [PubMed]
- 47. Zhang CB, Zhu P, Yang P, Cai JQ, Wang ZL, Li QB, Bao ZS, Zhang W, Jiang T. Identification of high risk anaplastic gliomas by a diagnostic and prognostic signature derived from mRNA expression profiling.Oncotarget.2015;6:36643-36651. https://doi.org/10.18632/oncotarget.5421 [PubMed]
- 48. Long J, Zhang L, Wan X, Lin J, Bai Y, Xu W, Xiong J, Zhao H. A four-gene-based prognostic model predicts overall survival in patients with hepatocellular carcinoma.J Cell Mol Med.2018; https://doi.org/10.1111/jcmm.13863 [PubMed]
- 49. Yan W, Zhang W, You G, Zhang J, Han L, Bao Z, Wang Y, Liu Y, Jiang C, Kang C, You Y, Jiang T. Molecular classification of gliomas based on whole genome gene expression: a systematic report of 225 samples from the Chinese Glioma Cooperative Group.Neuro-oncology.2012;14:1432-1440. https://doi.org/10.1093/neuonc/nos263 [PubMed]