Int J Biol Sci 2021; 17(3):702-711. doi:10.7150/ijbs.55552 This issue
Ferroptosis-Related Long Non-Coding RNA signature predicts the prognosis of Head and neck squamous cell carcinoma
1. Department of Critical Care Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, 430022, China.
2. Department of Otolaryngology Head and Neck Surgery, The Central Hospital of Wuhan, Tongji Medical College Huazhong University of Science and Technology, Wuhan, Hubei, China.
3. School of Public Health, Tongji Medical College, Huazhong University of Science and technology, Wuhan, China.
4. Department of Otorhinolaryngology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China.
5. Department of Infectious Diseases, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430022, China.
Tang Y, Li C, Zhang YJ, Wu ZH. Ferroptosis-Related Long Non-Coding RNA signature predicts the prognosis of Head and neck squamous cell carcinoma. Int J Biol Sci 2021; 17(3):702-711. doi:10.7150/ijbs.55552. Available from https://www.ijbs.com/v17p0702.htm
Background: Head and neck squamous cell carcinoma (HNSCC) are head and neck cancers. On the other hand, ferroptosis is a novel iron-dependent and ROS reliant type of cell death observed various disease conditions.
Method: We constructed a prognostic multilncRNA signature based on ferroptosis-related differentially expressed lncRNAs in HNSCC.
Results: We identified 25 differently expressed lncRNAs associated with prognosis of HNSCC. Kaplan-Meier analyses revealed the high-risk lncRNAs signature associated with poor prognosis of HNSCC. Moreover, the AUC of the lncRNAs signature was 0.782, underscoring their utility in prediction HNSCC prognosis. Indeed, our risk assessment model was superior to traditional clinicopathological features in predicting HNSCC prognosis. GSEA revealed the immune and tumor-related pathways in the low risk group individuals. Moreover, TCGA revealed T cell functions including cytolytic activity, HLA, regulation of inflammationp, co-stimulation, co-inhibition and coordination of type II INF response were significantly different between the low-risk and high-risk groups. Immune checkpoints such as PDCD-1 (PD-1), CTLA4 and LAG3, were also expressed differently between the two risk groups.
Conclusion: A novel ferroptosis-related lncRNAs signature impacts on the prognosis of HNSCC.
Keywords: Head and neck squamous cell carcinoma, ferroptosis, genes, lncRNAs, immune infiltration, data mining.
Head and neck squamous cell carcinoma (HNSCC) are neoplasms affecting different tissues and organs in the head and neck region including the tongue, mouth, nasopharynx, larynx and throat . Every year, there are over 655,000 new HNSCC cases (globally), resulting in 90,000 deaths . Tobacco smoking, alcohol consumption and human papillomavirus (HPV) infection are important risk factors for developing HNSCC . Moreover, Epstein-Barr virus (EBV) is strongly associated with nasopharyngeal carcinoma (NPC) . Although several treatment options such as surgery, radiotherapy and chemotherapy are available, the asymptomatic nature and lack of early detection of HNSCC result in less than 50% 5-year survival rate of individuals with the cancers. In addition, due to the same reasons, even after treatment, the risk of recurrence and metastasis of HNSCC remain high [5, 6]. Efficient clinical management of HNSCC has been constrained by suboptimal preclinical models or lack of accurate biomarkers for early diagnosis of the cancers. Currently, there are several anti-tumor drugs that induce apoptosis of cancer cells . Thus, it is particularly critical to explore other forms of cell death to overcome tumor cell resistance while at the same time uncovering new and efficient prognostic biomarkers for HNSCC.
Past decades have seen a rapid increase in research on tumor ferroptosis. This is an iron-dependent cell death reliant on intracellular accumulation of reactive oxygen species (ROS), distinct from apoptosis and autophagy . Dysregulation of iron metabolism is a risk factor for cancer and as well promotes tumor growth. Compared with normal cells, cancer cells exhibit iron addiction, which is the over dependence on iron for proliferation . Indeed, activation of ferroptosis pathways may override drug resistance experienced by current chemotherapeutic agents, opening a new therapeutic frontier for cancer treatment. Meanwhile, long non-coding RNAs (lncRNA) are a subset of RNA molecule about 200 nt, which regulate the expression of genes . Apart from gene regulation, lncRNA participates in various biological regulatory processes, including those implicated for occurrence, development and metastasis of tumors . At the moment, there are very few studies on ferroptosis-related lncRNAs. One recent study revealed that lncRNA LINC00618 is under-expressed in human leukemia cells. The mRNA promotes ferroptosis by increasing the levels of lipid ROS and iron, while decreasing the expression of SLC7A11. It further induces vincristine-induced ferroptosis and apoptosis of . In a related study, Lu et al.  found that in acute ischemic stroke, lncRNA PVT1 regulated ferroptosis via miR-214, which modulates the expression of TFR1 and TP53. Similarly, silencing lncRNA ZFAS1 attenuated ferroptosis and significantly modulated inflammation and lipid peroxidation . However, sequence-based studies that systematically evaluate ferroptosis-related lncRNA signature and its association with overall survival (OS) in HNSCC patients remain scanty. In this study, we first constructed a prognostic multi-lncRNA signature of differentially expressed ferroptosis-related lncRNA based on the Cancer Genome Atlas (TCGA) data. We then explored the role of ferroptosis-related mRNA, N6-methyladenosine (m6A) mRNA statues and immune responses in HNSCC prognosis.
RNA-sequence (42 normal and 487 tumor) data of 511 patients was extracted from the TCGA-HNSCC database. The clinical characteristics of the patients are shown in Table 1. The corresponding ferroptosis-related genes were downloaded from FerrDb , a web-based consortium that provided a comprehensive and up-to-date database for ferroptosis markers, their regulatory molecules and associated diseases. Overall, we identified 259 (Driver: 108; suppressor: 69; marker: 111) ferroptosis-related genes (Table S1). The relationship between the ferroptosis-related lncRNAs and HNSCC was assessed using Pearson correlation. The association was considered significant if the correlation coefficient |R2|>0.3 at P<0.001. The collected clinical-pathological data of the HNSCC patients included gender, age, stage, grade, TMN, survival status and survival time. Significant differential expression of ferroptosis-related lncRNAs was set at FDR<0.05 and |log2FC|≥1. First, we explored the function of both upregulated and downregulated ferroptosis-related differentially expressed genes (DEGs). We then used Gene ontology (GO) to evaluate the biological pathways associated with the DEGs. Further function analysis of biological processes (BP), molecular functions (MF) and cellular components (CC)) regulated by the differently expressed ferroptosis-related lncRNAs were analyzed based on Kyoto Encyclopedia of Genes and Genomes (KEGG) data using R software, ggplot2 package.
The clinical characteristics of patients in the TCGA dataset.
|Variable||Number of samples|
|Age at diagnosis|
Development of the ferroptosis-related lncRNAs prognostic signature
We used Lasso‐penalized Cox regression and Univariate Cox regression analyses to construct the ferroptosis-related lncRNAs signature, stratified based on risk score (Coefficient lncRNA1 × expression of lncRNA1) + (Coefficient lncRNA2 × expression of lncRNA2) + ⋯ + (Coefficient lncRNAn × expression lncRNAn). The associated risk score for each HNSCC patients were also evaluated. The RNAs were classified in either low-risk (<median number) or high-risk (≥ median number) group based on the median score.
The predictive nomogram
We performed Gene set enrichment analyses (GSEA) to define the lncRNAs signatures in the KEGG, which were then searched in the TCGA-HNSCC database. Statistical significance was set at P<0.05 and false discovery rate (FDR) q<0.25. A nomogram was constructed integrating the prognostic signatures, for predictive of 1, 3, and 5‐year OS of HNSCC patients.
Immunity analysis and gene expression
At the same time, the CIBERSORT [16, 17], ESTIMATE , MCPcounter , single-sample gene set enrichment analysis (ssGSEA)  and TIMER  algorithms were compared to assess cellular components or cell immune responses between high risk and low risk group based on ferroptosis-related lncRNAs signature. The differences in immune response under different algorithms were uncovered using a Heatmap. In addition, ssGSEA was used to quantify the tumor-infiltrating immune cell subgroups between the two groups as well as assessing their immune function. Potential immune checkpoint was also retrieved from previous literature.
Data was analyzed using Bioconductor packages in R software, version 4.0.2. Normally and non-normally distributed variables were analyzed using the unpaired student's t-test and the Wilcoxon test, respectively. Benjamini-Hochberg method was used to identify the differently expressed lncRNAs, based on FDR. The ssGSEA-normalized HNSCC DEGs were compared with a genome using "GSVA" (R-package). The sensitivity and specificity of the derive prognostic signatures for HNSCC in comparison to other clinicopathological was assessed using the operating characteristic curve (ROC) and decision curve analysis (DCA) . The relationship between ferroptosis-related lncRNAs and clinicopathological manifestations was evaluated using logistic regression analyses and a heatmap graph. The survival analysis of HNSCC patients based on the ferroptosis-related lncRNAs signature was assessed using the Kaplan-Meier survival analysis. For each analysis, statistical significance was set at P <0.05.
Enrichment Analysis of ferroptosis-related genes
We uncovered 65 ferroptosis-related DEGs (18 downregulated and 47 upregulated; Table S2). BP participated in the production of oxidative stress superoxide anion, homeostasis and metabolism of superoxides among others. MF mainly regulated production of nicotinamide adenine dinucleotide phosphate (NADPH) oxidase, molecular oxygen and binding of ferric iron.CC were mainly up-regulated in NADPH oxidase complex, oxidoreductase complex and apical plasma membrane synthesis pathways. KEGG based analysis revealed the over-expressed genes were mainly involved in ferroptosis, hypoxia-inducible factor (HIF)-1 signaling pathway, synthesis of microRNAs in cancer cells, IL-17 signaling pathway, central carbon metabolism in cancer cells, bladder cancer, human cytomegalovirus infection and mTOR signaling pathway (Figure 1 and Table S3).
The ferroptosis-based lncRNAs prognostic signature
We uncovered 935 ferroptosis-related lncRNAs (Table S4). Univariate COX analysis identified 75 significant ferroptosis-related lncRNAs, which were included in the multivariate COX analysis. Overall, 25 differently expressed lncRNAs (LINC01963, AL357033.4, LINC01980, AL132989.1, AATBC, ELF3-AS1, AC135050.6, AL161431.1, AC106820.3, AL451085.2, AC007991.2, AC136475.2, AC012467.2, AC144831.1, AC116914.2, AC008115.3, PSMA3-AS1, PCED1B-AS1, AL139158.2, EP300-AS1, AC104083.1, AL022328.2, AC012640.2, PAX8-AS1, and AL450992.2) were found to be independent prognosis predictors of HNSCC (Table S5). Thus, we calculated risk scores and constructed a prognostic signature for the lncRNAs.
Survival results and multivariate examination
Kaplan-Meier analyses revealed the expression of high-risk lncRNAs signatures corresponded with poorer survival (P<0.001, Figure 2A). Meanwhile, the AUC of the signature lncRNAs was 0.782, exhibiting superior performance than the traditional clinicopathological features in predicting the prognosis of HNSCC (Figure 2B, 2E). Using patient's risk survival status plot, we found the patient's risk score was inversely proportional to the survival of patients with HNSCC. Interestingly, our heatmap demonstrated that most of the novel lncRNAs uncovered in this study were negative correlation with our risk model and needs more researches to conduct (Figure 2C). The AUC predictive value of the novel IncRNAs signature for 1, 3, 5-year survival rate was 0.78, 0.83 and 0.71, respectively (Figures 2D). Univariate and multivariate COX analysis revealed that IncRNAs signature (HR: 1.38, 95CI: 1.25-1.51) as well as tumor stage (HR: 1.88, 95CI: 1.28-2.76) were independent prognosis factors of OS of HNSCC patients (Figure 3A, 3B). The relationship between lncRNA and mRNA is shown in Figure 3C. The heatmap for the association between ferroptosis-related lncRNAs prognostic signature and clinicopathological manifestations were also analyzed (Figure 4). The hybrid nomogram incorporating clinicopathological characteristics and the novel ferroptosis-related lncRNAs prognostic signature (Figure 5) was stable and accurate, thus may be applied in clinical management of HNSCC patients.
Gene set enrichment analyses
Gene set enrichment analyses (GSEA) revealed the majority of the novel ferroptosis-related lncRNAs prognostic signature regulated immune and tumor-related pathways such as primary immunodeficiency, natural killer cell mediated cytotoxicity, T/B cell receptor signaling pathway, FC gamma R mediated phagocytosis, homologous recombination, notch signaling pathway, acute myeloid leukemia, colorectal cancer, and jak stat signaling pathway (Figure 6) and (Table S6).
GO and KEGG analyses for ferroptosis-related differentially expressed genes. (A) GO and (B) KEGG.
Ferroptosis-related lncRNAs signature based on TCGA. (A) Kaplan-Meier curves result, (B). The AUC values of the risk factors, (C). Risk survival status plot, (D) The AUC of the for the prediction of 1, 3, 5-year survival rate of HNSCC and (E) The DCA of the risk factors.
Univariate and multivariate COX analysis for the expression of ferroptosis-related lncRNAs. (A). univariate, (B). multivariate, (C). The relationship between the novel lncRNA and mRNA expression.
Heatmap for ferroptosis-related lncRNAs prognostic signature and clinicopathological manifestations.
Immunity and gene expression
The heatmap of immune responses based on CIBERSORT, ESTIMATE, MCP counter, single-sample gene set enrichment analysis (ssGSEA) and TIMER algorithms is shown in Figure 7. Correlation analysis between immune cell subpopulations and related functions based on ssGSEA of TCGA-HNSCC data revealed that T cell functions including checkpoint (inhibition), cytolytic, HLA, regulation of inflammation, co-stimulation of, co-inhibition of and type II INF response were significantly different between the low-risk and high-risk groups Figure 8A. Given the importance of checkpoint inhibitor-based immunotherapies, we further explored the difference in the expression of immune checkpoints between the two groups. We found a substantial difference in the expression of PDCD-1 (PD-1), CTLA4, LAG3, BTLA among others, between the two groups of patients (Figure 8B). The comparison in the expression of m6A-related mRNA expression between high and low risk group suggested that the expression of RBM15, YTHDC1, and TDHDC2 were significance (Figure 9).
Ferroptosis can overcome resistance of malignant cells to chemotherapy and as well facilitate removal of defective cells. As such, it's potentially a novel approach for tumor treatment. In this study, we first identified novel ferroptosis-related prognostic lncRNAs signature based on the TCGA dataset. We then explored the roles of immune infiltrating cells in tumor microenvironment and immune checkpoint inhibitors in the prognosis of HNSCC. Findings of this study uncovered potential biomarker and therapeutic target in the ferroptosis signaling pathways.
A nomogram for both clinic-pathological factors and prognostic ferroptosis-related lncRNAs.
Gene enrichment analysis for ferroptosis-related lncRNAs based on TCGA.
Heatmap for immune responses based on CIBERSORT, ESTIMATE, MCPcounter, ssGSEA, and TIMER algorithms among high and low risk group.
(A). ssGSEA for the association between immune cell subpopulations and related functions (B). Expression of immune checkpoints among high and low HNSCC risk groups.
The expression of m6A-related genes between high and low HNSCC risk group.
Overall, out analyses uncovered 65 ferroptosis-related DEGs. KEGG analyses further revealed the genes mainly participated in hypoxia-inducible factor (HIF)-1 signaling pathway, expression of cancer related microRNAs, IL-17 signaling pathway, central carbon metabolism in cancer cells and mTOR signaling pathway. A recent study found that FG-4592 (an inhibitor of prolyl hydroxylase of HIF) pretreatment reduces ferroptosis in advanced stage of kidney injury by activating Nrf2 via the Akt/GSK-3β-mediated . Takahashi et al.  reported that cancer spheroids utilize the mammalian target of rapamycin (mTOR) for proliferation and the lipid peroxidase GPX4 for defense against ferroptosis. Overall, in this study, 25 differently expressed lncRNAs were found to be independent prognosis factors for HNSCC. A recent study found over-expression of LINC01963 in pancreatic carcinoma cells inhibited proliferation of pancretic cells via the miR-641/TMEFF2 pathway . Up-regulated expression of LINC01980 in esophageal squamous cell carcinoma confers poor prognosis and as well promotes epithelial-mesenchymal transition (EMT) via the miR-190a-5p and MYO5A pathways . In a related study, lncRNA AATBC was found to be over-expressed in NPC, corresponding to poor survival. Over-expression of lncRNA AATBC also promoted migration and invasion of NPC cells, mediated by miR-1237-3p . lncRNA ELF3-antisense RNA 1 (ELF3-AS1) promotes the proliferation of oral squamous cell carcinoma (OSCC) cells by disrupting glucose metabolism. Overall, over-expression of ELF3-AS1 conferred poor OS of OSCC patients [29, 30]. Similarly, one recent study reported that AL161431.1 target and binds miR-1252-5p, inhibiting MAPK signaling in endometrial carcinoma cells . PSMA3-AS1 knockdown suppresses the migration, proliferation and invasion of glioma cells in vivo and in vitro and represses proteasome inhibitors [32, 33]. Expression of PCED1B-AS1, closely associated with larger tumor size, higher grade, poor survival, rapid proliferation, glucose uptake and lactate release in cancer cells, was found to be upregulated in glioblastoma cells and tissues . Paired box 8 antisense RNA 1 (PAX8-AS1) polymorphisms rs4848320 and rs6726151 are risk factor for the development of childhood acute lymphoblastic leukemia, but decreases the risk of developing cervical cancer [35, 36]. On the other hand, one study reported thatdown-regulation of AL450992.2 was associated with autoimmune disorders . However, to date, there is no study on the role of ferroptosis-associated lncRNAs in the prognosis of HNSCC. Such findings may provide invaluable sight into the future control of cancers.
Here, the differently expressed ferroptosis-associated lncRNAs were stratified into two categories of high‐ and low‐risk to explore their potential roles in HNSCC. Ferroptosis combined with immune checkpoint inhibitors (ICIs) synergistically enhance antitumor activity, even in ICI-resistant types . There are currently very few studies exploring the relationship between ICI and ferroptosis. Increasing evidence suggests that microRNA (miRNA) and lncRNA are critical in mediating the regulation of ferroptosis. Nrf2 inhibits iron absorption, reducing production of ROS. As such, miRNA disrupt ferroptosis by regulating the expression of Nrf2 [39-41]. Meanwhile, miRNA participates in regulating transport, storage, utilization and absorption of iron. In recent years, new ferroptosis regulatory factors such as P53, ATF3/4, SLC7A11, ACSL4 and BECN1 pathway have been discovered. Intriguingly, lncRNA plays a critical role in regulating the expression of these factors .
Ferroptosis is a new form of cell death that can potentially provide a new approach in tumor treatment. However, many key issues such as the interconnection between ferroptosis and other cell deaths and host immunogenicity remain unsolved. Thus, this study explored ferroptosis biomarkers useful in predicting the prognosis of HNSCC, which can inform treatment modalities for the diseases. Nonetheless, our signature profiles need further validation using different cohorts. Given that our findings were not validated using clinical samples, the reliability of our results cannot be fully guaranteed. Also, based on the limited clinical data, our findings should be applied with caution. In general, the prognostic prediction model developed in this study needs further validation.
Specific ferroptosis associated lncRNAs can predict the prognosis of HNSCC.
HNSCC: Head and neck squamous cell carcinoma; ROS: reactive oxygen species; miRNA: microRNA; OS: overall survival (OS); DEGs: differentially expressed genes; TCGA: The Cancer Genome Atlas; GO: Gene ontology; BP: biological processes; MF: molecular function; CC: cellular components; KEGG: Kyoto Encyclopedia of Genes and Genomes; GSEA: Gene set enrichment analyses; FDR: false discovery rate; ssGSEA: single-sample gene set enrichment analysis (ssGSEA).
W.Z.H. designed and analyzed the research study; W.Z.H. wrote and revised the manuscript, T.Y.; L.C. and Z.Y.J. collected and analysis the data and all authors have read and approved the manuscript.
The authors have declared that no competing interest exists.
1. Fang XN, Yin M, Li H, Liang C, Xu C, Yang GW. et al. Comprehensive analysis of competitive endogenous RNAs network associated with head and neck squamous cell carcinoma. Sci Rep. 2018;8(1):10544
2. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M. et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359-86
3. Huang Y, Liu Z, Zhong L, Wen Y, Ye Q, Cao D. et al. Construction of an 11-microRNA-based signature and a prognostic nomogram to predict the overall survival of head and neck squamous cell carcinoma patients. BMC Genomics. 2020;21(1):691
4. Hitt M M, Allday M J, Hara T. EBV gene expression in an NPC-related tumour. Embo Journal. 1989;8(9):2639-2651
5. Posner M, Vermorken JB. Induction therapy in the modern era of combined-modality therapy for locally advanced head and neck cancer. Semin Oncol. 2008;35(3):221-8
6. Kumpitsch C, Moissl-Eichinger C, Pock J, Thurnher D, Wolf A. Preliminary insights into the impact of primary radiochemotherapy on the salivary microbiome in head and neck squamous cell carcinoma. Sci Rep. 2020;10(1):16582
7. Sung ES, Kim JY, Ahn YT, Lee IW, Choi SW, Jang HB. et al. An WG. Melatonin Exerts Anticancer Effects in Human Tongue Squamous Cell Carcinoma Cells by Promoting Autophagy. Anticancer Res. 2020;40(11):6295-6303
8. Huang L, McClatchy DB, Maher P, Liang Z, Diedrich JK, Soriano-Castell D. et al. Intracellular amyloid toxicity induces oxytosis/ferroptosis regulated cell death. Cell Death Dis. 2020;11(10):828
9. Manz DH, Blanchette NL, Paul BT, Torti FM, Torti SV. Iron and cancer: recent insights. Ann N Y Acad Sci. 2016;1368(1):149-61
10. Zhang T, Li K, Zhang ZL, Gao K, Lv CL. LncRNA Airsci increases the inflammatory response after spinal cord injury in rats through the nuclear factor kappa B signaling pathway. Neural Regen Res. 2020;16(4):772-777
11. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ. et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464(7291):1071-6
12. Wang Z, Chen X, Liu N, Shi Y, Liu Y, Ouyang L. et al. Nuclear Long Non-Coding RNA LINC00618 Accelerates Ferroptosis in a Manner Dependent upon Apoptosis. Mol Ther. 2020:S1525-0016(20)30481-0.
13. Lu J, Xu F, Lu H. LncRNA PVT1 regulates ferroptosis through miR-214-mediated TFR1 and p53. Life Sci. 2020;260:118305
14. Yang Y, Tai W, Lu N, Li T, Liu Y, Wu W. et al. lncRNA ZFAS1 promotes lung fibroblast-to-myofibroblast transition and ferroptosis via functioning as a ceRNA through miR-150-5p/SLC38A1 axis. Aging (Albany NY). 2020;12(10):9085-9102
15. Zhou N, Bao J. FerrDb: a manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database (Oxford). 2020;2020:baaa021
16. Newman AM, Liu CL, Green MR. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453-457
17. Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D. et al. Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint Blockade. Cell Rep. 2017;18(1):248-262
18. 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
19. Shi J, Jiang D, Yang S, Zhang X, Wang J, Liu Y. et al. LPAR1, Correlated With Immune Infiltrates, Is a Potential Prognostic Biomarker in Prostate Cancer. Front Oncol. 2020;10:846
20. Yi M, Nissley DV, McCormick F, Stephens RM. ssGSEA score-based Ras dependency indexes derived from gene expression data reveal potential Ras addiction mechanisms with possible clinical implications. Sci Rep. 2020;10(1):10258
21. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS. et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017;77(21):e108-e110
22. Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. 2006;26(6):565-74
23. Havel JJ, Chowell D, Chan TA. The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy. Nat Rev Cancer. 2019;19(3):133-150
24. Li X, Zou Y, Xing J, Fu YY, Wang KY, Wan PZ. et al. Pretreatment with Roxadustat (FG-4592) Attenuates Folic Acid-Induced Kidney Injury through Antiferroptosis via Akt/GSK-3β/Nrf2 Pathway. Oxid Med Cell Longev. 2020;2020:6286984
25. Takahashi N, Cho P, Selfors LM, Kuiken HJ, Kaul R, Fujiwara T. et al. 3D Culture Models with CRISPR Screens Reveal Hyperactive NRF2 as a Prerequisite for Spheroid Formation via Regulation of Proliferation and Ferroptosis. Mol Cell. 2020:S1097-2765(20)30693-6.
26. Li K, Han H, Gu W, Cao C, Zheng P. Long non-coding RNA LINC01963 inhibits progression of pancreatic carcinoma by targeting miR-641/TMEFF2. Biomed Pharmacother. 2020;129:110346
27. Liang X, Wu Z, Shen S, Niu Y, Guo Y, Liang J. et al. LINC01980 facilitates esophageal squamous cell carcinoma progression via regulation of miR-190a-5p/MYO5A pathway. Arch Biochem Biophys. 2020;686:108371
28. Tang T, Yang L, Cao Y, Wang M, Zhang S, Gong Z. et al. LncRNA AATBC regulates Pinin to promote metastasis in nasopharyngeal carcinoma. Mol Oncol. 2020;14(9):2251-2270
29. Chu H, Li Z, Gan Z, Yang Z, Wu Z, Rong M. LncRNA ELF3-AS1 is involved in the regulation of oral squamous cell carcinoma cell proliferation by reprogramming glucose metabolism. Onco Targets Ther. 2019;12:6857-6863
30. Mei JC, Yan G, Mei SQ. Diagnostic and Prognostic Potentials of Long Noncoding RNA ELF3-AS1 in Glioma Patients. Dis Markers. 2020;2020:8871746
31. Gu ZR, Liu W. The LncRNA AL161431.1 targets miR-1252-5p and facilitates cellular proliferation and migration via MAPK signaling in endometrial carcinoma. Eur Rev Med Pharmacol Sci. 2020;24(5):2294-2302
32. Zhou LL, Zhang M, Zhang YZ, Sun MF. Long non-coding RNA PSMA3-AS1 enhances cell proliferation, migration and invasion by regulating miR-302a-3p/RAB22A in glioma. Biosci Rep. 2020;40(9):BSR20191571
33. Xu H, Han H, Song S, Yi N, Qian C, Qiu Y. et al. Exosome-Transmitted PSMA3 and PSMA3-AS1 Promote Proteasome Inhibitor Resistance in Multiple Myeloma. Clin Cancer Res. 2019;25(6):1923-1935
34. Yao Z, Zhang Q, Guo F, Guo S, Yang B, Liu B. et al. Long Noncoding RNA PCED1B-AS1 Promotes the Warburg Effect and Tumorigenesis by Upregulating HIF-1α in Glioblastoma. Cell Transplant. 2020;29:963689720906777
35. Bahari G, Hashemi M, Naderi M, Sadeghi-Bojd S, Taheri M. Long non-coding RNA PAX8-AS1 polymorphisms increase the risk of childhood acute lymphoblastic leukemia. Biomed Rep. 2018;8(2):184-190
36. Han J, Zhou W, Jia M, Wen J, Jiang J, Shi J. et al. Expression quantitative trait loci in long non-coding RNA PAX8-AS1 are associated with decreased risk of cervical cancer. Mol Genet Genomics. 2016;291(4):1743-8
37. Teimuri S, Hosseini A, Rezaenasab A, Ghaedi K, Ghoveud E, Etemadifar M. et al. Integrative Analysis of lncRNAs in Th17 Cell Lineage to Discover New Potential Biomarkers and Therapeutic Targets in Autoimmune Diseases. Mol Ther Nucleic Acids. 2018;12:393-404
38. Tang R, Xu J, Zhang B, Liu J, Liang C, Hua J. et al. Ferroptosis, necroptosis, and pyroptosis in anticancer immunity. J Hematol Oncol. 2020;13(1):110
39. Kabaria S, Choi DC, Chaudhuri AD. MicroRNA- 7 activates Nrf2 pathway by targeting Keap1 expression. Free Radic Biol Med. 2015;89:548-556
40. Yang M, Yao Y, Eades G, Zhang Y, Zhou Q. MiR-28 regulates Nrf2 expression through a Keap1-independent mechanism. Breast Cancer Res Treat. 2011Oct;129(3):983-91
41. Xu D, Zhu H, Wang C, Zhu X, Liu G, Chen C. microRNA-455 targets cullin 3 to activate Nrf2 signaling and protect human osteoblasts from hydrogen peroxide. Oncotarget. 2017Jul22;8(35):59225-59234
42. Song X, Long D. Nrf2 and Ferroptosis: A New Research Direction for Neurodegenerative Diseases. Front Neurosci. 2020;14:267
Corresponding author: Zeng-Hong Wu, Department of Otorhinolaryngology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China (Email: D201981596edu.cn).