Qihang Yuan
Weizhi Zhang3†- 1Department of General Surgery, First Affiliated Hospital of Dalian Medical University, Dalian, Liaoning, China
- 2Laboratory of Integrative Medicine, First Affiliated Hospital of Dalian Medical University, Dalian, Liaoning, China
- 3Dalian No.24 High School, Dalian, Liaoning, China
Colon adenocarcinoma (COAD) is the primary factor responsible for cancer-related mortalities in western countries, and its development and progression are affected by altered sphingolipid metabolism. The current study aimed at investigating the effects of sphingolipid metabolism-related (SLP) genes on multiple human cancers, especially on COAD. We obtained 1287 SLP genes from the GeneCard and MsigDb databases along with the public transcriptome data and the related clinical information. The univariate Cox regression analysis suggested that 26 SLP genes were substantially related to the prognosis of COAD, and a majority of SLP genes served as the risk genes for the tumor, insinuating a potential pathogenic effect of SLP in COAD development. Pan-cancer characterization of SLP genes summarized their expression traits, mutation traits, and methylation levels. Subsequently, we focused on the thorough research of COAD. With the help of unsupervised clustering, 1008 COAD patients were successfully divided into two distinct subtypes (C1 and C2). C1 subtype is characterized by a poor prognosis, activation of SLP pathways, high expression of SLP genes, disordered carcinogenic pathways, and immune microenvironment. Based on the clusters of SLP, we developed and validated a novel prognostic model, consisting of ANO1, C2CD4A, EEF1A2, GRP, HEYL, IGF1, LAMA2, LSAMP, RBP1, and TCEAL2, to quantitatively evaluate the clinical outcomes of COAD. The Kaplain-Meier survival curves and ROC curves highlighted the accuracy of our SLP model in both internal and external cohorts. Compared to normal colon tissues, expression of C2CD4A was detected to be significantly higher in COAD; whereas, expression levels of EEF1A2, IGF1, and TCEAL2 were detected to be significantly lower in COAD. Overall, our research emphasized the pathogenic role of SLP in COAD and found that targeting SLP might help improve the clinical outcomes of COAD. The risk model based on SLP metabolism provided a new horizon for prognosis assessment and customized patient intervention.
Introduction
Colon adenocarcinoma (COAD) is a type of malignant tumor originating from the colon gland epithelium and the most prevalent pathological form of colon tumor. Colorectal cancer is the fourth most prevailing malignant tumor in China, and its prevalence is increasing (1). The clinical manifestations of COAD differ based on the tumor’s location and stage. Due to the unique nature of the signs of early colon cancer, some individuals are diagnosed at an advanced stage. Currently, the majority of patients with colon cancer may benefit from a comprehensive treatment regimen that includes extensive surgery, radiation, and chemotherapy (2). However, the unfavorable effects of surgery, radiation, and chemotherapy will also contribute to a reduction in the quality of life of individuals suffering from colon cancer, and the prognosis for those who cannot withstand surgical treatment, radiotherapy, and chemotherapy remains grim. In recent years, the in-depth investigation of tumor metabolic reprogramming and immunological microenvironment has provided evidence that both metabolic and immune heterogeneity plays a role in the onset and progression of colon cancer (3). Patients having diverse metabolic features may have a distinct immune microenvironment and prognosis. Accurate identification of the metabolic heterogeneity of colon cancer is a crucial problem for the advancement of precise oncology therapy.
According to previous biological studies on sphingolipid metabolism and its function, ceramide, ceramide-1-phosphate, glucosylceramide, lactose ceramide, galactosylceramide, sphingosine, sphingosine galactoside, and 1-phosphate-sphingosine are not only inactive precursors of sphingomyelin metabolism but also crucial effector molecules in cell signal transduction (4). It has been shown that sphingolipids and their metabolites play a role in several crucial processes of signal transduction, such as cell growth, differentiation, senescence, and programmed cell death, resulting in a wide range of cellular biological activities (5). Since ceramide was discovered to be the second messenger of lipids to promote apoptosis, multiple studies have demonstrated that aberrant sphingolipid metabolism is strongly linked to the onset and progression of cancer; however, the exact mechanism remains unclear (6).
This study aimed at exploring the possible role of sphingolipid metabolism in pan-cancer, particularly colon cancer, and to provide a new scheme for the prognosis assessment as well as tailored patient therapy. In addition, the influence of sphingolipid metabolism on the tumor immune milieu was investigated further, laying the groundwork for the investigation of the link between metabolic reprogramming and the tumor immune microenvironment. Considering the substantial impact of sphingolipid metabolism on the survival and treatment options of COAD patients, we developed a unique and reliable prognostic discriminant model based on the sphingolipid metabolism classification technique.
Material and methods
Acquisition of COAD datasets and SLP genes
The gene expression profiles and relevant clinicopathological data from COAD patients were retrieved from the TCGA (https://portal.gdc.cancer.gov/) and GEO (https://www.ncbi.nlm.nih.gov/) datasets. The TCGA dataset consisting of 473 COAD samples was used as the training dataset, while the GSE39582 cohort containing 566 COAD samples was used as the validation dataset. The “sva” package (7) in R was implemented for generating the relative gene expression matrices from the TCGA and GSE39582 cohorts. Until now, over 1,200 genes related to sphingolipid metabolism have been discovered on the “GeneCards (8)” and “MSigDB (9)” platforms. These SLP genes were selected for further expression profile analyses.
Identification of prognostic SLP genes
Univariate Cox regression analysis adjusted by the Benjamini & Hochberg (BH) technique was implemented to screen the prognostic SLP genes in both the TCGA and GSE39582 datasets (10). After the intersection, SLP genes with potential prognostic values in both datasets were obtained. In addition to the univariate Cox regression analysis, KM survival curves were also plotted for the purpose of verifying the prognostic performances of the above 26 SLP genes with the help of “survival” and “survminer” packages in R (11). Finally, the “cor” function in R was applied to explore the co-expression relationship between these 26 prognostic SLP genes in both the TCGA and GSE39582 datasets.
The pan-cancer landscape of SLP genes
The raw data of pan-cancer comprising the mRNA expression profiles, clinical information, single nucleotide variation (SNV) data, copy number variation (CNV) data, and methylation data were gathered from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/). The “limma” package in R was utilized for analyzing the differential expression of SLP genes in pan-cancer-related tumor tissues and paracancerous tissues.
TCGA database was used to collect the SNV data across 33 cancer types. The collected data include variant-type values: Missense_Mutation, Silent, 5’ Flank, 3’ UTR, RNA, In_Frame_Del, Nonsense_Mutation, SLPice_Site, Intron, 5’ UTR, In_Frame_Ins, Frame_Shift_Del, Nonstop_Mutation, 3’ Flank, Frame_Shift_Ins, and Translation_Start_Site. The Silent, Intron, IGR, 3’ UTR, 5’ UTR, 3’ Flank, and 5’ Flank were extracted for calculating the SNV percentage. SNV mutation frequency (percentage) of the coding region of each gene was computed using the formula: Number of Mutated Samples/Number of Cancer Samples. An SNV oncoplot plot was produced by maftools (12). Moreover, each gene’s CNV was elevated by incorporating the heterozygosity and homozygosity of amplification and deletion, in which over 5% was considered to be a high-frequency CNV. The R language was applied to visualize the results of SNV and CNV of SLP genes in pan-cancer. The methylation probe for the individual gene’s promoter was annotated using the R package “IlluminaHumanMethylation-450kanno.ilmn12.hg19” from Bioconductor. The Wilcoxon signed rank test was utilized for examining the differential methylation between all genes in the tumor as well as healthy samples, and a P-value cutoff of 0.05 was used for identifying the genes that were substantially hypo-methylated or hyper-methylated. Furthermore, single sample gene set enrichment analysis (ssGSEA) was employed for computing SLP scores in individual tumor samples so as to identify the pathways related to sphingolipid metabolism in several human malignancies (13). In addition, based on the SLP scores, the samples of every tumor type were classified into two groups, comprising the top 30% and the bottom 30%. Gene set enrichment analysis (GSEA) was then carried out similarly to our previous studies (14, 15).
SLP-based cluster analysis
A total of 1008 COAD patients were clustered into two subgroups following the SLP-related gene expression profiles using the R software package “ConsensusClusterPlus” (16). The highest number of clusters was 9, and 80% of the total samples were drawn 50 times, clusterAlg = “km”, distance = “euclidean”. With the aid of R software packages “survival” and “survminer”, the Kaplan-Meier survival analysis was then performed for the purpose of comparing the survival difference between the two groups. Moreover, the Wilcoxon test was conducted to observe the expression distribution of the SLP genes in both subtypes related to SLP (17). Furthermore, 9 typical SLP-associated pathways and 50 typical cancer-associated pathways were downloaded from the MSigDb platform (9). The “GSVA” package and the “Wilcox.test” function in R were employed to assess and compare the activities of the aforementioned pathways between the two subtypes (18).
Cluster-based analysis of tumor immune microenvironment
Sphingolipids are implicated in the interactions between cancer cells and immune system, as shown by a rising number of studies. In order to characterize the discrepancy in the immune response between C1 and C2 subtypes at the macroscopic level, the R program “estimate” was first used for calculating the immune score, stromal score, tumor purity, and estimate the score of individual COAD samples (19). The “estimate” algorithm was developed by Kosuke et al (20). Based on the expression data, the “estimate” method could estimate the matrix percentage (stromalscore) and immune score (immunescore) of tumour samples, which might be used to indicate the presence of matrix and immune cells. Adding the two fractions yields the estimatescore, which might be used to estimate the tumor’s purity. Subsequently, the R package “ggpubr” was utilized for visualizing the results. In order to uncover the infiltration abundance of each type of immunocytes, the CIBERSORT, MCPcounter, QUANTISEQ, XCELL, CIBERSORT-ABS, TIMER, and EPIC algorithms were employed to assess the immunological characteristics of both subtypes (21). Above 7 immunological algorithms were conducted with the help of TIMER2.0 platform, which was a public website using the immunedeconv method to assess the abundance of varied immunocyte infiltration. Subsequently, a heatmap was plotted to visualize the results at the microscopic level. Lastly, we studied the expression of common immune checkpoint genes (ICGs) across various subtypes for predicting the efficacy of immune checkpoint blockade therapy (22).
Cluster-based prediction of drug sensitivity
For estimating the sensitivity of chemotherapeutic agents, the R package pRRophetic was employed for determining the half-maximal inhibitory concentration (IC50) of samples in various groups via ridge regression (23). The pRRophetic package is a program developed by Prof. Paul Geeleher of the University of Minnesota (23, 24). The basic principle of this algorithm is based on the expression profile of the cell line and the corresponding IC50 information, and then establish a model through ridge regression, and then use this model to predict the chemotherapeutic response of clinical samples. Subsequently, the IC50 values across various subtypes were compared using the Wilcoxon test. A lower IC50 value indicated a better response to chemotherapy.
Cluster-based annotation of differentially expressed genes
The “limma” R package was utilized for screening DEGs across several clusters (|logFC| > 0.585, FDR < 0.05) (25). Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were carried out for selecting and visualizing the significantly enriched GO terms and KEGG pathways in DEGs (26–28).
DEGs-based development and verification of a novel SLP-associated prognostic panel
The 446 samples of COAD in the TCGA data set were sorted into training and validation cohorts. All the samples were placed back into the random grouping 100 times in advance to prevent the impact of random assignment bias on the stability of subsequent modeling, and the group sampling was done with a training cohort-to-validation cohort ratio of 1:1. Moreover, 224 samples were present in the final training cohort and 222 samples in the final internal validation cohort 1 (i.e. test1). All the 446 samples of COAD in the TCGA data set were rated as the internal validation cohort 2 (i.e., test2). Moreover, 562 COAD samples in the GSE39582 data set were rated as the external validation cohort (i.e., test 3). For DEGs between different molecular subtypes, the univariate Cox regression analysis was carried out using the coxph function of the survival package in R, and a P-value of < 0.05 was taken as the filtering threshold. Moreover, Lasso regression using R package glmnet was carried out for the purpose of eliminating the redundancy of prognostic genes for subsequently developing the prognostic model (29). Finally, multivariate Cox regression analysis and the “predict” function in R were employed for computing the risk scores of every sample of COAD. The samples in the training, test1, test2, and test3 were all allocated into high and low-risk subpopulations following the median risk score of the training cohort. KM survival curves and ROC curves were drawn for evaluating as well as calculating the prognostic values of SLP-APP in the aforementioned cohorts.
Differential expression analysis of SLP-APP genes
A multidimensional cancer genomics dataset called GEPIA incorporated mass data from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression project (GTEx) (30, 31). Boxplot could be employed to determine the expression level of a single gene in various cancer types. We identified the expression levels of SLP-APP genes in COAD as well as healthy colon tissues based on TCGA and GTEx data. Human Protein Atlas (HPA, https://www.proteinatlas.org/) was used for determining the protein expression in different tissues and organs of humans from the RNA and protein levels by using transcriptomics and proteomics techniques (32, 33). We identified the SLP-APP protein expression in both COAD and healthy colon tissues based on HPA data.
Functional analysis of SLP metabolism-associated genes in COAD based on the knock-down of HIF1α gene
On the GEO platform, we discovered knock-down sequencing data for the HIF1a gene after a search of more than 1200 sphingolipid metabolism genes previously compiled. Consequently, we investigated the function and possible mechanism of HIF1α in COAD. GSE155104 is an open cell line knock-down sequencing data set submitted by Prof.Glaus Garzon JF et al (34). The authors cultivated the mouse colon cancer cell line MC38 and silenced the HIF1α gene. The transcriptomes of control cells and knock-down cell lines were sequenced. We gathered and compiled these sequencing data, examined the differentially expressed genes of the two cell lines using the limma programme, and assessed the enrichment of the GO and KEGG ontologies.
Results
Identification of prognosis-related SLP genes
Using the “GeneCards” and “MSigDB” databases, 1287 SLP genes were detected (Table S1). It was found that 104 SLP genes were significantly linked to COAD prognosis using univariate Cox analysis in the TCGA cohort (Table S2). Similarly, 164 SLP genes were also significantly associated with the prognosis of COAD according to the univariate Cox analysis in the GSE39582 cohort (Table S3). After the intersection, 26 SLP genes were finally preserved as reliable SLP genes related to prognosis (Figure 1A). KM survival curves of 26 SLP genes were constructed to verify their prognostic performances (Figure 1B). The co-expression relationship of 26 genes was shown in Figure 1C. It was found that KCNE4 and CAVIN1 were positively correlated with most genes, such as NRP1 and S1PR3. These 26 SLP genes included LGALS4, FDFT1, FLOT1, RAB7A, CAVIN1, BIRC5, AGPAT1, TNFRSF11A, AGRN, SERPINE1, NRP1, MAPK12, S1PR3, FOXC1, ATP10A, NOS2, NOTCH4, MAPK11, IDUA, FABP4, KCNE4, ALOX12B, ISM1, NGFR, ALPP, PAQR9. As shown in Table 1 (35–56), the majority of genes may operate as proto-oncogenes and affect the cell proliferation, invasion, migration, and metastasis of COAD patients.
Figure 1 Prognostic performances of SLP genes in colon adenocarcinoma. (A) The results of univariate Cox regression analysis: Green represents 104 genes obtained by TCGA screening, red represents 164 genes obtained by geo screening, and 26 SLP genes obtained by intersection; (B) KM survival curves of 26 SLP genes; (C) Co-expression relationship between these 26 SLP genes in both TCGA and GEO cohorts.
Table 1 Contributions of 26 sphingolipid metabolism-related genes in COAD.
The pan-cancer landscape of SLP genes
To show the pan-cancer overview of the above-mentioned 26 SLP genes and reveal their potential biological function in cancers, we comprehensively investigated the mRNA expression traits, mutation landscapes, and methylation levels of these genes in multiple human cancers. Our results showed that BIRC5 was up-regulated, whereas FABP4 was down-regulated in most cancers (Figure 2A). Increased expression of BIRC5, SERPINE1, and PAQR9 was detected in COAD compared to that in the paracancerous tissues. Decreased expression of LGALS4, TNFRSF11A, FABP4, and NGFR was detected in COAD compared to that in the paracancerous tissues (Figure 2A). As depicted in Figure 2B, the CNVs of RAB7A, BIRC5, SERPINE1, and FABP4 were relatively high in most cancers. Significant SNV mutation of SLP genes (especially for NOTCH4 and ATP10A) was observed in COAD, LUAD, SKCM, and UCEC (Figures 2C, D). High methylation levels of PAQR9 were observed in most cancers, especially in the COAD (Figure 2E). Nevertheless, NOS2 exhibited a relatively low methylation level in most cancers. More importantly, these SLP genes were significantly linked to many typical immune-related pathways, including PD-1, IL-10, and chemokine signaling pathways (Figure 2F). Overall, these findings suggested a close association between SLP metabolism and tumor immune microenvironment.
Figure 2 Pan-cancer overview of SLP genes. (A) Expression levels of 26 SLP genes in pan-cancer: Orange represents upregulation, green represents downregulation, and the circle size shows the degree of difference; (B) Copy number variation status of 26 SLP genes in pan-cancer; (C) Single nucleotide variant frequency of 26 SLP genes in pan-cancer; (D) Single nucleotide variant types of 26 SLP genes in pan-cancer; (E) The methylation levels of 26 SLP genes in pan-cancer; (F) Relationship of 26 SLP genes with immune pathways of pan-cancer.
SLP-based cluster analysis
The 1008 samples of COAD could be classified into various clusters by following the consensus of mRNA expression of SLP genes. When the value of the clustering index “k” increased from 2 to 9, k = 2 was shown as the optimal point for attaining the largest differences (variations) between clusters (Figures 3A, B). Moreover, the interference between clusters was the least when k = 2 (Figures 3C, D). The COAD cohort was then sorted into cluster 1 (n = 255) and cluster 2 (n = 753) (Figure 3D). KM survival curves of the clusters indicated that C2-like COAD patients had higher rates of survival when compared to the C1-like COAD patients (p < 0.001, Figure 3E). The expression levels and pathway activities of SLP metabolism were substantially different between the two subtypes (Figure 3F). Most SLP genes were up-regulated in the C1 subtype. Meanwhile, the activities of most SLP pathways were also significantly increased in the C1 subtype. These findings indicated the activation of SLP metabolism might promote the progression and degree of COAD malignancy, thus contributing to the poor prognosis. Furthermore, the activities of typical cancer-related pathways exhibited remarkable differences between the C1 and C2 subtypes (Figure 3F). Specifically, the NOTCH signaling pathway, TGF-β signaling pathway, apoptosis, angiogenesis, and hypoxia pathways were significantly activated in the C1 subtype, which also might be responsible for the poor prognosis of the C1 subtype.
Figure 3 Consensus clustering and molecular characterization of colon adenocarcinoma. (A, B) Results of hierarchical clustering analysis show that the difference between clusters is the largest when k = 2; (C) When k = 2, the interference between clusters is minimal; (D) The cohort was classified into C1 and C2 clusters; (E) KM survival curve of two clusters: Purple represents C1 and green represents C2; (F) Heat map results of 9 typical SLP-related pathway activation, 26 SLP gene expression and 50 typical cancer-related pathway activation between the two subtypes. ** indicated p < 0.01, *** indicated p < 0.001, **** indicated p < 0.0001.
Comparison of the immune landscape between two clusters
First, the R “estimate” package was employed for investigating the immune-related scores between the two subtypes, and these scores including stromal, immune, and estimate scores were substantially higher in the C1 subgroup (Figure 4A). To elucidate the infiltration distribution of each type of immunocytes between C1 and C2 subtypes, multiple deconvolution algorithms, including TIMER, CIBERSOFT, MCPCOUNTER QUANTISEQ, and EPIC, were performed. Similarly, the results demonstrated a higher abundance of B cells, CD4+ T cells, CD8+ T cells, and macrophages in C1 than in C2 (Figure 4B). Moreover, the expression levels of common genes of immune checkpoint were then compared between both subgroups; and the results showed the up-regulation of most of these genes (including PDCD1LG2, TIGIT, TNFSF4, LAG3, CD86, CD40, and CD48) in C1 (Figure 4C).
Figure 4 Cluster-based tumor immune microenvironment analysis. (A) ESTIMATE algorithm-derived immune-related scores between the two clusters, including stromal score, immune score, tumor purity, and estimate score; (B) The distribution of infiltration of the cells involved in the immune system between the clusters using a variety of immune algorithms, such as TIMER, CIBERSOFT, CIBERSOFT-ABS, QUANTISEQ, MCPCOUNTER, XCELL, and EPIC; (C) The common immune checkpoint genes’ expression levels in both the subgroups. * indicated p < 0.05, ** indicated p < 0.01, *** indicated p < 0.001, ns indicated not significance.
Comparison of the drug sensitivity between the subgroups
While taking into account the significant contribution of molecularly targeted therapy to the improvement of prognosis among individuals suffering from COAD, we evaluated the expression profile characteristics of various clusters by using the R package “pRRophetic” so as to find sensitive targeted therapeutic agents for C1 and C2 subtypes. The findings of our study revealed that the C1 subtype might gain benefit from Bexarotene, Dasatinib, DMOG, Embelin, Imatinib, Nilotinib, Pazopanib, and Sunitinib; however, the C2 subtype might benefit from AKT.inhibitor.VIII, BIBW2992, BIRB.0796, GW.441756, Metformin, Methotrexate, Pyrimethamine, and Sorafenib (Figure 5).
Figure 5 Cluster-based targeted-drug sensitivity analysis.
Identification of differentially expressed genes and functional enrichment analysis
To explore the variations of biological functions between the clusters, 1381 DEGs between these clusters were detected. GO analyses indicated that DEGs were primarily enriched in extracellular matrix structural constituents, glycosaminoglycan binding, collagen binding, integrin binding, growth factor binding, and regulation of cell adhesion (Figure 6A). KEGG analysis demonstrated a primary involvement of DEGs in complement as well as coagulation cascades, focal adhesion, ECM-receptor interaction, as well as in chemokine and PI3K-Akt signaling pathways (Figure 6B).
Figure 6 Functional enrichment analysis of DEGs between the subgroups. (A) The gene ontology enrichment results; (B) The KEGG enrichment results.
Identification and verification of SLP-associated prognostic panel
446 and 562 COAD samples retrieved from the TCGA database and the GEO database, respectively (all containing entire mRNA expression and survival data) were utilized for sorting the samples into various cohorts. The current study consisted of four cohorts: the training cohort (224 TCGA samples, making up 50% of the total), test1 cohort (222 TCGA samples, making up 50% of the total), test2 (all 446 TCGA samples), and test3 cohort (562 GEO samples). Particularly during the validation of SLP-APP, the test1 and test2 cohorts were employed for the internal, whereas, the test3 cohort was employed for the external validation.
The risk score model was initially constructed as per the training cohort. Univariate Cox regression analysis carried out on every DEG indicated that 55 of the 1381 DEGs were potential prognostic indicators (Tables S4, S5). Subsequently, LASSO regression analysis was carried out for eliminating collinearity among the 55 candidate genes and to prevent over-fitting of the prognostic model (Figures 7A, B). Finally, ten genes (i.e., TCEAL2, GRP, EEF1A2, LSAMP, C2CD4A, RBP1, ANO1, HEYL, IGF1, and LAMA2) were utilized in a multivariate Cox proportional hazards regression analysis for establishing a novel SLP-APP (Figure 7C). The risk scores on the basis of the risk model were calculated by the “predict” function in R, and the samples were categorized into high- and low-risk groups having the median risk score as the cut-off value. The levels of the 10 signature genes were then examined in the two groups (Figures 8A–D). The samples’ risk scores were ranked from low to high, and a scatter plot was used to visualize the survival status. According to the findings, an increase in risk scores led to a rise in patient fatalities and a decrease in their duration of survival (Figures 8E–H). The survival curves of training, test1, test2, and test3 cohorts displayed a remarkably lower length of survival in the group with a high risk as compared to the group with a low risk (Figures 9A–D). Moreover, the ROC curve indicated that, in the training dataset, AUC values for one-, two- and three-year OS were 0.818, 0.878, and 0.854, respectively (Figure 9E). In the test1 dataset, AUC values for one-, two- and three-year OS were 0.646, 0.660, and 0.603, respectively (Figure 9F). In the test2 dataset, AUC values for one-, two- and three-year OS were 0.730, 0.784, and 0.749, respectively (Figure 9G). In the test3 dataset, AUC values for one-, two- and three-year OS were 0.615, 0.641, and 0.607, respectively (Figure 9H). Overall, these findings verified that the 10-gene signature showed a satisfactory prognostic prediction effect.
Figure 7 The process of LASSO-Cox regression analysis. (A, B) Lasso regression was performed to filter the candidate genes. (C) Multivariate Cox regression analysis aids in developing the prognostic model. * indicated p < 0.05, ** indicated p < 0.01, *** indicated p < 0.001
Figure 8 Distribution of model gene expression and patients’ survival between the subgroups with high and low risk. The model gene’s expression levels between subpopulations with high as well as low risk in (A) train, (B) test1, (C) test2, and (D) test3 cohorts. The survival status of colon adenocarcinoma between subpopulations of high and low risk in (E) train, (F) test1, (G) test2, and (H) test3 cohorts.
Figure 9 Assessment of the predictive ability of the SLP cluster-based prognostic model. Survival curve plots in (A) train, (B) test1, (C) test2, and (D) test3 cohorts. Receiver operating characteristic curves in (E) train, (F) test1, (G) test2, and (H) test3 cohorts.
Differential expression analysis of SLP-APP genes
The GEPIA2 database was utilized for comparing the mRNA expression levels of SLP-APP genes between cancer and paracancerous tissues. In the TCGA cohort, the transcriptomic levels of C2CD4A and GRP were substantially higher in COAD samples than in paracancerous tissues; nevertheless, the transcriptomic levels of EEF1A2, IGF1, and TCEAL2 were substantially lower in COAD samples than in paracancerous tissues (Figure 10A). In the TCGA and GTEx cohorts, the transcriptomic level of C2CD4A was substantially higher in COAD samples as compared to the paracancerous tissues; however, the transcriptomic levels of EEF1A2, IGF1, LAMA2, LSAMP, and TCEAL2 were significantly lower in COAD samples than paracancerous tissues (Figure 10B). The protein expression levels and cell location of SLP-APP genes in both cancer and adjacent tissues were observed using the HPA database (Figure 11).
Figure 10 mRNA expression levels of model genes between colon adenocarcinoma and paracancerous/normal tissues using (A) TCGA cohort and (B) TCGA+GTEx cohorts. * indicated p < 0.05.
Figure 11 Protein expression levels and cellular location of model genes between colon adenocarcinoma and paracancerous/normal tissues using the HPA database.
Functional analysis of SLP metabolism-associated genes in COAD based on the knock-down of HIF1α gene
We obtained the raw data of transcriptome sequencing after knocking down HIF1α gene, including 16941 gene expression of 3 paired MC38 cell lines. Compared with control cell lines, knock-down of HIF1α gene in MC38 cell lines showed extremely distinct molecular characteristics. A total of 1074 DEGs were identified, with 660 up-regulated genes and 414 down-regulated genes (Figure 12A). The top 20 DEGs were shown in the heatmap form (Figure 12B). GO and KEGG analysis indicated that these DEGs were mainly enriched in the following terms: regulation of oxygen levels, NAD metabolic process, pyruvate metabolic process, cell−cell junction, GTPase regulator activity, MAPK signaling pathway, HIF−1 signaling pathway, FoxO signaling pathway, and glycolysis (Figures 12C, D). Overall, these findings validated the regulatory role of HIF1α in metabolism reprogramming of COAD.
Figure 12 Functional validation of HIF1a gene in MC38 cell lines of mouse colon adenocarcinoma. (A) The volcano plot represents the differentially expressed genes between HIF1a-KD (knock down) and control cell lines. (B) The top 20 differentially expressed genes were shown in the heatmap. (C) GO anlysis of differentially expressed genes. (D) KEGG analysis of differentially expressed genes.
Discussion
Colon adenocarcinoma is the primary cause of cancer-related mortalities in humans, and its onset and progression are influenced by the presence and activity of sphingolipid components in tumors (57, 58). According to reports, sphingolipid is a bioactive molecule that plays a crucial regulatory function in the majority of body cells, including cell development control, proliferation, and programmed cell death, as well as the regulation of exosome production (59, 60). As demonstrated in the studies, the majority of the enzymes engaged in sphingolipid metabolism are responsible for the development of colon cancer (61, 62). Accumulating evidence has suggested that ceramide and sphingosine 1-phosphate (S1P) served as the crucial contributor in the occurrence and progression of COAD. The process of SLP metabolism can be roughly divided into the following steps: sphingomyelin-ceramides-sphinganine-S1P-hexadecenal. Some SLP-metabolizing enzymes have been reported to be dysregulated in human colon cancers (63). Majority of the enzymes have the potential to increase the ratio of sphingosine-1-phosphate (S1P) to ceramide, which in turn promote colon cancer cell survival, proliferation and cancer progression (63, 64). Colon cancer is associated with a reduction in alkaline-sphingomyelinase activity, which in turn decreases the hydrolysis of sphingomyelin and production of ceramide (65). The absence of ceramide in colon cancer may contribute to its growth because of its ability to suppress the proliferation of tumor cells (66). In addition, the phosphorylation of sphingosine by sphingosine kinases 1 (SK1) results in the formation of S1P, which promotes the growth of colon cancer cells. Studies on human colon cancer specimens have shown that the expression level of SK1 was significantly higher in COAD than normal mucosa (67). During metastasis, the expression of SK1 steadily rises (68). Collectively, disorders of enzyme activities contribute to the aberrant expression of substrates and products of SLP metabolism, and further promote the occurrence and progression of COAD. In addition, the aberrant gene expressions involved in sphingolipid metabolism may contribute to colon cancer growth and influence the treatment response of colon cancers as well. Additional research into SLP metabolism may offer a new therapeutic target for colon cancer (58).
Animal research and clinical observations have shown that the SLP metabolism of colon cancer cells has altered dramatically (66). To determine the extent to which SLP metabolism influences the prognosis of COAD patients and which genes influence the course of COAD, further research is required. In this work, we began by identifying the genes of SLP metabolism that exhibited a strong link to the prognosis of individuals suffering from COAD. Then, a pan-cancer study of these genes was conducted to highlight the significant impact of genes associated with SLP metabolism on a range of human malignancies. Based on the SLP metabolic genes, an unsupervised cluster analysis was carried out, and a prognostic risk score model was developed for estimating the probability of survival among patients with COAD.
Studies have shown that sphingolipids and enzymes that comprise the S1P pathway influence the development of ovarian cancer (69). Moreover, SLP metabolism affects the amount of estrogen in the body, thereby influencing the onset and progression of prostate cancer; however, the precise mechanism has not been elucidated. In our investigation, we chose sphingolipid metabolic genes for pan-cancer analysis, and the findings indicated that the BIRC5 protein was expressed at a high level in most of the malignancies, consistent with earlier research (38, 70, 71). BIRC5 is a tiny, multi-subtype protein whose role is to suppress apoptosis and promote cell growth (72). Additionally, reduced NGFR expression was observed in COAD. NGFR, a member of the tumor necrosis factor receptor family, is a multifunctional cell surface receptor that stimulates cell survival and differentiation during neural development, among other tasks (73). Yang and colleagues have shown that the NGFR expression is mainly down-regulated in colorectal cancer (55). NGFR inhibits tumor cell proliferation, migration, and invasion by inducing apoptosis and G1 phase arrest, which is useful in inhibiting the onset and advancement of tumors (55). More importantly, through the activation of S100A9, NGFR enhances the chemosensitivity of colorectal cancer cells by promoting the apoptotic and autophagic effects of 5-fluorouracil. Consistent with our result, the deletion of NGFR protein in colorectal cancer is linked to a poor OS in people suffering from colorectal cancer, which is an independent predictor of the prognosis of colorectal cancer (74).
Sphingolipids play a crucial role in the plasma membrane and have been shown to modulate the activity of surface receptors on immune cell surfaces (75). These metabolites have a role in the secretion of bioactive mediators such as S1P and ceramides, which control the crucial pathways necessary for the activation of immune cells and affect lymphocyte efflux and migration into the tumour microenvironment (TME) (76). Ceramide promoted the tumor-associated macrophages (TAMs) to differentiate into M1-like macrophages through activation of the protein kinase C pathway (77). Following this, M1-like macrophages induced tumor death via secreting IL-6, IL-10, IL-12, and TNF-α (75). Ceramide-mediated alterations of the rate of T helper 1 cell (Th1) and Th2 in host-TME led to TME skewed toward proinflammatory milieu, and then inhibit tumor growth (75). Moreover, C2-ceramide exerts its anti-tumor activity by increasing the percentage of CD8+T cells and producing perforin and granzyme B (77). However, S1P enhanced the expression of Bcl-2 in macrophages and encouraged the polarization of M2-like macrophages, which further facilitated tumor evasion (78, 79). Collectively, SLP metabolism is intimately associated with tumour immunity, and changing levels of sphingolipid metabolism are crucial determinants regulating TME. Similarly, our findings revealed that PD-1, IL-10, and chemokine signal pathways are among the several immune-related pathways that are substantially connected with these SLP genes. Overall, our findings indicate that individuals with colon adenocarcinoma have an aberrant SLP process and that the metabolic reprogramming of tumors is intimately tied to their immunological microenvironment. The pan-cancer description of the SLP genes provides a good framework for future research into additional malignancies.
Consequently, based on genes associated with sphingolipid metabolism, 1008 individuals with COAD were effectively classified into two metabolic subgroups (C1 and C2 subtypes). Patients with subtype C1 had an active sphingolipid metabolism pathway and relatively high levels of sphingolipid metabolism genes, but their prognoses were poor. These findings showed that the sphingolipid metabolism might have a detrimental role in the pathophysiology of colon adenocarcinoma, and that focused intervention of the sphingolipid metabolic pathway might improve the clinical outcomes of COAD patients. More importantly, we discovered that the NOTCH signaling pathway, TGF-β signaling pathway, apoptosis, angiogenesis, and hypoxia pathways were considerably active in the C1 subtype, which might be one of the causes of the subtype’s poor prognosis. The NOTCH signaling system is a highly conserved developmental pathway essential for apoptosis, tissue structure, and morphogenesis (80, 81). Increasing data suggest that this pathway is aberrant in several types of malignant tumors and may regulate oncogenes or tumor suppressors (82). Staudacher et al. revealed that the TGF-β signaling pathway is strongly linked to the incidence and progression of advanced colon cancer (83).
Considering that the process of metabolic reprogramming in tumors may be accompanied by an immune microenvironment imbalance, we investigated the possible variations in the immune microenvironment between C1 and C2 metabolic subtypes in depth. Intriguingly, the C1 subtype with active sphingolipid metabolism is coupled by strong infiltration of cells of the immune system and high expression levels of the immune checkpoints. Specifically, the infiltration abundance of B cells, CD4+T cells, CD8+T cells, and macrophages was greater in the C1 subgroup than that in the C2 subgroup, and immune checkpoint-associated genes (such as PDCD1LG2, TIGIT, TNFSF4, LAG3, CD86, CD40, and CD48) were up-regulated in C1 subgroup. Masugi et al. found that the PDCD1LG2 expression in colorectal cancer tumors was inversely connected with Crohn’s lymphoid response, indicating that the PDCD1LG2 might prevent the formation of tertiary lymphoid tissue in colorectal cancer (84). Furthermore, by studying the evolutionary genetic algorithms, Alderdice and colleagues found that IL2RB is a potential predictive biomarker of immune checkpoint therapy for colorectal cancer, which is consistent with our study (85). Immune checkpoint inhibitors are a new method of immunotherapy, which does not simply refer to improving the immunity of the body, but to improving the immune microenvironment around the tumor, thereby activating the activity of immune cells in the body to achieve the purpose of anti-tumor (86). The current immune checkpoint inhibitor treatment method has set off a boom in many tumors and has also achieved significant clinical progress, including lung cancer, lymphoma, melanoma, bowel cancer, and multiple other tumor treatments (87). The findings of this research provide a fresh approach to the future therapy of immunological checkpoints for colon cancer.
Finally, we assessed the ability of genes involved in sphingolipid metabolism to predict COAD survival rate. The findings demonstrated that patients with high risk had considerably lower survival rates than those with low risk. The AUC value of our prognostic model shows excellent diagnostic performance in the training set, internal verification set, and external verification set. In conclusion, the findings of this study revealed that the prognostic model based on sphingolipid metabolism might reliably predict the likelihood of survival for COAD patients. Our predictive diagnostic model comprises ten genes associated with sphingolipid metabolism, including TCEAL2, GRP, EEF1A2, LSAMP, C2CD4A, RBP1, ANO1, HEYL, IGF1, and LAMA2. Taglia et al. have reported that GRP/GRPR is responsible for promoting the binding of CD16+ and CD94+ natural killer cells following inducing expression of Hsp72 by signaling via focal adhesion kinase, thus contributing to tumor cell cytolysis (88). EEF1A2 served as an epithelial-mesenchymal transition-related gene that was found to be closely linked to the clinical outcomes of colon cancer, which was consistent with the outcomes of our study. Zhou et al. demonstrated that reduced expression of TCEAL2 may be associated with renal cell carcinoma carcinogenesis (89). Furthermore, Chang et al. demonstrated that LSMAP, a tumor suppressor, influences the incidence and progression of lung cancer by regulating the epithelial-mesenchymal transition pathway (90). C2CD4A is a possible diagnostic marker for colon cancer, and its expression is strongly associated with the tumor stage of colon cancer, which is almost compatible with our results (91). As an intracellular molecular chaperone, RBP1 is highly down-regulated in most malignancies in humans, such as breast cancer (92). Jiang et al. demonstrated that the expression of ANO1 is associated with advanced colon cancer lymphatic metastasis (93). The epidermal growth factor receptor/extracellular signal-regulated kinase signal pathway activation influences the onset and progression of colon cancer by upregulating the expression level of ANO1 in colon cancer (93). Weber et al. reported that the overexpression of HEYL might suppress the formation and migration of tumor cells following inhibiting the intravasation of metastasis-initiating cells (94). Gao et al. have reported that IGF1 was associated with the pivotal precursor to colorectal cancer (95). Lee et al. uncovered that the methylation level of LAMA2 had a crucial role in the onset and advancement of colorectal cancer (96). Although the functions of these hub genes in COAD need additional investigation, our analysis reveals that they are key prognostic variables and may be possible treatment targets.
There are still various significant limitations that need additional investigation. Initially, we simply investigated the prognostic significance of SLP genes in COAD. Due to the limitation of RNA-sequencing, this study presents genes that encode enzymes that are mainly transcription factors and growth factors that are indirectly involved in the metabolism of sphingolipids, such as EEF1A2, GRP, HESR, IGF1, LAMA2, etc. However, the expression of certain genes such as SMPD1, SMPD3, SGMS1, SGMS2, ASAH1, ASAH2, etc, that encode enzymes directly related to the SLP metabolism was not detected, and some genes were not associated with the clinical outcomes of COAD. More research based on new sequencing technologies are required to uncover the potential contributions of genes that encode enzymes directly related to the SLP metabolism in COAD. Secondly, even though our analysis confirmed the varied expression levels of genes related to the risk model at the mRNA as well as the protein levels, this research would become more comprehensive if in-depth biological experiments are further conducted. Consequently, more well-designed experimentation is required to confirm our results. Nevertheless, despite the aforementioned limitations, it is undeniable that this study made the first thorough investigation of SLP genes in COAD and summarization of the pan-cancer overview of SLP genes. Additionally, the cluster and risk score may be used independently to assess the prognosis in people suffering from COAD.
Conclusion
In our investigation, we examined 26 SLP genes strongly linked to the prognosis of COAD patients. Pan-cancer characterization of SLP genes highlighted their crucial role in the onset and progression of tumors. Individuals suffering from COAD were effectively sorted into two metabolic subgroups as per the features of SLP gene expression. Active sphingolipid metabolism subtypes were associated with the unbalanced local tumor immune microenvironment and poor clinical outcomes. Based on the metabolic classifier, we developed a unique and reliable prognostic model for patients with COAD, which gives a new perspective on the prognostic assessment of cancer patients.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
QY, WZ, and WS contributed to this study equally. All the authors participated in the design of study, data collection and processing, bioinformatics analysis, and writing and revising the manusript. All authors read and approved the final manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank Bullet Edits Limited for the linguistic editing of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.1045167/full#supplementary-material
References
1. Lin C, Ng HL, Pan W, Chen H, Zhang G, Bian Z, et al. Exploring different strategies for efficient delivery of colorectal cancer therapy. Int J Mol Sci (2015) 16(11):26936–52. doi: 10.3390/ijms161125995
2. Xu M, Chang J, Wang W, Wang X, Wang X, Weng W, et al. Classification of colon adenocarcinoma based on immunological characterizations: Implications for prognosis and immunotherapy. Front Immunol (2022) 13:934083. doi: 10.3389/fimmu.2022.934083
3. Xiao J, Wang X, Liu Y, Liu X, Yi J, Hu J. Lactate metabolism-associated lncRNA pairs: A prognostic signature to reveal the immunological landscape and mediate therapeutic response in patients with colon adenocarcinoma. Front Immunol (2022) 13:881359. doi: 10.3389/fimmu.2022.881359
4. Yamaji T, Hanada K. Sphingolipid metabolism and interorganellar transport: localization of sphingolipid enzymes and lipid transfer proteins. Traffic (2015) 16(2):101–22. doi: 10.1111/tra.12239
5. Hannun YA, Obeid LM. Sphingolipids and their metabolism in physiology and disease. Nat Rev Mol Cell Biol (2018) 19(3):175–91. doi: 10.1038/nrm.2017.107
6. Modrak DE, Gold DV, Goldenberg DM. Sphingolipid targets in cancer therapy. Mol Cancer Ther (2006) 5(2):200–8. doi: 10.1158/1535-7163.MCT-05-0420
7. Leek J, Johnson W, Parker H, Jaffe A, Storey JJB. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics (2012) 28(6):882–3. doi: 10.1093/bioinformatics/bts034
8. Safran M, Dalah I, Alexander J, Rosen N, Iny Stein T, Shmoish M, et al. GeneCards version 3: the human gene integrator. Database J Biol Database Curation (2010) 2010:baq020. doi: 10.1093/database/baq020
9. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics (2011) 27(12):1739–40. doi: 10.1093/bioinformatics/btr260
10. Green GH, Diggle PJ. On the operational characteristics of the benjamini and hochberg false discovery rate procedure. Stat Appl Genet Mol Biol (2007) 6:Article27. doi: 10.2202/1544-6115.1302
11. Rogula B, Lozano-Ortega G, Johnston KM. A method for reconstructing individual patient data from Kaplan-Meier survival curves that incorporate marked censoring times. MDM Policy Pract (2022) 7(1):23814683221077643. doi: 10.1177/23814683221077643
12. Mayakonda A, Lin D, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118
13. Jin Y, Wang Z, He D, Zhu Y, Chen X, Cao K. Identification of novel subtypes based on ssGSEA in immune-related prognostic signature for tongue squamous cell carcinoma. Cancer Med (2021) 10(23):8693–707. doi: 10.1002/cam4.4341
14. Ren J, Yuan Q, Liu J, Zhong L, Li H, Wu G, et al. Identifying the role of transient receptor potential channels (TRPs) in kidney renal clear cell carcinoma and their potential therapeutic significances using genomic and transcriptome analyses. BMC Med Genomics (2022) 15(1):156. doi: 10.1186/s12920-022-01312-x
15. Yuan Q, Deng D, Pan C, Ren J, Wei T, Wu Z, et 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. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics (2010) 26(12):1572–3. doi: 10.1093/bioinformatics/btq170
17. Bürkner PC, Doebler P, Holling H. Optimal design of the wilcoxon-Mann-Whitney-test. Biometrical J Biometrische Z (2017) 59(1):25–40. doi: 10.1002/bimj.201600022
18. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinf (2013) 14:7. doi: 10.1186/1471-2105-14-7
19. Wang J, Ren J, Liu J, Zhang L, Yuan Q, Dong B. Identification and verification of the ferroptosis- and pyroptosis-associated prognostic signature for low-grade glioma. Bosnian J basic Med Sci (2022) 22(5):728–50. doi: 10.17305/bjbms.2021.6888
20. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
21. Ren J, Wang A, Liu J, Yuan Q. Identification and validation of a novel redox-related lncRNA prognostic signature in lung adenocarcinoma. Bioengineered (2021) 12(1):4331–48. doi: 10.1080/21655979.2021.1951522
22. Miao Y, Yuan Q, Wang C, Feng X, Ren J, Wang C. Comprehensive characterization of RNA-binding proteins in colon adenocarcinoma identifies a novel prognostic signature for predicting clinical outcomes and immunotherapy responses based on machine learning. Combinatorial Chem High throughput screening (2022). doi: 10.2174/1386207325666220404125228
23. Geeleher P, Cox N, Huang RS. pRRophetic: an r package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PloS One (2014) 9(9):e107468. doi: 10.1371/journal.pone.0107468
24. Geeleher P, Cox NJ, Huang RS. Clinical drug response can be predicted using baseline gene expression levels and in vitro drug sensitivity in cell lines. Genome Biol (2014) 15(3):R47. doi: 10.1186/gb-2014-15-3-r47
25. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007
26. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res (2000) 28(1):27–30. doi: 10.1093/nar/28.1.27
27. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci Publ Protein Soc (2019) 28(11):1947–51. doi: 10.1002/pro.3715
28. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res (2021) 49(D1):D545–d51. doi: 10.1093/nar/gkaa970
29. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat software (2010) 33(1):1–22. doi: 10.18637/jss.v033.i01
30. Li C, Tang Z, Zhang W, Ye Z, Liu F. GEPIA2021: integrating multiple deconvolution-based analysis into GEPIA. Nucleic Acids Res (2021) 49(W1):W242–6. doi: 10.1093/nar/gkab418
31. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res (2017) 45(W1):W98–102. doi: 10.1093/nar/gkx247
32. Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. tissue-based map of the human proteome. Science (2015) 347(6220):1260419. doi: 10.1126/science.1260419
33. Thul PJ, Åkesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, et al. A subcellular map of the human proteome. Science (2017) 356(6340). doi: 10.1126/science.aal3321
34. Glaus Garzon JF, Pastrello C, Jurisica I, Hottiger MO, Wenger RH, Borsig L. Tumor cell endogenous HIF-1α activity induces aberrant angiogenesis and interacts with TRAF6 pathway required for colorectal cancer development. Neoplasia (2020) 22(12):745–58. doi: 10.1016/j.neo.2020.10.006
35. Kim SW, Park KC, Jeon SM, Ohn TB, Kim TI, Kim WH, et al. Abrogation of galectin-4 expression promotes tumorigenesis in colorectal cancer. Cell Oncol (2013) 36(2):169–78. doi: 10.1007/s13402-013-0124-x
36. Jiang H, Tang E, Chen Y, Liu H, Zhao Y, Lin M, et al. Squalene synthase predicts poor prognosis in stage I-III colon adenocarcinoma and synergizes squalene epoxidase to promote tumor progression. Cancer Sci (2022) 113(3):971–85. doi: 10.1111/cas.15248
37. Baig N, Li Z, Lu J, Chen H, Yu S, Li T, et al. Clinical significance and comparison of flotillin 1 expression in left and right colon cancer. Oncol Lett (2019) 18(2):997–1004. doi: 10.3892/ol.2019.10401
38. Wang R, Kang Y, Löhr CV, Fischer KA, Bradford CS, Johnson G, et al. Reciprocal regulation of BMF and BIRC5 (Survivin) linked to eomes overexpression in colorectal cancer. Cancer Lett (2016) 381(2):341–8. doi: 10.1016/j.canlet.2016.08.008
39. Karagiota A, Chachami G, Paraskeva E. Lipid metabolism in cancer: The role of acylglycerolphosphate acyltransferases (AGPATs). Cancers (2022) 14(1):228. doi: 10.3390/cancers14010228
40. Liang Q, Wang Y, Lu Y, Zhu Q, Xie W, Tang N, et al. RANK promotes colorectal cancer migration and invasion by activating the Ca(2+)-calcineurin/NFATC1-ACP5 axis. Cell Death Dis (2021) 12(4):336. doi: 10.1038/s41419-021-03642-7
41. Kim WT, Mun JY, Baek SW, Kim MH, Yang GE, Jeong MS, et al. Secretory SERPINE1 expression is increased by antiplatelet therapy, inducing MMP1 expression and increasing colon cancer metastasis. Int J Mol Sci (2022) 23(17):9596. doi: 10.3390/ijms23179596
42. Parikh AA, Fan F, Liu WB, Ahmad SA, Stoeltzing O, Reinmuth N, et al. Neuropilin-1 in human colon cancer: expression, regulation, and role in induction of angiogenesis. Am J Pathol (2004) 164(6):2139–51. doi: 10.1016/S0002-9440(10)63772-8
43. Yin N, Qi X, Tsai S, Lu Y, Basir Z, Oshima K, et al. p38γ MAPK is required for inflammation-associated colon tumorigenesis. Oncogene (2016) 35(8):1039–48. doi: 10.1038/onc.2015.158
44. Grbčić P, Sedić M. Sphingosine 1-phosphate signaling and metabolism in chemoprevention and chemoresistance in colon cancer. Molecules (2020) 25(10):2436. doi: 10.3390/molecules25102436
45. Li D, Li Q, Zhuo C, Xu Y, Cai S. Contribution of FOXC1 to the progression and metastasis and prognosis of human colon cancer. J Clin Oncol (2015) 33(3_suppl):636. doi: 10.1200/jco.2015.33.3_suppl.636
46. Castro ED, Mathias PPM, Batista WL, Sato AYS, Toledo MS, de Almeida VT, et al. Knockdown of the inducible nitric oxide synthase (NOS2) splicing variant S3 promotes autophagic cell death from nitrosative stress in SW480 human colon cancer cells. Cell Biol Int (2022) 46(1):158–69. doi: 10.1002/cbin.11717
47. Spiegel A, Hundley TR, Chen J, Gao J, Ouyang N, Liu X, et al. NO-donating aspirin inhibits both the expression and catalytic activity of inducible nitric oxide synthase in HT-29 human colon cancer cells. Biochem Pharmacol (2005) 70(7):993–1000. doi: 10.1016/j.bcp.2005.06.027
48. Scheurlen KM, Chariker JH, Kanaan Z, Littlefield AB, George JB, Seraphine C, et al. The NOTCH4-GATA4-IRG1 axis as a novel target in early-onset colorectal cancer. Cytokine Growth Factor Rev (2022) 67:25–34. doi: 10.1016/j.cytogfr.2022.06.002
49. Wu G, Chen Z, Li J, Ye F, Chen G, Fan Q, et al. NOTCH4 is a novel prognostic marker that correlates with colorectal cancer progression and prognosis. J Cancer (2018) 9(13):2374–9. doi: 10.7150/jca.26359
50. Zhang Z, Bu X, Yang J, Zhu S, He S, Zheng J, et al. NOTCH4 regulates colorectal cancer proliferation, invasiveness, and determines clinical outcome of patients. J Cell Physiol (2018) 233(10):6975–85. doi: 10.1002/jcp.26619
51. Dziegielewska B, Brautigan DL, Larner JM, Dziegielewski J. T-Type Ca2+ channel inhibition induces p53-dependent cell growth arrest and apoptosis through activation of p38-MAPK in colon cancer cells. Mol Cancer Res (2014) 12(3):348–58. doi: 10.1158/1541-7786.MCR-13-0485
52. Tian W, Zhang W, Zhang Y, Zhu T, Hua Y, Li H, et al. FABP4 promotes invasion and metastasis of colon cancer by regulating fatty acid transport. Cancer Cell Int (2020) 20:512. doi: 10.1186/s12935-020-01582-4
53. Liu D, Xiang X, Chen Y, Jiao Y, Wang L, Wang Y, et al. Identification KCNE4 and DNASE1L3 for prognostic utility in colorectal cancer: A bioinformatics analysis. Research Square (2021) doi: 10.21203/rs.3.rs-1039855/v1.
54. Wu Y, Liang X, Ni J, Zhao R, Shao S, Lu S, et al. Effect of ISM1 on the immune microenvironment and epithelial-mesenchymal transition in colorectal cancer. Front Cell Dev Biol (2021) 9:681240. doi: 10.3389/fcell.2021.681240
55. Yang Z, Chen H, Huo L, Yang Z, Bai Y, Fan X, et al. Epigenetic inactivation and tumor-suppressor behavior of NGFR in human colorectal cancer. Mol Cancer Res (2015) 13(1):107–19. doi: 10.1158/1541-7786.MCR-13-0247
56. Chen H, Huang J, Chen C, Jiang Y, Feng X, Liao Y, et al. NGFR increases the chemosensitivity of colorectal cancer cells by enhancing the apoptotic and autophagic effects of 5-fluorouracil via the activation of S100A9. Front Oncol (2021) 11:652081. doi: 10.3389/fonc.2021.652081
57. Chen W, Huang J, Xiong J, Fu P, Chen C, Liu Y, et al. Identification of a tumor microenvironment-related gene signature indicative of disease prognosis and treatment response in colon cancer. Oxid Med Cell Longevity (2021) 2021:6290261. doi: 10.1155/2021/6290261
58. Bataller M, Sánchez-García A, Garcia-Mayea Y, Mir C, Rodriguez I, ME LL. The role of sphingolipids metabolism in cancer drug resistance. Front Oncol (2021) 11:807636. doi: 10.3389/fonc.2021.807636
59. Lin M, Liao W, Dong M, Zhu R, Xiao J, Sun T, et al. Exosomal neutral sphingomyelinase 1 suppresses hepatocellular carcinoma via decreasing the ratio of sphingomyelin/ceramide. FEBS J (2018) 285(20):3835–48. doi: 10.1111/febs.14635
60. Cheng Q, Li X, Wang Y, Dong M, Zhan FH, Liu J. The ceramide pathway is involved in the survival, apoptosis and exosome functions of human multiple myeloma cells in vitro. Acta pharmacol Sin (2018) 39(4):561–8. doi: 10.1038/aps.2017.118
61. Wang D, Madunić K, Zhang T, Mayboroda OA, Lageveen-Kammeijer GSM, Wuhrer M. High diversity of glycosphingolipid glycans of colorectal cancer cell lines reflects the cellular differentiation phenotype. Mol Cell Proteomics (2022) 21(6):100239. doi: 10.1016/j.mcpro.2022.100239
62. Jennemann R, Volz M, Bestvater F, Schmidt C, Richter K, Kaden S, et al. Blockade of glycosphingolipid synthesis inhibits cell cycle and spheroid growth of colon cancer cells in vitro and experimental colon cancer incidence In vivo. Int J Mol Sci (2021) 22(19):10539. doi: 10.3390/ijms221910539
63. Machala M, Procházková J, Hofmanová J, Králiková L, Slavík J, Tylichová Z, et al. Colon cancer and perturbations of the sphingolipid metabolism. Int J Mol Sci (2019) 20(23):6051. doi: 10.3390/ijms20236051
64. Hait NC, Maiti A. The role of sphingosine-1-Phosphate and ceramide-1-Phosphate in inflammation and cancer. Mediators Inflamm (2017) 2017:4806541. doi: 10.1155/2017/4806541
65. Duan RD. Alkaline sphingomyelinase: an old enzyme with novel implications. Biochim Biophys Acta (2006) 1761(3):281–91. doi: 10.1016/j.bbalip.2006.03.007
66. García-Barros M, Coant N, Truman JP, Snider AJ, Hannun YA. Sphingolipids in colon cancer. Biochim Biophys Acta (2014) 1841(5):773–82. doi: 10.1016/j.bbalip.2013.09.007
67. Kawamori T, Osta W, Johnson KR, Pettus BJ, Bielawski J, Tanaka T, et al. Sphingosine kinase 1 is up-regulated in colon carcinogenesis. FASEB J (2006) 20(2):386–8. doi: 10.1096/fj.05-4331fje
68. Kawamori T, Kaneshiro T, Okumura M, Maalouf S, Uflacker A, Bielawski J, et al. Role for sphingosine kinase 1 in colon carcinogenesis. FASEB J (2009) 23(2):405–14. doi: 10.1096/fj.08-117572
69. Kreitzburg KM, van Waardenburg R, Yoon KJ. Sphingolipid metabolism and drug resistance in ovarian cancer. Cancer Drug resistance (2018) 1:181–97. doi: 10.20517/cdr.2018.06
70. Kondapuram SK, Coumar MS. Pan-cancer gene expression analysis: Identification of deregulated autophagy genes and drugs to target them. Gene (2022) 844:146821. doi: 10.1016/j.gene.2022.146821
71. Jaiswal PK, Goel A, Mittal RD. Survivin: A molecular biomarker in cancer. Indian J Med Res (2015) 141(4):389–97. doi: 10.4103/0971-5916.159250
72. Wang Y, Li X, Wang H, Zhang G. CircCAMSAP1 promotes non-small cell lung cancer proliferation and inhibits cell apoptosis by sponging miR-1182 and regulating BIRC5. Bioengineered (2022) 13(2):2428–39. doi: 10.1080/21655979.2021.2011639
73. Schor NF. The p75 neurotrophin receptor in human development and disease. Prog Neurobiol (2005) 77(3):201–14. doi: 10.1016/j.pneurobio.2005.10.006
74. Xiao Z, Nie K, Han T, Cheng L, Zhang Z, Peng W, et al. Development and validation of a TNF family-based signature for predicting prognosis, tumor immune characteristics, and immunotherapy response in colorectal cancer patients. J Immunol Res (2021) 2021:6439975. doi: 10.1155/2021/6439975
75. Li RZ, Wang XR, Wang J, Xie C, Wang XX, Pan HD, et al. The key role of sphingolipid metabolism in cancer: New therapeutic targets, diagnostic and prognostic values, and anti-tumor immunotherapy resistance. Front Oncol (2022) 12:941643. doi: 10.3389/fonc.2022.941643
76. Liu YY, Hill RA, Li YT. Ceramide glycosylation catalyzed by glucosylceramide synthase and cancer drug resistance. Adv Cancer Res (2013) 117:59–89. doi: 10.1016/B978-0-12-394274-6.00003-0
77. Morad SAF, Cabot MC. The onus of sphingolipid enzymes in cancer drug resistance. Adv Cancer Res (2018) 140:235–63. doi: 10.1016/bs.acr.2018.04.013
78. Ghosh S, Juin SK, Majumdar S. Cancer stem cells and ceramide signaling: the cutting edges of immunotherapy. Mol Biol Rep (2020) 47(10):8101–11. doi: 10.1007/s11033-020-05790-z
79. Oskouian B, Saba JD. Cancer treatment strategies targeting sphingolipid metabolism. Adv Exp Med Biol (2010) 688:185–205. doi: 10.1007/978-1-4419-6741-1_13
80. Li Z, Xin S, Yu S, Liang J, Zhang X. Prognostic signatures and therapeutic value based on the notch pathway in renal clear cell carcinoma. Oxid Med Cell Longevity (2022) 2022:1669664. doi: 10.1155/2022/1669664
81. Guruharsha KG, Kankel MW, Artavanis-Tsakonas S. The notch signalling system: recent insights into the complexity of a conserved pathway. Nat Rev Genet (2012) 13(9):654–66. doi: 10.1038/nrg3272
82. Yin L, Velazquez OC, Liu ZJ. Notch signaling: emerging molecular targets for cancer therapy. Biochem Pharmacol (2010) 80(5):690–701. doi: 10.1016/j.bcp.2010.03.026
83. Staudacher JJ, Bauer J, Jana A, Tian J, Carroll T, Mancinelli G, et al. Activin signaling is an essential component of the TGF-β induced pro-metastatic phenotype in colorectal cancer. Sci Rep (2017) 7(1):5569. doi: 10.1038/s41598-017-05907-8
84. Masugi Y, Nishihara R, Hamada T, Song M, da Silva A, Kosumi K, et al. Tumor PDCD1LG2 (PD-L2) expression and the lymphocytic reaction to colorectal cancer. Cancer Immunol Res (2017) 5(11):1046–55. doi: 10.1158/2326-6066.CIR-17-0122
85. Alderdice M, Craig SG, Humphries MP, Gilmore A, Johnston N, Bingham V, et al. Evolutionary genetic algorithm identifies IL2RB as a potential predictive biomarker for immune-checkpoint therapy in colorectal cancer. NAR Genomics Bioinf (2021) 3(2):lqab016. doi: 10.1093/nargab/lqab016
86. Topalian SL, Drake CG, Pardoll DM. Immune checkpoint blockade: a common denominator approach to cancer therapy. Cancer Cell (2015) 27(4):450–61. doi: 10.1016/j.ccell.2015.03.001
87. Drake CG, Lipson EJ, Brahmer JR. Breathing new life into immunotherapy: review of melanoma, lung and kidney cancer. Nat Rev Clin Oncol (2014) 11(1):24–37. doi: 10.1038/nrclinonc.2013.208
88. Taglia L, Matusiak D, Benya RV. GRP-induced up-regulation of Hsp72 promotes CD16+/94+ natural killer cell binding to colon cancer cells causing tumor cell cytolysis. Clin Exp Metastasis (2008) 25(4):451–63. doi: 10.1007/s10585-008-9151-9
89. Zhou Y, Zhang Y, Li W, Xu J, He X, Li X, et al. TCEAL2 as a tumor suppressor in renal cell carcinoma is associated with the good prognosis of patients. Cancer Manage Res (2020) 12:9589–97. doi: 10.2147/CMAR.S271647
90. Chang CY, Wu KL, Chang YY, Liu YW, Huang YC, Jian SF, et al. The downregulation of LSAMP expression promotes lung cancer progression and is associated with poor survival prognosis. J personalized Med (2021) 11(6). doi: 10.3390/jpm11060578
91. Lu J, Chen Q. Transcriptome-based identification of molecular markers related to the development and prognosis of colon cancer. Nucleosides nucleotides Nucleic Acids (2021) 40(11):1114–24. doi: 10.1080/15257770.2021.1975298
92. Yu J, Perri M, Jones JW, Pierzchalski K, Ceaicovscaia N, Cione E, et al. Altered RBP1 gene expression impacts epithelial cell retinoic acid, proliferation, and microenvironment. Cells (2022) 11(5):792. doi: 10.3390/cells11050792
93. Jiang Y, Cai Y, Shao W, Li F, Guan Z, Zhou Y, et al. MicroRNA−144 suppresses aggressive phenotypes of tumor cells by targeting ANO1 in colorectal cancer. Oncol Rep (2019) 41(4):2361–70. doi: 10.3892/or.2019.7025
94. Weber S, Koschade SE, Hoffmann CM, Dubash TD, Giessler KM, Dieter SM, et al. The notch target gene HEYL modulates metastasis forming capacity of colorectal cancer patient-derived spheroid cells in vivo. BMC Cancer (2019) 19(1):1181. doi: 10.1186/s12885-019-6396-4
95. Gao Y, Katki H, Graubard B, Pollak M, Martin M, Tao Y, et al. Serum IGF1, IGF2 and IGFBP3 and risk of advanced colorectal adenoma. Int J Cancer (2012) 131(2):E105–13. doi: 10.1002/ijc.26438
Keywords: colon adenocarcinoma, sphingolipid metabolism, pan-cancer analysis, hierarchical clustering, stratification model, prognostic biomarker
Citation: Yuan Q, Zhang W and Shang W (2022) Identification and validation of a prognostic risk-scoring model based on sphingolipid metabolism-associated cluster in colon adenocarcinoma. Front. Endocrinol. 13:1045167. doi: 10.3389/fendo.2022.1045167
Received: 15 September 2022; Accepted: 24 October 2022;
Published: 28 November 2022.
Edited by:
Rashika El Ridi, Cairo University, EgyptReviewed by:
Nicolas Coant, Stony Brook Medicine, United StatesAlice Vladimirovna Alessenko, Russian Academy of Sciences, Russia
Copyright © 2022 Yuan, Zhang and Shang. 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: Qihang Yuan, cWloYW5nZHlAMTYzLmNvbQ==
†These authors have contributed equally to this work