Abstract
Background: Ferroptosis is a potential target for cancer therapy, and lncRNAs can also affect ferroptosis by regulating related genes. The pathogenesis of clear cell renal cell carcinoma (ccRCC) regarding the regulation of ferroptosis by lncRNAs is still unknown.
Methods: We constructed a risk model based on data in ccRCC patients obtained from the TCGA database and validated the diagnostic and prognostic value of the model. In addition, immune function and immune checkpoint variability analysis validated the association of ferroptosis with ccRCC tumor immunity.
Results: The characteristics of ferroptosis-related lncRNAs (FRLs) were significantly correlated with the prognosis of ccRCC patients. The prognostic characteristics of FRLs were independent prognostic factors in ccRCC patients. Gene function in the high-risk group was associated with oxygen metabolic processes and immune pathways. Immune checkpoint variability analysis showed that HAVCR2, NRP1, and HHLA2 were upregulated in the low-risk group, while CD44, TNFRSF18, TNFSF14, TNFRSF8, CD276, and TNFRSF25 were upregulated in the high-risk group.
Conclusions: The prognostic characteristics of FRLs can effectively predict the prognosis of ccRCC patients and provide a new direction for clinical treatment.
Introduction
Clear cell renal cell carcinoma (ccRCC) is a common tumor of the urinary tract with a high mortality rate (Moch et al., 2016; Sung et al., 2021). Currently, immunotherapy is the latest treatment for ccRCC, and immune checkpoint inhibitors (ICIs) such as anti-PD-1/PDL1 inhibitors have been gradually used in clinical practice (Gonzalez et al., 2018; Postow et al., 2018; Chen et al., 2019; Khan et al., 2019; Zhang et al., 2019; Buonerba et al., 2020). Unfortunately, some patients are not sensitive to ICIs (Carosella et al., 2015; Tan et al., 2020). Therefore, it is critical to find biomarkers and potential target drugs for the progression and prognosis of ccRCC.
Ferroptosis is a novel form of cell death caused by iron-dependent oxidative damage (Stockwell et al., 2017). The pathological mechanism is the failure of glutathione peroxidase (GPX4), which leads to the accumulation of reactive oxygen species (ROS) on membrane (Dixon et al., 2012). Studies have shown that ferroptosis involved in tumors, degenerative diseases, and ischemia-reperfusion injury (Guiney et al., 2017; Hassannia et al., 2019; Li et al., 2019; Yan and Zhang, 2019). Interestingly, lncRNAs play a pivotal role in ferroptosis. Silencing lncRNA MEG8 can induce ferroptosis of hemangioma endothelial cells and inhibit cell proliferation (Mao et al., 2018). LncRNAOIP5-AS1 inhibited ferroptosis and promoted cell proliferation in prostate cancer cells exposed to cadmium (Lv et al., 2022). Ketamine can induce ferroptosis of hepatoma cells by targeting lncRNAPVT1/miR-214-3p/GPX4 (Wang et al., 2019; Bohosova et al., 2021). In glioma cells, upregulation of lncrNATMEM161B-AS1 promotes ferroptosis by sponging miR-27a-3p (Kapoor et al., 2021). In addition, lncRNARP11-89 inhibited ferroptosis and promoted the development of bladder cancer by sponging miR-129-5p (Yang et al., 2020).
The specific regulatory role of lncRNAs in ferroptosis remains to be further investigated. Moreover, how to regulate the aberrant expression of lncRNAs with ferroptosis is an urgent question to be investigated. Thus, we attempted to discover the regulatory role of lncRNAs in ferroptosis and provide the theoretical basis for its application in ccRCC.
Methods and Materials
Data Download and Study Design
The gene expression and clinical data were obtained from The Cancer Genome Atlas (https://www.cancer.gov/tcga), the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo), and FerrDb (http://www.zhounan.org/ferrdb/) databases. Immune-related clinical data were obtained from Timer (http://timer.comp-genomics.org). OS and Disease-Free Survival (DFS) data were obtained from the GEPIA database (http://gepia.cancer-pku.cn/).
Identification of FRGs and FRLs
A correlation coefficient filter criterion of 0.6 and a p-value filter criterion of 0.001 were used to screen for FRLs. Based on the expression data of FRLs with the corresponding survival data, prognosis-related FRLs were screened, and these FRLs would be involved in the prognostic model. The “limma” (Ritchie et al., 2015) packages were used to screen co-expression analysis of FRGs and lncRNA with p < 0.05 and log2FC > 2. The “limma” package was performed to identify FRGs (p < 0.05 and log2FC > 1) and identify expression level of FRG (Score Filter = 0.4, p-Value Filter = 0.001). The “limma” package was performed to identify FRLs (p < 0.05 and log2FC > 1).
Co-Expression Network Analysis and Nomogram
COX models were constructed for lncRNAs related to survival time, survival status, and prognosis of patients. Based on the model equations, the risk scores of patients were calculated and patients were divided into high-risk and low-risk groups according to the median risk scores. The FRGs and FRLs were inputted into the STRING online tool (https://string-db.org) to construct the PPI network. Then, Cytoscape software (Shannon et al., 2003) was used to analyze the PPI network. And Cytoscape software drew the figure of FRGs and FRLs co-expression network. The “regplot” and “survival” packages were utilized to perform the nomogram.
Cox Regression Analysis and Construction of a Proportional Hazards Model
The “limma” package was performed to identify the survival time. The “survival” and “glmnet” and “survminer” packages (Friedman et al., 2010) were utilized to perform univariate and multivariate Cox regression analyses. Through integrated Cox analysis, key FRLs were screened to construct the risk model. The risk curve was finished by the “pheatmap” treatment. The “survival”, “glmnet” and “survminer” packages were utilized to perform independent prognostic analysis.
Survival Analysis and Decision Curve Analysis
Based on the median gene expression/risk score, ccRCC patients are classified into high and low-risk groups. Then, survival curves of OS and DFS were drawn by the “survival” and “survminer” packages in R (v.4.1.1) and GraphPad software (version 8.0). A p < 0.05 was considered statistically significant. Moreover, the “survival”, “survminer”, and “timeROC” packages were used to generate a time-dependent ROC curve to evaluate the predictive value of the risk model. The “survival” and “ggDCA” packages drew the decision curve.
Gene Set Enrichment Analysis
All patients were divided into high and low-risk groups according to the median gene expression/risk score. GSEA (Subramanian et al., 2005) was performed to discover potential mechanisms and downstream signaling pathways.
Kyoto Encyclopedia of Genes and Genomes Pathway Enrichment Analysis and Gene Ontology Analysis
The “colorspace”, “stringi”, “clusterProfiler” (Wu T. et al., 2021), “DOSE”, “org.Hs.eg.db”, and “enrichplot” packages were used to conduct KEGG enrichment analysis of FRGs and FRLs. The “colorspace”, “stringi”, “clusterProfiler”, “DOSE”, “org.Hs.eg.db”, and “enrichplot” packages were used to conduct gene ontology analysis of FRGs and FRLs. The results were visualized by the “ggplot2” package in the R program. A p < 0.05 was selected as the cut-off point.
Immunohistochemistry
IHC staining data are obtained from the Human Protein Atlas website (HPA, https://www.proteinatlas.org), a database based on proteomic, transcriptomic, and systems biology data that can map tissues, cells, and organs. IHC staining of NRP1, HAVCR2, HHLA2CD44, TNFRSF18, TNFSF14, TNFRSF8, and CD276 in tumor tissues and normal tissues were obtained from HPA.
Statistics
All data analyses were performed using the R platform or GraphPad Prism 8.0. The FRLs of different groups were measured using Kaplan-Meier log-rank test method. The FRLs between high and low-risk groups were determined using the “limma” R package. p < 0.05 was considered statistically significant.
Results
Differential Expression Analysis and Functional Enrichment Analysis
We downloaded the transcriptomic data, along with the corresponding clinical data from the TCGA database. Meanwhile, we found all the genes related to ferroptosis from the FerrDb database. By taking intersections, we screened 230 FRGs differentially expressed in ccRCC. We screened out 230 differentially expressed genes (137 upregulated, 93 downregulated). Next, we screened out 76 ferroptosis-related genes (FRGs) (42 upregulated and 34 downregulated) and 1502 ccRCC-related lncRNAs (1265 upregulated and 237 downregulated), respectively.
Then, we performed GO and KEGG pathway enrichment analysis on FRGs. In biological processes (BP), FRGs are involved in chemical stress cell response, reactive oxygen species metabolism response, low oxygen response, oxygen level response, and oxidative stress cell response. Interestingly, in molecular function (MF), they are involved in ion binding, pyridoxal phosphate binding, and vitamin B6 binding processes (Figures 1A,B). The KEGG pathway analysis also revealed that FRGs were associated with multiple signaling pathways (such as HIF1 signaling, adipocytokine signaling, and PPAR signaling) and metabolic pathways (arachidonic acid metabolism, 2-oxocarboxylic acid metabolism, cysteine, and methionine metabolism) (Figures 1C,D).
FIGURE 1
Prognostic Signature
From the pre-acquired data, we screened ferroptosis-related lncRNAs (FRLs) with a correlation coefficient filter criterion of 0.6 and a p-value filter criterion of 0.001. We then analyzed the expression data of the screened FRLs with the corresponding survival data to select prognostically relevant FRLs, which will be involved in the prognostic model. We found that 67 FRLs were involved in the OS of ccRCC (Figure 2A). We took the obtained data of FRLs and performed COX model construction on patients’ survival time, survival status, and prognosis-related lncRNAs. We calculated the risk scores of the patients and divided them into high and low-risk groups with median value: risk score = (1.546 × LINC00894 expression level) + (2.5111 × AL139123.1 expression level) + (0.0119 × ASMTL−AS1 expression level) + (0.026 × AL157392.4 expression level) + (0.9667 × AL031714.1 expression level) + (−0.32554 × AC135050.3 expression level) + 0.644 × AP006621.2 expression level) + (3.271 × NARF−IT1 expression level) + (0.142 × YEATS2−AS1 expression level) + (0.375 × LINC02804 expression level) + (0.536 × AC024361.3 expression level) + (0.15983 × KIF1C−AS1 expression level) + (−0.5863 × PCED1B−AS1 expression level) + (1.2393 × UBE2Q1−AS1 expression level) + (−2.1882 × AL031705.1 expression level) + (1.4877 × AC005306.1 expression level) + (0.12394 × PTOV1−AS2 expression level) + (-0.0795 × AC114730.3 expression level) + (1.68734 × AC073487.1 expression level) + (−0.693 × AC104564.3 expression level) + (1.1236 × AC020907.4 expression level) + (1.81 × AC005387.2 expression level) + (−0.798 × AL513218.1 expression level) + (1.4485 × AC025766.1 expression level).
FIGURE 2
We randomly assign all samples in a 1:1 ratio and divided them into training and testing datasets to construct risk features. The clinical characteristics of all patients showed in Table 1. Subsequently, we developed risk curves for both groups based on the prognostic characteristics of the FRLs, and the risk curves showed the survival status and risk scores for each ccRCC sample (Figure 2B). The scatter plot and Kaplan-Meier survival curves showed significantly higher survival rates in the low-risk group than in the high-risk group (Figures 2C,D). These results show the excellent predictive performance of the risk score model in the training set. Similarly, the prognostic value of the risk signature is verified in both the validation set and the external independent test set (Figures 2E,F). Based on the risk score, we ranked the samples and extracted the expression data of FRLs with relevant clinical information to draw a clinical relevance heat map to show the expression levels of FRLs in high- and low-risk patients. The heat map showed that 23 of the 67 FRLs were the best candidates (Figure 2G).
TABLE 1
| Variables | Training Set | Validation Set | p-Value | ||
|---|---|---|---|---|---|
| (n = 268) | (n = 269) | ||||
| No. | % | No. | % | ||
| Age | 0.812 | ||||
| ≤65 | 174 | 64.9 | 178 | 66.2 | |
| >65 | 94 | 35.1 | 91 | 33.8 | |
| Stage | 0.875 | ||||
| I | 134 | 50.0 | 135 | 50.2 | |
| II | 29 | 10.8 | 28 | 10.4 | |
| III | 68 | 25.4 | 57 | 21.2 | |
| IV | 37 | 13.8 | 49 | 18.2 | |
| T stage | 0.926 | ||||
| T1 | 138 | 51.5 | 137 | 50.9 | |
| T2 | 30 | 11.2 | 39 | 14.5 | |
| T3 | 95 | 35.4 | 87 | 32.3 | |
| T4 | 5 | 1.9 | 6 | 2.3 | |
| N stage | 0.578 | ||||
| N0 | 134 | 50.0 | 146 | 54.3 | |
| N1 | 34 | 12.7 | 55 | 20.4 | |
| N2 | 20 | 7.5 | 45 | 16.7 | |
| N3 | 10 | 3.7 | 23 | 8.6 | |
| M stage | 0.688 | ||||
| M0 | 221 | 82.5 | 205 | 76.2 | |
| M1 | 47 | 17.5 | 64 | 23.8 | |
The association between risk score and patients’ clinical features in the training set.
Evaluation of Independent Prognostic Factors
Univariate analysis revealed that age (p < 0.001), stage (p < 0.001), and FRLs risk score (p < 0.001) were all independent prognostic factors for ccRCC patients, except gender (p = 0.748) (Figure 3A). Similarly, multivariate analysis showed a significant correlation between risk scores for FRLs and OS as well as other clinicopathological and prognostic characteristics except for gender (p = 0.694) (Figure 3B). In the training dataset, the ROC curves for years 1, 2, and 3 demonstrated that the FRLs risk score was strongly predictive: 1-year AUC (0.900), 2-year AUC (0.897), and 3-year AUC (0.912). Multivariate ROC curves showed that the FRLs risk score had better predictive performance than those clinicopathologic features, suggesting that the FRLs risk score is an independent predictor of survival in ccRCC patients (Figures 3C,D). The results of decision curve analysis (DCA) showed that the risk score of FRLs had better predictive power than other indicators (Figure 3E). In addition, In addition, similar results are also demonstrated in the validation set (Figures 3F–H).
FIGURE 3
Construction of Nomograph
We integrated the prognostic features and clinicopathologic factors (age, gender, grade, and stage) of FRLs to construct nomograms to accurately estimate the survival probability in ccRCC patients (Figure 3I). Then, we constructed a co-expression network using Cytoscape software (blue nodes represent FRLs and green nodes represent FRGs) (Figure 3J).
We also plotted a heat map on the clinical relevance of FRLs (horizontal coordinates indicate samples and vertical coordinates indicate FRLs). Among the different clinical characteristics of the high and low-risk groups, several clinical characteristics of grade, T, and M stage varied significantly in both groups (Figure 4A). We identified potential downstream signaling pathways of FRLs and found that it was associated with metabolism-related pathways by GSEA (Figures 4B,C).
FIGURE 4
Immune-Cell Correlation and Immune Function Analysis
We obtained immune-related clinical data from the Timer database and plotted the associated heat map (horizontal coordinates represent samples, vertical coordinates represent immune cells, and different colors represent predicted outcomes). The heat map showed that M1 macrophages, regulatory T cells, follicular helper T cells, and B cells, which are immune cells, had elevated expression in the high-risk group, while activated mast cells, neutrophils, and hematopoietic stem cells (HSC) had elevated expression in the low-risk group (Figure 5A). We then analyzed immune function variability and found that four immune functions in inflammation promotion, T cell co-inhibition, T cell co-stimulation, and type II IFN Response were actively involved in the high-risk group. These data suggest that ferroptosis is involved in ccRCC through the immune process pathway (Figure 5B).
FIGURE 5
Immunological Checkpoint Variability Analysis
We performed an immune checkpoint differential analysis and found that three genes, HAVCR2, NRP1, and HHLA2, were highly expressed in the low-risk group, and six genes, CD44, TNFRSF18, TNFSF14, TNFRSF8, CF276, and TNFRSF25, had increased expression in the high-risk group (Figure 5C). Then, we verified that the expression levels of CD44, CD276, HAVCR2, HHLA2, and NRP1 were statistically significant in ccRCC from the GEPIA database (Figure 5D). The expression levels of TNFSF14, TNFRSF18, NRP1, HAVCR2, and HHLA2 were significantly associated with the prognosis of ccRCC patients, and their OS and DFS were statistically significant (p < 0.05, Figure 6A). Meanwhile, we show the IHC staining of NRP1, HAVCR2, HHLA2CD44, TNFRSF18, TNFSF14, TNFRSF8, and CD276 in tumor tissues and normal tissues (Figure 6B).
FIGURE 6
Discussion
Ferroptosis is involved in the occurrence, development, and progression of tumors, which aroused our interest in studying its mechanism in ccRCC (Hambright et al., 2017; Stockwell et al., 2020; Ye et al., 2020). Tumor immunology, as an emerging research field, greatly complements and improves the research system of oncology (Galon and Bruni, 2020). After tumor cells escape the surveillance of the body’s immune system under the action of various factors, they can rapidly divide and proliferate in the body, accelerating tumor deterioration, which is immune escape (Dunn et al., 2002). Immune checkpoint blockade therapy based on programmed death receptors and their ligands enhances the aggressiveness of the host immune system against tumor cells by inhibiting the binding of them (Dunn et al., 2004; La-Beck et al., 2015; Postow et al., 2018; Khan et al., 2019). Non-specificity and drug resistance are the main problems faced by conventional treatments for ccRCC (Gai et al., 2020; Zou et al., 2021). And the role of lncRNA in cancer is receiving increasing attention from researchers (Inthagard et al., 2019; Resch et al., 2021).
We highlighted the lncRNAs associated with ferroptosis and ccRCC, which may yield promising results in cancer treatment. We validated the accuracy of the FRLs risk score in predicting the prognosis of ccRCC. Next, we tested the model with ROC curves and found that it applies to different situations. The prognostic characteristics of FRLs are an independent predictor of survival in patients with ccRCC and that they can be used as multiple indicators to diagnose or predict the onset or progression of ccRCC. We integrated the prognostic characteristics of FRLs with clinicopathological factors such as age, gender, grade, and stage and constructed a nomogram to predict the probability of survival in ccRCC patients.
LncRNAs regulate abnormal tumor lipid metabolism, thereby exerting oncogenic effects in tumorigenesis, affecting tumor cell proliferation, apoptosis, migration, invasion, and ferroptosis. LncRNA can also enable ferroptosis to produce an apoptosis-independent form of cell killing. Recent studies have found that the lncRNA SNHG12/SP1/CDCA3 axis promotes progression and sunitinib resistance in RCC, which provides a new therapeutic target for sunitinib-resistant RCC (Liu et al., 2020). In addition, lncRNA 00312 inhibits RCC proliferation and invasion and promotes apoptosis in RCC by inhibiting miR-34a-5p and overexpressing ASS1 (Zeng et al., 2020). SNHG17/miR-328-3p/H2AXaxis may be involved in RCC progression, providing a potential therapeutic target for RCC (Wu J. et al., 2021). Thus, we believe that targeting lncRNAs and combining both immune checkpoints and FRLs prognostic signature may create new opportunities for the treatment of ccRCC.
However, we still have several issues to resolve. First, identifying the most significant lncRNAs associated with ccRCC remains challenging. Our study lacks the necessary experimental validation related to lncRNA expression and other relevant experimental validation. In addition, lncRNAs are generally poorly conserved across species compared to protein-coding genes. Therapeutic strategies developed based on cellular and animal models are still a long way from clinical application and may require further research.
Conclusion
In conclusion, future research should reveal the relationship between each lncRNAs and ferroptosis and develop appropriate model systems to help the diagnosis, treatment, and prognosis of ccRCC. As the research continues, lncRNA-based therapeutic strategies will hopefully improve the prognosis of ccRCC patients.
Statements
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
YD and DL performed the statistical analyses and drafted the manuscript. YG and HZ supervised the statistical analyses. YX and YN participated in data analysis and interpretation and provided critical feedback. YD provided data for fine-mapping. All authors contributed to the article and approved the submitted version.
Funding
Supported by the National Natural Science Foundation of China (Nos. 81971371 and 82101671).
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/fgene.2022.912190/full#supplementary-material
References
1
BohosovaJ.KubickovaA.SlabyO. (2021). lncRNA PVT1 in the Pathogenesis and Clinical Management of Renal Cell Carcinoma. Biomolecules11 (5), 664. 10.3390/biom11050664
2
BuonerbaC.DolceP.IaccarinoS.ScafuriL.VerdeA.CostabileF.et al (2020). Outcomes Associated with First-Line Anti-PD-1/PD-L1 Agents vs. Sunitinib in Patients with Sarcomatoid Renal Cell Carcinoma: A Systematic Review and Meta-Analysis. Cancers12 (2), 408. 10.3390/cancers12020408
3
CarosellaE. D.PloussardG.LeMaoultJ.DesgrandchampsF. (2015). A Systematic Review of Immunotherapy in Urologic Cancer: Evolving Roles for Targeting of CTLA-4, PD-1/pd-L1, and HLA-G. Eur. Urol.68 (2), 267–279. 10.1016/j.eururo.2015.02.032
4
ChenL.PengT.LuoY.ZhouF.WangG.QianK.et al (2019). ACAT1 and Metabolism-Related Pathways Are Essential for the Progression of Clear Cell Renal Cell Carcinoma (ccRCC), as Determined by Co-expression Network Analysis. Front. Oncol.9, 957. 10.3389/fonc.2019.00957
5
DixonS. J.LembergK. M.LamprechtM. R.SkoutaR.ZaitsevE. M.GleasonC. E.et al (2012). Ferroptosis: an Iron-dependent Form of Nonapoptotic Cell Death. Cell.149 (5), 1060–1072. 10.1016/j.cell.2012.03.042
6
DunnG. P.BruceA. T.IkedaH.OldL. J.SchreiberR. D. (2002). Cancer Immunoediting: from Immunosurveillance to Tumor Escape. Nat. Immunol.3 (11), 991–998. 10.1038/ni1102-991
7
DunnG. P.OldL. J.SchreiberR. D. (2004). The Immunobiology of Cancer Immunosurveillance and Immunoediting. Immunity21 (2), 137–148. 10.1016/j.immuni.2004.07.017
8
FriedmanJ.HastieT.TibshiraniR. (2010). Regularization Paths for Generalized Linear Models via Coordinate Descent. J. Stat. Softw.33 (1), 1–22. 10.18637/jss.v033.i01
9
GaiC.LiuC.WuX.YuM.ZhengJ.ZhangW.et al (2020). MT1DP Loaded by Folate-Modified Liposomes Sensitizes Erastin-Induced Ferroptosis via Regulating miR-365a-3p/NRF2 axis in Non-small Cell Lung Cancer Cells. Cell. Death Dis.11 (9), 751. 10.1038/s41419-020-02939-3
10
GalonJ.BruniD. (2020). Tumor Immunology and Tumor Evolution: Intertwined Histories. Immunity52 (1), 55–81. 10.1016/j.immuni.2019.12.018
11
GonzalezH.HagerlingC.WerbZ. (2018). Roles of the Immune System in Cancer: from Tumor Initiation to Metastatic Progression. Genes Dev.32 (19-20), 1267–1284. 10.1101/gad.314617.118
12
GuineyS. J.AdlardP. A.BushA. I.FinkelsteinD. I.AytonS. (2017). Ferroptosis and Cell Death Mechanisms in Parkinson's Disease. Neurochem. Int.104, 34–48. 10.1016/j.neuint.2017.01.004
13
HambrightW. S.FonsecaR. S.ChenL.NaR.RanQ. (2017). Ablation of Ferroptosis Regulator Glutathione Peroxidase 4 in Forebrain Neurons Promotes Cognitive Impairment and Neurodegeneration. Redox Biol.12, 8–17. 10.1016/j.redox.2017.01.021
14
HassanniaB.VandenabeeleP.Vanden BergheT. (2019). Targeting Ferroptosis to Iron Out Cancer. Cancer Cell.35 (6), 830–849. 10.1016/j.ccell.2019.04.002
15
InthagardJ.EdwardsJ.RoseweirA. K. (2019). Immunotherapy: Enhancing the Efficacy of This Promising Therapeutic in Multiple Cancers. Clin. Sci. (Lond)133 (2), 181–193. 10.1042/CS20181003
16
KapoorP.ChowdhryA.BaggaD. K.BhargavaD.AishwaryaS. (2021). MicroRNAs in Oral Fluids (Saliva and Gingival Crevicular Fluid) as Biomarkers in Orthodontics: Systematic Review and Integrated Bioinformatic Analysis. Prog. Orthod.22 (1), 31. 10.1186/s40510-021-00377-1
17
KhanZ.HammerC.GuardinoE.ChandlerG. S.AlbertM. L. (2019). Mechanisms of Immune-Related Adverse Events Associated with Immune Checkpoint Blockade: Using Germline Genetics to Develop a Personalized Approach. Genome Med.11 (1), 39. 10.1186/s13073-019-0652-8
18
La-BeckN. M.JeanG. W.HuynhC.AlzghariS. K.LoweD. B. (2015). Immune Checkpoint Inhibitors: New Insights and Current Place in Cancer Therapy. Pharmacotherapy35 (10), 963–976. 10.1002/phar.1643
19
LiY.FengD.WangZ.ZhaoY.SunR.TianD.et al (2019). Ischemia-induced ACSL4 Activation Contributes to Ferroptosis-Mediated Tissue Injury in Intestinal Ischemia/reperfusion. Cell. Death Differ.26 (11), 2284–2299. 10.1038/s41418-019-0299-4
20
LiuY.ChengG.HuangZ.BaoL.LiuJ.WangC.et al (2020). Long Noncoding RNA SNHG12 Promotes Tumour Progression and Sunitinib Resistance by Upregulating CDCA3 in Renal Cell Carcinoma. Cell. Death Dis.11 (7), 515. 10.1038/s41419-020-2713-8
21
LvT.JinY.MiaoY.XuT.JiaF.FengH.et al (2022). LncRNA PVT1 Promotes Tumorigenesis of Glioblastoma by Recruiting COPS5 to Deubiquitinate and Stabilize TRIM24. Mol. Ther. - Nucleic Acids27, 109–121. 10.1016/j.omtn.2021.11.012
22
MaoC.WangX.LiuY.WangM.YanB.JiangY.et al (2018). A G3BP1-Interacting lncRNA Promotes Ferroptosis and Apoptosis in Cancer via Nuclear Sequestration of P53. Cancer Res.78 (13), 3484–3496. 10.1158/0008-5472.CAN-17-3454
23
MochH.CubillaA. L.HumphreyP. A.ReuterV. E.UlbrightT. M. (2016). The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs-Part A: Renal, Penile, and Testicular Tumours. Eur. Urol.70 (1), 93–105. 10.1016/j.eururo.2016.02.029
24
PostowM. A.SidlowR.HellmannM. D. (2018). Immune-Related Adverse Events Associated with Immune Checkpoint Blockade. N. Engl. J. Med.378 (2), 158–168. 10.1056/NEJMra1703481
25
ReschI.BruchbacherA.FrankeJ.FajkovicH.RemziM.ShariatS. F.et al (2021). Outcome of Immune Checkpoint Inhibitors in Metastatic Renal Cell Carcinoma across Different Treatment Lines. ESMO Open6 (4), 100122. 10.1016/j.esmoop.2021.100122
26
RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res.43 (7), e47. 10.1093/nar/gkv007
27
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res.13 (11), 2498–2504. 10.1101/gr.1239303
28
StockwellB. R.Friedmann AngeliJ. P.BayirH.BushA. I.ConradM.DixonS. J.et al (2017). Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell.171 (2), 273–285. 10.1016/j.cell.2017.09.021
29
StockwellB. R.JiangX.GuW. (2020). Emerging Mechanisms and Disease Relevance of Ferroptosis. Trends Cell. Biol.30 (6), 478–490. 10.1016/j.tcb.2020.02.009
30
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. U.S.A.102 (43), 15545–15550. 10.1073/pnas.0506580102
31
SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.et al (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin.71 (3), 209–249. 10.3322/caac.21660
32
TanG.XuanZ.LiZ.HuangS.ChenG.WuY.et al (2020). The Critical Role of BAP1 Mutation in the Prognosis and Treatment Selection of Kidney Renal Clear Cell Carcinoma. Transl. Androl. Urol.9 (4), 1725–1734. 10.21037/tau-20-1079
33
WangM.MaoC.OuyangL.LiuY.LaiW.LiuN.et al (2019). Long Noncoding RNA LINC00336 Inhibits Ferroptosis in Lung Cancer by Functioning as a Competing Endogenous RNA. Cell. Death Differ.26 (11), 2329–2343. 10.1038/s41418-019-0304-y
34
WuJ.DongG.LiuT.ZhangS.SunL.LiangW. (2021). LncRNA SNHG17 Promotes Tumor Progression and Predicts Poor Survival in Human Renal Cell Carcinoma via Sponging miR-328-3p. Aging13 (17), 21232–21250. 10.18632/aging.203440
35
WuT.HuE.XuS.ChenM.GuoP.DaiZ.et al (2021). clusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innovation2 (3), 100141. 10.1016/j.xinn.2021.100141
36
YanN.ZhangJ. (2019). Iron Metabolism, Ferroptosis, and the Links with Alzheimer's Disease. Front. Neurosci.13, 1443. 10.3389/fnins.2019.01443
37
YangY.TaiW.LuN.LiT.LiuY.WuW.et al (2020). lncRNA ZFAS1 Promotes Lung Fibroblast-To-Myofibroblast Transition and Ferroptosis via Functioning as a ceRNA through miR-150-5p/SLC38A1 axis. Aging12 (10), 9085–9102. 10.18632/aging.103176
38
YeZ.LiuW.ZhuoQ.HuQ.LiuM.SunQ.et al (2020). Ferroptosis: Final Destination for Cancer?Cell. Prolif.53 (3), e12761. 10.1111/cpr.12761
39
ZengJ.LiY.WangY.XieG.FengQ.YangY.et al (2020). lncRNA 00312 Attenuates Cell Proliferation and Invasion and Promotes Apoptosis in Renal Cell Carcinoma via miR-34a-5p/ASS1 Axis. Oxid. Med. Cell. Longev.2020, 1–16. 10.1155/2020/5737289
40
ZhangB.WuQ.WangZ.XuR.HuX.SunY.et al (2019). The Promising Novel Biomarkers and Candidate Small Molecule Drugs in Kidney Renal Clear Cell Carcinoma: Evidence from Bioinformatics Analysis of High‐throughput Data. Mol. Genet. Genomic Med.7 (5), e607. 10.1002/mgg3.607
41
ZouJ.WangL.TangH.LiuX.PengF.PengC. (2021). Ferroptosis in Non-small Cell Lung Cancer: Progression and Therapeutic Potential on it. Int. J. Mol. Sci.22 (24), 13335. 10.3390/ijms222413335
Summary
Keywords
ccRCC, ferroptosis, lncRNA, immune checkpoints, prognosis signatures
Citation
Dong Y, Liu D, Zhou H, Gao Y, Nueraihemaiti Y and Xu Y (2022) A Prognostic Signature for Clear Cell Renal Cell Carcinoma Based on Ferroptosis-Related lncRNAs and Immune Checkpoints. Front. Genet. 13:912190. doi: 10.3389/fgene.2022.912190
Received
04 April 2022
Accepted
09 May 2022
Published
24 May 2022
Volume
13 - 2022
Edited by
Yong Zhang, Capital Medical University, China
Reviewed by
Bosen You, Harbin Medical University, China
Huiying Yang, Sir Run Run Shaw Hospital, China
Updates
Copyright
© 2022 Dong, Liu, Zhou, Gao, Nueraihemaiti and Xu.
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: Yunfei Xu, xuyunfeibb@sina.com
†These authors have contributed equally to this work
‡ORCID ID: Ding Liu, orcid.org/0000-0002-9381-1069
This article was submitted to Cancer Genetics and Oncogenomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.