Abstract
Background:
The purpose of this study is to utilize single-cell sequencing data to explore glioma heterogeneity and identify key biomarkers associated with glioblastoma multiforme (GBM) relapse using machine learning.
Methods:
Single-cell sequencing and transcriptome data for gliomas were obtained from the GEO (GSE159416, GSE159605, and GSE186057) and TCGA databases. A prognostic model based on differentiation-related genes (DRGs) was constructed using weighted correlation network analysis, univariate Cox regression, and LASSO analysis. Key genes were identified using LASSO and SVM-RFE, with intersecting genes selected as the final set of key genes. Further analyses examined immune infiltration patterns and functional pathways. Importantly, we analyzed the relationship between prognostic-related genes and ubiquitination, and further characterized the characteristics of ubiquitination-related prognostic genes. In addition, we performed CCK-8 assays, colony formation, Transwell invasion assays, apoptosis assays to determine the role of ETV4 in glioma.
Results:
Examination of single-cell RNA-seq data from the GEO database revealed three distinct cell differentiation stages in glioma tissues. Marker genes for each of these cell states were combined to form DRGs. A 16-gene DRG signature was developed for predicting the survival of glioma patients. Machine learning identified four important genes with high AUCs in both training and test sets. Notably, 13 out of 16 genes in the DRG signature are ubiquitin-related, highlighting the involvement of ubiquitination in GBM. Moreover, we reported that inhibition of ETV4 attenuates cell proliferation and invasion in glioma cells.
Conclusion:
Our prognostic model, based on the differentiation-related gene signatures, may be valuable for predicting prognosis and immunotherapy response in glioma patients. Characterizing these ubiquitination-associated features may elucidate the molecular mechanisms driving GBM progression and offer novel insights for its diagnosis and treatment. Additionally, machine learning identified four biomarkers with potential for aiding in the diagnosis and treatment of GBM.
1 Introduction
Glioblastoma Multiforme (GBM) is the most common and aggressive primary brain tumor, characterized by invasive growth and resistance to conventional treatments (1, 2). The five-year survival rate for patients with GBM is only 7%, highlighting the poor prognosis associated with this disease (3). The World Health Organization (WHO) classifies gliomas into four grades, with high-grade gliomas (WHO grade 3 to 4) often exhibiting poor prognoses (4). For patients with GBM (WHO grade 4), median survival ranges from 14.6 to 17 months (5).
The complexity of GBM treatment is complicated by tumor heterogeneity, arising from the diverse types of GBM cells and a complex tumor microenvironment (6–8). This heterogeneity exists not only between individuals but also within individual tumors (9, 10). Single-cell sequencing (SCS) has been reported to explore the mechanisms of disease development (11, 12). SCS allows researchers to examine interactions between various types of GBM cells and neoplastic cells at a more detailed level (13). Unlike single-cell sequencing, next-generation sequencing (NGS), which is commonly used, analyzes the entire cell population and is unable to capture cellular heterogeneity (14). SCS amplifies and sequences the genome or transcriptome at the single-cell level, providing information on single nucleotide variations (SNVs), gene copy number variations (CNVs), single-cell genome structure variations, gene expression, gene fusions, alternative splicing in the single-cell transcriptome, and DNA methylation in the single-cell epigenome (15). SCS enables the study of genetic characteristics in diseases and biological processes at the single-cell level, including early embryonic development, tumorigenesis mechanism, tumor heterogeneity and evolution, as well as circulating tumor cells (CTCs) and clonal evolution (16). Meanwhile machine learning brings new opportunities for identifying biomarkers (17). By efficiently filtering out irrelevant features, machine learning is highly suitable as a tool for pre-screening features (18, 19).
Protein post-translational modifications (PTMs) are covalent and enzymatic alterations that occur during or after biosynthesis, modulating protein properties and functions. Among them, non-histone PTMs, including acetylation, lactylation, methylation, ubiquitination, phosphorylation, and SUMOylation, have been reported to be closely associated with cancer progression (20). Abnormalities in PTMs have been observed to influence cancer cell proliferation, migration, and invasion. Ubiquitination, a crucial PTM, governs protein stability and a wide range of cellular processes by attachment of ubiquitin molecules (21). Emerging evidence suggests that ubiquitination plays a pivotal role in mediating resistance to cancer immunotherapy (22). For instance, the ubiquitin-conjugating enzyme E2S has been reported to reduce the sensitivity of GBM cells to temozolomide by upregulating PGAM1 through interaction with OTUB2 (23). In addition, centromere protein U facilitates temozolomide resistance by mediating the ubiquitination and degradation of RPS3 in GBM (24). RNF8-mediated ubiquitination of KRT80 has been shown to drive glucose metabolic reprogramming and GBM progression (25). Moreover, the tryptophan-metabolizing enzyme IL4I1 inhibits ferroptosis in GBM by decreasing the ubiquitination of Nrf2 via I3P (26). Therefore, ubiquitination is critically involved in GBM development and progression.
This study aims to analyze glioma single-cell RNA sequencing (scRNA-seq) data to explore heterogeneity of GBMs and identify differentiation-related genes (DRGs) for prognostic prediction using bulk RNA-seq data. Additionally, support vector machines recursive feature elimination (SVM-RFE) and least absolute shrinkage and selection operator (LASSO) machine learning algorithms were employed to select key biomarkers associated with GBMs and enhance predictive accuracy for GBM prognosis. Importantly, this study also investigated transcription factors, protein-protein interaction networks, enriched pathways and biological functions, differential expression patterns, and survival associations of ubiquitin-related prognostic genes.
2 Methods
2.1 Data collection
Single-cell sequencing data were obtained from the GEO database (Gene Expression Omnbius) under accession GSE159416 (27) (https://www.ncbi.nlm.nih.gov/geo/), which includes single-cell RNA-seq data from 18 glioma patients. Additionally, clinical data and gene expression data were sourced from the TCGA (The Cancer Genome Atlas, https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga). Array-based expression profiles were collected from GEO datasets GSE159605 and GSE186057. A total of 1041 samples were included in this study.
2.2 Quality control for SCS data processing
ScRNA-seq data were preprocessed using “Seurat” and “Monocle” R packages. The PercentageFeatureSet function was used to calculate the number of mitochondrial genes, and cells with less than 500 or more than 5000 were removed based on quality control results. The LogNormalize method was used to normalize scRNA-seq data, and the ‘vst’ selection method was used to identify the top 1000 highly variable genes. Principal component analysis (PCA) was used to reduce dimensions for glioma cells. Using the t-distributed stochastic neighbor embedding (tSNE) approach, the top 10 principal components (PCs) with significant values were chosen for clustering. The ‘limma’ package was used to find marker genes in each cluster, selecting those with an adjusted p value of less than 0.05 and a |Log2 fold change (FC)| > 2 value. ScRNA-seq data were automatically annotated using the “SingleR” tool. Data from the “celldex” package’s major human cell atlas were utilized as reference data.
2.3 Pathway and pseudotime analysis
For astrocyte and tissue stem cells, pseudotime and trajectory analyses were carried out using the “Monocle” package, with the distribution of cells along each branch representing a single differentiation condition. |log2 (FC)|> 2 and adjusted p- values < 0.05 were used to identify DEGs in cells that were in separate development states. GO (Gene ontology) analysis, which categorizes gene functions, was applied to differentially expressed transcripts to provide a structural functional description. In the transcriptome project, GO functional analysis gives a categorical annotation of GO function for differentially expressed transcripts. ClusterProfiler, org.Hs.eg.db, enrichplot, and ggplot2 packages were used to analyze the KEGG (Kyoto Encyclopedia of Genes and Genomes, https://www.kegg.jp/kegg/kegg1.html) pathways and GO for DRGs.
2.4 Consensus clustering and prognosis analysis
Consensus clustering, an unsupervised clustering method, was used to distinguish samples into several subtypes according to different omics datasets. Data were classified using the “ConsensusClusterPlus” R package. Kaplan-Meier analysis was conducted to compare overall survival across different clusters, and a prognostic modle was developed using Cox regression analysis.
2.5 Ubiquitination−related prognosis genes
A list of Ubiquitination-related genes (URGs) was obtained by querying the human gene database GeneCards (https://www.genecards.org/) using the keyword “Ubiquitination” (28). The “venn” R package was then utilized to identify the intersection between URGs list and glioma-associated prognostic genes (GAPGs). Transcription factor analysis was conducted using the transcriptional regulatory relationships unraveled by sentence-based text mining (TRRUST, version 2, https://www.grnpedia.org/trrust/Network_search_form.php), which employs sentence-based text mining to uncover transcriptional regulatory relationships (29). Protein-protein interaction (PPI) networks were constructed using the STRING database (Version 12.0, https://cn.string-db.org/) (30). Functional annotation and enrichment analyses of gene lists were performed using Metascape (https://metascape.org/gp/index.html#/main/step1) (31). Additionally, differential gene expression and survival analyses were carried out using the GEPIA database (http://gepia.cancer-pku.cn/), which integrates tumor and normal tissue data from TCGA and GTEx (32).
2.5 Machine learning
Distinctive genes were screened using SVM-RFE (33) and LASSO (34). Key genes were individually filtered by LASSO and SVM-RFE, and the intersecting genes were then chosen as the final key genes. LASSO reduced model complexity by selectively incorporating variables, optimizing performance while controlling complexity to avoid overfitting. SVM-RFE provides high-accuracy feature screening and uses incremental regularization to prevent overfitting. Both algorithms were employed to screen for prognostic genes, and the intersecting genes from both methods were chosen as the final prognostic markers.
2.6 WGCNA
Weighted correlation network analysis (WGCNA) is a systems biology approach that identifies highly correlated gene sets and potential biomarkers by analyzing the association between gene modules and phenotype genes or therapeutic targets. WGCNA classifies co-expressed genes into modules, facilitating investigation of module-phenotype relationships. In this study, WGCNA was used to identify genes linked to the LLPS phenotype. The pickSoftThreshold function from the “WGCNA” R package was used to determine the optimal soft field value.
2.7 Immune inflation and immune checkpoint analysis
The “cibersort” package (35) was used to quantify related infiltration and activity levels for 22 immune cell types based on published gene signature lists across all tumors and normal samples. Immune cells types in this study included activated dendritic cells (DCs), neutrophils, mast cells, eosinophils, macrophages, and components of adaptive immunity such as B cells, T cells, central memory T cells. A heatmap was created to represent the findings of the investigation on the immunological checkpoint differences between high- and low- risk groups.
2.7 Cell culture and transfection
The human glioma cell lines LN229 and U251 were obtained from the Shanghai Cell Bank (Shanghai, China). Cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1% penicillin-streptomycin. Cells were maintained in a humidified incubator at 37°C with 5% CO2.
Short hairpin RNA (shRNA) targeting ETV4 gene and shRNA control (shNC) were purchased from GenePharma (Shanghai, China) company and transfected into LN229 and U251 cells using Lipofectamine 3000 (Invitrogen) according to the manufacturer’s protocol. shETV4-1: GCT GGA TGA CCC AAC AAA T; shETV4-2: CCC TGT GTA CAT ATA AAT GAA. Knockdown efficiency was validated by Western blotting.
2.8 Western blotting analysis
Cells were lysed in RIPA buffer (Beyotime, China) and protein concentrations were determined using the BCA Protein Assay Kit. Protein were separated by SDS-PAGE and transferred onto PVDF membranes. Membranes were blocked with 5% non-fat milk for 1 hour at room temperature, then incubated with anti-ETV4 antibody (CST #65763) overnight at 4°C. Protein bands were visualized using an ECL detection kit as described previously.
2.9 CCK-8 assays and colony formation
Cell viability was assessed using the Cell Counting Kit-8 (CCK-8) according to the manufacturer’s instructions. Briefly, transfected cells were seeded into 96-well plates. After the specified durations, 10 µL of CCK-8 reagent was added to each well. After cells were incubated at 37°C for 2 hours. Absorbance was measured at 450 nm by a microplate reader. Transfected cells were seeded into 6-well plates and cultured in complete medium for 14 days, allowing colonies to form. Then, cells were fixed with 4% paraformaldehyde for 15 minutes and stained with 0.1% crystal violet for 30 minutes at room temperature. After washing with PBS, colonies were imaged and counted.
2.10 Cell apoptosis assays
Cell apoptosis was evaluated using an Annexin V-FITC/PI Apoptosis Detection Kit according to the manufacturer’s instructions. Briefly, cells were harvested, washed, and resuspended in binding buffer. Then, 100 μL of the cell suspension was incubated with 5 μL Annexin V-FITC and 5 μL propidium iodide (PI) for 15 minutes at room temperature in the dark. After incubation, samples were analyzed using a flow cytometer. Cells were classified as viable, early apoptotic, late apoptotic, and necrotic.
2.11 Cell invasion assays
Cell invasion was assessed using Transwell chambers (Corning, USA) coated with Matrigel. Briefly, Cells were suspended in 200 μL serum-free medium and seeded into the upper chamber. The lower chamber was filled with medium containing 10% FBS. After incubation at 37°C for 24 hours, non-invading cells on the upper surface were removed with a cotton swab, and the invading cells on the lower surface were fixed with 4% paraformaldehyde and stained with Calcein AM. Invaded cells were imaged and counted under a microscope.
2.12 Statistical analysis
All statistical analyses were calculated using R software version 4.1.3 and the corresponding packages. Data were analyzed using student t-tests and one-way ANOVA for in vitro experiments. p < 0.05 was considered statistically significant.
3 Results
3.1 Single-cell sequencing analysis
The ranges of single-cell RNA counts, as well as RNA counts per cell, demonstrated a high level of sample quality control (Figure 1A). The 2000 most variable genes and the top 10 genes identified across all samples are displayed (Figure 1B). The number of identified genes and sequencing death showed a strong positive association (R = 0.92, Figure 1C). Next, we performed PCA on the normalized cell data (Figure 1D). Scatterplots were generated to depict the relationship between PC scores and selected genes, where larger absolute score values indicate stronger correlations. For each PC, the top 20 genes were chosen based on their coefficients. In PC_1, genes SOX2 and GPM6B exhibited significant correlations with score values (Figure 1E).
Figure 1
A gene expression heatmap was generated to visualize the expression levels of the top 20 selected genes for each PC, where darker colors indicate higher expression levels (Figure 1F). To investigate the relationship and clustering patterns among all cells, the t-distributed stochastic neighbor embedding (t-SNE) algorithm was employed. This algorithm effectively reduces high-dimensional data to a two-dimensional space. This analysis identified 14 distinct clusters, with cells closer together in the t-SNE plot indicating greater similarity (Figure 1G). To validate the optimal clustering results and identify potential feature genes for each subgroup, the top 20 genes for each cellular subgroup were generated, including ANLN, RRM2, PDIA2, CPS1, LGALS2, TMED3, and AQP5.
TMED3 exhibited high expression across all subclusters, indicating a ubiquitous expression pattern within the analyzed dataset (Figure 1H). Pseudo-temporal analysis was used to analyze cellular motion trajectories across the 14 subclusters, inferring the temporal progression and cell movement dynamics within each subcluster (Figure 1I). Cell type annotation was performed on cells distributed into 14 clusters based on marker gene expression. Additionally, pseudo-temporal analysis on the annotated cells inferred their temporal order or developmental trajectory (Figure 1J). A total of 17,846 astrocytes and stem cells underwent pseudotime and differentiation trajectory analysis (Figures 2A, B), revealing 3 branches of cells with diverse differentiation patterns. Cells in state 1 were inferred as the initial cell type, subsequently differentiating into various states (Figures 2C, D).
Figure 2
3.2 Consensus clustering
Consensus clustering was performed on cell populations to classify cells. The results suggested that the optimal clustering occurred at K = 2. The delta area plot indicated the best K value by identifying the point with the least increase in stability. Figure 2 showed the sample clusters at each K value with a heatmap, which qualitatively assesses unstable clusters and samples (Figure 2E).
3.3 KEGG and GO analysis
KEGG pathway analysis indicated that GBM development is associated with pathways including Coronavirus disease (COVID-19), ribosome, human papillomavirus infection, focal adhesion, proteoglycans in cancer, shear stress and atherosclerosism, ECM-receptor interaction, relaxin signaling pathway, amoebiasis and proteasome. GO analysis suggested that the molecular functions of the genes may relate to cytoplasmic translation, axonogenesis, axon development, glial cell differentiation and gliogenesis (Figure 3).
Figure 3
3.4 Immune infiltration and immune checkpoint analysis
Prognostic analysis revealed no significant differences among clusters 1, 2, and 3 (p = 0.833) (Figure 4A). However, immune checkpoints PDC1LG2, CD274, JAK2, HAVCR2, CD86, ICOSLG, YTHDF1, CD40, PVR, and TNFSF9 showed significant differential expression across clusters (Figure 4B). Glioma prognosis was correlated to immune checkpoints LGALS9, PVR, TNFSF9 and ICOSLG (Figure 4C). Cluster C2 samples had the highest immune, stromal, and ESTIMATE scores, but the lowest tumor purity (p < 0.001), showing the highest quantity of immune and stromal cells. In contrast, C1 had the highest tumor purity and the least amount of immune and stromal cells (Figure 4D). Moreover, WGCNA clustering showed co-expression patterns and functional modules (Figure 4E).
Figure 4
3.5 Weighted gene co-expression network analysis
WGCNA was employed to identify co-expressed gene modules and explore association between the gene network, phenotype, and core genes (36). Analysis was feasible with the soft domain value set to 6 (Figure 5A). WGCNA organized genes into modules, ultimately categorizing glioma gene expression into five distinct modules (Figure 5B). Mean connectivity remained constant when the soft threshold increased. All genes were grouped into 15 modules (Figure 5C), and the ME yellow module showed the strongest correlation with the survival time (futime) (Figure 5D, P < 0.05). Differential expression analysis on the key genes from the key modules identified two downregulated genes, while the remaining genes showed no differential expression (Figure 5E). The heatmap of differential expression is displayed in Figure 5E.
Figure 5
3.6 Prognosis model of GBM
A prognostic model for glioma was constructed using multiple Cox regression analysis of 532 genes screened from WGCNA. The LASSO algorithm was used to prevent overfitting. Univariate Cox regression analysis narrowed down 39 genes, followed by multiple distant linear regression (Figure 5F), resulting in 16 glioma-associated prognostic genes for the model (Table 1). Data were randomly split into training and test datasets and further divided into high-risk and low-risk groups based on the median risk score. In the training dataset, low-risk group patients had a significantly better prognosis than those in the high-risk group. However, in the test group, no significant difference in survival prognosis was observed between the two groups. The AUC values for the training dataset were 0.853 at one year, 0.930 at three years, and 0.968 at five years, which implies an improvement in predictive ability over time. For the test dataset, AUC values were 0.566 at one year, 0.591 at three years, and 0.586 at five years (Figure 5G). The 12 genes in the model were put into the single-cell data to assess their enrichment. Results demonstrated differential expression of the 12 genes across enriched clusters. Notably, NBP20 was primarily expressed in clusters 5, 12 and 13. RPL13 and CDK5 were significantly expressed in clusters 7 and 6, respectively. Additionally, the gene expression of LDHA, ALDOA, VKORC1, TUBB6, and PLAUR was enhanced across clusters (Figure 6).
Table 1
| Gene | HR | Z | P value | Lower | Upper |
|---|---|---|---|---|---|
| NBPF20 | 0.295241339188521 | 1.34345054682984 | 0.948486233498846 | 1.90288409892836 | 0.0964716449193646 |
| CDK5RAP2 | -0.398591676511448 | 0.671264738561516 | 0.514408353931362 | 0.875950683522888 | 0.00333199439858064 |
| RHEB | 0.264384784067291 | 1.30262933089252 | 1.08445898377724 | 1.56469096488212 | 0.00469979295937277 |
| RPL13AP5 | -0.131214541821952 | 0.877029594678408 | 0.74572031523489 | 1.03146031324022 | 0.112820108845984 |
| CAMK2N1 | 0.196446090232239 | 1.21706970716434 | 1.06321647782965 | 1.3931863387979 | 0.00438646096147386 |
| CPE | 0.16664125524179 | 1.18133039319568 | 0.998408106815502 | 1.39776659300079 | 0.0522093735949419 |
| NNAT | 0.16628013446608 | 1.18090386726605 | 1.06192535913427 | 1.31321277124486 | 0.00214871222075014 |
| LDHA | -0.331448737693059 | 0.717882955609934 | 0.597480014991583 | 0.862549248551042 | 0.000402327844752551 |
| ALDOA | 0.288032331830181 | 1.33380042763517 | 1.11458506722804 | 1.59613082309112 | 0.00166568621885795 |
| VKORC1 | 0.358234035314785 | 1.43080043912984 | 1.15793313183257 | 1.76796901335245 | 0.000905945183210735 |
| TUBB6 | 0.165766785316175 | 1.18029780684328 | 1.00330350190752 | 1.38851594775701 | 0.0455272061046904 |
| PLAUR | 0.214402745243588 | 1.23912160464437 | 1.01796544943357 | 1.50832462138159 | 0.0325589704722943 |
| RGS16 | 0.177090760257149 | 1.19373943230846 | 1.04858203072293 | 1.35899127631024 | 0.00742620195389444 |
| STC1 | 0.096837195341802 | 1.10168100022131 | 0.981281348695747 | 1.23685325096803 | 0.101014379923527 |
| MT2A | -0.0824804353428026 | 0.920829453393163 | 0.834120439835026 | 1.01655209696582 | 0.102129572069046 |
| CST3 | 0.199474018574925 | 1.22076049192171 | 1.05553281189418 | 1.41185206356837 | 0.0071806675559779 |
16 glioma-associated prognostic genes.
Figure 6
3.7 Ubiquitination−related prognostic genes in GBM
A total of 16 glioma-associated prognostic genes and 16,384 ubiquitination-related genes (URGs) were obtained from the GeneCards database. The intersection of these two gene sets identified 13 ubiquitination-related glioma-associated prognostic genes (UR-GAPGs): LDHA, TUBB6, ALDOA, RHEB, MT2A, CST3, CDK5RAP2, VKORC1, NNAT, CAMK2N1, PLAUR, CPE, and STC1 (Figure 7A). Transcription factor analysis revealed that three of these genes—ALDOA, LDHA, and PLAUR—are co-regulated by transcription factors JUN, HIF1A, SP1, and ATF1 (Figure 7B). For the PPI network, a minimum interaction confidence score of 0.15 (low confidence) was set, and disconnected nodes were excluded, resulting in a network of 10 genes (Figure 7C).
Figure 7
3.8 Functional characteristics of ubiquitination-related prognostic genes
Using the Metascape web portal, enrichment analysis of the 13 UR-GAPGs revealed that these genes are involved in pathways such as clear cell renal carcinoma, cellular homeostasis, and brain development (Figure 7D). Gene Ontology (GO) biological process analysis highlighted major roles in homeostatic and developmental processes (Figure 7E). Furthermore, a COVID-related enrichment analysis based on Blanco-Melo A549-ACE2-ruxolitinib RNA-seq data was conducted (Figure 8A). The DisGeNET database confirmed associations between these 13 genes and inflammatory processes, emphasizing the relevance of inflammation in their functional roles (Figure 8B). Analysis using the Cell Type Signatures database identified ZHONG PFC C1 OPC as the cell type most enriched for these genes (Figure 8C). Among transcription factor targets, CEBPB 01 emerged as the most significantly enriched (Figure 8D).
Figure 8
3.9 Differential expression and survival analysis of ubiquitination-related prognostic genes
Differential expression analysis using the GEPIA database, which integrates data from TCGA and GTEx, included 163 GBM tumor samples and 207 normal tissue samples. Box plots showed that LDHA, TUBB6, RHEB, MT2A, VKORC1, NNAT, CAMK2N1, PLAUR, and STC1 were significantly differentially expressed between tumor and normal tissues (P < 0.05) (Figures 9A–I). For survival analysis, patients were divided into high- and low-risk groups based on the average expression levels of the 13 UR-GAPGs. Kaplan–Meier analysis demonstrated that NNAT, PLAUR, and STC1 were significantly associated with overall survival (P < 0.05), with lower expression levels corresponding to improved survival outcomes (Figures 9J–L).
Figure 9
3.10 Machine learning
The study included two GEO datasets, GSE159605 and GSE186057, which were combined for differential analysis. The six most significantly different genes were visualized in a volcano plot and heatmap (Figures 10A, B). SVM-RFE, effective when there are more predictors than observations (33, 37), and LASSO were used to screen genes. The results of LASSO and SVM-RFE were displayed in Figures 7C, D. LASSO identified four key genes: ETV4, CACNA2D3, HIST1H3B, and HSPA1A (Figure 10C), while SVM-RFE identified seven genes: ETV4, HSPA1A, CACNA2D3, HIST1H3B, PTPRR, BDNF, and HIST1H3G (Figure 10D). The intersecting genes include ETV4, CACNA2D3, HIST1H3B and HSPA1A (Figure 10E).The combined data were divided into a training set and a test set to evaluate the sensitivity and specificity of selected biomarkers. In the training set, the AUCs of HSPA1A, HIST1H3B, ET4V and CACNA2D3 were 0.852, 0.861, 0.904 and 0.841, respectively (Figure 10F). Testing the biomarkers on the test set yielded AUCs of 1, 1, 1 and 0.900, indicating high accuracy (Figure 11A). The immune infiltration patterns for 22 cell types in the control group and treatment group are shown in Figure 11B.
Figure 10
Figure 11
3.11 Inhibition of ETV4 attenuates cell proliferation and invasion
Evidence revealed that HSPA1A promotes tumor cell proliferation and invasion, contributing to the progression and recurrence of GBM (38). CACNA2D3 suppressed cell proliferation, migration and invasion in glioma (39). The prognostic impact of HIST1H3B/C mutations in diffuse midline gliomas varies depending on patient age (40). However, the role of ETV4 in glioma development and progression remains largely uncharacterized. To explore whether ETV4 regulates cell viability in glioma cells, shETV4 transfection was conducted in LN229 and U251 cells. Western blotting results demonstrated downregulation of ETV4 in shETV4-transfected cells (Figure 12A). CCK-8 assays were performed to evaluate cell viability in glioma cells after shETV4 transfection. Cell viability was decreased in shETV4-transfected cells compared with control group (Figure 12B). Moreover, colony formation results revealed that shETV4 transfection decreased colony formation ability in glioma cells (Figures 12C, D). Additionally, cell apoptosis was increased in shETV4-transfected glioma cells (Figures 13A, B). Strikingly, Transwell invasion assays showed that shETV4 transfection reduced cell invasive ability in glioma cells (Figures 13C, D). Hence, ETV4 downregulation inhibits cell viability, proliferation and invasion in glioma cells.
Figure 12
Figure 13
4 Discussion
The advent of single-cell sequencing has enabled high-throughput analyses of the genome, transcriptome, and epigenome at the single-cell level, exposing gene structure and expression profiles within individual cells and reflecting cellular heterogeneity (41–43). This approach has allowed researchers to explore putative monocyte oncogenes in gliomas, assess monocyte functional states across various phases, and examine intercellular communication channels using published single-cell sequencing data. Chen et al. employed a machine learning-based approach and identified UPP1 as a critical oncogene involved in tumorigenesis and immune evasion in gliomas (44). Similarly, Hu et al. integrated scRNA-seq with machine learning techniques to uncover the role of ALPK1 in shaping tumor immune heterogeneity and regulating the TGF-β signaling pathway in glioma (45). Integrating scRNA-seq with machine learning, Yang et al. revealed the heterogeneity of GBM-associated neutrophils and developed a prognostic model based on VEGFA-expressing neutrophils (46). In this study, single-cell sequencing data of GBM (GSE159416) was used to analyze three distinct cell differentiation states in glioma tissues, with most of the 14 enriched clusters being associated with fibroblasts. Marker genes for each of these cell states were combined to form DRGs. To predict the prognosis of glioma patients, a prognostic model based on 16 DRGs was developed.
The AUC values of our model were higher than those of other glioma prognostic models, suggesting greater accuracy in predicting the prognosis of glioma patients. By means of LASSO and SVM-RFE machine learning, four disease-associated genes, including ETV4, CACNA2D3, HSPA1A and HIST1H3B, were identified. These genes demonstrated high sensitivity and specificity in both training and test groups, with AUC values exceeding 0.9. ETV4 is a member of the polyoma enhancer activator 3 (PEA3) family, which plays an important role in cell growth, invasion and metastasis (47). During embryogenesis, ETV4 is widely expressed in various tissues, promoting the morphogenesis of epithelial-derived organs such as the kidneys (48), lungs (49), and breasts (50). However, in adults, ETV4 is rarely expressed in normal tissues, mainly appearing in tumor tissues such as those in breast, gastric, prostate, colon and ovarian cancers. Numerous studies have shown that high expression of ETV4 in prostate cancer (51), breast cancer (52), colorectal cancer (53), pancreatic cancer (54), and cholangiocarcinoma (55) often correlates with poor prognosis. Notably, one study showed that ETV4 expression increases with glioma grade progression (56). ERK kinase promoted phosphorylation of ETV4, leading to blockade of ETV4 ubiquitination and degradation in colorectal cancer (57). PTK6 induced phosphorylation of ETV4 and increased nuclear translocation of ETV4, leading to enhanced metastasis in bladder cancer (58). We found that inhibition of ETV4 attenuates cell proliferation and invasion in glioma cells.
CACNA2D3 is an auxiliary member of the α-2/δ subunit triple family with voltage-dependent calcium channel complexes and plays a key role in tumor suppression. CACNA2D3 is downregulated in gliomas and functions as a tumor suppressor (39). This downregulation in glioma cells and high-grade glioma tissues is associated with increased methylation (39). Additionally, a case report indicates that HIST1H3B K27 mutation is associated with glioma development of (59–61). Genomic analysis of HIST1H3B mutations may aid in timely glioma diagnosis, supporting surgical and clinical management of these patients (61). In gliomas, heat shock protein 70 (Hsp70, HSPA1A) is overexpressed in the cytoplasm (62). Overexpression of HSPA1A significantly enhances cell proliferation, with cellular immunofluorescence revealing its primary localization in the cytoplasm, where it promotes tumor cell proliferation (62).
For ubiquitin-related biomarkers, LDHA-mediated metabolic reprogramming has been shown to promote cardiomyocyte proliferation by reducing reactive oxygen species (ROS) and inducing M2 macrophage polarization (63). In clinical tissue samples, PTMs were quantitatively assessed, highlighting the pivotal role of ALDOA K330 ubiquitination/acetylation in tumor progression (64). RHEB ubiquitination plays a regulatory role in growth factor-induced activation of mTORC1 (65), and the Cullin3-Rbx1-KLHL9 E3 ubiquitin ligase complex has been shown to mediate RHEB ubiquitination, thereby facilitating amino acid-induced mTORC1 activation (66). Additionally, Malvidin has been found to reduce trauma-induced heterotopic ossification of tendons in rats by promoting RHEB degradation via the ubiquitin–proteasome pathway (67). In GSH-depleted RAW 264.7 cells, hydrogen peroxide triggers Beclin 1-independent autophagic cell death through suppression of the mTOR pathway via ubiquitination and degradation of RHEB (68). Moreover, STC1, a ubiquitin-related gene, is significantly upregulated in lung adenocarcinoma and is associated with poor prognosis, suggesting its potential as a biomarker for prognosis evaluation, tumor characterization, and therapeutic decision-making (69). However, despite these findings, it remains unclear whether these ubiquitin-related biomarkers play a role in glioblastoma (GBM). Their potential as novel therapeutic targets or biomarkers for precision medicine in GBM warrants further investigation.
5 Conclusions
In conclusion, our prognostic model based on differentiation-related gene signatures shows promise for predicting glioma prognosis and immunotherapy response. Furthermore, the characterization of ubiquitination-related features and machine learning–identified biomarkers provides novel insights into GBM progression, diagnosis, and treatment. Our study has several drawbacks. First, this research is based on bioinformatics analysis, and experimental validation is needed to further explore the underlying molecular mechanism. Our data is sourced from the TCGA and GEO databases. Therefore, validating the findings using independent clinical cohorts or third-party datasets is crucial to enhance the robustness and reliability of the results. There is a lack of experimental validation to support the transcriptomic predictions of immune infiltration, such as flow cytometry or immunohistochemical (IHC) staining of immune cell markers. All immune profiling in the study is based solely on transcriptomic data, with no spatial information to determine whether immune cells are excluded, infiltrating, or peripherally localized within the tumor microenvironment. Moreover, the clinical utility of immune signatures could be strengthened by correlating gene expression profiles with treatment responses in cohorts receiving immune checkpoint blockade therapy. In addition, the glioma cell lines LN229 and U251 both harbor p53 mutations, which may limit their ability to represent the broader biological heterogeneity of glioma. Orthotopic or subcutaneous GBM xenograft models should be employed to validate the in vitro findings and assess their relevance in vivo. Lastly, the functions of CACNA2D3, HIST1H3B, and HSPA1A in GBM should be validated through both in vitro and in vivo experiments. Collectively, we have constructed a novel GBM-related model to assess patient prognosis and identified four new signatures for diagnostic prediction, which may prove beneficial for future treatment strategies in GBM patients.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
YY: Data curation, Formal Analysis, Investigation, Methodology, Resources, Software, Writing – original draft. ZC: Data curation, Formal Analysis, Investigation, Methodology, Resources, Software, Writing – original draft. QZ: Data curation, Formal Analysis, Investigation, Methodology, Resources, Software, Writing – original draft. GW: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Validation, Visualization, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by Jiaxing Public Welfare Research Project (2024AY30019) and Zhejiang Basic Public Welfare Research Project (LGF19H160013).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.
References
1
Sanchez GilJRabkinSD. An armed oncolytic virus for GBM destruction. Nat Cancer. (2022) 3:1274–6. doi: 10.1038/s43018-022-00457-z
2
LlaciLRubinJMitraR. EPCO-22. single cell analysis to understand sex differences in GBM outcome. Neuro-Oncology. (2022) 24:vii120–1. doi: 10.1093/neuonc/noac209.457
3
MillerKDOstromQTKruchkoCPatilNTihanTCioffiGet al. Brain and other central nervous system tumor statistics, 2021. CA Cancer J Clin. (2021) 71:381–406. doi: 10.3322/caac.21693
4
NesselhutJMarxDCillienNChangRYBrockmannW-PMartinMet al. Comparison of early versus late onset of cellular immunotherapy in glioblastoma multiforme WHO IV. J Clin Oncol. (2017) 35:e13531–1. doi: 10.1200/JCO.2017.35.15_suppl.e13531
5
StuppRLukasRVHegiME. Improving survival in molecularly selected glioblastoma. Lancet. (2019) 393:615–7. doi: 10.1016/S0140-6736(18)33211-2
6
BeckerAPSellsBEHaqueSJChakravartiA. Tumor heterogeneity in glioblastomas: from light microscopy to molecular pathology. Cancers. (2021) 13:761. doi: 10.3390/cancers13040761
7
KumariSGuptaRAmbastaRKKumarP. Multiple therapeutic approaches of glioblastoma multiforme: From terminal to therapy. Biochim Biophys Acta Rev Cancer. (2023) 1878(4):188913. doi: 10.1016/j.bbcan.2023.188913
8
van SolingeTSNielandLChioccaEABroekmanMLD. Advances in local therapy for glioblastoma - taking the fight to the tumour. Nat Rev Neurol. (2022) 18:221–36. doi: 10.1038/s41582-022-00621-0
9
Cruz Da SilvaEFoppoloSLhermitteBIngremeauMJustinianoHKleinLet al. Bioimaging nucleic-acid aptamers with different specificities in human glioblastoma tissues highlights tumoral heterogeneity. Pharmaceutics. (2022) 14:1980. doi: 10.3390/pharmaceutics14101980
10
VenkataramaniVSchneiderMGiordanoFAKunerTWickWHerrlingerUet al. Disconnecting multicellular networks in brain tumours. Nat Rev Cancer. (2022) 22:481–91. doi: 10.1038/s41568-022-00475-0
11
LinPGanYBHeJLinSEXuJKChangLet al. Advancing skeletal health and disease research with single-cell RNA sequencing. Mil Med Res. (2024) 11:33. doi: 10.1186/s40779-024-00538-3
12
FengDCZhuWZWangJLiDXShiXXiongQet al. The implications of single-cell RNA-seq analysis in prostate cancer: unraveling tumor heterogeneity, therapeutic implications and pathways towards personalized therapy. Mil Med Res. (2024) 11:21. doi: 10.1186/s40779-024-00526-7
13
Degl’InnocentiAdi LeoNCiofaniG. Genetic hallmarks and heterogeneity of glioblastoma in the single-cell omics era. Adv Ther. (2020) 3:1900152. doi: 10.1002/adtp.201900152
14
HuangRHWangLXHeJGaoW. Application and prospects of single cell sequencing in tumors. biomark Res. (2021) 9:88. doi: 10.1186/s40364-021-00336-2
15
XueRLiRBaiF. Single cell sequencing: technique, application, and future development. Sci Bull. (2015) 60:33–42. doi: 10.1007/s11434-014-0634-6
16
ZhouW-MYanY-YGuoQ-RJiHWangHXuT-Tet al. Microfluidics applications for high-throughput single cell sequencing. J Nanobiotechnol. (2021) 19:312. doi: 10.1186/s12951-021-01045-6
17
LedesmaDSymesSRichardsS. Advancements within modern machine learning methodology: impacts and prospects in biomarker discovery. Curr Med Chem. (2021) 28:6512–31. doi: 10.2174/0929867328666210208111821
18
FedikNZubatyukRKulichenkoMLubbersNSmithJSNebgenBet al. Extending machine learning beyond interatomic potentials for predicting molecular properties. Nat Rev Chem. (2022) 6:653–72. doi: 10.1038/s41570-022-00416-3
19
KleppeASkredeOJDe RaedtSLiestolKKerrDJDanielsenHE. Designing deep learning studies in cancer diagnostics. Nat Rev Cancer. (2021) 21:199–211. doi: 10.1038/s41568-020-00327-9
20
LiZZhuTWuYYuYZangYYuLet al. Functions and mechanisms of non-histone post-translational modifications in cancer progression. Cell Death Discov. (2025) 11:125. doi: 10.1038/s41420-025-02410-2
21
MadhukarGHaqueMAKhanSKimJJDanishuddin. E3 ubiquitin ligases and their therapeutic potential in disease Management. Biochem Pharmacol. (2025) 236:116875. doi: 10.1016/j.bcp.2025.116875
22
HongZLiuFZhangZ. Ubiquitin modification in the regulation of tumor immunotherapy resistance mechanisms and potential therapeutic targets. Exp Hematol Oncol. (2024) 13:91. doi: 10.1186/s40164-024-00552-0
23
XuLWangBGangZHanZWangALiuQet al. Ubiquitin-conjugating enzyme E2S decreases the sensitivity of glioblastoma cells to temozolomide by upregulating PGAM1 via the interaction with OTUB2. Int J Biol Macromol. (2025) 302:140583. doi: 10.1016/j.ijbiomac.2025.140583
24
SunJZhaoWZhangLWuSXueSCaoHet al. Centromere protein U mediates the ubiquitination and degradation of RPS3 to facilitate temozolomide resistance in glioblastoma. Drug Resist Update. (2025) 80:101214. doi: 10.1016/j.drup.2025.101214
25
LiuCHeWZhaoHWangSQianZ. KRT80, regulated by RNF8-mediated ubiquitination, contributes to glucose metabolic reprogramming and progression of glioblastoma. Neurochem Res. (2025) 50:128. doi: 10.1007/s11064-025-04380-4
26
XuYHongYYanTSunQYuanFLiangSet al. Tryptophan metabolic enzyme IL4I1 inhibits ferroptosis by decreasing ubiquitination of nrf2 via I3P in glioblastoma. Cell Prolif. (2025) 58:e13816. doi: 10.1111/cpr.13816
27
CloughEBarrettT. “The gene expression omnibus database”. In: Statistical genomics. Methods in molecular biology, vol 1418. New York: Springer (2016). p. 93–110.
28
StelzerGRosenNPlaschkesIZimmermanSTwikMFishilevichSet al. The geneCards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinf. (2016) 54:31–3. doi: 10.1002/cpbi.5
29
HanHChoJWLeeSYunAKimHBaeDet al. TRRUST v2: an expanded reference database of human and mouse transcriptional regulatory interactions. Nucleic Acids Res. (2018) 46:D380–6. doi: 10.1093/nar/gkx1013
30
SzklarczykDKirschRKoutrouliMNastouKMehryaryFHachilifRet al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. (2023) 51:D638–46. doi: 10.1093/nar/gkac1000
31
ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6
32
TangZLiCKangBGaoGLiCZhangZ. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. (2017) 45:W98–W102. doi: 10.1093/nar/gkx247
33
SanzHValimCVegasEOllerJMReverterF. SVM-RFE: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinf. (2018) 19:1–18. doi: 10.1186/s12859-018-2451-4
34
ParkTCasellaG. The bayesian lasso. J Am Stat Assoc. (2008) 103:681–6. doi: 10.1198/016214508000000337
35
ChenBKhodadoustMSLiuCLNewmanAMAlizadehAA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. (2018) 1711:243–59.
36
ZhangBHorvathS. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. (2005) 4. doi: 10.2202/1544-6115.1128
37
DuanK-BRajapakseJCWangHAzuajeF. Multiple SVM-RFE for gene selection in cancer classification with expression data. IEEE Trans Nanobiosci. (2005) 4:228–34. doi: 10.1109/TNB.2005.853657
38
ZhaoNZhangJZhaoLFuXZhaoQChaoMet al. Long Noncoding RNA NONHSAT079852.2 Contributes to GBM Recurrence by Functioning as a ceRNA for has-mir-10401-3p to Facilitate HSPA1A Upregulation. Front Oncol. (2021) 11:636632. doi: 10.3389/fonc.2021.636632
39
JinYCuiDRenJWangKZengTGaoL. CACNA2D3 is downregulated in gliomas and functions as a tumor suppressor. Mol Carcinogenesis. (2017) 56:945–59. doi: 10.1002/mc.22548
40
VuongHGNgoTNMLeHTDunnIF. The prognostic significance of HIST1H3B/C and H3F3A K27M mutations in diffuse midline gliomas is influenced by patient age. J Neurooncol. (2022) 158:405–12. doi: 10.1007/s11060-022-04027-2
41
WangJSongY. Single cell sequencing: a distinct new field. Clin Trans Med. (2017) 6:e10. doi: 10.1186/s40169-017-0139-4
42
TangLHuangZPMeiHHuY. Insights gained from single-cell analysis of chimeric antigen receptor T-cell immunotherapy in cancer. Mil Med Res. (2023) 10:52. doi: 10.1186/s40779-023-00486-4
43
BaiYMYangFLuoPXieLLChenJHGuanYDet al. Single-cell transcriptomic dissection of the cellular and molecular events underlying the triclosan-induced liver fibrosis in mice. Mil Med Res. (2023) 10:7. doi: 10.1186/s40779-023-00441-3
44
ChenZLiuCZhangCXiaYPengJMiaoCet al. Machine learning-based discovery of UPP1 as a key oncogene in tumorigenesis and immune escape in gliomas. Front Immunol. (2024) 15:1475206. doi: 10.3389/fimmu.2024.1475206
45
HuYQinSDengR. Impact of glioma metabolism-related gene ALPK1 on tumor immune heterogeneity and the regulation of the TGF-beta pathway. Front Immunol. (2024) 15:1512491. doi: 10.3389/fimmu.2024.1512491
46
YangYLiuZWangZFuXLiZLiJet al. Large-scale bulk and single-cell RNA sequencing combined with machine learning reveals glioblastoma-associated neutrophil heterogeneity and establishes a VEGFA(+) neutrophil prognostic model. Biol Direct. (2025) 20:45. doi: 10.1186/s13062-025-00640-z
47
GaoXJiangMChuYHanYJinYZhangWet al. ETV4 promotes pancreatic ductal adenocarcinoma metastasis through activation of the CXCL13/CXCR5 signaling axis. Cancer Lett. (2022) 524:42–56. doi: 10.1016/j.canlet.2021.09.026
48
LuBCCebrianCChiXKuureSKuoRBatesCMet al. Etv4 and Etv5 are required downstream of GDNF and Ret for kidney branching morphogenesis. Nat Genet. (2009) 41:1295–302. doi: 10.1038/ng.476
49
ChengTZhangZChengYZhangJTangJTanZet al. ETV4 promotes proliferation and invasion of lung adenocarcinoma by transcriptionally upregulating MSI2. Biochem Biophys Res Commun. (2019) 516:278–84. doi: 10.1016/j.bbrc.2019.06.115
50
DumortierMLadamFDamourIVacherSBiècheIMarchandNet al. ETV4 transcription factor and MMP13 metalloprotease are interplaying actors of breast tumorigenesis. Breast Cancer Res. (2018) 20:73. doi: 10.1186/s13058-018-0992-0
51
HellwegeJNStallingsSTorstensonESCarrollRBorthwickKMBrilliantMHet al. Heritability and genome-wide association study of benign prostatic hyperplasia (BPH) in the eMERGE network. Sci Rep. (2019) 9:6077. doi: 10.1038/s41598-019-42427-z
52
YuanZYDaiTWangSSPengRJLiXHQinTet al. Overexpression of ETV4 protein in triple-negative breast cancer is associated with a higher risk of distant metastasis. Onco Targets Ther. (2014) 7:1733–42. doi: 10.2147/ott.S66692
53
FonsecaASRamãoABürgerMCde SouzaJESZanetteDLde MolfettaGAet al. ETV4 plays a role on the primary events during the adenoma-adenocarcinoma progression in colorectal cancer. BMC Cancer. (2021) 21:207. doi: 10.1186/s12885-021-07857-x
54
YoshiyaSItohSYoshizumiTYugawaKKuriharaTToshimaTet al. Impact of capicua on pancreatic cancer progression. Ann Surg Oncol. (2021) 28:3198–207. doi: 10.1245/s10434-020-09339-z
55
SingsuksawatEThuwajitCCharngkaewKThuwajitP. Increased ETV4 expression correlates with estrogen-enhanced proliferation and invasiveness of cholangiocarcinoma cells. Cancer Cell Int. (2018) 18:25. doi: 10.1186/s12935-018-0525-z
56
BabalYKKandemirBKurnazIA. Gene regulatory network of ETS domain transcription factors in different stages of glioma. J Personalized Med. (2021) 11:138. doi: 10.3390/jpm11020138
57
XiaoJYangSShenPWangYSunHJiFet al. Phosphorylation of ETV4 at Ser73 by ERK kinase could block ETV4 ubiquitination degradation in colorectal cancer. Biochem Biophys Res Commun. (2017) 486:1062–8. doi: 10.1016/j.bbrc.2017.03.163
58
ZhangQLiuSWangHXiaoKLuJChenSet al. ETV4 mediated tumor-associated neutrophil infiltration facilitates lymphangiogenesis and lymphatic metastasis of bladder cancer. Adv Sci (Weinh). (2023) 10:e2205613. doi: 10.1002/advs.202205613
59
GuidiM. P14.43 Histone H3F3A and HIST1H3B K27M mutations in pediatric high-grade gliomas: the Florentine experience. Neuro-Oncology. (2019) 21:iii76–7. doi: 10.1093/neuonc/noz126.278
60
LópezGYVan ZiffleJOnoderaCGrenertJPYehIBastianBCet al. The genetic landscape of gliomas arising after therapeutic radiation. Acta Neuropathol. (2019) 137:139–50. doi: 10.1007/s00401-018-1906-z
61
Martínez-RicarteFMayorRMartínez-SáezERubio-PérezCPinedaECorderoEet al. Molecular diagnosis of diffuse gliomas through sequencing of cell-free circulating tumor DNA from cerebrospinal fluid. Clin Cancer Res. (2018) 24:2812–9. doi: 10.1158/1078-0432.Ccr-17-3800
62
ThorsteinsdottirJStanglSFuPGuoKAlbrechtVEigenbrodSet al. Overexpression of cytosolic, plasma membrane bound and extracellular heat shock protein 70 (Hsp70) in primary glioblastomas. J Neuro-Oncol. (2017) 135:443–52. doi: 10.1007/s11060-017-2600-z
63
ChenYWuGLiMHesseMMaYChenWet al. LDHA-mediated metabolic reprogramming promoted cardiomyocyte proliferation by alleviating ROS and inducing M2 macrophage polarization. Redox Biol. (2022) 56:102446. doi: 10.1016/j.redox.2022.102446
64
LinZPGanGXuXWenCDingXChenXYet al. Comprehensive PTM profiling with SCASP-PTM uncovers mechanisms of p62 degradation and ALDOA-mediated tumor progression. Cell Rep. (2025) 44:115500. doi: 10.1016/j.celrep.2025.115500
65
DengLChenLZhaoLXuYPengXWangXet al. Ubiquitination of Rheb governs growth factor-induced mTORC1 activation. Cell Res. (2019) 29:136–50. doi: 10.1038/s41422-018-0120-9
66
YaoYHongSYoshidaSSwaroopVCurtinBInokiK. The Cullin3-Rbx1-KLHL9 E3 ubiquitin ligase complex ubiquitinates Rheb and supports amino acid-induced mTORC1 activation. Cell Rep. (2025) 44:115101. doi: 10.1016/j.celrep.2024.115101
67
JiangHDingYLinXTianQLiuYHeHet al. Malvidin attenuates trauma-induced heterotopic ossification of tendon in rats by targeting Rheb for degradation via the ubiquitin-proteasome pathway. J Cell Mol Med. (2024) 28:e18349. doi: 10.1111/jcmm.18349
68
SeoGKimSKByunYJOhEJeongSWChaeGTet al. Hydrogen peroxide induces Beclin 1-independent autophagic cell death by suppressing the mTOR pathway via promoting the ubiquitination and degradation of Rheb in GSH-depleted RAW 264.7 cells. Free Radic Res. (2011) 45:389–99. doi: 10.3109/10715762.2010.535530
69
SunDDuanXLiNQiaoOHouYMaZet al. Construction of ubiquitination-related risk model for predicting prognosis in lung adenocarcinoma. Sci Rep. (2025) 15:11787. doi: 10.1038/s41598-025-92177-4
Summary
Keywords
Glioma, ScRNA-seq, prognosis, biomarker, machine learning
Citation
Yan Y, Chu Z, Zhong Q and Wang G (2025) Single-cell sequencing combined with machine learning to identify glioma biomarkers and therapeutic targets. Front. Oncol. 15:1629102. doi: 10.3389/fonc.2025.1629102
Received
15 May 2025
Accepted
07 July 2025
Published
24 July 2025
Volume
15 - 2025
Edited by
Xiangpeng Dai, Jilin University, China
Reviewed by
Jia Liu, Xi’an Jiaotong University, China
Lixia Wang, Soochow University, China
Xuelin Zhong, Lenox Hill Hospital, United States
Updates
Copyright
© 2025 Yan, Chu, Zhong and Wang.
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: Genghuan Wang, 13967380141@163.com
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.