Int J Med Sci 2024; 21(6):1103-1116. doi:10.7150/ijms.91446 This issue Cite
Research Paper
1. Department of General Surgery, Cancer Hospital of Dalian University of Technology, Cancer Hospital of China Medical University, Liaoning Cancer Hospital & Institute, Shenyang, Liaoning 110042, China.
2. Liaoning University of Traditional Chinese Medicine, Shenyang, Liaoning 110847, China.
3. Department of Colorectal Surgery, Shengjing Hospital of China Medical University, Shenyang, Liaoning 110004, China.
4. Department of Endoscopy, Cancer Hospital of Dalian University of Technology, Cancer Hospital of China Medical University, Liaoning Cancer Hospital & Institute, Shenyang, Liaoning 110042, China.
5. Department of Gynecology, People's Hospital of Liaoning Province, Shenyang, Liaoning 110016, China.
6. Department of Breast Surgery, Panjin Central Hospital, Panjin, Liaoning 124010, China.
7. Liaoning Provincial Key Laboratory of Precision Medicine for Malignant Tumors, Shenyang, Liaoning 110042, China.
8. Markey Cancer Center, University of Kentucky, Lexington, KY 40536, United States of America.
# These authors contributed equally to this work.
Background: Colorectal cancer (CRC) has a high morbidity and mortality. Ferroptosis is a phenomenon in which metabolism and cell death are closely related. The role of ferroptosis-related genes in the progression of CRC is still not clear. Therefore, we screened and validated the ferroptosis-related genes which could determine the prevalence, risk and prognosis of patients with CRC.
Methods: We firstly screened differentially expressed ferroptosis-related genes by The Cancer Genome Atlas (TCGA) database. Then, these genes were used to construct a risk-score model using the least absolute shrinkage and selection operator (LASSO) regression algorithm. The function and prognosis of the ferroptosis-related genes were confirmed using multi-omics analysis. The gene expression results were validated using publicly available databases and qPCR. We also used publicly available data and ferroptosis-related genes to construct a prognostic prediction nomogram.
Results: A total of 24 differential expressed genes associated with ferroptosis were screened in this study. A three-gene risk score model was then established based on these 24 genes and GPX3, CDKN2A and SLC7A11 were selected. The significant prognostic value of this novel three-gene signature was also assessed. Furthermore, we conducted RT-qPCR analysis on cell lines and tissues, and validated the high expression of CDKN2A, GPX3 and low expression of SLC7A11 in CRC cells. The observed mRNA expression of GPX3, CDKN2A and SLC7A11 was consistent with the predicted outcomes. Besides, eight variables including selected ferroptosis related genes were included to establish the prognostic prediction nomogram for patients with CRC. The calibration plots showed favorable consistency between the prediction of the nomogram and actual observations. Also, the time-dependent AUC (>0.7) indicated satisfactory discriminative ability of the nomogram.
Conclusions: The present study constructed and validated a novel ferroptosis-related three-gene risk score signature and a prognostic prediction nomogram for patients with CRC. Also, we screened and validated the ferroptosis-related genes GPX3, CDKN2A, and SLC7A11 which could serve as novel biomarkers for patients with CRC.
Keywords: Colorectal Cancer, Ferroptosis, Risk model, Nomogram, Prognosis
In 2020, colorectal cancer (CRC) was ranked third in incidence and second in mortality among the most prevalent cancers [1]. The incidence of CRC is on the rise in countries with an increasing human development index [2]. Despite improvements in diagnosis and treatment, some patients with CRC may experience disease progression, leading to less effective therapy and poorer survival outcomes [3-5]. Therefore, there is an urgent need to identify the novel prognostic indicators and explore the responding regulatory mechanisms for patients with CRC.
Iron has significant biological importance in the development of CRC and other solid cancer [6, 7]. Unlike normal cells, many kinds of tumor cells exhibit a heightened dependence on iron for accelerated growth, rendering them vulnerable to iron depletion. Besides, tumor cells also have a high proliferation rate owing to their addiction to iron [8-10]. Ferroptosis is a novel form of cell death represented by dysregulation of iron and related lipid metabolism [11, 12]. Dysregulation of ferroptosis has been shown to contribute to the proliferation, invasion, and metastasis of solid cancer cells, including CRC [13-15]. As for the specific relationships between ferroptosis and CRC, and the importance of ferroptosis in the development of CRC, it may mainly lie on the function of ferroptosis-related genes. Ferroptosis-related genes significantly regulate molecular changes and physiological activities of CRC cells [16, 17]. Abnormal expression of ferroptosis-related genes may affect the occurrence, development and even the treatment sensitivity of CRC. It has been shown that decreased expression of ferroptosis-related gene Glutathione Peroxidase 4 (GPX4) may increase CRC cells' sensitivity to oxidative stress and decreased CRC cells' tolerance to lipid peroxides, which may finally lead to the death of CRC cells [18]. Meanwhile, it has also been found that the down-regulation of ferroptosis-related gene Metallothionein-1G resulted in a significant inhibition of CRC cells' proliferation [19]. Other studies even found that the enhanced expression and activity of ferroptosis-related genes in CRC cells were also crucial for many kinds of anti-tumor drugs to activate iron-dependent cell death and inhibit CRC development. For example, the use of ferroptosis inducers can significantly reverse Oxaliplatin resistance in CRC cells [20]. All these above findings proved that ferroptosis may significantly affect the malignant development of CRC. Ferroptosis-related genes may serve as novel biomarkers and therapeutic targets for the patients with CRC.
Our study aimed to investigate and validate ferroptosis-related genes which could serve as novel prognostic biomarkers for patients with CRC by conducting a comprehensive bioinformatics analysis on The Cancer Genome Atlas (TCGA) datasets, integrating gene expression levels, nomograms, immune infiltration, spatial and single-cell location, DNA methylation, clinical data and Quantitative RT-PCR (RT-qPCR) analysis. In addition, gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were also used to explore the potential pathways and mechanisms associated with selected ferroptosis-related genes in the various groups. In conclusion, we screened and validated the differentiated expressed ferroptosis-related genes GPX3, CDKN2A, and SLC7A11, which may serve as novel prognostic biomarkers for patients with CRC. Furthermore, this novel three-gene risk score signature and nomogram could be helpful for monitoring the survival outcomes of patients with CRC.
Ferroptosis-related genes were retrieved using MSigDB version 7.1 [21, 22]. This study incorporated previously reported genes identified through literature review [23, 24]. A gene set containing 404 ferroptosis-related genes was obtained by eliminating overlapping genes.
Gene expression data for 647 CRC samples from TCGA and 51 normal samples from the GTEx project were downloaded from the UCSC Xena. The dataset included samples from TCGA, TARGET, and GTEx. The R package “TCGA biolinks” was used to retrieve clinical data pertaining to patients with CRC. Gene expression data were collected from patients. Patients with prognostic information-deficient data were excluded from the analysis. Finally, 619 patients were included in TCGA dataset. Ethics committee approval was not required because all data were publicly accessible.
Initially, we screened 1765 differentially expressed genes (DEGs) identified in TCGA gene expression matrix. The “DESeq2” package in R software (v 3.6.3) was used to analyze and identify DEGs between TCGA-CRC samples with normal tissues [25-28]. DEGs were filtered using an adjusted P-value threshold of < 0.05 and an absolute log2-fold change >1. Twenty-four ferroptosis-related genes were selected for downstream analysis. Metascape (https://metascape.org/gp/index.html#/main/step1) was used for functional enrichment analysis of the selected DEGs [29].
We used Cox proportional hazards regression and least absolute shrinkage and selection operator (LASSO) regression analysis to select 24 differentially expressed ferroptosis-related genes from the TCGA dataset. LASSO regression analysis with 10-fold cross-validation was used for establishing an ferroptosis-related genes-based scoring model to predict the patient's prognosis. The risk score was calculated as follows:
Risk Score = β1 × exp1 + β2 × exp2 + … + βi × exp i
where β = the coefficient of ferroptosis-related genes, and exp = the expression levels of ferroptosis-related genes.
Normalized gene expression levels were calculated using the “edgeR” package. We also calculated a risk score for each available patient with CRC and generated distribution and receiver operating characteristic (ROC) curves using “time ROC” software [30]. The patients were stratified into high- and low-risk categories based on their median risk scores. Kaplan-Meier (KM) survival analysis was used to determine the performance of the ferroptosis-related genes-based scoring model in predicting prognosis. The difference in overall survival (OS) was compared in patients with CRC in the two subgroups.
FDR correction was applied to the significant enrichment P-values. Pre-ranked GSEA was used to confirm the biological functions and pathways of each module. The Molecular Signature Database v5.1 predefined gene sets were used. The enrichment score was calculated based on the rank order of genes determined by the random forest accuracy. GSVA was conducted using the GSVA package in R software.
The human CRC cell lines SW480, SW620, and HCT116 and the human colonic epithelial cell line NCM460 were procured from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). The cells were cultured in L-15 or RPMI-1640 medium (Gibco, Carlsbad, USA) supplemented with 10% fetal bovine serum (FBS, Gibco) and 1% penicillin-streptomycin (HyClone, Logan, USA) at 37°C in 5% CO2.
RT-qPCR was performed to assess the expression level of GPX3, CDKN2A and SLC7A11 mRNA in CRC cell lines and tissue samples. Total RNA was isolated from the cells or tissue samples using TRIzol reagent (Invitrogen). Complementary DNA (cDNA) was synthesized using the PrimeScript RT Reagent Kit (TaKaRa, Japan). RT-qPCR was performed on an ABI 7500 H T system using SYBR Green PCR master mix (Applied Biosystems, CA, USA) and primers. The specific mRNA expression levels were quantified using the 2-ΔΔCT method. β-actin served as the control for normalization. The primer sequences for RT-qPCR were shown in Table 1.
The sequences of primers.
Gene name | Sequences | |
---|---|---|
GPX3 | Forward | 5'-CAAGAGAAGTCGAAG ATGGACTG-3' |
Reverse | 5'-GGGGATGTACTCCTCCCCAT-3' | |
CDKN2A | Forward | 5'-GGGGTCGGGTAGAGGAGG-3' |
Reverse | 5'- GCCCATCATCATGACCTGGA-3' | |
SLC7A11 | Forward | 5'-GGAGAAGGAATTCCAGGTCAT-3' |
Reverse | 5'- AGCAATACAAGGAAGCCTTAGGT-3' | |
β-actin | Forward | 5'- GTGGACATCCGCAAAGAC-3' |
Reverse | 5'- AAAGGGTGTAACGCAACTA-3' |
The code for transcriptomic, spatial, and statistical analyses is available in GitHub (https://github.com/lrb0533/CRC-Spatial). RNA-seq data for single-cell studies were collected from public databases (www.10xgenomics.com).
The Colon Adenocarcinoma (COAD) single-cell dataset, Genomic Spatial Event (GSE), comprising GSM6061704-11773, GSM6061704-15564, and GSM6061704-8810, was obtained from the Gene Expression Omnibus (GEO) database. The dataset comprised 11 samples. Quality control was subsequently performed on the data. We selected cells with <10% mitochondrial genes, >200 total genes, and genes with expression levels between 200 and 6000 expressed in at least three cells. The number of highly variable genes was set to 2000. SCT correction was applied to these 11 samples for integration. The tSNE method was used to reduce data dimensionality by setting the “DIMS” parameter to 20. Cell clustering was performed using the “KNN” method with a resolution of 0.6. These cells were subsequently annotated using various cell surface markers. Finally, the percentage of ferroptosis-related genes in each cell was determined by using the “Percentage Feature Set” function to import the relevant genes.
Univariate and multivariate Cox regression analyses were conducted to analyze the clinicopathological parameters, including T stage, N stage, M stage, age, CEA level, CDKN2A expression, GPX3 expression, and SLC7A11 expression, to determine the independent predictive ability of the risk score system. The clinical characteristics of these patients with CRC were obtained from TCGA-COAD dataset (Table 2). The rms package was used to integrate independent prognostic factors and construct a nomogram for predicting the 1-, 3-, and 5-year OS. ROC analysis and calibration were used to evaluate the discriminative ability of the nomogram [31].
Tumor IMmune Estimation Resource (TIMER) enables automated visualization and analysis of immune cell infiltration in all TCGA tumors [32, 33]. The TIMER algorithm uses six discrete immune cell subsets to estimate infiltration. The cellular composition includes B cells, CD4+ T cells, CD8+ T cells, macrophages, neutrophils, and dendritic cells. This study aimed to analyze the results of infiltration estimation and compare immune cell subsets in low- and high-risk groups (http://timer.cistrome.org/) [34]. The MethSurv database was used to analyze the DNA methylation sites of ferroptosis-related genes in TCGA database (https://biit.cs.ut.ee/methsurv/).
The miRNAs and TFs of the three selected ferroptosis-related genes were predicted using TargetScan Human 7.2 and GeneCards. Results were visualized using the Cytoscape software (v3.8.0).
Statistical analyses were performed using GraphPad Prism and the R software. Kaplan-Meier survival analyses were performed using log-rank tests. When applicable, hazard ratios (HRs) and 95% confidence intervals (CIs) were reported. Student's t-test and independent sample t-test were used for two-group comparisons. A two-tailed P value of <0.05 was considered statistically significant without explicit annotation.
Clinical characteristics of patients with CRC from TCGA-COAD dataset.
Characteristic | Levels | Overall |
---|---|---|
n | 619 | |
T stage, n (%) | T1 | 20 (3.2%) |
T2 | 105 (17%) | |
T3 | 422 (68.4%) | |
T4 | 70 (11.3%) | |
N stage, n (%) | N0 | 351 (57%) |
N1 | 150 (24.4%) | |
N2 | 115 (18.7%) | |
M stage, n (%) | M0 | 459 (84.1%) |
M1 | 87 (15.9%) | |
Pathologic stage, n (%) | Stage I | 105 (17.5%) |
Stage II | 227 (37.9%) | |
Stage III | 179 (29.9%) | |
Stage IV | 88 (14.7%) | |
Gender, n (%) | Female | 289 (46.7%) |
Male | 330 (53.3%) | |
Race, n (%) | Asian | 12 (3.3%) |
Black or African American | 65 (17.6%) | |
White | 292 (79.1%) | |
Age, n (%) | ≤65 | 269 (43.5%) |
>65 | 350 (56.5%) | |
BMI, n (%) | <25 | 98 (32.2%) |
>=25 | 206 (67.8%) | |
Residual tumor, n (%) | R0 | 450 (91.5%) |
R1 | 6 (1.2%) | |
R2 | 36 (7.3%) | |
CEA level, n (%) | ≤5 | 252 (63.5%) |
>5 | 145 (36.5%) | |
Perineural invasion, n (%) | No | 172 (74.1%) |
Yes | 60 (25.9%) | |
Lymphatic invasion, n (%) | No | 331 (59.3%) |
Yes | 227 (40.7%) | |
History of colon polyps, n (%) | No | 364 (68.4%) |
Yes | 168 (31.6%) | |
Colon polyps present, n (%) | No | 207 (69.5%) |
Yes | 91 (30.5%) | |
Neoplasm type, n (%) | Colon adenocarcinoma | 454 (73.3%) |
Rectum adenocarcinoma | 165 (26.7%) | |
OS event, n (%) | Alive | 492 (79.5%) |
Dead | 127 (20.5%) |
This study was approved by the Ethics Committee of Cancer Hospital of China Medical University. Written informed consent was obtained from all participants, including healthy volunteers and patients who provided clinical specimens.
We firstly conducted an integrated analysis to investigate ferroptosis-related DEGs in CRC, the schematic overview of the whole study was shown in Figure 1A. A comprehensive analysis of 404 ferroptosis-related genes was conducted using the GeneCard database and literature review. We identified 1765 DEGs in 647 tumor and 51 normal samples from the TCGA COAD database. Volcano and heat maps were generated to visualize the DEGs (Figure 1B-C). Both datasets were plotted using Venn diagrams and 24 ferroptosis-related DEGs were selected (Figure 1D). The Log2-fold change value and adjusted P-value of these 24 genes were also shown in Table 3. Figure 1E displayed the heatmap of these 24 ferroptosis-related DEGs based on the TCGA datasets We then sequenced 24 ferroptosis-related DEGs with varying degrees of correlation (Figure 1F) and conducted GO and KEGG enrichment analyses to examine the functions of the selected genes. The study found these 24 DEGs were mainly enriched in processes related to homocysteine and sulfur amino acid metabolism, cellular response to oxidative stress, and the Vitamin D receptor pathway (Figure 1G-J). Subsequently, an interaction network of these 24 DEGs was established using the STRING database (Figure 1K). Furthermore, the overall correlations of the 24 DEGs were also analyzed (Figure 1L).
The selected 24 ferroptosis-related DEGs.
Gene name | Log2-fold change value | Adjusted P-value |
---|---|---|
ALB | 2.770574661 | 1.07E-58 |
CA9 | 2.308236488 | 4.31E-47 |
CBS | 2.109997134 | 2.35E-06 |
CDKN2A | 2.217311987 | 3.65E-15 |
CP | 2.238971828 | 9.95E-21 |
DPEP1 | -2.33747343 | 4.18E-15 |
ETV4 | -2.41166948 | 2.06E-30 |
GDPD5 | 4.140235566 | 2.76E-25 |
GPX3 | 3.880235305 | 5.23E-46 |
HCAR1 | 3.682487512 | 5.68E-63 |
HILPDA | 2.066355368 | 5.35E-39 |
KRT16 | 5.849761786 | 3.60E-72 |
KRT6B | 5.060435957 | 1.89E-65 |
LCN2 | 5.121558867 | 2.23E-10 |
MT1G | 5.725881403 | 3.80E-36 |
MYCN | 6.123645098 | 5.03E-58 |
NGB | 2.543142061 | 8.06E-79 |
NOX4 | 2.826449302 | 4.88E-53 |
SCD | -3.29715582 | 1.97E-41 |
SLC7A11 | -2.0894971 | 7.07E-18 |
TFAP2A | -3.63533677 | 4.23E-42 |
TUBB4A | 4.082463138 | 6.25E-17 |
TYRO3 | -5.58165687 | 9.02E-19 |
UCHL1 | 3.548301809 | 1.31E-35 |
We performed LASSO, COX regression analysis on these 24 ferroptosis-related DEGs. Finally, a total of three ferroptosis-related genes were identified. The formula for calculating the risk score: Risk score= (0.0853 × CDKN2A expression level) + (0.0745 × GPX3 expression level) + (-0.0413 × SLC7A11 expression level). The gene set was optimized using regression coefficient calculations (Figure 2A-B). The survival curve distribution of the risk model was analyzed. The coefficients were determined through 10-fold cross-validation, and GPX3, CDKN2A, and SLC7A11were identified as the most predictive variables based on minimum criteria.
Patients categorized as high risk exhibited a statistically significant decrease in survival rate compared to those categorized as low risk (Figure 2C). Moreover, Kaplan-Meier analysis and log-rank testing revealed a significant reduction in OS time for high-risk patients compared to low-risk patients (HR=1.85, P<0.001, Figure 2D). In addition, the ROC curves of this signature exhibited AUC values of 0.627, 0.627, and 0.620 at 1, 3, and 5 years, respectively (Figure 2E). GSEA analysis was conducted to investigate the potential influence of selected 3 ferroptosis-related gene expression levels on the transcriptomic profiles of both high-risk and low-risk groups of CRC. The results showed that pathways related to Fatty Acid Metabolism, IL6_JAK_STAT3 signaling, IL2_STAT5 signaling, WNT_β-catenin Signaling, Apical Junction, myogenesis, KARS Signaling, and interferon-gamma response were enriched in the high-risk population (Supplementary figure 1).
We conducted RT-qPCR analysis on three colorectal cancer cell lines (SW480, SW620, and HCT116) and one normal human colonic epithelial cell line (NCM460) to validate the aforementioned hypothesis. The results revealed that the observed mRNA expression of GPX3、CDKN2A and SLC7A11 in CRC cell lines was almost consistent with the aforementioned bioinformation results (Figures 3A-C). The expression of GPX3 and CDKN2A mRNA was significantly higher in CRC cell lines than in normal human colonic epithelial cell line. However, SLC7A11 had opposite results. These results were subsequently validated in tissues from patients with colorectal cancer. We collected 4 pairs CRC and normal tissues samples, and validated the mRNA expression of these three selected genes. As a result, we found that the expression of CDKN2A mRNA and GPX3 mRNA was significantly higher in CRC tissues than in normal tissues, the expression of SLC7A11 mRNA was significantly lower than in normal tissues (Figure 3D-F).
Identification and functional enrichment analysis of dysregulated ferroptosis-related genes in CRC. A: Schematic overview of the whole study. B-C: Differential gene volcano, heatmap of CRC based on TGCA data. D: Venn diagram representing intersections of CRC DEGs and ferroptosis-related genes. E: Heatmap of the expression levels of 24 ferroptosis-related DEGs based on TCGA data. F: Map of ranking of differences of 24 ferroptosis-related DEGs. G-I: Enriched Gene Ontology terms and KEGG pathways of 24 ferroptosis-related DEGs. J: Enriched Gene Ontology terms and KEGG pathways of 24 ferroptosis-related DEGs according to the Metascape. K: The correlation network of 24 ferroptosis-related DEGs according to the STRING database. L: The correlations heatmap of 24 ferroptosis-related DEGs.
LASSO regression and risk score calculation. A: Coefficient value of 24 ferroptosis-related DEGs. B: Partial likelihood deviance of 24 ferroptosis-related DEGs. C: Risk score and survival time distributions, and heatmaps of gene-expression levels of the ferroptosis-related signature based on the TCGA data. D: Kaplan-Meier analysis suggests the survival outcome in the high-risk and low-risk groups in the TCGA cohort. E: ROC curves and AUC values of the risk score model for predicting the 1-, 3-, and 5-year OS times in the TCGA cohort.
The diagnostic values of GPX3, CDKN2A, and SLC7A11 for death were evaluated using ROC curve analysis, and we found all three genes exhibited satisfactory diagnostic potential for death (Figure 4A-D). Furthermore, we analyzed the correlations between the immune cell fraction of each sample, as precomputed in the single-molecule and cloud data, to compare the correlations of the three genes with immune cells. The five immune cells that exhibited the strongest correlation with these three genes are listed. GPX3 was more closely related to IDC, macrophages, DC, and mast cells; CDKN2A was more closely related to NK CD56bright cells, IDC, NK cells, TReg, Cytotoxic cells, and TH1 cells. SLC7A11 was found to be more correlated with TH2 cells, TCM, T helper cells, macrophages, and Tgd (Figure 4B-D). Finally, approximately 50 immune checkpoint-related genes were assessed in relation to GPX3, CDKN2A, and SLC7A11 to investigate their involvement in CRC (Figure 4E-G). Heat map analysis of single gene co-expression revealed a strong correlations between GPX3 and ADORA2A, BTLA, CD200, CD200R1, and CD244; CDKN2A have higher correlations with CD276, CD27, CD70, HAVCR2, IDO1; SLC7A11 have higher correlations with CD160, CD27, CD44, CD70, CD80.
Spatial transcriptome and single-cell sequencing analyses were used to evaluate the expression pattern and localization of GPX3, SLC7A11, and CDKN2A in CRC. The single-cell correlation technique enables the identification of cell types and tissue heterogeneity. Spatial transcriptomics enables the simultaneous acquisition of gene expression and spatial information, but not at the single-cell level. Integrating single-cell detection with spatial transcriptome detection enhances the precision of 3D state characterization, thereby expanding potential avenues for future investigations. First, spatial transcriptome analysis was conducted using a publicly available database. Seventeen clusters were identified using UMAP analysis (Figure 5A-B). Feature plots and spatial feature plots demonstrated the high expression of CDKN2A and SLC7A11 in CRC tissues. GPX3 was not found in the public database. The spatial expression and distribution of CDKN2A and SLC7A1 in CRC tissues were mapped using this database (Figure 5C, D).
Subsequently, single-cell sequencing was performed to determine the expression of these three genes across distinct cell types. Seven primary cell clusters were identified using UMAP analysis (Figure 5E). No clusters were exclusively derived from a single individual. An additional dot plot is shown in Figure 5F. The seven single-cell subcohorts were identified as “Epithelial cells,” “iPS cells,” “Fibroblasts,” “B cells,” “T cells,” “Marcrophage,” and “Endothelial cells,” in sequential order. The results indicated that SLC7A11 was predominantly enriched in subpopulations of Fibroblasts, Epithelial cells, and iPS cells. CDKN2A was predominantly enriched in subpopulations of Fibroblasts and Epithelial cells. However, the distribution of GPX3 expression was diffuse, and the level of expression was low. Furthermore, we evaluated the density profiles of the three genes using single-cell sequencing (Figure 5G-I). Pseudotime analysis was used to depict gene expression patterns during the initial stages of CRC development (Figure 5J). Based on pseudo-temporal data, we investigated the correlation between the expression of the three genes and the early stage of CRC, suggesting that the expression of SLC7A11 and CDKN2A may be significantly associated with the development and onset of early-stage CRC (Figure 5K).
To facilitate the clinical practice, we further converted the complex mathematical model into a nomogram. A nomogram was constructed to predict the 1-, 3-, and 5-year survival of patients with CRC using factors such as T stage, N stage, M stage, age, CEA level, CDKN2A expression, GPX3 expression, and SLC7A11 expression (Figure 6A). Calibration curves were created to predict the clinical outcomes for the 1-, 3-, and 5-year periods with reasonable accuracy (Figure 6B-D). The predictions made by the model were close to the actual outcomes. Subsequently, the ROC curves and area under the curve (AUC) values of the 1-, 3-, and 5-year OS were predicted (Figure 6E-G). The ROC curves are significantly above the diagonal dashed lines, with AUC values above 0.7. The above results showed that this nomogram had significant high prognostic prediction ability for the survival outcomes of patients with CRC.
The validation of CDKN2A, GPX3 and SLC7A11 mRNA expression in patients with CRC. A-C: The validation of CDKN2A, GPX3 and SLC7A11 mRNA expression in several CRC and normal cell lines. D-E: The validation of CDKN2A, GPX3 and SLC7A11 mRNA expression in CRC and normal tissues. (Statistical significances were calculated using independent sample t-test. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001).
By the way, we analyzed the DNA methylation levels of GPX3, CDKN2A, and SLC7A11 in CRC. DNA methylation heatmaps for the CpG sites in these three genes were also plotted. The CpG sites of CDKN2A (cg18849169, cg26891370, cg04201367, cg11109721, and cg22005145), GPX3 (cg12840719 and cg04026675), and SLC7A11 (cg06623625, cg02734904, cg06206831, cg21877274, cg01309945, cg04474257, and cg24869834) exhibited elevated methylation levels in CRC (Supplementary figure 2A-C). Methylation analysis of the GPX3, CDKN2A, and SLC7A11 promoter regions, along with the corresponding expression data, was further conducted using publicly available datasets from the UALCAN database. The results revealed that the methylation levels of CDKN2A and GPX3 promoters were significantly higher in CRC tissues than in normal tissues. In contrast, SLC7A11 exhibited the opposite trend (Supplementary figure 2D-F). We lastly predicted miRNAs and TFs for three genes using TargetScan Human 7.2 and GeneCards databases. Based on these 28 TFs and three miRNAs, the TFs-mRNA-miRNA network was constructed, as shown in Supplementary figure 2G. We may explore the deep regulation mechanism of the three selected ferroptosis-related genes in the development of CRC based on this result.
Associations between clinicopathologic features, tumor immunity and three selected ferroptosis-related genes in the TCGA-COAD dataset. A: ROC and AUC of CDKN2A, GPX3 and SLC7A11 to diagnose death for patients with CRC in the TCGA-COAD database. B-D: The relationships between GPX3, CDKN2A and SLC7A11 expression and infiltrating immune cells in the TCGA-COAD database. E-G: The relationships between GPX3, CDKN2A and SLC7A11 expression and CRC immune checkpoint genes of TCGA-COAD database.
Single cell localization of three selected ferroptosis-related genes and spatial transcriptome validation. A-B: Seventeen clusters were identified by t-SNE and uniform manifold approximation and projection (UMAP). C-D: Feature plots and spatial feature plots were utilized to illustrate the distribution expression of CDKN2A and SLC7A11. E: Seven major cell clusters were identified via UMAP. F: Dot plot of cells proportion in the respective cluster expressing selected marker genes (dot size), and average expression (color scale). G-I: UMAP plot of GPX3, SLC7A11 and CDKN2A expression across all cell clusters. J: Pseudo-time trajectory and modules of genes whose expression varied with pseudo-time. K: Pseudo-time analysis was used to plot the early-stage trajectories of CRC.
Ferroptosis is primarily induced by the abnormal metabolism of amino acids, lipids, or iron, resulting in the excessive accumulation of L-ROS. This phenomenon can occur in both normal and tumor cells. The intracellular levels of iron and GSH in cancerous cells compared to those in normal cells can be utilized for targeted therapy based on ferroptosis [35, 36]. Ferroptosis has demonstrated potential in the development of CRC and even in CRC therapy [37-39]. However, the actual regulation mechanisms are still not very clear. Therefore, identifying the novel driving genes of ferroptosis can significantly contribute to the exploration of CRC development and the clinical treatment of CRC. In this study, we collected gene expression and clinicopathological data from TCGA dataset. Initially, 24 DEGs related to ferroptosis were identified. A prognostic model was developed using LASSO and Cox regression analyses, which identified three genes, GPX3, CDKN2A, and SLC7A11 as novel potential prognostic biomarkers.
Prognostic nomogram for the 1-, 3-, and 5-year survival times of patients with CRC. A: Independent risk factors screened by multivariate Cox regression in the TCGA cohort were integrated into the nomogram model. B-D: Calibration curves of the nomogram for predicting 1-, 3-, and 5-year OS in the TCGA cohort. E-G: ROC curves and AUC values of the nomogram for predicting 1-, 3-, and 5-year OS.
CDKN2A mainly encodes two proteins: p14 and p16. p16 functions as a tumor suppressor by inhibiting cell cycle progression from the G1 to the S phase, thereby slowing down cell division. p14, a splice variant of CDKN2A, is significantly associated with cancer poor prognosis [40]. CDKN2A has multifaceted involvement in cancer, encompassing its regulatory function and effects on different cancer types. However, there are few studies which explore the role of ferroptosis-related CDKN2A in the progression of CRC. Specifically, Xing X et al. found that the level of methylation of CDKN2A gene was an effective predictor of poor survival outcomes of patients with CRC [41]. Also, Shi WK et al. found that CDKN2A has a pro-tumor effect in CRC. In CRC cell lines, silence the expression of CDKN2A can significantly inhibit the proliferation, induce cell apoptosis, arrest cell cycle and affect the process of epithelial-mesenchymal transition [42]. Furthermore, Dong Y et al. found that high expression of CDKN2A can not only contribute to a poor survival outcome of patients with CRC, but also may guide PD-1 mediated immunotherapy for patients with CRC [43]. In short, our data validated a significant correlation between the contributing role of CDKN2A and the development of CRC. However, the specifical molecular mechanisms of CDKN2A in the malignant progression of CRC are still unclear and need further investigation. More validated experiments are required to be performed in the future.
GPX3 could catalyze the reduction of organic hydrogen peroxide and hydrogen peroxide (H2O2) through the action of glutathione. Research indicates that the dysfunction of GPX3 expression in tumor cells is linked to a poor prognosis for patients with cancer and chemotherapeutic resistance [44], including endometrial adenocarcinoma [45], lung cancer [46], gastric cancer [47], prostate cancer [48], cervical cancer [49], and thyroid cancer [50]. Among these kinds of solid cancer, GPX3 can not only play a cancer-promoting role, but also play a cancer-suppressing effect. However, there few researches about the role of GPX3 in the development of CRC. Pelosof L et al. found that GPX3 promoter methylation may predict platinum sensitivity in patients with CRC [51]. Haug U et al. has also found that genetic variability in GPX3 may contribute to risk of rectal cancer but not of colon cancer [52]. In conclusion, our research is the first study to explore the expression and prognostic effect of GPX3 in patients with CRC. However, the specific function and molecular mechanism of GPX3 in the development of CRC still need to be further explored.
Many studies have shown that SLC7A11 is one of the key regulators of ferroptosis [11, 53-55]. Its main biological role is to control the transmembrane exchange of cystine and glutamic acid, and to mediate the entry of cystine into cells, thus providing raw material for intracellular GSH synthesis. Then, GSH may convert lipid peroxidation into non-toxic fatty alcohols, which protect cells from oxidative stress [56, 57]. The expression of SLC7A11 is dysregulated in many kinds of tumor cells, such as lung cancer, liver cancer, breast cancer and CRC, and plays an important role in the occurrence and development of tumor. Inhibition of SLCA11 expression can affect the growth and drug sensitivity of many kinds of tumors. However, in different condition, the effect of SLC7A11 on the development of cancers may be different [58].
Some researchers found that SLC7A11 was overexpressed in various types of tumors, and the silence of SLC7A11 expression may increase the level of intracellular reactive oxygen species (ROS) and induce death of tumor cells [59, 60]. However, other researchers observed that inhibition of SLC7A11 expression may contribute to the growth and proliferation of tumor cells. They found that tumor cells having high level of SLC7A11 expression were more sensitive to the level of glucose and glutamine restrictions than these cells with low level of SLC7A11 expression, which may due to the high nutrient dependence of tumor cells. Meanwhile, inhibition of SLC7A11 expression may contribute to tumor cell's adaptation to a hypoxic tumor environment, as a result, facilitating tumor development [61]. The down-regulated expression of SLC7A11 can also increase the expression of multidrug-resistant proteins which can result in drug resistance [62]. Just like our study, we found that the overall expression of SLC7A11 was lower in CRC than in normal tissues. Low expression of SLC7A11 may contribute to a worse survival outcome for patients with CRC. CRC cells with SLC7A11 low expression may more adaptive to the hypoxic tumor microenvironment and have more malignant phenotypes. However, there were still other studies which found that the inhibition of SLC7A11 expression may promote ferroptosis of CRC cells, then inhibit the malignant development of CRC and even CRC stem cells [63-67]. From what we have discussed above, we can conclude that the actual action mechanism of SLC7A11 in the malignant development of CRC may be very complex, we need perform more experiments to make it clear in the future.
In this study, a risk-score signature and a prognostic prediction nomogram were also constructed to accurately predict the 1-, 3-, and 5-year OS rates. The established prediction model based on Lasso-Cox regression still had large predictive advantages. To validate our results again, OS survival curves were plotted and found that patients with different risk had significantly different outcomes. Furthermore, the calibration plot of the nomogram showed favorable consistency between the prediction of the nomogram and actual observations. Also, the time-dependent AUC (>0.7) indicated satisfactory discriminative ability of the nomogram. Just like we have stated, ferroptosis is a prominent subject in tumor research. However, the regulatory mechanisms linking tumor immunity and ferroptosis remain unclear. In our study, KEGG and GO analyses of DEGs in high-risk and low-risk groups reveal immune-related biological pathways and functions. High-risk patients demonstrate increased expression of immune responses and tumor progression pathways. A potential correlation between ferroptosis and tumor immune function can be deduced. Immunological checkpoint analysis revealed a positive correlation between risk score and immune checkpoint protein expression levels. To further investigate the expression of ferroptosis-related genes, we used spatial transcriptome validation and single-cell localization to analyze the expression patterns of the three genes in CRC. Our findings indicate that CDKN2A and SLC7A11 are primarily expressed in fibroblast subpopulations. Pseudotime analysis predicts that the CDKN2A and SLC7A11 are closely related to the early stage of CRC. Collectively, we proposed that these three ferroptosis-related genes have significant implications for the development, prognosis, and clinical relevance of CRC, which provides a promising avenue for investigating the relationship between ferroptosis and CRC.
However, this study has some limitations. First, the primary focus of this study was the predictive methodology. The relevant gene expression in CRC cell lines and clinical specimens has been validated, but the related functional analysis requires further in vitro and in vivo assays. Second, further analysis is required to explore the conclusions drawn from single-cell and spatial transcriptome analyses. Third, further investigation is required to explore the action mechanisms of GPX3, CDKN2A, and SLC7A11 in the development of CRC, including their upstream transcription factors and interacting proteins.
Supplementary figures.
This work was supported by the National Natural Science Cultivation Foundation of China of Liaoning Cancer Hospital (2021-ZLLH-18), Doctoral Start-up Foundation of Liaoning Province (2023-BS-048) and “5310” Talent Strategy Project Funding Plan of Liaoning Cancer Hospital & Institute (To Dr. Xudong Zhu).
The authors have declared that no competing interest exists.
1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7-30
2. Siegel RL, Miller KD, Goding Sauer A, Fedewa SA, Butterly LF, Anderson JC. et al. Colorectal cancer statistics, 2020. CA Cancer J Clin. 2020;70:145-64
3. Morazan-Fernandez D, Mora J, Molina-Mora JA. In Silico Pipeline to Identify Tumor-Specific Antigens for Cancer Immunotherapy Using Exome Sequencing Data. Phenomics. 2023;3:130-7
4. Biller LH, Schrag D. Diagnosis and Treatment of Metastatic Colorectal Cancer: A Review. JAMA. 2021;325:669-85
5. Gao B, Li X, Li S, Wang S, Wu J, Li J. Pan-cancer analysis identifies RNA helicase DDX1 as a prognostic marker. Phenomics. 2022;2:33-49
6. Phipps O, Brookes MJ, Al-Hassi HO. Iron deficiency, immunology, and colorectal cancer. Nutr Rev. 2021;79:88-97
7. Koleini N, Shapiro JS, Geier J, Ardehali H. Ironing out mechanisms of iron homeostasis and disorders of iron deficiency. J Clin Invest. 2021 131
8. Qu L, He X, Tang Q, Fan X, Liu J, Lin A. Iron metabolism, ferroptosis, and lncRNA in cancer: knowns and unknowns. J Zhejiang Univ Sci B. 2022;23:844-62
9. Khan A, Singh P, Srivastava A. Iron: Key player in cancer and cell cycle?. J Trace Elem Med Biol. 2020;62:126582
10. Wang Y, Yu L, Ding J, Chen Y. Iron Metabolism in Cancer. Int J Mol Sci. 2018 20
11. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE. et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149:1060-72
12. Liu S, Yue M, Lu Y, Wang Y, Luo S, Liu X. et al. Advancing the frontiers of colorectal cancer treatment: harnessing ferroptosis regulation. Apoptosis. 2023
13. Dixon SJ, Stockwell BR. The role of iron and reactive oxygen species in cell death. Nat Chem Biol. 2014;10:9-17
14. Xue X, Ramakrishnan SK, Weisz K, Triner D, Xie L, Attili D. et al. Iron Uptake via DMT1 Integrates Cell Cycle with JAK-STAT3 Signaling to Promote Colorectal Tumorigenesis. Cell Metab. 2016;24:447-61
15. Zhang Q, Deng T, Zhang H, Zuo D, Zhu Q, Bai M. et al. Adipocyte-Derived Exosomal MTTP Suppresses Ferroptosis and Promotes Chemoresistance in Colorectal Cancer. Adv Sci (Weinh). 2022;9:e2203357
16. Chen Y, Li H. Prognostic and Predictive Models for Left- and Right- Colorectal Cancer Patients: A Bioinformatics Analysis Based on Ferroptosis-Related Genes. Front Oncol. 2022;12:833834
17. Zhu X, Li S. Ferroptosis, Necroptosis, and Pyroptosis in Gastrointestinal Cancers: The Chief Culprits of Tumor Progression and Drug Resistance. Adv Sci (Weinh). 2023;10:e2300824
18. Sui X, Zhang R, Liu S, Duan T, Zhai L, Zhang M. et al. RSL3 Drives Ferroptosis Through GPX4 Inactivation and ROS Production in Colorectal Cancer. Front Pharmacol. 2018;9:1371
19. Peng B, Peng J, Kang F, Zhang W, Peng E, He Q. Ferroptosis-Related Gene MT1G as a Novel Biomarker Correlated With Prognosis and Immune Infiltration in Colorectal Cancer. Front Cell Dev Biol. 2022;10:881447
20. Yang C, Zhang Y, Lin S, Liu Y, Li W. Suppressing the KIF20A/NUAK1/Nrf2/GPX4 signaling pathway induces ferroptosis and enhances the sensitivity of colorectal cancer to oxaliplatin. Aging (Albany NY). 2021;13:13515-34
21. Subramaniana A. et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. 2005.
22. Liberzon A, Birger C, Thorvaldsdottir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1:417-25
23. Zhang S, Chang W, Wu H, Wang YH, Gong YW, Zhao YL. et al. Pan-cancer analysis of iron metabolic landscape across the Cancer Genome Atlas. J Cell Physiol. 2020;235:1013-24
24. Mou Y, Zhang Y, Wu J, Hu B, Zhang C, Duan C. et al. The Landscape of Iron Metabolism-Related and Methylated Genes in the Prognosis Prediction of Clear Cell Renal Cell Carcinoma. Front Oncol. 2020;10:788
25. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550
26. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139-40
27. McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288-97
28. 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:e47
29. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O. et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523
30. Blanche P, Dartigues JF, Jacqmin-Gadda H. Estimating and comparing time-dependent areas under receiver operating characteristic curves for censored event times with competing risks. Stat Med. 2013;32:5381-97
31. Alba AC, Agoritsas T, Walsh M, Hanna S, Iorio A, Devereaux PJ. et al. Discrimination and Calibration of Clinical Prediction Models: Users' Guides to the Medical Literature. JAMA. 2017;318:1377-84
32. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q. et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48:W509-W14
33. Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH. et al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019;35:i436-i45
34. Lee OJ, Schneider-Stock R, McChesney PA, Kuester D, Roessner A, Vieth M. et al. Hypermethylation and loss of expression of glutathione peroxidase-3 in Barrett's tumorigenesis. Neoplasia. 2005;7:854-61
35. Mentoor I, Engelbrecht AM, van Jaarsveld PJ, Nell T. Chemoresistance: Intricate Interplay Between Breast Tumor Cells and Adipocytes in the Tumor Microenvironment. Front Endocrinol (Lausanne). 2018;9:758
36. Stockwell BR, Friedmann Angeli JP, Bayir H, Bush AI, Conrad M, Dixon SJ. et al. Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell. 2017;171:273-85
37. Bozatzi P, Sapkota GP. The FAM83 family of proteins: from pseudo-PLDs to anchors for CK1 isoforms. Biochem Soc Trans. 2018;46:761-71
38. Hassannia B, Vandenabeele P, Vanden Berghe T. Targeting Ferroptosis to Iron Out Cancer. Cancer Cell. 2019;35:830-49
39. Wei R, Zhao Y, Wang J, Yang X, Li S, Wang Y. et al. Tagitinin C induces ferroptosis through PERK-Nrf2-HO-1 signaling pathway in colorectal cancer cells. Int J Biol Sci. 2021;17:2703-17
40. Kingsbury JM, Shamaprasad N, Billmyre RB, Heitman J, Cardenas ME. Cancer-associated isocitrate dehydrogenase mutations induce mitochondrial DNA instability. Hum Mol Genet. 2016;25:3524-38
41. Xing X, Cai W, Shi H, Wang Y, Li M, Jiao J. et al. The prognostic value of CDKN2A hypermethylation in colorectal cancer: a meta-analysis. Br J Cancer. 2013;108:2542-8
42. Shi WK, Li YH, Bai XS, Lin GL. The Cell Cycle-Associated Protein CDKN2A May Promotes Colorectal Cancer Cell Metastasis by Inducing Epithelial-Mesenchymal Transition. Front Oncol. 2022;12:834235
43. Dong Y, Zheng M, Wang X, Yu C, Qin T, Shen X. High expression of CDKN2A is associated with poor prognosis in colorectal cancer and may guide PD-1-mediated immunotherapy. BMC Cancer. 2023;23:1097
44. Chen B, Rao X, House MG, Nephew KP, Cullen KJ, Guo Z. GPx3 promoter hypermethylation is a frequent event in human cancer and is associated with tumorigenesis and chemotherapy response. Cancer Lett. 2011;309:37-45
45. Falck E, Karlsson S, Carlsson J, Helenius G, Karlsson M, Klinga-Levan K. Loss of glutathione peroxidase 3 expression is correlated with epigenetic mechanisms in endometrial adenocarcinoma. Cancer Cell Int. 2010;10:46
46. An BC, Choi YD, Oh IJ, Kim JH, Park JI, Lee SW. GPx3-mediated redox signaling arrests the cell cycle and acts as a tumor suppressor in lung cancer cell lines. PLoS One. 2018;13:e0204170
47. Cai M, Sikong Y, Wang Q, Zhu S, Pang F, Cui X. Gpx3 prevents migration and invasion in gastric cancer by targeting NFsmall ka, CyrillicB/Wnt5a/JNK signaling. Int J Clin Exp Pathol. 2019;12:1194-203
48. Yu YP, Yu G, Tseng G, Cieply K, Nelson J, Defrances M. et al. Glutathione peroxidase 3, deleted or methylated in prostate cancer, suppresses prostate cancer growth and metastasis. Cancer Res. 2007;67:8043-50
49. Zhang X, Zheng Z, Yingji S, Kim H, Jin R, Renshu L. et al. Downregulation of glutathione peroxidase 3 is associated with lymph node metastasis and prognosis in cervical cancer. Oncol Rep. 2014;31:2587-92
50. Zhao H, Li J, Li X, Han C, Zhang Y, Zheng L. et al. Silencing GPX3 Expression Promotes Tumor Metastasis in Human Thyroid Cancer. Curr Protein Pept Sci. 2015;16:316-21
51. Pelosof L, Yerram S, Armstrong T, Chu N, Danilova L, Yanagisawa B. et al. GPX3 promoter methylation predicts platinum sensitivity in colorectal cancer. Epigenetics. 2017;12:540-50
52. Haug U, Poole EM, Xiao L, Curtin K, Duggan D, Hsu L. et al. Glutathione peroxidase tagSNPs: associations with rectal cancer but not with colon cancer. Genes Chromosomes Cancer. 2012;51:598-605
53. Koppula P, Zhuang L, Gan B. Cystine transporter SLC7A11/xCT in cancer: ferroptosis, nutrient dependency, and cancer therapy. Protein Cell. 2021;12:599-620
54. Chen X, Li J, Kang R, Klionsky DJ, Tang D. Ferroptosis: machinery and regulation. Autophagy. 2021;17:2054-81
55. Zhang Y, Shi J, Liu X, Feng L, Gong Z, Koppula P. et al. BAP1 links metabolic regulation of ferroptosis to tumour suppression. Nat Cell Biol. 2018;20:1181-92
56. He F, Zhang P, Liu J, Wang R, Kaufman RJ, Yaden BC. et al. ATF4 suppresses hepatocarcinogenesis by inducing SLC7A11 (xCT) to block stress-related ferroptosis. J Hepatol. 2023;79:362-77
57. Yan Y, Teng H, Hang Q, Kondiparthi L, Lei G, Horbath A. et al. SLC7A11 expression level dictates differential responses to oxidative stress in cancer cells. Nat Commun. 2023;14:3673
58. Li S, Lu Z, Sun R, Guo S, Gao F, Cao B. et al. The Role of SLC7A11 in Cancer: Friend or Foe?. Cancers (Basel). 2022 14
59. Balza E, Castellani P, Delfino L, Truini M, Rubartelli A. The pharmacologic inhibition of the xc- antioxidant system improves the antitumor efficacy of COX inhibitors in the in vivo model of 3-MCA tumorigenesis. Carcinogenesis. 2013;34:620-6
60. Cramer SL, Saha A, Liu J, Tadi S, Tiziani S, Yan W. et al. Systemic depletion of L-cyst(e)ine with cyst(e)inase increases reactive oxygen species and suppresses tumor growth. Nat Med. 2017;23:120-7
61. Briggs KJ, Koivunen P, Cao S, Backus KM, Olenchock BA, Patel H. et al. Paracrine Induction of HIF by Glutamate in Breast Cancer: EglN1 Senses Cysteine. Cell. 2016;166:126-39
62. Li SJ, Cao B, Lu ZY, Sun RB, Guo SH, Xie Y. et al. Cystine supplementation rebalances the redox homeostasis of microenvironment in non-small cell lung cancer cells and reverses their resistance to docetaxel. Acta Pharmacol Sin. 2021;42:2132-43
63. Martino E, Balestrieri A, Aragona F, Bifulco G, Mele L, Campanile G. et al. MiR-148a-3p Promotes Colorectal Cancer Cell Ferroptosis by Targeting SLC7A11. Cancers (Basel). 2023 15
64. Lei S, Chen C, Han F, Deng J, Huang D, Qian L. et al. AMER1 deficiency promotes the distant metastasis of colorectal cancer by inhibiting SLC7A11- and FTL-mediated ferroptosis. Cell Rep. 2023;42:113110
65. Saini KK, Chaturvedi P, Sinha A, Singh MP, Khan MA, Verma A. et al. Loss of PERK function promotes ferroptosis by downregulating SLC7A11 (System Xc(-)) in colorectal cancer. Redox Biol. 2023;65:102833
66. Guo S, Zhao W, Zhang W, Li S, Teng G, Liu L. Vitamin D Promotes Ferroptosis in Colorectal Cancer Stem Cells via SLC7A11 Downregulation. Oxid Med Cell Longev. 2023;2023:4772134
67. Elrebehy MA, Abdelghany TM, Elshafey MM, Gomaa MH, Doghish AS. miR-509-5p promotes colorectal cancer cell ferroptosis by targeting SLC7A11. Pathol Res Pract. 2023;247:154557
Corresponding authors: Xudong Zhu, Department of General Surgery, Cancer Hospital of Dalian University of Technology, Cancer Hospital of China Medical University, Liaoning Cancer Hospital & Institute, Liaoning Provincial Key Laboratory of Precision Medicine for Malignant Tumors, Shenyang, Liaoning 110042, P.R. China; Markey Cancer Center, University of Kentucky, Lexington, KY 40536, United States of America. Telephone: +86 13354204706; E-mail: xdzhuedu.cn, XudongZhuedu. Tong Zhu, Department of Breast Surgery, Panjin Central Hospital, Panjin, Liaoning 124010, P.R. China. Telephone: +86 15909835580; E-mail: jimmiezhunet. Guohua Zhao, Department of General Surgery, Cancer Hospital of Dalian University of Technology, Cancer Hospital of China Medical University, Liaoning Cancer Hospital & Institute, Shenyang, Liaoning 110042, P.R. China. Telephone: +86 18900917009; E-mail: zhaoguohua1975com.
Received 2023-10-23
Accepted 2024-3-25
Published 2024-4-22