Abstract
Background:
The diagnosis of oligodendroglioma based on the latest World Health Organization Classification of Tumors of the Central Nervous System (WHO CNS 5) criteria requires the codeletion of chromosome arms 1p and 19q and isocitrate dehydrogenase gene (IDH) mutation (mut). Previously identified prognostic indicators may not be completely suitable for patients with oligodendroglioma based on the new diagnostic criteria. To find potential prognostic indicators for oligodendroglioma, we analyzed the expression of mRNAs of oligodendrogliomas in Chinese Glioma Genome Atlas (CGGA).
Methods:
We collected 165 CGGA oligodendroglioma mRNA-sequence datasets and divided them into two cohorts. Patients in the two cohorts were further classified into long-survival and short-survival subgroups. The most predictive mRNAs were filtered out of differentially expressed mRNAs (DE mRNAs) between long-survival and short-survival patients in the training cohort by least absolute shrinkage and selection operator (LASSO), and risk scores of patients were calculated. Univariate and multivariate analyses were performed to screen factors associated with survival and establish the prognostic model. qRT-PCR was used to validate the expression differences of mRNAs.
Results:
A total of 88 DE mRNAs were identified between the long-survival and the short-survival groups in the training cohort. Seven RNAs were selected to calculate risk scores. Univariate analysis showed that risk level, age, and primary-or-recurrent status (PRS) type were statistically correlated with survival and were used as factors to establish a prognostic model for patients with oligodendroglioma. The model showed an optimal predictive accuracy with a C-index of 0.912 (95% CI, 0.679–0.981) and harbored a good agreement between the predictions and observations in both training and validation cohorts.
Conclusion:
We established a prognostic model based on mRNA-sequence data for patients with oligodendroglioma. The predictive ability of this model was validated in a validation cohort, which demonstrated optimal accuracy. The 7 mRNAs included in the model would help predict the prognosis of patients and guide personalized treatment.
Introduction
Oligodendroglioma is a subtype of glioma with a relatively favorable prognosis compared with those of other entities (1–3). In contrast with the previous diagnostic criteria, which relied only on histopathological evidence (4), the integrated diagnosis of oligodendroglioma based on the World Health Organization Classification of Tumors of the Central Nervous System (WHO CNS 5) criteria requires the codeletion of chromosome arms 1p and 19q and isocitrate dehydrogenase gene (IDH) mutation (mut) (5). The WHO CNS 5 proposes an integrated diagnosis based on the consideration of pathological features and molecular profiles, such as IDH, 1p/19q, and TERT molecular markers. According to the new diagnostic criteria, the variant of oligoastrocytomas has been removed, along with the term “anaplastic” (5), which indicates the consensus that the histopathological classification is not adequate to predict the prognosis of tumors has been established (6). These revised diagnostic criteria for diffuse glioma are more reliable in predicting prognosis than classic histopathological methods (7).
According to prior reports, the 5 year survival rate of oligodendroglioma almost reaches 90% (8). However, we also observed that some patients showed a shorter survival time, irrespective of the treatment they received or the WHO grade. Moreover, the disparity in survival also suggests that a more objective and optimal stratification of oligodendrogliomas that is not merely dependent on histopathological evidence is urgently needed.
To identify potential prognostic factors and establish a prognostic estimation model for oligodendrogliomas in the context of the WHO CNS 5 classification, we analyzed the data from the China Glioma Genome Atlas (CGGA) database and the information of patients in the Department of Neurosurgical Oncology, Beijing Tiantan Hospital, Capital Medical University.
Materials and methods
CGGA data
RNA sequencing and corresponding clinical data were retrieved from the CGGA database (DataSet ID: mRNAseq_693 and mRNAseq_325) (9–13). A total of 165 glioma datasets with IDH mutation and 1p/19q co-deletion were selected and divided into two cohorts (109 in the training cohort and 56 in the validation cohort). The flowchart of the study is shown in Supplementary Figure 1. The corresponding clinical data of the patients with oligodendroglioma are summarized in Table 1.
Table 1
| Variables | Training cohort (n = 109) | Validation cohort (n = 56) | P-value |
|---|---|---|---|
| Survival status | 0.9161 | ||
| Alive | 77 | 40 | |
| Dead | 32 | 16 | |
| PRS status | 0.5614 | ||
| Primary | 65 | 36 | |
| Recurrent | 44 | 20 | |
| Gender | 0.0709 | ||
| Male | 50 | 34 | |
| Female | 59 | 22 | |
| Age | 0.2618 | ||
| ≤ 60 | 103 | 55 | |
| >60 | 6 | 1 | |
| Adjuvant therapy | 0.3429 | ||
| Radiotherapy + Chemotherapy | 52 | 26 | |
| Chemotherapy only | 17 | 25 | |
| Radiotherapy only | 30 | 3 | |
| None | 10 | 2 |
Corresponding clinical data of the patients in two cohorts.
Quantitative reverse transcription polymerase chain reaction using clinical samples
Tumor tissues of the patient collected from the Department of Neurosurgical Oncology of Beijing Tiantan Hospital were used for qRT–PCR. We chose patients with IDH mut and 1p/19q-codeleted oligodendroglioma and obtained their cryopreserved tumor samples along with the corresponding clinical data. Total RNA was isolated from tumor tissues. M-MLV reverse transcriptase (TaKaRa Bio, Japan) was used to synthesize cDNA. qRT–PCR was performed using a Roche LightCycler® 480II system (Roche, Switzerland). GAPDH was used as an internal reference to quantify the mRNAs. qRT–PCR assays were performed in triplicate. The primer and probe sequences of the genes used for qRT–PCR analysis are listed in Supplementary Table 1. Using patients with a survival time of more than 5 years as the control group, cycle threshold (CT) values were converted to relative expression to validate the expression differences.
Screening of differentially expressed mRNAs and calculating risk score
Since many low-grade glioma studies have used 5-year survival as an indicator of long-term survival (14, 15), we adopted 1,825 days (≈5 years) as the cutoff point to classify the length of patient survival time. The DESeq2 package was used to filter the differentially expressed mRNAs (DE mRNAs) between the long survival (>1,825 days) and short survival patients (≤ 1,825 days) in the training cohort with the criteria of an |log2-fold-change| ≥ 1.0 and a false discovery rate of the p < 0.001. The package made use of empirical Bayes techniques to estimate priors for log-fold-change and dispersion and to calculate posterior estimates for these quantities (16). The least absolute shrinkage and selection operator (LASSO) regression analysis was used to select the most useful predictive mRNAs from the mRNA-sequence and survival data (17). The risk score was calculated for each patient and was used to divide patients into two different groups (high-risk and low-risk). Survival analysis conducted with the DE mRNA accompanied by the corresponding clinical data was performed using the survival package. The ROCR, glmnet, and caret packages were used to assess the accuracy of the prediction of survival by risk level. Receiver operating characteristic (ROC) curves were generated to validate the accuracy of the risk score for predicting patients' OS.
Establishment of a prognostic model
Clinical data were combined with the risk scores of patients. A multivariate Cox proportional hazards model was constructed using the significant parameters identified by univariate analysis of all factors in the combined data. The multivariable Cox regression analysis included the following factors: age (>60 or not), primary-or-recurrent status (PRS) type, and risk level. A prognostic model was generated to predict the survival time of patients suffering from oligodendroglioma. To provide the clinician with an intuitive method for assessing the possibility of survival at a specific time, we established a prognostic nomogram model prepared by binary logistics regression analysis in the training cohort.
Apparent performance of the prognostic model in the training cohort
The Hosmer–Lemeshow test was used to verify the accuracy of the prognostic model. A calibration curve was plotted to visualize the results, and the ability of the model to predict a prognosis of death was evaluated by ROC. The model was subjected to bootstrapping validation (1,000 bootstrap resamples) to plot the calibration curves. Harrell's concordance index (C-index) was used to calculate the degree of discrimination between the predicted value and the observed value.
Independent validation of the prognostic model
The performance of the prognostic model was tested in a validation cohort. The Cox regression formula established in the training cohort was applied to all patients in the validation cohort. The C-index and calibration curve were used to evaluate the accuracy of our prediction model.
Statistical analysis
Statistical analysis was conducted using R software (version 4.1.3. http://www.r-project.org). The packages in R that were used in this study are reported in Supplementary Table 2. The reported statistical significance levels were all obtained from two-sided tests, and a p < 0.05 was considered statistically significant.
Results
Identification of DE mRNAs
A total of 165 oligodendroglioma datasets from the CGGA mRNA-sequencing database were selected, including 109 datasets in the training cohort and 56 datasets in the validation cohort. The baseline of the two cohorts showed no statistically significant differences (Table 1). After normalization, 88 mRNAs were identified to be significantly differentially expressed (|log2-fold-change| ≥ 1.0 and false discovery rate of p < 0.001) between the long-survival and short-survival patients in the training cohort. The expression levels of the DE mRNAs are presented in Supplementary Figure 2.
LASSO gene selection and risk score calculation
In total, 88 DE mRNAs were reduced to 7 mRNAs (LOXL2, APLN, TROAP, NUF2, CTHRC1, SLC12A5, and ANGPT2) in the LASSO regression model using the minimum criteria for risk factors (Figures 1A,B). These mRNAs were features with non-zero coefficients in the LASSO logistic regression model. This approach can be used not only to select mRNAs based on the strength of their univariate association with survival but also to calculate a risk score and perform risk stratification. The distribution of risk scores is shown in Figures 1C,D. The distribution showed that as the risk score increased, the survival time of patients decreased, and the mortality rate increased. The area under the ROC curve (AUC) was 0.679, indicating that the risk score had acceptable prediction accuracy for survival.
Figure 1
Establishment of the individualized prognostic model
The indicators with p < 0.05 in multivariate analysis (age, PRS status, and risk level) were selected as components of the established prognostic model. Age > 60 years, recurrent tumors, and high-risk level predicted a worse prognosis in patients with oligodendroglioma. With these parameters, a prognostic model was established. A nomogram was generated to present the results (Figure 2A). The ROC curve indicated that the predictive ability of the model was improved after age and PRS status was included (AUC: 0.738 vs. 0.679) (Figure 2B). Compared with the prognostic model merely based on baseline data, the predictive value of the combined model showed higher accuracy (Figure 3). At the same time, based on the prognostic model, the DynNom package and the shinyPredict package were used to generate scripts to individually predict the prognosis of patients with oligodendroglioma (Figures 4A,B).
Figure 2
Figure 3
Figure 4
Evaluation of the prognostic model in the training cohort
The calibration curves of the prognostic model for the probability of survival at 3, 5, and 8 years showed similarities in the observed and predicted values in the training cohort (Figure 4C). The Hosmer–Lemeshow test suggested that there was no departure from perfect fit (p = 0.9997; Table 2). The prognostic model yielded a C-index of 0.912 (95% CI, 0.679–0.981) in the training cohort, which implied that the model had good accuracy for predicting prognosis.
Table 2
| Interval | Alive | Dead | Predicted alive | Predicted dead |
|---|---|---|---|---|
| [0.118, 0.125) | 38 | 6 | 38.829471 | 5.170529 |
| [0.125, 0.153) | 15 | 2 | 14.392965 | 2.607035 |
| [0.153, 0.422) | 6 | 3 | 5.205934 | 3.794066 |
| [0.422, 0.498) | 16 | 17 | 16.571630 | 16.428370 |
| [0.498, 0.889] | 2 | 4 | 2.000000 | 4.000000 |
The predicted value in the Hosmer–Lemeshow test is close to the observed value.
Independent validation of the prognostic model
The calibration curves of the prognosis model showed good agreement between the predictions and observations in the validation cohort (Figure 4D). The Hosmer–Lemeshow test yielded a non-significant p-value (p = 0.8921), and the C-index of the nomogram for the prediction of survival was 0.778 (95% CI, 0.769–0.787). The prognosis model predicted a good consistency between the probability of death among patients and the actual percentage of death in the observed population.
Validation of RNA expression differences between tumor tissues of patients
Samples from the long-survival and short-survival patients were included for qRT–PCR to calculate the relative expression levels. We sought to validate the differential expression of the seven target genes (LOXL2, APLN, SLC12A5, TROAP, NUF2, ANGPT2, and CTHRC1) among these samples (Figure 5A). Unfortunately, although the differences in the mean relative gene expression between groups were consistent with expectations, except for APLN and ANGPT2, there were no statistically significant differences in gene expression between the two groups (Figure 5B). Due to the limitations of preservative methods and the long survival time of patients with oligodendroglioma, a suitable number of samples was difficult to obtain, which resulted in an insufficient sample size.
Figure 5
Discussion
Using the clinical and RNA-seq information contained in the CGGA database, we identified 7 mRNAs that could be used to predict prognosis. Then, we established a prognosis model for patients with oligodendroglioma based on seven selected mRNAs, which showed high predictive accuracy. Accurate prognosis prediction is an important component of the individualized treatment of tumors (18). Several studies aimed to find more accurate strategies to predict the prognosis of patients with oligodendroglioma (19, 20). Cao et al. analyzed the information of 4,568 patients with oligodendroglioma and established a prognostic model based on the Surveillance, Epidemiology, and End Results (SEER) database. Final results showed that radiotherapy, age, tumor location, grade, and surgical resection were independent prognostic factors of oligodendroglioma. However, they only focused on oligodendrogliomas diagnosed by histopathological criteria, and the molecular pathology factors have been neglected (19).
In addition to assessments of pathology physiological changes, consideration of the patients in daily life was also a necessary part of estimating the survival of patients with cancer (21). The survival analysis of high-grade oligodendroglioma by Liu et al. was more comprehensive and more detailed in the collection of clinical information (20). Family circumstances were also included in the survival analysis. This prognosis model was consistent with the International Classification of Functioning, Disability, and Health guidelines. Their model introduced new metrics for the establishment of prognostic models. The disadvantage of this model was that it lacked molecular pathological information for a more accurate diagnosis of patients.
Unlike previous studies, our study is one of the rare prognostic studies of oligodendroglioma combined with clinical and molecular data. In our study, the AUC value of the ROC curve of the prognostic model that integrated clinical and molecular data was higher than that of the prognostic model that contained only molecular data, and the C-index value of our prognosis model was higher than that of the Cao L study, which contained only clinical information. These results demonstrated that the integration of clinical and molecular factors could predict prognosis more accurately. In the future, more clinical information can be collected from patients for larger and more detailed prognostic analysis to establish a more accurate individualized prognostic prediction model.
Biomarkers, as prognostic indicators, are always of great interest to researchers. Some studies regarding the molecular markers associated with the survival of patients with oligodendroglioma have aimed to stratify oligodendroglioma, with an attempt to identify possible therapeutic targets and improve the survival of patients (22–26). Several previous studies that aimed to identify prognostic indicators of oligodendroglioma proposed the use of CDKN2A/B, PTEN, NOTCH1, and other biomarkers as classification criteria to reclassify oligodendroglioma, but there is no consistent conclusion (22–25). These prognostic indicators that were assessed in previous studies were discussed based on histology-confirmed oligodendrogliomas rather than an integrated diagnosed oligodendroglioma. For oligodendroglioma with IDH mutation and 1p/19q codeletion, our previous work showed that patients with oligodendroglioma with 1q/19p copolysomy had a worse prognosis (26).
In this study, we found that 7 DE mRNAs, namely, LOXL2, APLN, TROAP, NUF2, CTHRC1, SLC12A5, and ANGPT2, showed the greatest impact on prognosis. These mRNAs may play a vital role in the tumorigenesis of oligodendroglioma. Previous studies have shown that LOXL2, a member of the lysyl oxidase (LOX) family, not only promotes glioma cell proliferation, migration, and invasion and induces the epithelial-to-mesenchymal transition (EMT) process but also reduces the sensitivity of glioma cells to temozolomide (TMZ) (27). APLN is activated by VEGF signaling and hypoxia-responsive elements in the APLN promoter, stimulates angiogenic sprouting, and plays a necessary and sufficient role in tumor angiogenesis (28). TROAP activates the Wnt/β-catenin pathway and upregulates the expression of its downstream targets to play a tumor-promoting role (29). ANGPT2 activates angiogenesis through VEGFA, normalizes tumor blood vessels, and promotes the malignant transformation of glioblastoma (30, 31). NUF2 has a potential role in glioma growth and TMZ resistance (32). The CTHRC1 gene contributes to tissue repair in vascular remodeling in response to injury by limiting collagen matrix deposition and promoting cell migration (33). In addition to the mRNAs with increased expression mentioned above, which are related to cell differentiation and proliferation, the decreased expression of SLC12A5 aroused our interest. SLC12A5 (encoding the KCC2 protein) acts to stabilize nerve cell potential, and its reduced expression correlates with the development of epilepsy (34–39). Consistent with our findings, Yang and Gao et al. found that the expression of SLC12A5 in patients with shorter survival time was significantly lower than that in patients with longer survival (40). On the one hand, this finding may be explained by the fact that epilepsy is associated with IDH1 mutations in low-grade gliomas (41), and IDH mutations predict better glioma prognosis. On the other hand, seizures may cause some patients to seek treatment more actively, which may result in a better prognosis (42, 43). These DE mRNAs can be used not only to predict the prognosis of patients with oligodendroglioma but also to provide directions for potential therapeutic targets of oligodendroglioma. ISL2, a nuclear and chromatin-associated transcription factor (44), regulates the transcription of ANGPT2 by binding to the ANGPT2 promoter. When ISL2 expression was found to be reduced, oligodendroglioma cell proliferation was reduced. TMZ combined with anti-ISL2 therapy may be effective in oligodendroglioma in vitro and tumor-bearing animal models (30). For this possible treatment, more clinical studies are required to study its safety and efficiency.
In our study, we established a prognosis model for patients with oligodendroglioma. The model might be used to assign individualized treatment plans for patients with oligodendroglioma. Although excellent predictive models were identified, some limitations still existed in our study. Constrained by limited conditions, the study only contained data from CGGA and lacked verification performed using other databases. Most of the cases in the CGGA database are patients treated in our hospital, so there may have been some bias in the selection of patients included in the analysis. At the same time, our study was a retrospective study, and the treatments that patients received were not strictly consistent, which may decrease the power of our model. Additionally, studies of molecular mechanisms are needed in the future to explain the differential expression of mRNAs among patients with oligodendroglioma.
Conclusion
We established an individualized prognostic model for patients with 1p/19q codeleted and IDH mutant oligodendroglioma based on mRNA seq data. The model would help predict the prognosis of patients and guide personalized treatment.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by Medical Ethics Committee of Tiantan Hospital, Beijing Tiantan Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
QZ and SL: conception and design. QZ, MingxL, XZhan, XZhao, YC, and MingL: collection and assembly of data. QZ, HL, SS, XR, and CY: data analysis and interpretation. QZ: cell biological experiments and manuscript writing. All authors: final approval of manuscript.
Funding
This study was supported by the Capital Funds for Health Improvement and Research (2020-2-1075).
Acknowledgments
The authors sincerely thank the patients and their families for their participation in this study. We are grateful to CGGA for its open database for researchers.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.1074593/full#supplementary-material
Supplementary Figure 1Flowchart of the study.
Supplementary Figure 2Heatmap of differentially expressed mRNAs (DE mRNAs).
Supplementary Table 1Information of primer sequences of the genes used for qRT–PCR analysis.
Supplementary Table 2R packages used in the study.
References
1.
TakahashiYNakamuraHMakinoKHideTMutaDKamadaHet al. Prognostic value of isocitrate dehydrogenase 1, O6-methylguanine-DNA methyltransferase promoter methylation, and 1p19q co-deletion in Japanese malignant glioma patients. World J Surg Oncol. (2013) 11:284. 10.1186/1477-7819-11-284
2.
TanakaKSasayamaTMizukawaKTakataKSulaimanNSNishiharaMet al. Combined IDH1 mutation and MGMT methylation status on long-term survival of patients with cerebral low-grade glioma. Clin Neurol Neurosurg. (2015) 138:37–44. 10.1016/j.clineuro.2015.07.019
3.
FranceschiEMuraADe BiaseDTalliniGPessionAFoschiniMPet al. The role of clinical and molecular factors in low-grade gliomas: what is their impact on survival?Future Oncol. (2018) 14:1559–67. 10.2217/fon-2017-0634
4.
EngelhardHH. Current diagnosis and treatment of oligodendroglioma. Neurosurg Focus. (2002) 12:E2. 10.3171/foc.2002.12.2.3
5.
LouisDNPerryAWesselingPBratDJCreeIAFigarella-BrangerDet al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol. (2021) 23:1231–51. 10.1093/neuonc/noab106
6.
ChenRJLuMYWangJWilliamsonDRodigSJLindemanNIet al. Pathomic fusion: an integrated framework for fusing histopathology and genomic features for cancer diagnosis and prognosis. IEEE Trans Med Imaging. (2022) 41:757–70. 10.1109/TMI.2020.3021387
7.
PareiraESShibuyaMOharaKNakagawaYKanazawaTKamamotoDet al. The oligodendroglial histological features are not independently predictive of patient prognosis in lower-grade gliomas. Brain Tumor Pathol. (2022) 39:79–87. 10.1007/s10014-022-00426-5
8.
JaeckleKABallmanKVvan den BentMGianniniCGalanisEBrownPDet al. CODEL: phase III study of RT, RT + TMZ, or TMZ for newly diagnosed 1p/19q codeleted oligodendroglioma. Analysis from the initial study design. Neuro Oncol. (2021) 23:457–67. 10.1093/neuonc/noaa168
9.
ZhaoZZhangKNWangQLiGZengFZhangYet al. Chinese glioma genome atlas (CGGA): a comprehensive resource with functional genomic data from Chinese glioma patients. Genomics Proteomics Bioinformatics. (2021) 19:1–12. 10.1016/j.gpb.2020.10.005
10.
ZhaoZMengFWangWWangZZhangCJiangT. Comprehensive RNA-seq transcriptomic profiling in the malignant progression of gliomas. Sci Data. (2017) 4:170024. 10.1038/sdata.2017.24
11.
WangYQianTYouGPengXChenCYouYet al. Localizing seizure-susceptible brain regions associated with low-grade gliomas using voxel-based lesion-symptom mapping. Neuro Oncol. (2015) 17:282–8. 10.1093/neuonc/nou130
12.
LiuXLiYQianZSunZXuKWangKet al. A radiomic signature as a non-invasive predictor of progression-free survival in patients with lower-grade gliomas. Neuroimage Clin. (2018) 20:1070–7. 10.1016/j.nicl.2018.10.014
13.
BaoZSChenHMYangMYZhangCBYuKYeWLet al. RNA-seq of 272 gliomas revealed a novel, recurrent PTPRZ1-MET fusion transcript in secondary glioblastomas. Genome Res. (2014) 24:1765–73. 10.1101/gr.165126.113
14.
ChintagumpalaMEckelSPKrailoMMorrisMAdesinaAPackerRet al. A pilot study using carboplatin, vincristine, and temozolomide in children with progressive/symptomatic low-grade glioma: a children's oncology group studydagger. Neuro Oncol. (2015) 17:1132–8. 10.1093/neuonc/nov057
15.
ArmstrongGTConklinHMHuangSSrivastavaDSanfordREllisonDWet al. Survival and long-term health and cognitive outcomes after low-grade glioma. Neuro Oncol. (2011) 13:223–34. 10.1093/neuonc/noq178
16.
LoveMIHuberWAndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. 10.1186/s13059-014-0550-8
17.
SauerbreiWRoystonPBinderH. Selection of important variables and determination of functional form for continuous predictors in multivariable model building. Stat Med. (2007) 26:5512–28. 10.1002/sim.3148
18.
WeiRLZhangLWLiJGYangFDXueYKWeiXT. Behavior-Oriented nomogram for the stratification of lower-grade gliomas to improve individualized treatment. Front Oncol. (2020) 10:538133. 10.3389/fonc.2020.538133
19.
CaoLRongPZhuGXuAChenS. Clinical characteristics and overall survival prognostic nomogram for oligodendroglioma: a surveillance, epidemiology, and end results population-based analysis. World Neurosurg. (2021) 151:e810–20. 10.1016/j.wneu.2021.04.122
20.
LiuSLiuXXiaoYChenSZhuangW. Prognostic factors associated with survival in patients with anaplastic oligodendroglioma. PLoS ONE. (2019) 14:e211513. 10.1371/journal.pone.0211513
21.
TschiesnerULinseisenECoenenMRogersSHarreusUBerghausAet al. Evaluating sequelae after head and neck cancer from the patient perspective with the help of the international classification of functioning, disability and health. Eur Arch Otorhinolaryngol. (2009) 266:425–36. 10.1007/s00405-008-0764-z
22.
WijnengaMFrenchPJDubbinkHJDinjensWAtmodimedjoPNKrosJMet al. Prognostic relevance of mutations and copy number alterations assessed with targeted next generation sequencing in IDH mutant grade II glioma. J Neurooncol. (2018) 139:349–57. 10.1007/s11060-018-2867-8
23.
HalaniSHYousefiSVelazquezVJRossiMRZhaoZAmrollahiFet al. Multi-faceted computational assessment of risk and progression in oligodendroglioma implicates NOTCH and PI3K pathways. NPJ Precis Oncol. (2018) 2:24. 10.1038/s41698-018-0067-9
24.
AppayRDehaisCMaurageCAAlentornACarpentierCColinCet al. CDKN2A homozygous deletion is a strong adverse prognosis factor in diffuse malignant IDH-mutant gliomas. Neuro Oncol. (2019) 21:1519–28. 10.1093/neuonc/noz124
25.
AokiKNakamuraHSuzukiHMatsuoKKataokaKShimamuraTet al. Prognostic relevance of genetic alterations in diffuse lower-grade gliomas. Neuro Oncol. (2018) 20:66–77. 10.1093/neuonc/nox132
26.
JiangHRenXZhangZZengWWangJLinS. Polysomy of chromosomes 1 and 19: an underestimated prognostic factor in oligodendroglial tumors. J Neurooncol. (2014) 120:131–8. 10.1007/s11060-014-1526-y
27.
ZhangQYangLGuanGChengPChengWWuA. LOXL2 upregulation in gliomas drives tumorigenicity by activating autophagy to promote TMZ resistance and trigger EMT. Front Oncol. (2020) 10:569584. 10.3389/fonc.2020.569584
28.
KalinREGlassR. APLN/APLNR signaling controls key pathological parameters of glioblastoma. Cancers. (2021) 13:3899. 10.3390/cancers13153899
29.
ZhaoZQWuXJChengYHZhouYFMaXMZhangJet al. TROAP regulates cell cycle and promotes tumor progression through Wnt/beta-Catenin signaling pathway in glioma cells. CNS Neurosci Ther. (2021) 27:1064–76. 10.1111/cns.13688
30.
QiLWangZYShaoXRLiMChenSNLiuXQet al. ISL2 modulates angiogenesis through transcriptional regulation of ANGPT2 to promote cell proliferation and malignant transformation in oligodendroglioma. Oncogene. (2020) 39:5964–78. 10.1038/s41388-020-01411-y
31.
Di TacchioMMacasJWeissenbergerJSommerKBahrOSteinbachJPet al. Tumor vessel normalization, immunostimulatory reprogramming, and improved survival in glioblastoma with combined inhibition of PD-1, angiopoietin-2, and VEGF. Cancer Immunol Res. (2019) 7:1910–27. 10.1158/2326-6066.CIR-18-0865
32.
GuoLWuZ. FOXM1-mediated NUF2 expression confers temozolomide resistance to human glioma cells by regulating autophagy via the PI3K/AKT/mTOR signaling pathway. Neuropathology. (2022) 42:430–46. 10.1111/neup.12824
33.
MeiPJBaiJMiaoFAChenCZhuYSLiZLet al. CTHRC1 mediates multiple pathways regulating cell invasion, migration and adhesion in glioma. Int J Clin Exp Pathol. (2017) 10:9318–29.
34.
StodbergTMctagueARuizAJHirataHZhenJLongPet al. Mutations in SLC12A5 in epilepsy of infancy with migrating focal seizures. Nat Commun. (2015) 6:8038. 10.1038/ncomms9038
35.
SaitsuHWatanabeMAkitaTOhbaCSugaiKOngWPet al. Impaired neuronal KCC2 function by biallelic SLC12A5 mutations in migrating focal seizures and severe developmental delay. Sci Rep. (2016) 6:30072. 10.1038/srep30072
36.
SaitoTIshiiASugaiKSasakiMHiroseS. A de novo missense mutation in SLC12A5 found in a compound heterozygote patient with epilepsy of infancy with migrating focal seizures. Clin Genet. (2017) 92:654–8. 10.1111/cge.13049
37.
PuskarjovMSejaPHeronSEWilliamsTCAhmadFIonaXet al. A variant of KCC2 from patients with febrile seizures impairs neuronal Cl-extrusion and dendritic spine formation. EMBO Rep. (2014) 15:723–9. 10.1002/embr.201438749
38.
KahleKTMernerNDFriedelPSilayevaLLiangBKhannaAet al. Genetically encoded impairment of neuronal KCC2 cotransporter function in human idiopathic generalized epilepsy. EMBO Rep. (2014) 15:766–74. 10.15252/embr.201438840
39.
TillASzalaiRHegyiMKovesdiEBukiGHadzsievKet al. [A rare form of ion channel gene mutation identified as underlying cause of generalized epilepsy]. Orv Hetil. (2019) 160:835–8. 10.1556/650.2019.31404
40.
YangZGongWZhangTGaoH. Molecular features of glioma determined and validated using combined TCGA and GTEx data analyses. Front Oncol. (2021) 11:729137. 10.3389/fonc.2021.729137
41.
LiYShanXWuZWangYLingMFanX. IDH1 mutation is associated with a higher preoperative seizure incidence in low-grade glioma: a systematic review and meta-analysis. Seizure. (2018) 55:76–82. 10.1016/j.seizure.2018.01.011
42.
StuppRHegiMEMasonWPvan den BentMJTaphoornMJJanzerRCet al. Effects of radiotherapy with concomitant and adjuvant temozolomide vs. radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-Year analysis of the EORTC-NCIC trial. Lancet Oncol. (2009) 10:459–66. 10.1016/S1470-2045(09)70025-7
43.
OhbaSHiroseY. Biological significance of mutant isocitrate dehydrogenase 1 and 2 in gliomagenesis. Neurol Med Chir. (2016) 56:170–9. 10.2176/nmc.ra.2015-0322
44.
OzturkHCingozHTufanTYangJAdairSJTummalaKSet al. ISL2 is a putative tumor suppressor whose epigenetic silencing reprograms the metabolism of pancreatic cancer. Dev Cell. (2022) 57:1331–46. 10.1016/j.devcel.2022.04.014
Summary
Keywords
oligodendroglioma, 1p/19q codeletion, prognostic model, WHO CNS 5, mRNA-sequence
Citation
Zhu Q, Shen S, Yang C, Li M, Zhang X, Li H, Zhao X, Li M, Cui Y, Ren X and Lin S (2022) A prognostic estimation model based on mRNA-sequence data for patients with oligodendroglioma. Front. Neurol. 13:1074593. doi: 10.3389/fneur.2022.1074593
Received
19 October 2022
Accepted
21 November 2022
Published
14 December 2022
Volume
13 - 2022
Edited by
Quan Cheng, Xiangya Hospital, Central South University, China
Reviewed by
Shufa Zheng, First Affiliated Hospital of Fujian Medical University, China; Qianquan Ma, Peking University Third Hospital, China
Updates
Copyright
© 2022 Zhu, Shen, Yang, Li, Zhang, Li, Zhao, Li, Cui, Ren and Lin.
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: Song Lin linsong2005@126.com
This article was submitted to Neuro-Oncology and Neurosurgical Oncology, a section of the journal Frontiers in Neurology
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.