Abstract
Background:
Laryngeal cancer (LCA) is the second most common type of head and neck malignancy, characterized by high recurrence rates and poor overall survival (OS). However, progress in curing LCA through molecular-targeted diagnostics and therapies is slow and limited. The occurrence and progression of cancer are closely associated with metabolic reprogramming. Therefore, this study aimed to identify metabolism-related LCA subtypes through a comprehensive analysis of transcriptomic, mutational, methylation, and single-cell RNA sequencing, in hopes of finding factors which influences the prognosis of LCA.
Methods:
First, to identify metabolism-related LCA subtypes, data from 114 patients with LCA from The Cancer Genome Atlas (TCGA) dataset were collected for an unsupervised clustering analysis, which focused on the expression characteristics of survival-related metabolic genes. Subsequently, prognostic and diagnostic models have been developed using machine learning techniques. Specifically, the prognostic model utilized the least absolute shrinkage and selection operator (LASSO) Cox regression, whereas the diagnostic model was built using the Random Forest (RF) algorithm. Furthermore, to ensure the reproducibility, the results of the subtypes and models were validated using three independent bulk RNA datasets and a scRNA-seq dataset.
Results:
Two robust subtypes were identified and independently validated. Each subtype has a distinct prognostic outcomes and molecular features. Specifically, the LCA1 subtype exhibited better prognosis, enriched metabolic pathways, and higher mutation frequencies. Notably, significant damaging mutations in the methyltransferases NSD1 were observed in this subtype. In contrast, the LCA2 subtype was associated with poorer prognosis, higher immune infiltration, and elevated methylation levels. Moreover, in LCA2 tumors, higher levels of T cell/APC co-inhibition and inhibitory checkpoints were observed. In addition, the diagnostic model demonstrated strong performance, achieving an area under the curve (AUC) values of 1.000 in the training group and 0.947 in the validation group. The prognostic model effectively predicted patient outcomes, with the RiskScore emerging as an independent prognostic factor.
Conclusion:
This study offers new perspectives for patient stratification and presents opportunities for therapeutic development in LCA. Furthermore, we explored the potentials of several key tumor markers for both diagnosis and prognosis prediction.
1 Introduction
As the second most common head and neck malignancy, Laryngeal cancer (LCA) accounted for over 184,615 new cases and approximately 99,840 deaths globally in 2020 (). Although medical diagnostic and treatment methods are continuously evolving, the diagnosis of LCA still relies primarily on endoscopic and pathological examinations (). LCA treatment typically involves surgery combined with radiotherapy and chemotherapy (). Cetuximab, Pembrolizumab, and Nivolumab have been used in patients with LCA, particularly in cases of recurrence and metastasis (–). Endoscopic techniques such as transoral laser microsurgery (TLM) and transoral robotic surgery (TORS) are increasingly being employed for the effective removal of suitable early-stage LCA (). However, owing to the concealed nature of the laryngeal site and the lack of early diagnostic methods, LCA is always diagnosed at advanced stages (III-IV) (). This late diagnosis limits the patients’ therapeutic options and significantly affects their quality of life. Numerous studies have increasingly focused on exploring the molecular biomarkers of LCA, however, despite the identification of some potential molecular targets such as TLSs, NOTCH1, and BMP2, new effective therapeutic targets for LCA have not yet been confirmed (–). Hence, it is crucial to identify early diagnostic markers and explore the molecular mechanisms which affect the prognosis of LCA.
The metabolic processes of malignant tumors were different from those of normal tissues, because they require higher and faster amounts of materials to support the proliferation of tumor cells. Consequently, tumor cells tend to undergo metabolic reprogramming to meet these elevated demands (, ). Because of such differences between malignant tumors and normal tissues, it is possible to explore targeted therapies against the metabolic dependence of tumor cells. According to some scholars, several anti-tumor pathways targeting metabolic enzymes have been identified, including: 5-fluorouracil, Capecitabine, Pemetrexed, and Raltitrexed targeting to Thymidylate synthase (TS); Methotrexate, Pemetrexed targeting to Dihydrofolate reductase (DHFR); Pemetrexed targeting to Glycinamide ribonucleotide formyltransferase (GARFT), et al. (). In addition, many studies have employed multi-omics strategies for metabolic subtype classification. Yuan et al. revealed the metabolic heterogeneity associated with HER2 in gastric cancer in response to immunotherapy and neoadjuvant chemotherapy (). Similarly, Li et al. employed a multi-omics analysis to reveal TAM2-related glycolysis and pyruvate metabolism remodeling in pancreatic cancer (). Hepatocellular carcinoma was classified into three metabolic subclasses: C1, C2, and C3 (). Similarly, uterine corpus endometrial carcinoma is classified into two types, C1 and C2 (). However, the metabolism-related molecular characteristics of LCA remained unreported.
In this study, we utilized unsupervised clustering analysis to categorize LCA into two distinct subtypes based on metabolism-related gene expression profiles for the first time. The reproducibility of the two subtypes was confirmed using two independent datasets. Moreover, each subtype features by different somatic alterations, immune infiltration profiles, DNA methylation patterns, metabolic features, and clinical outcomes. We further developed diagnostic and prognostic models by employing the Random Forest (RF) and least absolute shrinkage and selection operator (LASSO) Cox regression methods.
2 Materials and methods
2.1 Data source
The data used in this study was primarily obtained from two public databases. RNA-seq data, clinical information, DNA methylation data, and DNA mutation data from 114 LCA samples and 12 matched normal mucosal samples were downloaded from The Cancer Genome Atlas (TCGA, available at http://cancergenome.nih.gov). This dataset was employed to identify of LCA subtypes, and the construction of diagnostic and prognostic models. Additionally, a single-cell RNA sequencing dataset (GSE252490) () and three bulk RNA sequencing datasets (GSE130605 (), GSE27020 (), GSE142083 () were acquired from National Center for Biotechnology Information’s Gene Expression Omnibus (GEO, accessible at https://www.ncbi.nlm.nih.gov/geo). These datasets were used to validate the reliability and generalizability of the findings. Specifically, GSE130605, which included 50 LCA tissue samples and 50 matched normal tissue samples, was used to validate the metabolic subtype results. GSE27020, containing 109 LCA samples, was used for both metabolic subtype validation and prognostic risk model validation. GSE142083, comprising 53 LCA samples and 53 matched normal samples, was employed to validate the diagnostic prediction models. GSE252490 consists of three LCA samples with lymph node metastasis and was used to validate metabolic subtypes, diagnostic prediction models, and prognostic risk assessment models at the single-cell level.
2.2 Discovery and confirmation of subtypes related to metabolism
The previous 2752 genes associated with metabolism were collected for further analysis (). First, by selecting genes with high median absolute deviation values (MAD > 0.5) in the expression profiles of the TCGA dataset, we obtained 1534 genes. Using COX regression analysis of survival in the TCGA cohort with the R package “survival”, genes significantly associated with prognosis were identified (P < 0.05). These genes across all patients were then utilized for consensus clustering analysis with the R package “ConsensusClusterPlus” (, ). To identify robust clusters, the cumulative distribution function (CDF) and consensus heatmap were applied to choose the optimal K. Consequently, the TCGA dataset samples were divided into two subtypes, designated as LCA1 and LCA2. Next, by using the R package “DESeq2” (), differential gene expression analysis was conducted between the two subtypes. Genes with an absolute |log2 fold change (FC)| > 0.5 and an adjusted P-value < 0.05 were considered potential candidate genes. Finally, the same consensus clustering analysis was performed on datasets GSE130605 and GSE27020 using the same candidate genes to validate metabolism-associated subtypes in other cohorts.
2.3 Immunoinfiltration analysis
Single-sample GSEA (ssGSEA) was applied to estimate the relative fraction of 13 immune-related functions and 23 immune cells between the two subtypes by using the “Gene Set Variation Analysis (GSVA)” package (). We assessed the infiltration of immune components, stromal components, and tumor components using the ESTIMATE algorithm, which provided the ImmuneScore, StromalScore, ESTIMATEScore, and TumorPurity (). Additionally, we compared the expression levels of several immune checkpoint genes between the two subtypes ().
2.4 Calculation of gene signatures related to metabolism
A total of 114 metabolism-relevant pathways were gathered from previously published research (, ). The R package “GSVA” () was used to quantify the enrichment degree of metabolism-relevant signatures for each sample. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) () enrichment analyses were applied to the highly expressed metabolism-relevant genes of the two subtypes, respectively.
2.5 Characterization of mutations
Mutation annotation format (MAF) files for LCA were downloaded from the TCGA and processed by the “maftools” R package. Besides, mutation information was also analyzed by using the “maftools” R package () between the two subtypes.
2.6 Differences in DNA methylation between the subtypes
To explore the epigenetic differences between the two subtypes, methylation data for LCA was downloaded from the TCGA database. Differentially methylated CpG sites were identified by using the “limma” R package (). Then, 1471 CpG probes were considered as the most differentially methylated sites between the two subtypes (Absolute |log2 fold change (FC)| > 0.25 and Adjusted P-value < 0.01). In order to better explore the epigenetic differences, GO enrichment analysis () was applied to the highly methylated genes of the two subtypes, respectively.
2.7 Identification of metabolism-associated diagnosis model for LCA
To explore a robust diagnostic model, we initially identified differentially expressed genes between tumor and normal samples from the TCGA dataset by using the R package “DESeq2” (). By applying the stringent criteria of absolute |log2 fold change (FC)| > 2 and an adjusted P-value < 0.01, we successfully identified 293 differentially expressed genes. Subsequently, using the R package “randomForest” (), we developed a diagnostic model to distinguish tumor tissues from normal tissues. A Nomogram () was generated to facilitate clinical application by using the top six significant genes which were identified in the random forest model, by employing the R package “rmda”. Then, the R package “ROC” was used to evaluate the prediction accuracy of the Nomogram. Furthermore, decision curve analysis () was conducted to assess the clinical utility of the diagnostic model by using the R package “rmda”. For external validation, we used an independent GEO dataset (GSE142083, n = 106) to validate the performance and robustness of the diagnostic model.
2.8 Identification of metabolism-associated genes prognosis model
To construct a metabolism-related prognosis model, the 56 final metabolism-related genes were incorporated into LASSO Cox regression analyses. Kaplan-Meier (KM) curves were constructed to compare the survival differences between the high- and low-risk score groups (). To estimate the model’s predictive accuracy, the package “ROC” was used to assess the area under the curve (AUC) values of ROC curves for survival state and 1-, 3-, and 5-year survival rates. The RiskScore and other clinical parameters predictions of the 5-year survival rate in the TCGA dataset were also assessed by using ROC curves. Univariate and multivariate Cox regression analyses of clinicopathological features were further conducted to evaluate whether the RiskScore was an independent predictive factor. The GEO dataset (GSE27020, n = 109) was used for validation. The Wilcoxon test was used to analyze the differences in prognostic signature genes’ expression and RiskScore distribution between the two subtypes. ssGSEA was applied to estimate the relative fraction of 13 immune-related functions and 23 immune cells by using the “GSVA” package (). OncoPredict was used to estimate differences in drug sensitivity between the two groups using the R package “oncoPredict” (, ).
2.9 Single-cell analysis further validated the subtypes, diagnostic models and prognostic models
To analyze the single-cell RNA sequencing from the GEO database (GSE252490) for laryngeal squamous cell carcinoma (LSCC), we used the packages including “Seurat” and “SingleR” () to carry on standard quality control procedures including “PercentageFeatureSet”, “NormalizeData”, “HarmonyIntegration”, “RunPCA”, “FindNeighbors”, “FindClusters”, “RunUMAP”, and “FindAllMarkers” functions. Besides, we used the “SingleR” package and the existing markers from the previous studies () to distinguish the cell types. Copy number variation (CNV) scores for epithelial cells were calculated using the “CopyKat” package (), with B cells as the reference. To understand the biological roles of marker genes within each cluster of malignant epithelial cells, we utilized the “clusterProfiler” and “org.Hs.eg.db” software packages. Additionally, we investigated the developmental paths of various malignant epithelial cell categories using the Monocle2 algorithm (). Furthermore, we employed the “AddModuleScore” () to compute scores related to subtypes, diagnostic model, and prognostic model, which were based on the mean expression levels of genes.
2.10 Statistical analysis
All statistical analyses were performed using R software (version 4.3.2 and version 4.4.0). The Wilcoxon test was used to analyze the differences among groups. Survival data were analyzed using the KM curves. In particular, P < 0.05 was considered to be statistically significant.
3 Results
3.1 Two distinct tumor subtypes in LCA were identified by Consensus clustering
To investigate the metabolic heterogeneity in LCA, we collected 2752 metabolism-related genes from previous studies for clustering analysis (). The workflow of this study is illustrated in Figure 1. Initially, we adopted four study methods: MAD > 0.5, univariate COX regression analysis of survival, consensus clustering analysis, and differential expression analysis of the two subtypes. This process yielded 56 genes and two metabolism-related subtypes including LCA1 and LCA2 (Figure 2A, Supplementary Figure S1). Principal component analysis (PCA) was performed to assess the composition of the two subtypes and confirm the stable expression differences between the two subtypes (Figure 2B). Survival analysis revealed that the prognosis of LCA1 was significantly better than that of LCA2 (P<0.001) (Figure 2C), whereas the clinical characteristics of the two subtypes did not show any significant differences (Supplementary Table S1). Finally, we validated the reproducibility of the two subtypes in the two independent datasets including GSE130605 and GSE27020. We identified similar metabolism-related molecular subtypes of LCA in both the GSE130605 (Figures 3A–C; Supplementary Figure S2), and the GSE27020 datasets (Figures 3D–F, Supplementary Figure S3) by performing consensus clustering analysis of the final candidate genes and conducting PCA on the entire transcriptome data.
Figure 1
Figure 2
Figure 3
3.2 Metabolic profiling of the two subtypes showed significant differences
We further investigated the metabolic characteristics of the two subtypes. Initially, the GSVA method was used to calculate the scores of 114 metabolic related pathways’ scores. LCA1 exhibited 44 specific metabolic signatures, while LCA2 had nine (Figure 4A; Supplementary Table S2). Then, this study applied differentially expressed gene (DEGs) analysis to analyze the 2752 metabolism-related genes, and 143 highly expressed genes in LCA1 and 194 in LCA2 were identified (Absolute |log2 fold change (FC)| > 1 and Adjusted P-value < 0.05). Both KEGG pathway and GO functional enrichment analyses were used to analyze the highly expressed genes in the two subtypes, respectively (Figures 4B, C; Supplementary Figure S4A). Finally, the results of the enrichment analyses indicated that LCA1 had higher levels of metabolism processes compared with LCA2, which were related to amino acids, lipids, and vitamins. In contrast, LCA2 tumors showed few enrichments of metabolic signatures such as glycosylation, which indicated their lower metabolic activities.
Figure 4
3.3 LCA2 showed a higher degree of immune infiltration
It is well established that the amino acid metabolism reprogramming in tumors plays a critical role in controlling the differentiation and function of immune cells (). To explore the differences in immune infiltration between LCA1 and LCA2, we calculated the ImmuneScore, StromalScore, ESTIMATEScore, and TumorPurity by using the ESTIMATE algorithm. Consequently, LCA2 had a higher StromalScore and ESTIMATEScore, whereas LCA1 had a higher TumorPurity (Figure 5A). Then, the result of “ssGSEA” analysis showed higher scores for the three immune-related functions and 11 immune cells in LCA2 (Figures 5B, C). Additionally, because T cell and natural killer cell inhibition are known to be important mechanisms for immune escape in cancer, we continued to explore the different expression levels of various immune checkpoint genes between the two subtypes. Consequently, several checkpoint genes showed higher expression levels in LCA2, which indicated that a higher level of immunosuppression was observed in this group of LCA2 (Figure 5D).
Figure 5
3.4 Nonsense mutations of NSD1 and genome-wide hypomethylation occurred in the LCA1
The tumoral genomic landscape has been shown to correlated with antitumor immunity. To investigate whether differences exist in somatic mutation frequencies and patterns between the two subtypes, we compared the somatic variations of the two subtypes and found a higher tumor mutation burden in LCA1. In addition, the genes with high mutation frequencies in each subtype are shown in Figures 6A, B. Consequently, we found that 11 genes were shared between the two subtypes, whereas each subtype had nine unique high-frequency mutated genes. Notably, the mutation frequency of the histone methyltransferase NSD1(69%) in LCA1 was much higher than that in LCA2. However, in terms of variant and SNV classification, the two metabolism-related subtypes showed similar trends (Figures 6C, D).
Figure 6
The histone methyltransferase NSD1 is strongly associated with DNA methylation. To further investigate the methylation differences between the two subtypes, we examined the different methylated CpG sites in the TCGA cohort. A total of 1471 CpG probes were significantly differentially methylated between the subtypes (Absolute |log2 fold change (FC)| > 0.25 and Adjusted P-value < 0.01) (Figure 7A). Genes with methylation changes in LCA1 were highly enriched in calcium-dependent cysteine-type endopeptidase activity, negative regulation of neuronal differentiation, Fc-gamma receptor signaling pathway, and aggresome. Differentially methylated genes in LCA2 were highly enriched in hormone secretion, hormone transport, regulation of the G protein-coupled receptor signaling pathway, regulation of hormone secretion, and DNA-binding transcription activator activity. (Figure 7B).
Figure 7
3.5 The diagnostic panel, which was derived from machine learning, exhibited an excellent performance
To develop a robust diagnostic panel, we first identified 273 metabolic genes with significant differential expression between tumor and normal samples through DEGs analysis in the TCGA cohort (Absolute |log2 fold change (FC)| > 2 and Adjusted P-value < 0.01) (Figure 8A, Supplementary Table S3). Then, machine learning was used to develop a model to predict the clinical status in this study. The random forest model was trained using minimum error regression trees, and the six most essential genes were selected to discriminate the tumor from normal samples (Figures 8B, C). Next, a nomogram was developed for clinical use using logistic regression analysis (Figure 8D). To further investigate the characteristics of our diagnostic model, we examined the distribution of the model genes in the TCGA samples and the independent dataset GSE142083. The heatmap demonstrated substantial expression differences in the model genes between the tumor and normal samples (Figures 8E, F). In order to validate the performance of the diagnostic panel, we calculated the predicted scores of each sample. ROC analysis showed an AUCs of 1 for TCGA and 0.947 for GSE142083, indicating that the model had an excellent diagnostic predictive value (Figures 8G, J). We also observed significant differences in the predictive scores between tumor and normal samples (Figures 8I, L). Decision Curve Analysis (DCA) plots also confirmed the predictive superiority of our diagnostic model in both TCGA and GSE142083 datasets (Figures 8H, K).
Figure 8
3.6 The metabolic prognostic model accurately predicted the outcomes of patients with LCA
As a precise prognosis could enable precise intervention in tumors and help in the clinical treatment of patients, we also attempted to develop a machine-learning-derived prognostic model. First, using the LASSO regression algorithm and multivariate COX regression analysis, we selected ten essential genes, with six acting as protective factors (GPT, CHRNB1, CYP2D6, PLCH2, ABCF2, and EPHX2) and four as risk factors (HS3ST2, TPCN2, TRPC1, and PHYHD1) (Figures 9A–C). The RiskScore formula is shown in the Additional file1, RiskScore formular. Based on the median of the RiskScore, the TCGA samples were classified into high- and low-risk groups. Next, the result of KM analysis revealed significantly longer survival time in the low-risk group (Figure 9G). To assess the accuracy of our model, we performed ROC analysis, yielding AUC values of 0.802, 0.912, and 0.920 for the 1-, 3-, and 5-year survival rates, respectively (Figure 9D). The results of the comparative analysis indicated that the RiskScore had the highest predictive ability for 5-year survival compared with other clinical parameters (Figure 9F). Moreover, the clinical characteristics did not show significant differences between the two groups, implying that the prognostic model was actually an independent factor for prognosis prediction (Supplementary Table S4). In addition, we used the independent dataset GSE27020 from GEO to validate the model. However, owing to the lack of expression data for TPCN2 and PHYHD1, the RiskScores for GSE27020 were not perfectly verified for repeatability. However, the results of the KM analysis still showed significant differences in disease-free survival (DFS) between the groups (Figure 9H). Furthermore, ROC analysis for GSE27020 yielded AUC values of 0.649, 0.694, and 0.776 for the 1-, 3-, and 5-year progression-free survival rates, respectively (Figure 9E).
Figure 9
To better investigate the characteristics of our prognostic model, we examined the distribution of the RiskScore in TCGA samples (Figures 10A–D). Additionally, univariate and multivariate Cox regression analyses were performed to assess the prognostic significance of the variable factors, including sex, age, stage, Tumor size (T), Nodal status (N), and RiskScore. As a result, Sex, N, and RiskScore were identified as significant prognostic factors (Figure 10B, Supplementary Figure S4B). Notably, the distribution of the RiskScore between the two metabolic subtypes showed significant differences in both the TCGA and GSE27020 datasets (Figures 10E, F).
Figure 10
We also explored the differences in immune cell infiltration between the high- and low-risk groups. Several conclusions were drawn. Firstly, the result of ssGSEA analysis indicated that the scores of 5 immune cell scores and 1 immune-related function were higher in the high-risk group (Supplementary Figures S5A–C). Furthermore, correlation analysis between immune cells and the RiskScore was conducted using TIMER, XCELL, ABS, CIBERSORT-ABS, QUANTISEQ, MCPCOUNTER, EPIC, and CIBERSORT, with significant results (P < 0.05) displayed (Supplementary Figure S5D). Finally, OncoPredict was utilized to forecast drug sensitivity in the LCA. Differences in drug sensitivity differences with P < 0.001 between the two groups were presented in Supplementary Figure S5E.
3.7 Subtype, diagnostic, and prognostic related genes were analyzed in the single-cell sequencing
To explore gene expression related to metabolic subtypes, diagnostic and prognostic models at the single-cell sequencing level, we analyzed a single-cell dataset comprising three samples, GSE142083. After a series of quality control procedures, 16,452 cells were selected for further analysis. Initially, using an appropriate resolution, the cells were divided into 14 clusters and annotated with five cell types (Figure 11A; Supplementary Figure S6A). Marker genes for the five cell types are shown in (Figures 11B, D). The most differentially expressed genes in each cell type are shown in Figure 11C. Next, epithelial-derived cells were extracted for further analysis. Epithelial-derived cells were classified into six clusters using an appropriate resolution (Figure 12B). Duplicates were identified and removed from the six clusters (Supplementary Figures S6B, C). Furthermore, the “CopyKat” algorithm was then applied to each sample individually, which classified epithelial-derived cells into malignant and non-malignant categories (Figures 12A, C). Then, we examined the distribution of diagnostic model-related genes in epithelial-derived cells (Figure 12D) and scored these cells by using the “AddModuleScore” function. Notably, the diagnostic scores for malignant and non-malignant epithelial cells were significantly different (Figure 12E).
Figure 11
Figure 12
In addition, malignant epithelial cells (6451 cells) were then extracted and classified into five clusters using an appropriate resolution (Figure 13A). Highly expressed genes in the five clusters were analyzed by using Gene Ontology (GO) enrichment (Figure 13B). Next, the five clusters were categorized into high (G3 and G4) and low (G0, G1, and G2) variant groups based on CNV scores, which were obtained from the “CopyKat” algorithm (Figure 13C). Furthermore, we also investigated the distribution of genes related to prognostic risk models in malignant epithelial cells and scored the malignant epithelial cells to obtain the prognostic scores by using the ‘AddModuleScore’ function (Supplementary Figure S6D). The “AddModuleScore” function was used to score genes associated with subtyping (Supplementary Table S5). Cells were scored based on the characteristic genes that were highly expressed in the two subtypes, and subsequently were divided into LCA1 (G2 and G3) and LCA2 (G0, G1, and G4) groups (Figure 13D). The relationships between CNV groups and prognostic scores was also investigated. Significant differences in the prognostic scores were observed between the high- and low-variant groups (Figure 13F). Additionally, notable discrepancies in prognostic scores were observed between the metabolic subtypes (Figure 13G). Finally, a pseudotime trajectory analysis using “Monocle 2” was performed to explore the underlying evolution of epithelial cells with diverse CNV scores and prognostic scores (Figure 13E).
Figure 13
4 Discussion
Currently, the classification of LCA mainly relies on the sites of its occurrence and the types of pathology, but remains inadequate for treatment guidance (). In addition, the lack of diagnostic and therapeutic biomarkers has led to delayed late-stage diagnosis and unimproved 5-year survival rates in LCA. To enrich the classification method, we classified LCA into two subtypes based on the 56 metabolism-related genes. And we further explored the differences in clinical and metabolic characteristics, immune infiltration, DNA mutations, DNA methylation, transcriptional data, and single-cell data between the two subtypes. In addition, we established highly reliable diagnostic and prognostic models using machine learning methods. These findings were validated to ensure robustness by using three independent bulk RNA datasets and a single-cell RNA dataset. In conclusion, our findings extended the molecular subtyping of LCA and deepen our understanding of metabolic heterogeneity within this tumor.
There is a consensus that the metabolic reprogramming of tumors leads to distinct metabolic profiles compared to normal tissues (, ), however, the metabolism of different tumors and different sections of the same tumor is also heterogeneous (). Moreover, the metabolic preferences of tumors change dynamically during cancers progression (). Tumor metabolic heterogeneity is a prominent feature of tumor development and has an important impact on the effectiveness of treatment and prognosis (). A previous study demonstrated that metabolic heterogeneity in lung cancer tissues may lead to tumor progression and drug resistance (). In our study, the enrichment of metabolic pathways between the LCA subtypes also showed heterogeneity; LCA1 was enriched in the amino acid, lipid, and vitamin pathways, whereas LCA2 was mainly enriched in glycosylation. Glycosylation has been shown to promote tumor invasion (), suggesting that the therapies targeting glycosylation-related metabolic pathways could be explored for LCA2. Meanwhile, the degree of metabolic pathway abnormality was significantly higher in the LCA2 group. Previous studies have shown that a high degree of metabolic abnormality usually leads to a shorter survival time or resistance to anti-tumor treatments, including chemotherapy, radiotherapy, targeted therapies, and immunotherapy ().
Although immunotherapy has been widely adopted for cancer treatment (50), progress in LCA immunotherapy has been limited. Some drugs, including cetuximab, pembrolizumab, and nivolumab, have been just used in combination therapy after surgery for advanced LCA (). The analysis of immune infiltration across the subtypes of LCA1 and LCA2 in our study revealed that the infiltration of immune cells, ESTIMATE scores, and certain immune functions (APC_co_inhibition, Parainflammation, T_cell_co_inhibition) were significantly higher in LCA2. In addition, multiple immune checkpoint genes (HAVCR2, LAIR1, PDCD1LG2, BTLA, CD274) were expressed at higher levels in LCA2, suggesting that immunotherapy targeting these inhibitory sites (51) may be effective in this group.
The two subtypes also displayed different patterns of DNA changes. Notably, as a transcriptional regulatory protein with histone methyltransferase activity, NSD1 exhibits a markedly higher mutation frequency in LCA1. Moreover, NSD1 is primarily involved in processes such as methylation, transcriptional regulation, DNA binding, nuclear receptor binding, and histone modification. As nonsense mutations predominated in NSD1 variants within LCA1, the expression of NSD1 appeared to be a significantly reduced in LCA1 (Supplementary Figure S4C). This may account for the reduced methylation level observed in LCA1. Similarly, a previous study indicated that inactivating mutations in NSD1 and NSD2 in LCA are associated with a favorable prognosis (52). Furthermore, inactivating mutations in NSD1 can inhibit tumor growth, decrease methylation, and reduce immune infiltration (53, 54). All these observations are aligned with our observations. Hence, NSD1 has a great potential as a key molecule in the classification, prognosis estimation, and targeted treatment of LCA, warranting further explorations.
Indeed, these two isoforms really showed significant differences in their methylation levels. Methylation plays a role in gene silencing, X-chromosome inactivation, and genome stability (55, 56). Moreover, immune infiltration and function are also closely related to methylation levels (57). Up to now, methylation inhibitors are used to treat certain hematological malignancies. These DNA methylation inhibitors can mediate alterations in immune cells functions associated with acquired immunity (58). Typically, cancers exhibit genome-wide hypomethylation and site-specific hypermethylation, and the interactions between genome-wide hypomethylation and site-specific hypermethylation have something to do with the epigenetic and metabolic reprogramming (59). In the present study, LCA1 exhibited lower levels of methylation, lower immune infiltration, and higher metabolic enrichment, which may be significantly associated with high-frequency nonsense mutations of NSD1. Furthermore, the results of the enrichment analysis of methylated sites also revealed the significant differences between the two subtypes: LCA2 methylated sites were enriched in DNA transcription and hormone-related pathways, particularly in insulin secretion. Hormone secretion and transport have been better explored in some hormone-dependent tumors, such as breast and prostate cancers, and hormone therapies have been employed in their treatment (60, 61). Meanwhile, insulin-related traits are also thought to be strongly associated with cancer development (62, 63). Therefore, the relationship between LCA and insulin-related features requires further investigation.
Patients with LCA often experience a poor quality of life due to the late-stage diagnosis (). Thus, there is an urgent need to develop non-invasive early screening methods to enhance the precision of medical strategies. Over the past several years, researchers have utilized machine learning algorithms to uncover the hidden relationships between multi-omics data and disease conditions, as well as to develop predictive models (64, 65). In our study, the diagnostic model we had established by machine learning with six signature genes could effectively distinguished tumor samples from normal samples. This model was validated using the independent dataset GSE142083 to ensure its generalizability. The results of the ROC analysis and DCA further confirmed the predictive superiority of the diagnostic model, demonstrating its high clinical utility. On the other hand, the prognostic model with 10 signature genes established by machine learning could effectively predict the prognosis of LCA and was validated using an independent dataset. Previous studies, such as those using GPT and SMS as prognostic risk models (AUC=0.748, 0.823, and 0.781 for 1-, 3-, and 5-year survival, respectively) (66), and a study using TMEM2, DACT1, STMN2, GPR173 as prognostic risk models (AUC=0.814, 0.859, and 0.782 for 1-, 3-, and 5-year survival, respectively) (67) also established risk models for LCA. Compared to their models, our model (AUC=0.802, 0.912, and 0.920 for 1-, 3-, and 5-year survival, respectively) showed relatively better predictive accuracy. In the future studies, clinical participants and the improved machine learning algorithms should be utilized to further validate and refine the two models for clinical transformation.
Single-cell sequencing analysis has opened a new era in the exploration of the tumor cell heterogeneity and the tumor microenvironment infiltration, and has been applied to explore cell heterogeneity, immune microenvironment, and drug resistance mechanisms of various types of malignancies (68, 69). In this study, we used scRNA-seq analyses to validate the expression of genes related to metabolic subtypes as well as diagnostic, and prognostic models. The “copycat” was used to distinguish malignant from non-malignant epithelial cells and to obtain CNV scores. Diagnostic scores were calculated for all epithelial-derived cells and were significantly higher for malignant epithelial-derived cells. Malignant epithelial cells were classified into LCA1 and LCA2 groups based on the results of the metabolic subtype-related gene analysis. Prognostic risk scores were calculated for malignant cells. It was found that significantly higher prognostic scores were observed in the LCA2 group. These findings were consistent with the results observed in the bulk RNA datasets, which further confirmed the robustness of the metabolic subtypes as well as the diagnostic and prognostic models.
This study had several limitations. First, the large number of subtype-determining genes, along with their associated high costs, poses a challenge for clinical applicability. Second, although we included a substantial number of patients from both the microarray and RNA-seq platforms, which suggests that our conclusions may be highly reliable and robust without platform bias, the outcomes still require further validation in a prospective study owing to the retrospective nature of the current research. Third, additional experimental studies and validation of the signatures of clinical specimens are necessary for future researches.
5 Conclusions
In conclusion, metabolic reprogramming significantly influences tumor growth, progression, and the tumor microenvironment. Despite growing attention to therapeutic approaches targeting tumor metabolism, the diversity and heterogeneity of tumor metabolism pose challenges to metabolic therapy. A comprehensive exploration of metabolic pathways and mechanisms across various tumors is essential to advance metabolic therapies. In this study, we subtyped LCA based on metabolic factors and examined immune infiltration, DNA mutations, methylation, and single-cell data, which laid the foundation for future subclassifications in diagnosis and treatment. Using machine learning, we developed metabolism-related diagnostic and prognostic models with excellent performances. Our work enhances the understanding of patients with LCA classifications, facilitating the development of LCA early detection and shedding light on the precision treatments of LCA.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. The data can be found at TCGA (https://portal.gdc.cancer.gov/, accession information: TCGA-HNSC) and GEO (https://www.ncbi.nlm.nih.gov/geo, accession numbers: GSE252490, GSE130605, GSE27020, GSE142083).
Author contributions
DZ: Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. XP: Data curation, Software, Writing – review & editing. XD: Data curation, Software, Writing – review & editing. CL: Data curation, Software, Writing – review & editing. SL: Conceptualization, Supervision, 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 Research Development Project of the Affiliated Hospital of North Sichuan Medical College (NO. 2024GC014).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1512502/full#supplementary-material
Abbreviations
KEGG, Kyoto Encyclopedia of Genes and Genomes; OS, Overall survival; DFS, Disease-free survival; LCA, Laryngeal cancer; LSCC, Laryngeal squamous cell carcinoma; ssGSEA, Single sample gene set enrichment analysis; ROC, operating characteristic curve; LASSO, Least absolute shrinkage and selection operator cox regression analysis; RF, Random Forest algorithm; TCGA, The Cancer Genome Atlas; GEO, Gene Expression Omnibus; PCA, Principal component analysis; CDF, Cumulative distribution function; AUC, Area under the curve; MAF, Mutation annotation format files; DEGs, Differentially expressed genes CNV, Copy number variation; DCA, Decision curve; GO, Gene Ontology; CIBERSORT, Cell-type Identification using Estimating Relative Subsets of RNA Transcripts; TMB, Tumor mutation burden; ESTIMATE, Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data.
References
1
SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660
2
HanJSumerBD. The changing demographics and treatment of larynx cancer. Ann Surg Oncol. (2021) 28:6927–8. doi: 10.1245/s10434-021-10657-z
3
NociniRMolteniGMattiuzziCLippiG. Updates on larynx cancer epidemiology. Chin J Cancer Res. (2020) 32:18–25. doi: 10.21147/j.issn.1000-9604.2020.01.03
4
BurtnessBHarringtonKJGreilRSoulièresDTaharaMde CastroGJr.et al. Pembrolizumab alone or with chemotherapy versus cetuximab with chemotherapy for recurrent or metastatic squamous cell carcinoma of the head and neck (KEYNOTE-048): a randomised, open-label, phase 3 study. Lancet. (2019) 394:1915–28. doi: 10.1016/S0140-6736(19)32591-7
5
FerrisRLBlumenscheinGJr.FayetteJGuigayJColevasADLicitraLet al. Nivolumab for recurrent squamous-cell carcinoma of the head and neck. N Engl J Med. (2016) 375:1856–67. doi: 10.1056/NEJMoa1602252
6
VermorkenJBMesiaRRiveraFRemenarEKaweckiARotteySet al. Platinum-based chemotherapy plus cetuximab in head and neck cancer. N Engl J Med. (2008) 359:1116–27. doi: 10.1056/NEJMoa0802656
7
KoroulakisAAgarwalM. Laryngeal cancer. In: StatPearls. StatPearls Publishing LLC, Treasure Island (FL (2024).
8
SteuerCEEl-DeiryMParksJRHigginsKASabaNF. An update on larynx cancer. CA Cancer J Clin. (2017) 67:31–50. doi: 10.3322/caac.21386
9
ZhuQZhangXLuFMiaoSZhangCLiuZet al. RUNX1-BMP2 promotes vasculogenic mimicry in laryngeal squamous cell carcinoma via activation of the PI3K-AKT signaling pathway. Cell Commun Signal. (2024) 22:227. doi: 10.1186/s12964-024-01605-x
10
LiangHZhangZGuanZZhengSLouJLiuWet al. Follicle-like tertiary lymphoid structures: A potential biomarker for prognosis and immunotherapy response in patients with laryngeal squamous cell carcinoma. Front Immunol. (2023) 14:1096220. doi: 10.3389/fimmu.2023.1096220
11
GongXYChenHBZhangLQChenDSLiWChenDHet al. NOTCH1 mutation associates with impaired immune response and decreased relapse-free survival in patients with resected T1-2N0 laryngeal cancer. Front Immunol. (2022) 13:920253. doi: 10.3389/fimmu.2022.920253
12
YuanYLiHPuWChenLGuoDJiangHet al. Cancer metabolism and tumor microenvironment: fostering each other? Sci China Life Sci. (2022) 65:236–79. doi: 10.1007/s11427-021-1999-2
13
WardPSThompsonCB. Metabolic reprogramming: a cancer hallmark even warburg did not anticipate. Cancer Cell. (2012) 21:297–308. doi: 10.1016/j.ccr.2012.02.014
14
SchmidtDRPatelRKirschDGLewisCAVander HeidenMGLocasaleJW. Metabolomics in cancer research and emerging applications in clinical oncology. CA Cancer J Clin. (2021) 71:333–58. doi: 10.3322/caac.21670
15
YuanQDengDPanCRenJWeiTWuZet al. Integration of transcriptomics, proteomics, and metabolomics data to reveal HER2-associated metabolic heterogeneity in gastric cancer with response to immunotherapy and neoadjuvant chemotherapy. Front Immunol. (2022) 13:951137. doi: 10.3389/fimmu.2022.951137
16
LiXDuYJiangWDongSLiWTangHet al. Integrated transcriptomics, proteomics and metabolomics-based analysis uncover TAM2-associated glycolysis and pyruvate metabolic remodeling in pancreatic cancer. Front Immunol. (2023) 14:1170223. doi: 10.3389/fimmu.2023.1170223
17
YangCHuangXLiuZQinWWangC. Metabolism-associated molecular classification of hepatocellular carcinoma. Mol Oncol. (2020) 14:896–913. doi: 10.1002/1878-0261.12639
18
ZhaoMLiW. Metabolism-associated molecular classification of uterine corpus endometrial carcinoma. Front Genet. (2023) 14:955466. doi: 10.3389/fgene.2023.955466
19
SunYChenSLuYXuZFuWYanW. Single-cell transcriptomic analyses of tumor microenvironment and molecular reprograming landscape of metastatic laryngeal squamous cell carcinoma. Commun Biol. (2024) 7:63. doi: 10.1038/s42003-024-05765-x
20
ZhengXGaoWZhangZXueXMijitiMGuoQet al. Identification of a seven-lncRNAs panel that serves as a prognosis predictor and contributes to the Malignant progression of laryngeal squamous cell carcinoma. Front Oncol. (2023) 13:1106249. doi: 10.3389/fonc.2023.1106249
21
FountzilasEKotoulaVAngouridakisNKarasmanisIWirtzRMEleftherakiAGet al. Identification and validation of a multigene predictor of recurrence in primary laryngeal cancer. PloS One. (2013) 8:e70429. doi: 10.1371/journal.pone.0070429
22
GaoWZhangYLuoHNiuMZhengXHuWet al. Targeting SKA3 suppresses the proliferation and chemoresistance of laryngeal squamous cell carcinoma via impairing PLK1-AKT axis-mediated glycolysis. Cell Death Dis. (2020) 11:919. doi: 10.1038/s41419-020-03104-6
23
PossematoRMarksKMShaulYDPacoldMEKimDBirsoyKet al. Functional genomics reveal that the serine synthesis pathway is essential in breast cancer. Nature. (2011) 476:346–50. doi: 10.1038/nature10350
24
WuFYinYYFanWHZhaiYYuMCWangDet al. Immunological profiles of human oligodendrogliomas define two distinct molecular subtypes. EBioMedicine. (2023) 87:104410. doi: 10.1016/j.ebiom.2022.104410
25
WilkersonMDHayesDN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. (2010) 26:1572–3. doi: 10.1093/bioinformatics/btq170
26
LoveMIHuberWAndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
27
HänzelmannSCasteloRGuinneyJ. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinf. (2013) 14:7. doi: 10.1186/1471-2105-14-7
28
YoshiharaKShahmoradgoliMMartínezEVegesnaRKimHTorres-GarciaWet al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. (2013) 4:2612. doi: 10.1038/ncomms3612
29
PostowMASidlowRHellmannMD. Immune-related adverse events associated with immune checkpoint blockade. N Engl J Med. (2018) 378:158–68. doi: 10.1056/NEJMra1703481
30
ChenLZhangYHWangSZhangYHuangTCaiYD. Prediction and analysis of essential genes using the enrichments of gene ontology and KEGG pathways. PloS One. (2017) 12:e0184129. doi: 10.1371/journal.pone.0184129
31
MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. (2018) 28:1747–56. doi: 10.1101/gr.239244.118
32
RitchieMEPhipsonBWuDHuYLawCWShiWet al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
33
SavargivMMasoumiBKeyvanpourMR. A new random forest algorithm based on learning automata. Comput Intell Neurosci. (2021) 2021:5572781. doi: 10.1155/2021/5572781
34
ParkSY. Nomogram: An analogue tool to deliver digital knowledge. J Thorac Cardiovasc Surg. (2018) 155:1793. doi: 10.1016/j.jtcvs.2017.12.107
35
Van CalsterBWynantsLVerbeekJFMVerbakelJYChristodoulouEVickersAJet al. Reporting and interpreting decision curve analysis: A guide for investigators. Eur Urol. (2018) 74:796–804. doi: 10.1016/j.eururo.2018.08.038
36
JiangMGuXXuYWangJ. Metabolism-associated molecular classification and prognosis signature of head and neck squamous cell carcinoma. Heliyon. (2024) 10:e27587. doi: 10.1016/j.heliyon.2024.e27587
37
LiuZDingMQiuPPanKGuoQ. Natural killer cell-related prognostic risk model predicts prognosis and treatment outcomes in triple-negative breast cancer. Front Immunol. (2023) 14:1200282. doi: 10.3389/fimmu.2023.1200282
38
JiangPGuSPanDFuJSahuAHuXet al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. (2018) 24:1550–8. doi: 10.1038/s41591-018-0136-1
39
PatelAPTiroshITrombettaJJShalekAKGillespieSMWakimotoHet al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science. (2014) 344:1396–401. doi: 10.1126/science.1254257
40
GaoRBaiSHendersonYCLinYSchalckAYanYet al. Delineating copy number and clonal substructure in human tumors from single-cell transcriptomes. Nat Biotechnol. (2021) 39:599–608. doi: 10.1038/s41587-020-00795-2
41
TrapnellCCacchiarelliDGrimsbyJPokharelPLiSMorseMet al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. (2014) 32:381–6. doi: 10.1038/nbt.2859
42
WeiWLiuYShenYYangTDongYHanZet al. In situ tissue profile of rat trigeminal nerve in trigeminal neuralgia using spatial transcriptome sequencing. Int J Surg. (2024) 110:1463–74. doi: 10.1097/JS9.0000000000001110
43
YangLChuZLiuMZouQLiJLiuQet al. Amino acid metabolism in immune cells: essential regulators of the effector functions, and promising opportunities to enhance cancer immunotherapy. J Hematol Oncol. (2023) 16:59. doi: 10.1186/s13045-023-01453-1
44
Martinez-OutschoornUEPeiris-PagésMPestellRGSotgiaFLisantiMP. Cancer metabolism: a therapeutic perspective. Nat Rev Clin Oncol. (2017) 14:11–31. doi: 10.1038/nrclinonc.2016.60
45
TorresanoLNuevo-TapiolesCSantacatterinaFCuezvaJM. Metabolic reprogramming and disease progression in cancer patients. Biochim Biophys Acta Mol Basis Dis. (2020) 1866:165721. doi: 10.1016/j.bbadis.2020.165721
46
FaubertBSolmonsonADeBerardinisRJ. Metabolic reprogramming and cancer progression. Science. (2020) 368(6487):eaaw5473. doi: 10.1126/science.aaw5473
47
HensleyCTFaubertBYuanQLev-CohainNJinEKimJet al. Metabolic heterogeneity in human lung tumors. Cell. (2016) 164:681–94. doi: 10.1016/j.cell.2015.12.034
48
TikhonovASmoldovskayaOFeyzkhanovaGKushlinskiiNRubinaA. Glycan-specific antibodies as potential cancer biomarkers: a focus on microarray applications. Clin Chem Lab Med. (2020) 58:1611–22. doi: 10.1515/cclm-2019-1161
49
ZhuLZhuXWuY. Effects of glucose metabolism, lipid metabolism, and glutamine metabolism on tumor microenvironment and clinical implications. Biomolecules. (2022) 12(4):580. doi: 10.3390/biom12040580
50
Del PaggioJC. Immunotherapy: Cancer immunotherapy and the value of cure. Nat Rev Clin Oncol. (2018) 15:268–70. doi: 10.1038/nrclinonc.2018.27
51
PanCLiuHRobinsESongWLiuDLiZet al. Next-generation immuno-oncology agents: current momentum shifts in cancer immunotherapy. J Hematol Oncol. (2020) 13:29. doi: 10.1186/s13045-020-00862-w
52
PeriSIzumchenkoESchubertADSlifkerMJRuthKSerebriiskiiIGet al. NSD1- and NSD2-damaging mutations define a subset of laryngeal tumors with favorable prognosis. Nat Commun. (2017) 8:1772. doi: 10.1038/s41467-017-01877-7
53
TopchuIBychkovIGurselDMakhovPBoumberY. NSD1 supports cell growth and regulates autophagy in HPV-negative head and neck squamous cell carcinoma. Cell Death Discov. (2024) 10:75. doi: 10.1038/s41420-024-01842-6
54
LiYGoldbergEMChenXXuXMcGuireJTLeuzziGet al. Histone methylation antagonism drives tumor immune evasion in squamous cell carcinomas. Mol Cell. (2022) 82:3901–18.e7. doi: 10.1016/j.molcel.2022.09.007
55
Morales-NebredaLMcLaffertyFSSingerBD. DNA methylation as a transcriptional regulator of the immune system. Transl Res. (2019) 204:1–18. doi: 10.1016/j.trsl.2018.08.001
56
MooreLDLeTFanG. DNA methylation and its basic function. Neuropsychopharmacology. (2013) 38:23–38. doi: 10.1038/npp.2012.112
57
ZhangMWFujiwaraKCheXZhengSZhengL. DNA methylation in the tumor microenvironment. J Zhejiang Univ Sci B. (2017) 18:365–72. doi: 10.1631/jzus.B1600579
58
JonesPAOhtaniHChakravarthyADe CarvalhoDD. Epigenetic therapy in immune-oncology. Nat Rev Cancer. (2019) 19:151–61. doi: 10.1038/s41568-019-0109-9
59
SunLZhangHGaoP. Metabolic reprogramming and epigenetic modifications on the path to cancer. Protein Cell. (2022) 13:877–919. doi: 10.1007/s13238-021-00846-7
60
MoisandAMadéryMBoyerTDomblidesCBlayeCLarmonierN. Hormone receptor signaling and breast cancer resistance to anti-tumor immunity. Int J Mol Sci. (2023) 24(20):15048. doi: 10.3390/ijms242015048
61
ZhangZKarthausWRLeeYSGaoVRWuCRussoJWet al. Tumor microenvironment-derived NRG1 promotes antiandrogen resistance in prostate cancer. Cancer Cell. (2020) 38:279–96.e9. doi: 10.1016/j.ccell.2020.06.005
62
BrentaGDi FermoF. Thyroid cancer and insulin resistance. Rev Endocr Metab Disord. (2024) 25:19–34. doi: 10.1007/s11154-023-09849-7
63
WangXSunJLiJCaiLChenQWangYet al. Bidirectional Mendelian randomization study of insulin-related traits and risk of ovarian cancer. Front Endocrinol (Lausanne). (2023) 14:1131767. doi: 10.3389/fendo.2023.1131767
64
WangSWangRHuDZhangCCaoPHuangJ. Machine learning reveals diverse cell death patterns in lung adenocarcinoma prognosis and therapy. NPJ Precis Oncol. (2024) 8:49. doi: 10.1038/s41698-024-00538-5
65
ChenYWangBZhaoYShaoXWangMMaFet al. Metabolomic machine learning predictor for diagnosis and prognosis of gastric cancer. Nat Commun. (2024) 15:1657. doi: 10.1038/s41467-024-46043-y
66
ShenYHuangQZhangYHsuehCYZhouL. A novel signature derived from metabolism-related genes GPT and SMS to predict prognosis of laryngeal squamous cell carcinoma. Cancer Cell Int. (2022) 22:226. doi: 10.1186/s12935-022-02647-2
67
LiCZhuYShiS. Effective prognostic risk model with cuproptosis-related genes in laryngeal cancer. Braz J Otorhinolaryngol. (2024) 90:101384. doi: 10.1016/j.bjorl.2023.101384
68
TanZChenXZuoJFuSWangHWangJ. Comprehensive analysis of scRNA-Seq and bulk RNA-Seq reveals dynamic changes in the tumor immune microenvironment of bladder cancer and establishes a prognostic model. J Transl Med. (2023) 21:223. doi: 10.1186/s12967-023-04056-z
69
MoncadaRBarkleyDWagnerFChiodinMDevlinJCBaronMet al. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol. (2020) 38:333–42. doi: 10.1038/s41587-019-0392-8
Summary
Keywords
laryngeal cancer, molecular subtype, metabolic profiles, prognosis model, diagnostic model, ScRNA-seq
Citation
Zheng D, Pu X, Deng X, Liu C and Li S (2025) Metabolic profiles in laryngeal cancer defined two distinct molecular subtypes with divergent prognoses. Front. Immunol. 16:1512502. doi: 10.3389/fimmu.2025.1512502
Received
16 October 2024
Accepted
29 April 2025
Published
22 May 2025
Volume
16 - 2025
Edited by
Khosrow Kashfi, City University of New York, United States
Reviewed by
Vishnu Modur, Medpace, United States
Abhijeet R Patil, Teva Pharmaceuticals, United States
Updates
Copyright
© 2025 Zheng, Pu, Deng, Liu and Li.
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: SiJun Li, 1006038473@qq.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.