- 1Department of Neurosurgery, Center for Malignant Brain Tumors, National Glioma MDT Alliance, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 2China Anti-Cancer Association Specialty Committee of Glioma, Peking Union Medical College Hospital, Beijing, China
- 3’4 + 4’ Medical Doctor Program, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 4Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
- 5Research Unit of Accurate Diagnosis, Treatment, and Translational Medicine of Brain Tumors (No.2019RU011), Chinese Academy of Medical Sciences, Beijing, China
- 6Eight-year Medical Doctor Program, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Introduction: Glioblastoma (GBM), the most lethal primary brain malignancy, is divided into histological (hist-GBM) and molecular (mol-GBM) subtypes according to the 2021 World Health Organization classification of central nervous system tumors. This study aimed to characterize the clinical, radiological, molecular, and survival features of GBM under the current classification scheme and explore survival determinants.
Methods: We re-examined the genetic alterations of IDH-wildtype diffuse gliomas at our institute from 2011 to 2022, and enrolled GBMs for analysis after re-classification. Univariable and multivariable analyses were used to identify survival determinants.
Results: Among 209 IDH-wildtype gliomas, 191 were GBMs, including 146 hist-GBMs (76%) and 45 mol-GBMs (24%). Patients with mol-GBMs were younger, less likely to develop preoperative motor dysfunction, and more likely to develop epilepsy than hist-GBMs. Mol-GBMs exhibited lower radiographic incidences of contrast enhancement and intratumoral necrosis. Common molecular features included copy-number changes in chromosomes 1, 7, 9, 10, and 19, as well as alterations in EGFR, TERT, CDKN2A/B, and PTEN, with distinct patterns observed between the two subtypes. The median overall survival (mOS) of GMB was 12.6 months. Mol-GBMs had a higher mOS than hist-GBMs, although not statistically significant (15.6 vs. 11.4 months, p=0.17). Older age, male sex, tumor involvement of deep brain structure or functional area, and genetic alterations in CDK4, CDK6, CIC, FGFR3, KMT5B, and MYB were predictors for a worse prognosis, while MGMT promoter methylation, maximal tumor resection, and treatment based on the Stupp protocol were predictive for better survival.
Conclusion: The definition of GBM and its clinical, radiological, molecular, and prognostic characteristics have been altered under the current classification.
1 Introduction
Glioblastoma (GBM), the most common and lethal primary malignant brain tumor worldwide, is associated with a dismal prognosis (1, 2). Despite continuous efforts to improve glioblastoma-related outcomes, the median overall survival (mOS) is only 16 months, and the real-world 5-year survival rate is 6.8% (3). Moreover, glioblastomas invariably recur as aggressive, therapy-resistant relapses, highlighting the need for more effective treatment strategies (4, 5). Before the publication of the fifth edition of the World Health Organization (WHO) classification of central nervous system (CNS) tumors, known as the “WHO CNS5 classification” in June 2021, GBM classification primarily relied on histology. The traditional hallmarks of GBM were astrocytomas with the features of central pseudo-palisading necrosis and abnormal microvascular proliferation under the microscope (6).
Advancements in molecular testing have improved our understanding of the molecular categorization of gliomas, leading to a debate regarding whether GBM should be divided into isocitrate dehydrogenase (IDH)-mutant and IDH-wildtype subcategories (7, 8). IDH-mutant gliomas differ fundamentally from IDH-wildtype gliomas in terms of metabolism, epigenetics, biological behavior, aggressive infiltration, vulnerable population, and response to therapy (9–12). The Consortium to Inform Molecular and Practical Approaches to CNS Tumor Taxonomy (cIMPACT-NOW) proposed that IDH-wildtype diffuse gliomas that contain specific-altered genes strongly suggest a poor prognosis similar to GBM, even if the tumor exhibits low-grade histological features (11). Hence, in the WHO CNS5 classification scheme, molecular information and histological features have been incorporated for an integrated diagnosis that has substantially changed the overall classification of gliomas (13–15). Notably, one of the most significant changes in adult-type diffuse gliomas is the reclassification of ‘GBM, IDH-mutant, WHO grade 4’ as ‘Astrocytoma, IDH-mutant, WHO grade 4’ (10). GBM are now defined as diffuse, astrocytic gliomas that are IDH-wildtype and H3-wildtype, exhibiting one or more of the following histological or genetic features: microvascular proliferation, necrosis, telomerase reverse transcriptase (TERT) promoter mutation, epidermal growth factor receptor (EGFR) gene amplification, a concomitant gain of chromosome 7 and loss of chromosome 10 (+7/-10) (CNS WHO grade 4) (14, 16, 17). Histological diagnosis of GBM was determined by the presence of necrosis or microvascular proliferation, which is referred to as histological GBM (hist-GBM). IDH-wildtype diffuse astrocytic tumors without the histological features of GBM, which would have otherwise been classified as grade 2 or 3, are considered as molecular GBM (mol-GBM, WHO grade 4) if they harbor any of the following molecular abnormalities: TERT promoter mutation, EGFR amplification, or chromosomal + 7/−10 copy changes. An emerging concern is whether the results from previous studies fit GBMs as defined by the current classification system. It is also unclear whether the low-histological-grade mol-GBMs share the same clinical, radiological, pathological, molecular, and prognostic characteristics as hist-GBMs. Accumulating prospective biomarker alterations for GBM have been confirmed to have major consequences on the management of GBM, such as O-6-methylguanine-DNA methyltransferase (MGMT) promoter methylation, phosphatase and tension homolog (PTEN) loss, and histone methyltransferase KMT5B alteration (alias SUV420H1) (18–20). However, these gene alterations were not included in the current classification system for subtype categorization of gliomas due to the limited availability of evidence-based data in real-world settings(16, 21).
The present study aimed to comprehensively describe the clinical presentations, radiological features, pathological characteristics, and molecular alterations of GBMs according to the updated WHO CNS5 classification, as well as to assess patient survival and clinical determinants, and compare the aforementioned biological characteristics between hist-GBMs and mol-GBMs. To achieve these objectives, we collected and analyzed data from GBM patients treated at our institute over 11 years. In addition, we utilized public datasets to validate our findings, with the ultimate goal of establishing a more solid foundation for clinical recognition and decision-making of hist-GBMs and mol-GBMs.
2 Methods
2.1 Study participants
Patients with IDH-wildtype and H3-wildtype diffuse gliomas who underwent surgical tumor resection at Peking Union Medical College Hospital Neurosurgery between January 2011 and January 2022 were retrospectively screened. Patients’ tumors were defined as GBM (WHO grade 4), anaplastic astrocytomas (WHO grade 3), and diffuse astrocytomas (WHO grade 2) by the WHO CNS4 classification, while were recently re-classified as GBM, diffuse gliomas (NOS), or certain types of pediatric-type diffuse gliomas using the WHO CNS5 classification. In this analysis, only patients with GBMs as defined by the WHO CNS5 classification were enrolled and further analyzed.
The present study was approved by the institutional review board at Peking Union Medical College Hospital. Written informed consent was obtained from every subject.
2.2 Definition of histological and molecular glioblastomas
Hist-GBMs were defined as IDH-wildtype and H3-wildtype high-grade diffuse gliomas with microvascular proliferation and/or intratumoral necrosis. Mol-GBMs were defined as histological WHO grade 2-3 IDH-wildtype and H3-wildtype diffuse gliomas with TERT promoter mutation, EGFR amplification, or +7/-10 chromosome variations. When neither of the aforementioned criteria was met, the tumor was classified as diffuse gliomas (NEC) or pediatric-type diffuse glioma, including MYB- or MYBL1-altered low-grade diffuse astrocytoma and H3-wildtype and IDH-wildtype pediatric-type high-grade diffuse glioma.
2.3 Data acquisition
Clinical information, radiological characteristics, and histological data were extracted from the electronic medical record (Supplementary Table 1). Sixty molecular markers of interest, including gene alterations and chromosome copy number variations, were analyzed using the next-generation sequencing method, the PCR-based assay, and the fluorescence in situ hybridization. These markers were chosen to differentiate subtypes of diffuse gliomas or to predict patient prognosis.
2.4 Targeted sequencing
Genomic DNA from formalin-fixed paraffin-embedded (FFPE) tissue was isolated using the DNA FFPE Tissue Kit (Qiagen #56404). A custom glioma next-generation sequencing (NGS) panel was designed to investigate detailed molecular alterations in these patients. The 60-gene NGS panel contains some interested well-known genes and recently described genes that were associated with the diagnosis, grading, and treatment responses of gliomas, including IDH1, IDH2, H3F3A, HIST1H3B, HIST1H3C, EGFR, TERT, CDKN2A/B, CDK4, CDK6, CIC, FGFR3, MYB, MYBL1, KMT5B, MGMT, PTEN, chromosome copy number variations, etc (Supplementary Table 1). According to the manufacturer’s specifications, The DNA library was constructed by DNA repair (NEBNext® FFPE DNA Repair Mix, #M6630S), DNA fragmentation and end-repair (TIANSeq Fragment/Repair/Tailing Module, #NG301), and hybridization (KAPA HyperCapture Reagent Kit, #09075828001). Among all of the 209 FFPE tissues that had been acquired and sent to NGS test, only 78 cases had adequate tissue and passed the quality control and were analyzed for the molecular part in this article. After the quality assessment, the library was sequenced with the Illumina platform and the mean read depth was >500X. GRCH38 was used as the reference genome and subsequent analyses.
2.5 Validation dataset
RNA-sequencing and clinicopathological data of GBMs were downloaded from the Chinese Glioma Genome Atlas (CGGA; http://www.cgga.org.cn/; Dataset ID: mRNAseq_693, mRNAseq_325) to obtain a validation dataset. A total of 280 patients who were diagnosed with hist-GBMs according to the WHO CNS5 classification were enrolled. Information included were gender, age, OS, and data of molecular markers. The optimal cutoff values of the expression level of each parameter were calculated using the “surv_cutpoint” function of the “survminer” R package, and the patients were consequently divided into the high and low expression groups for validation.
2.6 Statistical analyses
Clinical, radiological, and pathological data expressed as categorical variables were shown as numbers and percentages, and continuous variables were presented as the means ± standard deviations or medians plus interquartile range according to data distribution. Comparisons of categorical variables were performed using the chi-squared test. Student’s t-test was used to assess the differences between normally distributed continuous variables, while the Mann-Whitney U test was used for variables that failed the normality test. For most parameters, all patients were included in the analysis. However, only the patients with available and complete data were enrolled for analysis of several variables. p<0.05 indicated statistical significance. The Sankey diagram was used to visualize the changes in tumor classifications using the WHO CNS4 and the WHO CNS5 classification systems. The waterfall chart was performed to illustrate the GBM samples’ molecular alterations and IHC expressions. Median OS and the 95% confidence interval (95%CI) were calculated for hist-GBMs and mol-GBMs and patient subgroups with distinct clinical, radiological, or molecular features. Kaplan-Meier curves were drawn to illustrate the OS, and a log-rank p<0.05 indicated a significant difference in survival. Comprehensive sub-analyses were performed by using both the uni- and multi-variable hazard ratio (HR) models to evaluate the associations between biological variables and OS. The above analyses were performed for all glioblastomas, hist-GBMs, and mol-GBMs. To further validate the study’s results, all variables significantly associated with patient OS were included in the validation set analysis, where univariate subgroup analysis was performed for each factor using the RNA-sequencing data. SPSS Statistics (version 26.0, IBM, USA) was used to conduct the data analyses, and RStudio (PBC & Certified B Corp.®, USA) was used to generate graphs.
3 Results
3.1 Changes in the classification of IDH-wildtype diffuse gliomas
Among 209 IDH-wildtype diffuse gliomas (146 GBM, 38 anaplastic astrocytomas, and 25 diffuse astrocytomas) according to the WHO CNS4 classification, 191 (91%) were reclassified as GBM (146 hist-GBMs and 45 mol-GBMs), 11 (5%) as H3-wildtype and IDH-wildtype pediatric-type high-grade diffuse gliomas, and 7 (3%) as MYB/MYBL1-altered pediatric-type low-grade diffuse astrocytomas based on the WHO CNS5 classification system (Figure 1). Out of 38 IDH-wildtype anaplastic astrocytomas (71%) and 25 diffuse astrocytomas (72%), 27 (71%) and 18 (72%) respectively were reclassified as mol-GBMs due to the presence of key genetic features. The number of molGBMs with one, two, three genetic alterations were 11 (24.4%), 18 (40.0%), and 16 (35.6%), respectively. The remaining IDH-wildtype anaplastic astrocytomas (29%) and diffuse astrocytomas (28%) were reclassified as pediatric-type diffuse high-grade and low-grade gliomas, respectively.
Figure 1 Changes of subtypes of IDH-wildtype diffuse gliomas from the WHO CNS4 to the WHO CNS5 classification. Each rectangle represents a subtype of IDH-wildtype diffuse gliomas. The blocks on the left include glioblastomas (WHO grade 4), anaplastic astrocytomas (WHO grade 3), and diffuse astrocytoma (WHO grade 2) classified using the WHO CNS4 classification system. The blocks on the right include glioblastomas, H3-wildtype and IDH-wildtype diffuse pediatric-type high-grade gliomas, and MYB- or MYBL1-altered diffuse pediatric-type astrocytomas classified using the WHO CNS5 classification system.
3.2 Clinical presentations
The average age of GBM patients was 55.5 years, with males comprising 52.4% of all samples (Table 1). At diagnosis, 67.5% of patients exhibited neurological symptoms, with motor dysfunction (37.2%) and aphasia (21.5%) being the most common. Intracranial hypertension was found in 48.2% of patients, while 20.9% had a history of epilepsy. The median duration of the disease before admission was 5 weeks, and the median Karnofsky performance scale was 80. The majority of patients underwent gross total resection (61.8%; 36.1% for those invading deep brain structures), and 64.1% received the standard Stupp treatment protocol following surgery.
Table 1 Clinical, radiological, and pathological features of patients with histological or molecular glioblastomas.
Compared with hist-GBMs, patients with mol-GBMs were diagnosed at a younger age (50.8 ± 15.8 vs. 56.9 ± 14.1 years, p=0.015), had a lower incidence of motor dysfunction (17.8% vs. 43.2%, p=0.002), a higher incidence of epilepsy (33.3% vs. 17.1%, p=0.019) before admission, and a lower rate of gross total resection (44.4% vs. 67.1%, p=0.006).
3.3 Radiological features
The frontal (46.3%) and temporal (42.0%) lobes were the most frequently affected regions, and 52.5% of tumors invaded more than one lobe (Table 1). 22.2% of included GBMs invaded deep structures (thalamus or callosum), and two were originated from the subtentorial region. Functional brain regions were involved in 55.6% of patients, with the motor cortex and tract being the most affected (30.9%). GBM appeared hypo-intensive (63.6%) or mixed (35.2%) on T1-weighed images and hyper-intensive (52.5%) or mixed (46.9%) on T2-weighed images. The average maximal tumor diameter was 4.1 ± 1.7 cm. Most tumors exhibited contrast-enhanced (92.0%) and intratumoral necrosis (80.9%).
Compared with hist-GBMs, mol-GBMs were less likely to have contrast-enhancement (78.8% vs. 95.3%, p=0.006) and intratumoral necrosis (63.6% vs. 85.3%, p=0.005), and the median diameter of necrosis inside mol-GBMs was smaller (1.0 vs. 2.3cm, p=0.002).
3.4 Under-the-microscope pathological characteristics
The median ki-67 index was 30% for all GBMs (Table 1). Hist-GBMs had a higher median ki-67 index compared to mol-GBMs (30% vs. 10%, p<0.001). Among mol-GBMs, 60% were upgraded from WHO grade 3 IDH-wildtype anaplastic astrocytomas, and the other 40% originated from WHO grade 2 IDH-wildtype diffuse astrocytomas. No differences in immunohistochemical expressions were found between hist-GBMs and mol-GBMs.
3.5 Genetic alterations and IHC expression
Mutation was considered the main type of genetic alteration in TERT promotor, TOP3A, and BRAF; while deletions were predominantly observed in chr-10p/10q/9p, CDKN2A/B, PTEN, PTPN11, and FGFR1/2/4. Amplification was most commonly seen in chr-7p/q, EGFR, CDK6/4, MYB/MYBL1, MET, PIK3CA, NTRK2, PDGFRA, PEG3, KIT, NOTCH1, PPM1D, RB1, and KRAS (Figure 2). Immunohistochemistry results indicated that GFAP, S-100, Olig-2, and ATRX were positive in more than 90% of glioblastomas (Table 1).
Figure 2 Molecular alterations and immunohistochemical presentations of histological and molecular glioblastomas. Molecular alterations and immunohistochemical presentations of histological and molecular glioblastomas. Each column represents an individual patient, and the classification of tumors is displayed at the bottom. Each row indicates a genetic or immunohistochemical parameter, and these parameters are listed from top to bottom based on the frequency of genetic alterations and immunohistochemical expressions. Genetic mutation is shown as green, deletion is shown as blue, amplification is shown as red, ambiguous is shown as yellow, and positive immunohistochemical staining is shown as brown. The frequencies of mutation, deletion, and amplification of each gene are shown in the right histogram. Differences in molecular alterations and immunohistochemical expressions are compared, and an asterisk (*) to the left indicates a p<0.05.
Compared with mol-GBMs, hist-GBMs had higher rates of genetic alterations in PTEN, TOP3A, CDK4, MYB, KIT, KRAS, and NTRK3 (Supplementary Table 2).
3.6 Overall survival of glioblastomas as defined by the WHO CNS5 classification system
Our results showed that the median OS for all GBMs was 12.6 months (95%CI: 11.1-15.8); it was 15.6 months (95%CI: 12.5-21.9) for mol-GBMs and 11.4 months (95%CI: 10.7-15.8) for hist-GBMs (Figure 3). However, despite that mol-GBMs showed better survival, no significant difference was found between the two subtypes (HR=0.76, 95%CI: 0.52-1.12; p=0.17). Additionally, the median OS of patients with pediatric-type diffuse gliomas (the other IDH-wildtype gliomas in this study) was 55.1 months, significantly longer than that of patients with mol-GBMs (p=0.006).
Figure 3 Overall survival of patients with glioblastomas using the current classification scheme and its clinical relevance. (A): Kaplan–Meier curves showing the overall survival (OS) of molecular glioblastomas vs. histological glioblastomas (median OS: 11.4 months vs. 15.6 months, hazard ratio [HR]: 0.76, P=0.17), and the OS curve for all glioblastomas (median OS: 12.6 months). In all glioblastomas (B-I), age ≥60 years, males, postoperative treatment without Stupp protocol, MGMT promotor unmethylation, alteration of CIC, KMT5B, CDK4, FGFR3 were correlated with a shorter OS. In histological glioblastomas (J-U), age ≥60 years, postoperative treatment without standardized Stupp protocol, having chromosome 7 gain and chromosome 10 loss, alteration of CDK6, TOP3A, PTEN, MYBL1, KRAS, MYB, FGFR2, CIC, MGMT, promotor unmethylation were correlated with a shorter OS. In molecular glioblastomas (V–AA), male sex, ki-67≥30%, alteration of TP53, CDK4, FGFR3, and KMT5B were correlated with a shorter OS.
Patients aged ≥60 years (11.4 vs. 14.3, months; p=0.008), male patients (11.3 vs. 15.2, months; p=0.044), and those who did not undergo the standard Stupp treatment protocol (9.9 vs. 18.0, months; p=0.012) had a shorter OS than the others. MGMT promotor unmethylation, FGFR3 alteration, CIC alteration, KMT5B alteration, and CDK4 alteration were molecular features that predicted a shorter OS in GBM (p<0.05).
3.7 Clinical correlations of OS in histological and molecular glioblastomas
In hist-GBMs, older age (≥60 years), not receiving the Stupp treatment protocol, +7/-10 chromosome copy-number changes, MGMT unmethylation, and alterations in CDK6, TOP3A, PTEN, MYB, MYBL1, KRAS, FGFR2, and CIC were significantly correlated with a shorter OS. In mol-GBMs, male sex, ki-67 ≥30%, and alterations in TP53, CDK4, FGFR3, and KMT5B were significantly correlated with a shorter OS (Figure 3). No other molecular correlations with OS were found (Supplementary Figures 1–3).
Multivariable regression analyses showed that older age, male sex, tumor involvement of functional brain areas and deep brain structures, and alterations in CDK4, CDK6, and KMT5B were predictive of a decreased OS, while maximal surgical resection as compared to biopsy, treatment with the Stupp protocol, and MGMT promoter methylation was predictive of a prolonged OS (Table 2). CIC and MYB alterations were predictive of a decreased OS in hist-GBMs. Presentation of motor dysfunction before surgery, a higher ki-67 index, TP53 alteration, and FGFR3 alteration were predictive of a decreased OS in mol-GBMs.
Table 2 Independent clinical, radiological, and pathological factors associated with overall survival of patients with glioblastomas.
3.8 Validation of the clinical correlations of OS using the CGGA dataset
In the RNA-sequencing validation set, high expression in CDK6, CIC, CDK4, TOP3A, and MYB and low expression in TP53 and FGFR3 were all significantly correlated with a shorter OS in hist-GBMs (Figure 4).
Figure 4 Clinical relevance of overall survival of patients with histological glioblastomas using the validation set. This figure included 15 factors showing the influence of gene expression and other clinical parameters on the OS (A-O), including CDK6, TP53, FGFR3, CIC, CDK4, TOP3A, MYB, FGFR2, age, KRAS, MYBL1, KMT5B, PTEN, MGMT methylation, and sex.
4 Discussion
Since the release of the 5th edition of the WHO classification of CNS tumors, the categorization of adult-type gliomas has been altered drastically. Among primary CNS tumors, GBM has undergone the greatest changes. In the present study, over 70% of previously defined grade 2-3 IDH-wildtype diffuse gliomas were re-classified as grade 4 GBMs due to specific molecule alterations, while the remaining 30% of them were re-classified into other subtypes. Furthermore, the previously recognized IDH-mutant GBM have been re-classified as astrocytoma, WHO grade 4. Overall, the composition of GBMs nowadays differs greatly from earlier ones in that all IDH-mutant astrocytomas are excluded, and a large number of WHO grade 2-3 IDH-wildtype gliomas with specific molecular features are involved (14). Consequently, previous conclusions on GBM’s clinical, radiographic, molecular, and prognostic features might no longer be compatible with the current clinical practice. Additionally, since mol-GBMs has become a newly defined group of tumors, increasing uncertainty about its diagnostic and therapeutic principles has emerged (22). Hence, it is imperative and crucial to reassess the characteristics of GBMs in the context of the updated classification. In this study, we conducted genetic analyses on IDH-wildtype diffuse gliomas resected at our institute. A total of 191 GBMs defined by the WHO CNS5 classification were identified and included in the study. The study also monitored the survival rates of the patients and investigated potential correlations between biological parameters and their prognoses. To examine the molecular changes in glioma patients, we developed and utilized a specialized NGS panel to compare molecular and histological characteristics of GBMs. The objective of this study was to provide valuable insights for the clinical treatment of this novel subset of mol-GBMs.
Based on this new classification, current GBMs were diagnosed at a younger age and were more predominant in male patients. Also, more than two-thirds were accompanied by neurological deficits, particularly motor dysfunction, which differed slightly from the prior studies (23, 24). Older age, especially in hist-GBMs, and male sex, particularly in mol-GBMs, were found to be independently associated with patient survival. Preoperative motor dysfunction also served as a negative prognostic factor in mol-GBMs. Individuals with mol-GBMs were diagnosed 6 years earlier than those with hist-GBMs. Patients with mol-GBMs had a lower incidence of motor dysfunction and higher incidence of preoperative seizures than those with hist-GBMs, which might be due to the aggressive behavior of the latter, tending to compress and invade surrounding tissue rapidly, thus speeding up the onset of neurological symptoms while delaying the development of epilepsy (25, 26). It is well known that IDH mutation is more common in low-grade gliomas (LGGs) with a better prognosis than high-grade gliomas and that epilepsy is significantly correlated with LGGs (27). Although mol-GBMs are IDH-wildtype, they are histologically classified as relatively-low-grade. Certain molecular features, other than IDH mutation, may correlate with seizures in mol-GBMs; however, further investigation is needed. In our analysis, a higher number of patients with hist-GBMs than those with mol-GBMs received gross total resection. This might be due to the lack of angiogenic presentations and intratumoral necrosis, making it challenging to clearly recognize the tumor border of mol-GBMs during surgery. Another explanation might be that the postoperative radiological evaluation of the residual tumor/edema is difficult for mol-GBMs, since many of these tumors do not show radiologically contrast-enhancement, increasing the possibility of misinterpreting gross total resection as subtotal resection. For both subtypes, the greater the extent of tumor resection, the longer the OS. Although the Stupp treatment protocol was predictive of a better prognosis for hist-GBMs, our result showed it was not a significant protective factor for survival in mol-GBMs. The optimal treatment strategy for this new subtype of GBMs remains to be determined and requires further clinical trials.
The frontal and temporal lobes are the most commonly involved areas in GBM, they are the largest brain lobes (28). In the present study, tumor involvement of functional brain regions was found to be an independent predictor of a worse prognosis, while the motor cortex/tract was the most affected functional area, consistent with existing literature (29). Involvement of the insular lobe, thalamus, and callosum also predicted a shortened OS, especially in hist-GBMs. The pathological findings contradict the imaging findings of central nervous system necrosis. While we observed imaging enhancement and necrosis in the mol-GBMs, their pathological findings did not show necrosis. It has been noted that the imaging results did not perfectly match the pathology findings. The presentations of T1WI and T2WI of GBMs using the current classification shared many similarities with the previous classification (30). Previous studies have found that contrast enhancement in mol-GBMs is a marker of aggressiveness (17). In our study, we found that the rates of contrast enhancement and intratumoral necrosis were lower, and the maximal diameter of intratumoral necrosis was shorter in mol-GBMs, which might be due to their relatively low histological grade.
The identification of specific biomarkers plays a crucial role in determining the prognosis of a medical condition and selecting the most suitable therapeutic approach (31). The occurrence of IDH mutation is an early event in the development of glioma and has significant consequences for glioma progression and treatment response. The IDH gene has emerged as a promising molecular biomarker for predicting the response to TMZ treatment and has been validated as a standalone favorable prognostic indicator in patients with GBMs (32). Histone H3s mutations may drive pediatric GBM (33). Copy-number change in certain chromosomes, or chromosomal instability, is one of the hallmarks of GBMs (34). In our cohort, changes in chromosomes, especially chr-1, 7, 9, and 10, occurred in over 90% of GBMs. Besides, alterations of EGFR, TERT, CDK6, PTEN, BRAF, FGFR2, and CDKN2A/B were found in over 70% of GBMs. The current classification for GBM is mainly based on IDH, EGFR, TERT, and chromosome copy-number changes. Giulia et al. have underscored the significance of isolated TERT promoter mutation in glioma and its correlation with improved prognosis, amidst the plethora of molecules under investigation (35). As our understanding of molecular pathogenesis evolves, additional biomarkers may be added to improve these tumors’ classification, survival prediction, and targeted therapy.
The rates of alterations of PTEN, CDK4, MYB, TOP3A, KIT, KRAS, and NTRK3 were higher in hist-GBMs than in mol-GBMs. Previous studies have reported that PTEN alteration is associated with therapeutic resistance in GBMs (36). CDK4/6, which regulates the cell cycle, are an important factor for the tumorigenesis of GBMs (37). MYB is a transcription factor from the myeloblastosis family linked to the development and progression of GBMs (38). TOP3A levels are higher in GBMs compared to LGGs (39). KIT is elevated at the time of both the first diagnosis and the recurrence of GBMs (40). KRAS alteration is identified as an oncogene for GBMs, and its over-expression has a crucial role in glioma cell growth and proliferation (41). NTRK1-3 genes have been identified as potential driver mutation partners in GBMs (42). Approximately 12% of mol-GBM in this cohort showed no nuclear ATRX staining. Although this is not identical to our concept, ATRX expression was reported in other publications, showing that not all GBMs expressed ATRX (43). Previous studies reported that ATRX-deficient GBM cells show enhanced sensitivity to irradiation(44). The differences in molecular alterations between hist-GBMs and mol-GBMs revealed that this might be the key biological mechanism behind their clinical disparities.
Few studies assessed the real-world survival of GBM patients based on the WHO CNS5 classification (17, 45, 46). In this cohort, the median OS of GBMs was 12.6 months, which is shorter than 16 months calculated using the WHO CNS4 classification (47). The shorter OS may be partly due to the exclusion of IDH-mutant astrocytomas with relatively better survival and the introduction of mol-GBMs that are more infiltrative and less responsive to TMZ chemotherapy. In addition, whether mol-GBMs with a relatively low histological grade have the same OS as hist-GBMs with necrosis and/or microvascular proliferation are less known. In our cohort, the median OS of mol-GBMs was slightly higher than that of hist-GBMs, but the difference was not statistically significant. This finding is consistent with the reports by Ramos-Fresnedo (17) and Grogan (45) but not Ostrom (46), and might be explained by the small sample size in this study. However, since the OS curve of mol-GBMs with only 34 patients remained above the hist-GBMs curve, the non-significant OS difference might be false-negative if the sample size was increased.
The WHO CNS5 classification places greater emphasis on molecular markers than ever before. Beyond IDH mutation, EGFR amplification, TERT promotor mutation, and +7/-10 chromosome copy-number variations, it is important to find additional molecular changes associated with the prognosis of GBMs that could assist in distinguishing further subgroups for a more individualized treatment (22). In addition to MGMT promoter methylation being linked to positive outcomes, alterations of CDK4/6, CIC, FGFR3, MYB, TP53, and KMT5B were strongly related to a worse OS in mol-GBMs, hist-GBMs, or both, which is consistent with the previous studies (20, 48–50). These findings provided a valuable hint for further exploration of these molecular alterations as potential categorization markers or therapeutic targets in GBMs.
Subsequently, we attempted to utilize the TCGA, CGGA, and GEO datasets to confirm our findings. However, since the release of the WHO CNS5 classification was too recent, we could not find any available dataset of DNA sequencing for validation, so we finally used the RNA-sequencing dataset from the CGGA to access the associations between the expression levels of genes and the prognosis of GBMs. The results revealed that the expression levels of CIC, CDK4, CDK6, TP53, FGFR3, TOP3A, and MYB might impact the survival of patients with hist-GBMs, thus validating our results from the perspective of the expression of gene alterations.
The current investigation exhibits certain constraints. The screening of IDH-wildtype gliomas during a certain period was incomplete due to unsuccessful DNA extraction from some FFPE samples, which experienced DNA degradation over an extended period. Consequently, a selection bias was inevitable. Furthermore, it is noteworthy that the present study was conducted solely on a cohort of Asian patients, thereby constraining the generalizability of the findings to a broader international population. To enhance the validity of the findings of this study, it is recommended that future collaborations across multiple international centers be pursued, to broaden the scope of sample collection. Additionally, the outcomes were not verified utilizing DNA-sequencing information because of the inadequacy of present public datasets.
5 Conclusions
According to the WHO CNS5 classification, GBMs are a new comprehensive entity that encompasses both hist-GBMs (histologically grade 4) and mol-GBMs (histologically grade 2-3). Clinicians should update their diagnostic and therapeutic approaches to this lethal brain malignancy, as the clinical, radiological, molecular, and prognostic characteristics of GBMs have significantly changed after the release of the current WHO classification. Mol-GBMs, a newly identified GBM subtype, lacks sufficient research and clinical awareness, despite having distinct genetic features that may lead to different clinical manifestations. To explore additional targeted therapies, it is necessary to incorporate more molecular features into the current WHO classification.
Data availability statement
The data presented in the study are deposited in the figshare repository (DOI: 10.6084/m9.figshare.23593203). This data can be found here: https://figshare.com/articles/dataset/Data_for_the_article_Histological_and_molecular_glioblastoma_IDH-wildtype_A_real-world_landscape_using_the_2021_WHO_classification_of_central_nervous_system_tumors_published_in_Frontiers_in_Oncology/23593203.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Review Committee of Peking Union Medical College Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Conceptualization: XG and WC. Methodology: XG, LG, YL, and ZZ. Formal analysis and investigation: XG, LG, YL, ZZ, WC, YNW, YKW, HX, YS, DL, YX, TY, JL, JW, TL, HW, XZ, QL, SJ, TQ, SG, HL, YG, CW, KZ, YW, and WM. Writing—original draft preparation: XG, LG, YL, and ZZ. Writing—review and editing: YW and WM. Supervision: YW and WM. All authors provided meaningful contributions to the manuscript. All authors read and approved the final manuscript.
Funding
This work was funded by the National Natural Science Foundation of China (82203527), the Fundamental Research Funds for the Central Universities (3332022009), and the China Postdoctoral Science Foundation (2022M710451) for XG, and by the Beijing Municipal Natural Science Foundation (7202150) and the National High Level Hospital Clinical Research Funding (2022-PUMCH-A-019) for YW, and by the National High Level Hospital Clinical Research Funding (2022-PUMCH-B-113), the Tsinghua University-Peking Union Medical College Hospital Initiative Scientific Research Program (2019ZLH101) and the Beijing Municipal Natural Science Foundation (19JCZDJC64200[Z]) for WM.
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.
The reviewer XF declared a shared affiliation with the author ZZ to the handling editor at the time of review.
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/fonc.2023.1200815/full#supplementary-material
References
1. Ohgaki H, Kleihues P. The definition of primary and secondary glioblastoma. Clin Cancer Res (2013) 19(4):764–72. doi: 10.1158/1078-0432.CCR-12-3002
2. Omuro A, DeAngelis LM. Glioblastoma and other malignant gliomas: a clinical review. JAMA (2013) 310(17):1842–50. doi: 10.1001/jama.2013.280319
3. Chen F, Wendl MC, Wyczalkowski MA, Bailey MH, Li Y, Ding L. Moving pan-cancer studies from basic research toward the clinic. Nat Cancer (2021) 2(9):879–90. doi: 10.1038/s43018-021-00250-4
4. Bagley SJ, Kothari S, Rahman R, Lee EQ, Dunn GP, Galanis E, et al. Glioblastoma clinical trials: current landscape and opportunities for improvement. Clin Cancer Res (2022) 28(4):594–602. doi: 10.1158/1078-0432.CCR-21-2750
5. van Solinge TS, Nieland L, Chiocca EA, Broekman MLD. Advances in local therapy for glioblastoma - taking the fight to the tumour. Nat Rev Neurol (2022) 18(4):221–36. doi: 10.1038/s41582-022-00621-0
6. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 world health organization classification of tumors of the central nervous system: a summary. Acta Neuropathol (2016) 131(6):803–20. doi: 10.1007/s00401-016-1545-1
7. Eckel-Passow JE, Lachance DH, Molinaro AM, Walsh KM, Decker PA, Sicotte H, et al. Glioma groups based on 1p/19q, IDH, and TERT promoter mutations in tumors. N Engl J Med (2015) 372(26):2499–508. doi: 10.1056/NEJMoa1407279
8. Akyerli CB, Yuksel S, Can O, Erson-Omay EZ, Oktay Y, Cosgun E, et al. Use of telomerase promoter mutations to mark specific molecular subsets with reciprocal clinical behavior in IDH mutant and IDH wild-type diffuse gliomas. J Neurosurg (2018) 128(4):1102–14. doi: 10.3171/2016.11.JNS16973
9. Horbinski C. What do we know about IDH1/2 mutations so far, and how do we use it? Acta Neuropathol (2013) 125(5):621–36. doi: 10.1007/s00401-013-1106-9
10. Cancer Genome Atlas Research Network, Brat DJ, Verhaak RG, Aldape KD, Yung WK, Salama SR, et al. Comprehensive, integrative genomic analysis of diffuse lower-grade gliomas. N Engl J Med (2015) 372(26):2481–98. doi: 10.1056/NEJMoa1402121
11. Brat DJ, Aldape K, Colman H, Holland EC, Louis DN, Jenkins RB, et al. cIMPACT-NOW update 3: recommended diagnostic criteria for “Diffuse astrocytic glioma, IDH-wildtype, with molecular features of glioblastoma, WHO grade IV”. Acta Neuropathol (2018) 136(5):805–10. doi: 10.1007/s00401-018-1913-0
12. Brat DJ, Aldape K, Colman H, Figrarella-Branger D, Fuller GN, Giannini C, et al. cIMPACT-NOW update 5: recommended grading criteria and terminologies for IDH-mutant astrocytomas. Acta Neuropathol (2020) 139(3):603–8. doi: 10.1007/s00401-020-02127-9
13. Foote MB, Papadopoulos N, Diaz LA Jr.Genetic classification of gliomas: refining histopathology Cancer Cell (2015) 28(1):9–11 doi: 10.1016/j.ccell.2015.06.014
14. Louis DN, Perry A, Wesseling P, Brat DJ, Cree IA, Figarella-Branger D, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro Oncol (2021) 23(8):1231–51. doi: 10.1093/neuonc/noab106
15. Melhem JM, Detsky J, Lim-Fat MJ, Perry JR. Updates in IDH-wildtype glioblastoma. Neurotherapeutics. (2022) 19(6):1705–23. doi: 10.1007/s13311-022-01251-6
16. Horbinski C, Berger T, Packer RJ, Wen PY. Clinical implications of the 2021 edition of the WHO classification of central nervous system tumours. Nat Rev Neurol (2022) 18(9):515–29. doi: 10.1038/s41582-022-00679-w
17. Ramos-Fresnedo A, Pullen MW, Perez-Vega C, Domingo RA, Akinduro OO, Almeida JP, et al. The survival outcomes of molecular glioblastoma IDH-wildtype: a multicenter study. J Neurooncol (2022) 157(1):177–85. doi: 10.1007/s11060-022-03960-6
18. Hegi ME, Diserens AC, Gorlia T, Hamou MF, de Tribolet N, Weller M, et al. MGMT gene silencing and benefit from temozolomide in glioblastoma. N Engl J Med (2005) 352(10):997–1003. doi: 10.1056/NEJMoa043331
19. Benitez JA, Finlay D, Castanza A, Parisian AD, Ma J, Longobardi C, et al. PTEN deficiency leads to proteasome addiction: a novel vulnerability in glioblastoma. Neuro Oncol (2021) 23(7):1072–86. doi: 10.1093/neuonc/noab001
20. Lopez V, Tejedor JR, Carella A, Garcia MG, Santamarina-Ojeda P, Perez RF, et al. Epigenetic deregulation of the histone methyltransferase KMT5B contributes to malignant transformation in glioblastoma. Front Cell Dev Biol (2021) 9:671838. doi: 10.3389/fcell.2021.671838
21. Brito C, Azevedo A, Esteves S, Marques AR, Martins C, Costa I, et al. Clinical insights gained by refining the 2016 WHO classification of diffuse gliomas with: EGFR amplification, TERT mutations, PTEN deletion and MGMT methylation. BMC Cancer (2019) 19(1):968. doi: 10.1186/s12885-019-6177-0
22. Berger TR, Wen PY, Lang-Orsini M, Chukwueke UN. World health organization 2021 classification of central nervous system tumors and implications for therapy for adult-type gliomas: a review. JAMA Oncol (2022) 8(10):1493–501. doi: 10.1001/jamaoncol.2022.2844
23. Lee D, Riestenberg RA, Haskell-Mendoza A, Bloch O. Diffuse astrocytic glioma, IDH-wildtype, with molecular features of glioblastoma, WHO grade IV: a single-institution case series and review. J Neurooncol (2021) 152(1):89–98. doi: 10.1007/s11060-020-03677-4
24. McKinnon C, Nandhabalan M, Murray SA, Plaha P. Glioblastoma: clinical presentation, diagnosis, and management. BMJ (2021) 374:n1560. doi: 10.1136/bmj.n1560
25. Lee JW, Wen PY, Hurwitz S, Black P, Kesari S, Drappatz J, et al. Morphological characteristics of brain tumors causing seizures. Arch Neurol (2010) 67(3):336–42. doi: 10.1001/archneurol.2010.2
26. Chen DY, Chen CC, Crawford JR, Wang SG. Tumor-related epilepsy: epidemiology, pathogenesis and management. J Neurooncol (2018) 139(1):13–21. doi: 10.1007/s11060-018-2862-0
27. Duan WC, Wang L, Li K, Wang WW, Zhan YB, Zhang FJ, et al. IDH mutations but not TERTp mutations are associated with seizures in lower-grade gliomas. Medicine (2018) 97(50):e13675. doi: 10.1097/MD.0000000000013675
28. Babu R, Komisarow JM, Agarwal VJ, Rahimpour S, Iyer A, Britt D, et al. Glioblastoma in the elderly: the effect of aggressive and modern therapies on survival. J Neurosurg (2016) 124(4):998–1007. doi: 10.3171/2015.4.JNS142200
29. Roux A, Roca P, Edjlali M, Sato K, Zanello M, Dezamis E, et al. MRI Atlas of IDH wild-type supratentorial glioblastoma: probabilistic maps of phenotype, management, and outcomes. Radiology (2019) 293(3):633–43. doi: 10.1148/radiol.2019190491
30. Ceravolo I, Barchetti G, Biraschi F, Gerace C, Pampana E, Pingi A, et al. Early stage glioblastoma: retrospective multicentric analysis of clinical and radiological features. Radiol Med (2021) 126(11):1468–76. doi: 10.1007/s11547-021-01401-4
31. Sledzinska P, Bebyn MG, Furtak J, Kowalewski J, Lewandowska MA. Prognostic and predictive biomarkers in gliomas. Int J Mol Sci (2021) 22(19):10373. doi: 10.3390/ijms221910373
32. SongTao Q, Lei Y, Si G, YanQing D, HuiXia H, XueLin Z, et al. IDH mutations predict longer survival and response to temozolomide in secondary glioblastoma. Cancer Sci (2012) 103(2):269–73. doi: 10.1111/j.1349-7006.2011.02134.x
33. Korshunov A, Schrimpf D, Ryzhova M, Sturm D, Chavez L, Hovestadt V, et al. H3-/IDH-wild type pediatric glioblastoma is comprised of molecularly and prognostically distinct subtypes with associated oncogenic drivers. Acta Neuropathol (2017) 134(3):507–16. doi: 10.1007/s00401-017-1710-1
34. Cancer Genome Atlas Research Network. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature (2008) 455(7216):1061–8. doi: 10.1038/nature07385
35. Berzero G, Di Stefano AL, Ronchi S, Bielle F, Villa C, Guillerm E, et al. IDH-wildtype lower-grade diffuse gliomas: the importance of histological grade and molecular assessment for prognostic stratification. Neuro Oncol (2021) 23(6):955–66. doi: 10.1093/neuonc/noaa258
36. Benitez JA, Ma J, D’Antonio M, Boyer A, Camargo MF, Zanca C, et al. PTEN regulates glioblastoma oncogenesis through chromatin-associated complexes of DAXX and histone H3.3. Nat Commun (2017) 8:15223. doi: 10.1038/ncomms15223
37. Xu G, Li JY. CDK4, CDK6, cyclin D1, p16(INK4a) and EGFR expression in glioblastoma with a primitive neuronal component. J Neurooncol (2018) 136(3):445–52. doi: 10.1007/s11060-017-2674-7
38. Zhang X, Lv QL, Huang YT, Zhang LH, Zhou HH. Akt/FoxM1 signaling pathway-mediated upregulation of MYBL2 promotes progression of human glioma. J Exp Clin Cancer Res (2017) 36(1):105. doi: 10.1186/s13046-017-0573-6
39. Arivazhagan A, Kumar DM, Sagar V, Patric IR, Sridevi S, Thota B, et al. Higher topoisomerase 2 alpha gene transcript levels predict better prognosis in GBM patients receiving temozolomide chemotherapy: identification of temozolomide as a TOP2A inhibitor. J Neurooncol (2012) 107(2):289–97. doi: 10.1007/s11060-011-0758-3
40. Joensuu H, Puputti M, Sihto H, Tynninen O, Nupponen NN. Amplification of genes encoding KIT, PDGFRalpha and VEGFR2 receptor tyrosine kinases is frequent in glioblastoma multiforme. J Pathol (2005) 207(2):224–31. doi: 10.1002/path.1823
41. Li Y, Li Y, Ge P, Ma C. MiR-126 regulates the ERK pathway via targeting KRAS to inhibit the glioma cell proliferation and invasion. Mol Neurobiol (2017) 54(1):137–45. doi: 10.1007/s12035-015-9654-8
42. Zhou Z, He H, Wang K, Shi X, Wang Y, Su Y, et al. Granzyme a from cytotoxic lymphocytes cleaves GSDMB to trigger pyroptosis in target cells. Science (2020) 368(6494):eaaz7548. doi: 10.1126/science.aaz7548
43. Qin T, Mullan B, Ravindran R, Messinger D, Siada R, Cummings JR, et al. ATRX loss in glioma results in dysregulation of cell-cycle phase transition and ATM inhibitor radio-sensitization. Cell Rep (2022) 38(2):110216. doi: 10.1016/j.celrep.2021.110216
44. Koschmann C, Calinescu AA, Nunez FJ, Mackay A, Fazal-Salom J, Thomas D, et al. ATRX loss promotes tumor growth and impairs nonhomologous end joining DNA repair in glioma. Sci Transl Med (2016) 8(328):328ra328. doi: 10.1126/scitranslmed.aac8228
45. Grogan D, Bray DP, Cosgrove M, Boucher A, Erwood A, Linder DF, et al. Clinical and radiographic characteristics of diffuse astrocytic glioma, IDH-wildtype, with molecular features of glioblastoma: a single institution review. J Neurooncol (2022) 157(1):187–95. doi: 10.1007/s11060-022-03961-5
46. Ostrom QT, Shoaf ML, Cioffi G, Waite K, Kruchko C, Wen PY, et al. National-level overall survival patterns for molecularly-defined diffuse glioma types in the united states. Neuro Oncol (2022) 25(4):799–807. doi: 10.1093/neuonc/noac198
47. Stupp R, Taillibert S, Kanner A, Read W, Steinberg D, Lhermitte B, et al. Effect of tumor-treating fields plus maintenance temozolomide vs maintenance temozolomide alone on survival in patients with glioblastoma: a randomized clinical trial. JAMA (2017) 318(23):2306–16. doi: 10.1001/jama.2017.18718
48. Di Stefano AL, Fucci A, Frattini V, Labussiere M, Mokhtari K, Zoppoli P, et al. Detection, characterization, and inhibition of FGFR-TACC fusions in IDH wild-type glioma. Clin Cancer Res (2015) 21(14):3307–17. doi: 10.1158/1078-0432.CCR-14-2199
49. Bunda S, Heir P, Metcalf J, Li ASC, Agnihotri S, Pusch S, et al. CIC protein instability contributes to tumorigenesis in glioblastoma. Nat Commun (2019) 10(1):661. doi: 10.1038/s41467-018-08087-9
Keywords: glioblastoma, histological glioblastoma, molecular glioblastoma, molecular alteration, chromosome copy number variation
Citation: Guo X, Gu L, Li Y, Zheng Z, Chen W, Wang Y, Wang Y, Xing H, Shi Y, Liu D, Yang T, Xia Y, Li J, Wu J, Zhang K, Liang T, Wang H, Liu Q, Jin S, Qu T, Guo S, Li H, Wang Y and Ma W (2023) Histological and molecular glioblastoma, IDH-wildtype: a real-world landscape using the 2021 WHO classification of central nervous system tumors. Front. Oncol. 13:1200815. doi: 10.3389/fonc.2023.1200815
Received: 05 April 2023; Accepted: 19 June 2023;
Published: 06 July 2023.
Edited by:
Manuela Caroli, Fondazione Politecnico di Milano, ItalyReviewed by:
Xing Fan, Capital Medical University, ChinaDavid Meredith, Brigham and Women’s Hospital and Harvard Medical School, United States
Vani Santosh, National Institute of Mental Health and Neurosciences (NIMHANS), India
Copyright © 2023 Guo, Gu, Li, Zheng, Chen, Wang, Wang, Xing, Shi, Liu, Yang, Xia, Li, Wu, Zhang, Liang, Wang, Liu, Jin, Qu, Guo, Li, Wang and Ma. 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: Yu Wang, eXdhbmdAcHVtY2guY24=; Wenbin Ma, bWF3YjIwMDFAaG90bWFpbC5jb20=
†These authors have contributed equally to this work