Abstract

Objectives. To investigate the predictors of telomerase reverse transcriptase (TERT) promoter mutations in adults suffered from high-grade glioma (HGG) through radiomics analysis, develop a noninvasive approach to evaluate TERT promoter mutations. Methods. 126 adult patients with HGG (88 in the training cohort and 38 in the validation cohort) were retrospectively enrolled. Totally 5064 radiomics features were, respectively, extracted from three VOIs (necrosis, enhanced, and edema) in MRI. Firstly, an optimal radiomics signature (Radscore) was established based on LASSO regression. Secondly, univariate and multivariate logistic regression analyses were performed to investigate important potential variables as predictors of TERT promoter mutations. Besides, multiparameter models were established and evaluated. Eventually, an optimal model was visualized as radiomics nomogram for clinical evaluations. Results. 6 radiomics features were selected to build Radscore signature through LASSO regression. Among them, 5 were from necrotic VOIs and 1 was from enhanced ones. With univariate and multivariate analysis, necrotic volume percentages of core (CNV), Age, Cho/Cr, Lac, and Radscore were significantly higher in TERTm than in TERTw (). 4 models were built in our study. Compared with Model B (Age, Cho/Cr, Lac, and Radscore), Model A (Age, Cho/Cr, Lac, Radscore, and CNV) has a larger AUC in both training (0.955 vs. 0.917, ) and validation (0.889 vs. 0.868, ) cohorts. It also has higher performances in net reclassification improvement (NRI), integrated discrimination improvement (IDI), and decision curve analysis (DCA) evaluation. Conclusively, Model A was visualized as a radiomics nomogram. Calibration curve shows a good agreement between estimated and actual probabilities. Conclusions. Age, Cho/Cr, Lac, CNV, and Radscore are important indicators for TERT promoter mutation predictions in HGG. Tumor necrosis seems to be closely related to TERT promoter mutations. Radiomics nomogram based on multiparameter MRI and CNV has higher prediction accuracies.

1. Introduction

As the most common primary brain tumor in adults, glioma can be classified as low-grade (LGG, WHO I-II) and high-grade (HGG, WHO III-IV) according to the classification criteria of the World Health Organization (WHO) [1]. Compared with LGG, HGG is characterized by more vigorous cell growth, more tumor angiogenesis, higher heterogeneity, and worse prognosis, especially glioblastoma (GBM), with a median survival of 12-14 months for the patients receiving standard treatments [2, 3]. As indicated in previous studies, hyperproliferativeness provides a significant biological basis for the occurrence, migration, diffusion, invasion, postoperative recurrence, and drug resistance of HGG, where the stability and nonshrinkage of telomeres play a crucial role [4].

Telomeres, controlling the limited division of normal cells, are shortened with each division of normal cells, whereas they could be continuously elongated by telomerase in cancer cells [5]. Telomerase consisted of the RNA subunit and reverse transcriptase (TERT), maintaining the length of telomere by adding hexamer repeats at the end of chromosomes [6]. Therefore, TERT plays a crucial role in the cancerization. Recently, the recurrent mutations at two hotspots termed as C228T and C250T in the TERT promoter have been identified in gliomas [7]. Besides, the mutations have been considered as one of the major mechanisms of telomerase activation in gliomas [8]. Up to now, a large number of researches have explored the role and value of TERT in glioma, which demonstrated that 80% of the TERT promoters in the primary GBM have mutated [9]. The TERT genotype is not only a vital prognostic and predictive biomarker for glioma, especially for HGG [10, 11], but also a promising indicator for the sensitivity of GBM to radiotherapy and temozolomide [12]. Currently, the TERT genotype has been prevailingly determined by sequencing tumor samples, which can only be obtained postoperatively. Accurate assessment of the TERT genotype before surgery can direct the exploitation of therapeutic strategies. Therefore, the establishment of noninvasive technology to identify the TERT genotype of tumors is urgently needed.

Radiomics has attracted increased attention in recent years due to the representation of medical images containing information on the pathophysiology and prognosis of diseases [13]. A variety of quantitative radiological features could be extracted from medical images to reveal the information on tumors [14]. Radiomics feature has been applied as the noninvasive alternative to identify the genomic and proteomic changes in tumors, which also broadly utilized in tumor diagnosis, prognosis prediction, treatment selection, gene prediction, and so on [1518]. However, the use of radiomics analysis to predict the mutant status of the TERT promoter has not been widely reported. This study aimed to investigate the predictors of TERT promoter mutations in HGG through radiomics analysis and develop a noninvasive approach to evaluate of TERT promoter mutations.

2. Materials and Methods

The institutional review board has approved this retrospective study, and the requirements for patient informed consent were waived for the anonymity of data.

2.1. Patients

Based on the inclusion and exclusion criteria as shown in Figure 1, 126 patients were finally enrolled in our study. All the subjects were randomly divided into the training and the validation cohorts by computer sampling at a ratio of 7 : 3. The TERT promoter mutations were determined by capillary electrophoresis. Clinical and pathological data were obtained via reviewing electronic medical records.

2.2. Data Acquisition of MRI

All preoperative MRI was performed on 3.0 T MR scanners (Magneto Trio, Siemens, Germany) with an eight-channel head coil. Contrast-enhanced T1-weighted (CE-T1w), T2-weighted imaging fluid-attenuated inversion recovery (T2flair), T1-weighted (T1w), T2-weighted (T2w), and magnetic resonance spectroscopy (MRS) sequences were applied for the following analyses. The acquisition parameters were summarized in Supplementary S1.

2.3. Preprocessing, Segmentation, and Feature Extraction of Images

The schemes of image preprocessing consisted of coregistration, reslice, and normalization. Firstly, the T2flair, T2w, and T1w images were coregistered to the corresponding CE-T1w images on the basis of affine transformation through the Linear Image Registration Tool (FLIRT) of Functional MRI of the Brain (FMRIB) Software Library (FSL) of the Oxford Center. Subsequently, the resolution of each modality was uniformly resampled to . Finally, the intensity of the images of each modality was normalized according to the Collewet normalization algorithm () to correct the effects of different acquisition protocols [19]. The success of coregistration and normalization was visually verified by two authors (T.H.A. and W.H., with 10 and 3 years of experience in brain MRI research, respectively).

Tumor segmentation was conducted on 3D slicer (version 4.10.1). As illustrated in Figures 2 and 3 heterogeneous regions (the enhanced lesion inclusive of necrosis (tumor core), enhanced, and necrosis) were plotted by the semiautomatic segment editor module. Detailed procedures and parameters for tumor separation were listed in Supplementary S2 to facilitate reproducibility. The volume of interest (VOI) was taken charge by T.H.A., who was trained and supervised by a board-certified radiologist specialized in neurooncology (W.G.Y. with 25-year experience). All the VOIs were eventually registered to the MNI152 standard space for obtaining the location information by FSL-FLIRT.

The radiomics features were extracted by PyRadiomics (version 2.1.2), a flexible open-source platform capable of obtaining a large panel of engineered features from medical images [20]. The extracted features included the “first-order statistics (First-order), Gray Level Cooccurrence Matrix (GLCM), Gray Level Dependence Matrix (GLDM), Gray Level Run Length Matrix (GLRLM), Gray Level Size Zone Matrix (GLSZM), Neighbouring Gray Tone Difference Matrix (NGTDM), Shap3D and Shap2D.” All the above features were displayed in the Supplementary S3 to facilitate the application of our findings. Only the features of T2flair and CE-T1w sequences were acquired, which were considered as the best series of HGG studies in various former studies [21]. Therefore, we obtained the features of necrosis, enhanced, and edema VOIs on the above two sequences, respectively.

2.4. Intraobserver and Interobserver Agreement

Interobserver and intraobserver agreement of VOI-based radiomics feature in 30 randomly chosen patients were evaluated in our study. Interobserver agreement of feature extraction by two authors (T.H.A. and W.H.) was initially analyzed. Meanwhile, to assess the intraobserver agreement, one of the researchers (T.H.A.) repeated the extraction twice in two weeks following the equivalent protocol. These radiomics features extracted from the VOIs were evaluated by the intraclass correlation coefficient (ICC). The score of ICC greater than 0.85 was considered a satisfactory agreement.

2.5. Dimensionality Reduction and Radiomics Features Selection

To reduce the dimension and decrease the redundant information, two steps were scheduled to select the features. Firstly, all features were analyzed by the independent samples test or Mann-Whiney test in the training cohort, and variables with were chosen as the potentially important parameters. Additionally, the least absolute shrinkage and selection operator (LASSO) feature selection algorithm was subsequently conducted for dimensionality reduction and feature selection. Nonzero coefficients chosen by LASSO as the optimal features were utilized to establish the Radscore formula, applying to calculate the Radscore for each patient to predict the TERT promoter mutations.

2.6. The Selection of Clinical and Radiological Characteristics

The clinical characteristics of age, gender, and tumor grade were identified. The radiological features were assessed by T.H.A., containing tumor location, the volume of the necrotic or cystic part of the tumor (NeV), the volume of the solid portion of the tumor (EnV), the volume of edema around tumor (EdV), the percentage of necrotic volume in core volume of the tumor (CNV, %) and that in overall (ONV, %), and MRS features based on enhanced region. CNV was represented as NeV/(NeV+EnV), and ONV was represented as NeV/(NeV+EnV+EdV), whose definition can be more intuitively shown in Figures 2(d) and 2(e). Univariate and multivariate logistic regression analyses were performed to determine potential important variables. The features with were selected as the elements for the establishment of the model.

2.7. Development of Radiomics Model and Evaluation of the Performance

Models were established based on the logistic regression with forwarding stepwise selection. To determine the best model, the performances of predictive models were assessed in validation cohorts. The discrimination was evaluated by the receiver operating characteristic (ROC) curve, net reclassification improvement (NRI), and integrated discrimination improvement (IDI), and the clinical availability was appraised by decision curve analysis (DCA) [2224]. To offer an individualized and easy-to-operate tool for noninvasive prediction of the TERT genotype, the optimal model was visualized as radiomics nomogram.

2.8. Statistical Analysis

All statistical analyses were performed according to the (version 3.4.3). Data were presented as for continuous variables and as frequency (%) for categorical variables. Detailed statistical steps and packages were listed in Supplementary S4.

3. Results

3.1. Clinical Characteristics of the Patients

126 patients who suffered from HGG were divided into the training cohort (88, 70%) and the validation cohort (38, 30%). There was no significant difference in terms of age (), gender (), grade (), or location () between the training and validation cohorts (Table 1).

3.2. Feature Extraction and Dimensionality Reduction

Eventually, 1688 features were extracted from each VOI, and a total of 5064 radiomics features were generated from every patient. Totally, 1230 features () were screened out as the potentially crucial variables at the first step of dimensionality reduction. Applying the optimal regulation weight λ () for the LASSO algorithm, 6 nonzero coefficient features were finally chosen (Figure 3). Among them, 5 were from necrotic VOIs and 1 was from enhanced ones. Three, two, and one radiomics features were selected from CE-T1w, T2flair, and original shape class, respectively (Table 2). In the end, the Radscore formula was established according to the coefficients of the six features obtained previously (Supplementary S5).

3.3. Interobserver and Intraobserver Agreement

For the intraobserver agreement, 209 features were unqualified (ICC <0.85), accounting for 4.13% of all variables (209/5064) (Figure 4(a)). For the interobserver agreement, 299 features were unqualified, occupying 5.91% of the variables (299/5064) (Figure 4(b)). Furthermore, none of the unqualified features consisted of the six features selected by LASSO regression (ICC were all higher than 0.85, Table 2).

3.4. Radiomics Model Construction and Validation

As indicated in Table 3, Radscore, Age, CNV, Cho/Cr, and Lac were significantly higher in TERTm than in TERTw () in the training cohort by univariate and multivariate analysis.

4 models were built in this study, Age, Cho/Cr, Lac, and Radscore were taken into consideration in Model B, while Model A further added the feature of CNV based on Model B. Model C was composed of a single Radscore feature, while Model D consisted of a single CNV feature. (In fact, we had built a total of 15 models using the above 5 features (Supplementary S6)). As described in Table 4, Model A and Model B performed better, and in order to compare the differences between the two, we conducted a further evaluation. The difference of AUCs in the two models was statistically significant analyzed by the DeLong test () (Figures 5(a) and 5(b)). Additionally, NRI was 0.187 (95% CI, 0.022-0.351) in the training cohort, 0.261 (95% CI, 0.045-0.478) in the validation cohort. Besides, the IDI was 0.146 (95% CI, 0.071-0.221) in training cohort, 0.128 (95% CI, 0.033-0.222) in the validation cohort. The DCA curves illustrated that Model A offered a higher overall net benefit in contrast to Model B, indicating Model A was superior across nearly the entire range of pt values (Figure 5(c)).

3.5. The Visualization of Radiomics Nomogram

From the above, Model A was visualized as the radiomics nomogram (Figure 6(a)). The AUC of the nomogram was 0.951 (95% CI, 0.906-0.982) in the training cohort, 0.883 (95% CI, 0.757-0.953) in the validation cohort (Figure 6(b)). Moreover, the calibration curves of the nomograms in training and validation cohorts both indicated good agreement between the predictability of TERT promoter mutations and actual status, respectively (Figure 6(c)).

4. Discussions

In this study, we found that the radiomics method can well predict TERT promoter mutations. Radiogenomics was a new field for studying the relationship between radiological features and genomic data, which explored the relationship between radiological features and gene phenotypes by extracting quantitative information on a large number of radiological data features [13]. Studies predicting glioma genes based on radiomics mainly focused on the IDH, 1p/19q, MGMT, ATRX, and EGFR genes, which all agreed that radiomics was an excellent noninvasive method for predicting the genetic status of glioma [14, 25, 26]. However, the relevant researches using radiomics to predict TERT genotypes were relatively rare. Arita et al. [27] used radiomics methods to predict IDH and TERT genes in grade II/III gliomas, suggesting that conventional MRI-based radiomics could be a noninvasive diagnostic technique for molecular characterization of grade II/III gliomas. As Gillies states [13], “radiomics: images are data, not just pictures,” radiomics analysis will play an increasingly important role in clinical work.

The LASSO regression was adopted to decrease redundant radiomics features. By shrinking irrelevant variables to zero and only maintaining useful features, LASSO could effectively reduce the number of variables for model fitting, which have been demonstrated to be available for high dimensional data [18, 28]. Eventually, 6 radiomics features were selected from the multiparametric and multiregional MR images in our study. With the information on regional angiogenesis and the destruction of the blood-brain barrier, the CE-T1w sequence can well display the active and necrotic areas of the tumor. Moreover, the anatomical information of the tumor can be obtained from the T2flair sequence, such as peritumoral edema [29]. In our study, 5 out of the 6 radiomics features were associated with necrosis, including 3 GLSZM features, 1 GLDM feature, and 1 shape-based feature. GLSZM features are high-order ones, which play a significant role in measuring the local heterogeneity of tumor in respect of size and grayscale. In previous studies, it has been substantiated that the information reflected by GLSZM features is most similar to the details observed by doctors at the time of manual film reading [30]. It is considered to be an important radiomics feature for the prediction of GBM survival [3]. GLDM refers to the adjacent gray correlation feature of the matrix, which reflects the gray homogeneity of the local focus. These two groups of features are common in suggesting a significant difference in the volume and grayscale of the necrotic area between the TERTm cohort and the TERTw cohort. As a statistic of grayscale texture complexity, GLCM indicates the heterogeneity and malignancy of the enhancement region. Wavelet filters are effective in improving the texture features of images significantly. An analysis led us to find out that the radiomics features of necrosis, especially GLSZM and GLDM have a potential to be closely associated with TERT promoter mutations. For further verification, nevertheless, a larger sample size and multiple centers are still required.

Yamashita et al. [31] investigated 112 GBM patients with IDH wild type based on radiomics and revealed that TERT promoter mutations were correlated with the percentage of necrotic volume and age. Our further research demonstrated that CNV dramatically upregulated in TERTm than TERTw, rather than ONV, probably because the latter contained edema volume, which was more susceptible to steroids and antihydrophobic agents [32]. To verify the discovery, 4 models were established, and there was only one more factor (CNV) in Model A than Model B. The AUC of Model A in both training cohort and validation cohort were higher than that of Model B in the ROC analyses (0.955 vs 0.917, 0.889 vs 0.868, respectively). Performing NRI and IDI analyses in both training cohort and validation cohort, the results also illustrated that Model A had a remarkable improvement over Model B (NRI >0, IDI >0). Additionally, the area under the DCA curve of Model A was larger, and the clinical net benefit was better. Therefore, it was reasonable to believe that CNV was an independent predictor for TERT promoter mutations. Various previous researches provided a rational explanation that the levels of epidermal growth factor receptor amplification and interleukin 6 could be upregulated by TERT promoter mutations, inducing tumor angiogenesis and necrosis [33].

MRS is a noninvasive MR-based imaging technique that provides data on cellular metabolism. Plenty of published studies have confirmed that MRS can not only improve the diagnostic accuracy of glioma, grade the tumor, but also identify the radioactive necrosis and recurrence and predict survival rate [34]. There was a strong connection between the Cho/Cr and Lac and TERTm in our study. Cho peak, located at 3.22 ppm, was an indicator of myelination, cell metabolism, and glial hyperplasia. The increase of choline suggested an enhancement of cell membrane conversion, which was a sign of accelerated cell proliferation [35]. The concentration of creatine in brain tissue and the position of its peak in the spectral line were relatively consistent. The peak was located at 3.02 ppm, which was frequently taken as a control value. The increasing level of Cho/Cr indicated that the tumor cells were active in proliferation, which made the malignancy of the tumor predictable to some extent [36]. Lac peak, located at 1.32 ppm, was the product of anaerobic glycolysis and could not be measured in normal brain tissue. It was demonstrated in previous studies that Lac peak is more frequent to appear in HGG, which is supposedly attributed to the inhibition of aerobic respiration and cerebral tissue ischemia and hypoxia [36, 37]. According to our findings, the Lac peak appears more frequently in the TERTm group than in the TERTw group, which is potentially ascribed to the abnormal increase in cell membrane metabolism, energy depletion, anaerobic glycolysis, and the cell necrosis caused by the increased cell proliferation and mitosis of the tumor. Though Lip peak was verified to be associated with the activity of tumor necrosis [38], there was no predictive value demonstrated in our research, which may be due to the Lip peak being more sensitive to the placement of voxels, thus affecting the definitive results [34]. Therefore, it is necessary to conduct a further study on the Lip feature for predicting TERT promoter mutations. Age and gender were also taken into consideration for model development. From the results of univariate analyses, the probability of TERT mutation upregulated with every additional year of age (4% and 6% in training and validation cohorts, respectively). However, gender was not statistically significant in this study. You et al. [39] analyzed 887 gliomas for TERT promoter mutations based on histological and genetic backgrounds by DNA sequencing, which demonstrated that the frequency of TERT mutations increased with age. From what has been discussed above, TERTm in adults with HGG were more likely to be older, with visible Lac peaks, the higher Cho/Cr indexes, Radscore, and CNV (Figure 7).

Despite the promising results, some limitations of this study were also necessary to be further investigated. Firstly, functional magnetic resonance imaging series were not taken into account in this study, such as diffusion-weighted imaging, dynamic contrast-enhanced, and intravoxel incoherent motion. We believed that the addition of the advanced imaging features would improve the performance of the established radiomics nomogram which was expected to be tried in the next step. Secondly, it has been indicated in the literature that TERT played different roles in different IDH and ATRX phenotypes [38]. However, only the feasibility and methods to predict TERT promoter mutations were explored in our study, and no above-stratified analysis was performed, which was also the emphasis in our next research. Lastly, more advanced machine learning modeling methods was not used, such as support vector machines and neural networks, because the simplest logistic regression we used has yielded satisfactory results.

5. Conclusions

In conclusion, Age, Cho/Cr, Lac, CNV, and Radscore were important indicators for TERT promoter mutation predictions in HGG. Furthermore, tumor necrosis seems to be closely related to TERT promoter mutations. Radiomics nomogram based on multiparameter MRI and CNV could effectively predict the TERT promoter mutations in adults who suffered from HGG with high accuracy.

Data Availability

The data that support the findings of this study are available on request from the corresponding author (Guangyao Wu). The data are not publicly available because of the aforementioned data containing information that could compromise research participant privacy.

Conflicts of Interest

All the authors declare no conflict of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (grant no. 81227902) and the National Key Basic Research Program (grant no. 2016YFC1304702).

Supplementary Materials

Supplementary S1: The acquisition parameters of MRI. Supplementary S2: Procedures and parameters for tumor separation. Supplementary S3: Radiomics feature extraction. Supplementary S4: Statistical analysis. Supplementary S5: Radscore formula. Supplementary S6: Performance of the models for predicting TERT promoter mutations. (Supplementary Materials)