- Open Access
Potential prognostic value of a eight ferroptosis-related lncRNAs model and the correlative immune activity in oral squamous cell carcinoma
BMC Genomic Data volume 23, Article number: 80 (2022)
To investigate the prognostic value of ferroptosis-related long noncoding RNAs (lncRNAs) in oral squamous cell carcinoma (OSCC) and to construct a prognostic risk and immune activity model.
We obtained clinical and RNA-seq information on OSCC patient data in The Cancer Genome Atlas (TCGA) Genome Data Sharing (GDC) portal. Through a combination of a differential analysis, Pearson correlation analysis and Cox regression analysis, ferroptosis-related lncRNAs were identified, and a prognostic model was established based on these ferroptosis-related lncRNAs. The accuracy of the model was evaluated via analyses based on survival curves, receiver operating characteristic (ROC) curves, and clinical decision curve analysis (DCA). Univariate Cox and multivariate Cox regression analyses were performed to evaluate independent prognostic factors. Then, the infiltration and functional enrichment of immune cells in high- and low-risk groups were compared. Finally, certain small-molecule drugs that potentially target OSCC were predicted via use of the L1000FWD database.
The prognostic model included 8 ferroptosis-related lncRNAs (FIRRE, LINC01305, AC099850.3, AL512274.1, AC090246.1, MIAT, AC079921.2 and LINC00524). The area under the ROC curve (AUC) was 0.726. The DCA revealed that the risk score based on the prognostic model was a better prognostic indicator than other clinical indicators. The multivariate Cox regression analysis showed that the risk score was an independent prognostic factor for OSCC. There were differences in immune cell infiltration, immune functions, m6A-related gene expression levels, and signal pathway enrichment between the high- and low-risk groups. Subsequently, several small-molecule drugs were predicted for use against differentially expressed ferroptosis-related genes in OSCC.
We constructed a new prognostic model of OSCC based on ferroptosis-related lncRNAs. The model is valuable for prognostic prediction and immune evaluation, laying a foundation for the study of ferroptosis-related lncRNAs in OSCC.
Oral cancer ranks among the most prevalent malignant tumours in the head and neck. In 2020, more than 350,000 newly confirmed cases and 175,000 deaths from oral cancers were reported worldwide . Oral squamous cell carcinoma (OSCC) accounts for 90% of oral cancers . At present, many clinical guidelines clearly indicate that the diagnosis and treatment of OSCC cannot be generalized, and the use of comprehensive sequence therapy should be accompanied by individualized treatment . However, despite this guidance, a OSCC diagnosis is a poor prognosis, with a 5-year survival rate of approximately 60% . OSCC is also associated with a high cervical lymph node metastasis rate, leading to a worsened prognosis . Therefore, finding new predictors of survival and developing new detection methods for better clinical decision-making are essential.
Ferroptosis refers to an iron-dependent cell death process, and the morphological characteristics and biochemical markers of ferroptosis are significantly different from those of apoptosis, necrosis, and autophagy [6,7]. Although ferroptosis was first described in 2012 , a clearer understanding of ferroptosis-related mechanisms and functions have since led researchers to show that ferroptosis is inseparable from tumours. Recent research has revealed the association of ferroptosis with tumorigenesis and progression in, for example, bladder cancer , ovarian cancer  and breast cancer . In addition, ferroptosis plays a role in tumours by interacting with different components in the tumour microenvironment (TME). Tumour cells with reduced E-cadherin levels and loss of intercellular adhesion have been reported to be highly sensitive to ferroptosis [12,13], and cell density is an important factor in determining the susceptibility to ferroptosis regardless of the cell-specific phenotype . Most solid tumours are hypoxic, and hypoxia increases the level of carbonic anhydrase 9 (CA9). Studies have shown that elevated CA9 can reduce ferroptosis by controlling intracellular iron metabolism . Ferroptosis also affects tumour cell sensitivity to radiotherapy and can be used to overcome chemotherapy resistance [16,17]. In OSCC, certain ferroptosis-related genes, such as SLC7A11  and GPX4 , can impact the prognosis of patients by regulating ferroptosis in cancer cells. These findings suggest that developing ferroptosis-related treatment strategies is an emerging direction for OSCC treatment.
Long noncoding RNAs (lncRNAs) are RNAs with a transcript length between 200 and 100,000 nt and that do not encode proteins but participates in many physiological processes . To date, more than 1 × 106 lncRNAs have been reported in the human genome, and it has been indicated that disordered lncRNAs are closely connected to the occurrence and development of human cancers . LncRNAs can regulate biological behaviours such as tumour cell proliferation, apoptosis, invasion, and metastasis. Recently, the effects of lncRNAs on ferroptosis regulation have been studied by researchers. Studies have shown that lncRNAs, as dual regulators of ferroptosis, either participate in ferroptosis by inactivating certain miRNAs, as endogenous competing RNAs, or binding to certain enzymes to regulate ferroptosis and influence the biological activity of cancer cells . The most recent reports revealed an association of ferroptosis-related lncRNAs with the prognosis of various cancers, such as colon adenocarcinoma  and breast cancer . However, the role played by ferroptosis as well as its associated lncRNAs in OSCC remains unclear. Therefore, studying lncRNAs associated with OSCC and ferroptosis is crucial for understanding the mechanisms underlying OSCC.
Bioinformatics techniques constitute a new technological approach by effectively combining bioinformatics with medicine. Functional genomics based on bioinformatics is a rapidly developing field . The TCGA database includes complete genome-sequencing studies of a variety of tumours, providing great help for scientific research and discovery of new molecular targets in tumours. Many tumour biomarkers have been discovered and applied clinically, significantly leading to early diagnosis of tumours and increasing the overall survival rate [26,27]. Recently, a model containing 8 ferroptosis-related lncRNAs has been reported; however, the model exhibited low predictive power for OSCC, with an area under the curve (AUC) = 0.690 . Other scholars constructed a prognostic model containing 9 ferroptosis-related lncRNAs , which only used bioinformatics to explore the relationship between ferroptosis-related lncRNAs and the prognosis of head and neck squamous cell carcinoma patients. In addition, this model was not specific for OSCC and lacks relevant in vitro experimental validation. A new prognostic model of OSCC incorporating ferroptosis-related lncRNAs was developed using bioinformatics methods. The prognostic ability of this model was confirmed, and immune function was analysed via different methods. In addition, we investigated differentially expressed ferroptosis genes in the L1000FWD database, identifying small-molecule drugs that potentially target ferroptosis genes in OSCC.
Materials and methods
We obtained RNA sequencing (FPKM) and clinical information on OSCC from the TCGA (https://portal.gdc.cancer.gov/). Table 1 presents the clinical data for 338 samples. According to the FerrDb website (http://www.zhounan.org/ferrdb/) and previous research, 382 ferroptosis-related genes were identified, including ferroptosis-inducing genes, ferroptosis-suppressing genes and ferroptosis markers. The codes used in this study can be found on Github (https://github.com/qiulin961028/ferroptosis-related-lncRNAs), and Fig. 1 shows the flow chart.
Construction and validation of the prognostic model
Ferroptosis-related gene expression was determined for the samples, and Pearson correlation analysis was performed to identify ferroptosis-related lncRNAs (|correlation coefficient|> 0.4, p < 0.001). Then, we acquired lncRNAs that show prognostic promise in ferroptosis as determined through univariate Cox regression (p < 0.05). Before establishing the model, we constructed a network with ferroptosis-related mRNAs and lncRNAs, followed by visualization using Cytoscape. The prognostic risk model was further refined by multivariate Cox regression analysis, and the risk score for patients was calculated using Eq. (1):
Coefi is the risk regression coefficient for every ferroptosis-related lncRNA, and X represents the lncRNA expression level. Based on this model, patients’ risk scores were measured, and the patients were assigned to a low- or high-risk group in with the median risk score serving as the cut-off value.
Immediately afterward this analysis, the overall survival (OS) for patients with OSCC was compared between the two risk groups via a survival analysis. The accuracy of the prognostic model was evaluated on the basis of ROC curves. We thus identified factors that independently predicted prognosis via univariate and multifactorial Cox regression. Prognostic correlation line graphs including age, risk score, sex, tumour grade, and TN stage were plotted with the "RMS" package in R language software, and internal calibration curves were plotted for the line graphs. LNCipedia (https://lncipedia.org/) was used to retrieve ferroptosis-related lncRNA sequences, and the lncLocator database (http://www.csbio.sjtu.edu.cn/bioinf/lncLocator/) was used to identify lncRNA cellular compartment localization based on its sequence.
Immune cell infiltration prediction
To evaluate the degree of immune cell infiltration, we performed a ssGSEA to quantify subgroups of infiltrating immune cells in conjunction with the immune function of both groups. The underlying immune checkpoint and m6A genes were identified based on previous research, and gene expression differences between the two groups were examined.
Pathway enrichment analysis
Further, Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed with both groups. Using GSEA (4.1.1) software, the data were analysed, and enrichment maps were created.
Potential small molecule drug prediction
Differentially expressed ferroptosis-related genes were classified into up- or downregulated groups and imported into the L1000FWD website (https://maayanlab.cloud/L1000FWD/) to obtain permuted outcomes. Drug structures are shown on PubChem.ncbi.nlm.nih.gov.
Human OSCC cell lines WSU-HN6 and CAL-27 were used in this study. WSU-HN6 was obtained from Ninth People's Hospital, Shanghai Jiao Tong University School of Medicine (Shanghai, China), and CAL-27 cell line was purchased from American Type Culture Collection (ATCC, Manassas, USA). All cells were passaged and preserved in the Central Laboratory of Peking University Hospital of Stomatology and regularly tested to ensure mycoplasma negative. All cells were cultured in high glucose DMEM medium (Gibco, CA, USA) containing 10% fetal bovine serum (Gibco, CA, USA) and 1% penicillin/streptomycin solution at 37 °C and 5% CO2.
Total RNA was extracted from cells and tissues using Trizol. Cytoplasmic and nuclear RNA were isolated and purified using the Nuc-Cyto-Mem Preparation Kit (APPLYGEN) and Trizol according to the manufacturers' instructions. Then totol RNA reverse transcribed into cDNA using a Prime Script™ RT Kit. The cDNA template was subsequently amplified by real-time PCR (RT‒PCR) using SYBR Green qPCR Master Mix (ABclonal, Beijing, China). GAPDH and U6 wer used as the internal reference and mRNA relative expression was measured by the 2−ΔΔCT method. The primer sequences were shown in supplementary Table 1.
For gene expression levels, the Wilcoxon test and unpaired Student's t test were performed with data showing with a normal and a nonnormal distribution, respectively. We assessed OSCC patient survival by Kaplan‒Meier curves, and ROC analysis and DCA were performed with the "timeROC" and "ggDCA" software packages, respectively. Data analysis was performed using R software (4.1.1), with P < 0. 05 indicating a significant difference.
Data processing and discovery of ferroptosis-associated lncRNAs with prognostic significance
A total of 386 differentially expressed lncRNAs in OSCC were obtained by rank sum test (Fig. 2A). We obtained differentially expressed ferroptosis-related lncRNAs via correlation analysis, and eight ferroptosis- and prognosis-related lncRNAs were recognized via univariate Cox survival analysis: FIRRE, LINC01305, AC099850.3, AL512274.1, AC090246.1, MIAT, AC079921.2 and LINC00524 (Fig. 2B). A correlation network between ferroptosis genes and these prognosis-related lncRNAs was constructed and visualized by Cytoscape (Fig. 2C). Among these lncRNAs, AC099850.3, LINC01305, and AL512274.1 were coexpressed with a relatively higher number of ferroptosis genes.
Prognostic model establishment and verification
A prognostic risk model was established on the basis of Cox regression analysis; then, we determined risk scores for all cases for the expression levels of risk regression coefficients and ferroptosis-related lncRNAs. Risk score = [FIRRE expression × (0.66714)] + [LINC01305 expression × (0.78751)] + [AC099850.3 expression × (0.029993)] + [AL512274.1 expression × (0.05794)] + [AC090246.1 expression × (0.541331)] + [MIAT expression × (0.24386)] + [AC079921.2 expression × (0.75098)] + LINC00524 expression × (0.105386)].
The survival analysis results revealed an the obviously lower OS rate in the high-risk group compared with that in the low-risk group (p < 0.001) (Fig. 3A). The ROC curve showed 1-, 2- and 3-year area under the curve (AUC) values of 0.726, 0.677, and 0.687, respectively (Fig. 3B), suggesting that the risk model showed good performance for predicting patient prognosis.
The risk score independently predicts OSCC prognosis
Univariate Cox analysis was performed on the basis of patients’ clinical features. The findings revealed that age, risk score, stage, and tumour grade were differed greatly and that these characteristics were risk factors for OSCC (Fig. 4A). However, another multifactorial Cox analysis revealed that the risk score may independently predict OSCC prognosis (Fig. 4B) (HR = 1.444, 95% CI = 1.207–1.728).
Ranking of patients according to risk scores to analyse their survival status revealed a lower survival status and higher death likelihood for high-risk patients (Fig. 4C and D). The differential expression profiles for the eight lncRNAs between the two groups are displayed in a heatmap (Fig. 4E), which shows that FIRRE, AC099850.3, and AC090246.1 expression was obviously increased in the high-risk group, whereas that of LINC01305, AL512274.1, MIAT, and AC079921.2 was significantly decreased. Therefore, the risk model’s accuracy in predicting the prognosis of OSCC patients was confirmed.
Relationship of clinicopathological features with the risk model
To assess the difference in prognosis predicted by the risk model and analysis clinicopathological features, ROC curves of clinical features and risk scores were drawn. As shown in Fig. 5A, the risk model’s AUC value exceeded that of other clinical indicators (AUC = 0.726, 1 year). We then plotted a DCA curve, which indicated that the risk score was a better prognostic factor than other clinical indicators (Fig. 5B). Immediately afterwards, we evaluated the relationship between clinical indicators and risk values for each patient, and the results were plotted in a heatmap (Fig. 5C), which showed a significant difference in the T stage of OSCC of both groups (p < 0.05). Subsequently, we constructed a nomogram including age, sex, stage, grade, risk score, TN stage and other prognostic factors with the nomogram’s internal calibration curves. Then, we selected an OSCC patient and used the patient’s data for scoring. Based on the score, the probability of this patient’s surviving less than 1, 3 and 5 years was predicted (the probability of survival less than 1, 3 and 5 years was 8.33, 21.8 and 28%, respectively), and personalized treatment was determined to be an option (Fig. 5D). In addition, the results also showed that the nomogram correction curves at 1, 3 and 5 years were very close to the ideal line, which indicated that the nomogram exhibited high accuracy in predicting the survival rate of the patient at 1, 3 and 5 years (Fig. S1).
Moreover, considering that the cellular localization of lncRNAs determines the underlying mechanisms, we analysed the subcellular localization of the eight lncRNAs via lncLocator. As shown in Fig. 6A-H, AC099850.3, AC090246.1, MIAT, AC079921.2 and LINC00524 were mainly located in the cytoplasm, the other two lncRNAs (LINC01305 and AL512274.1) were mainly distributed in the cytosol, and FIRRE was mainly located in the nucleus. Subsequently, the results of in vitro experiments were consistent with the predicted results of the database. In two OSCC cell lines, FIRRE, LINC01305 and AL512274.1 were localized in the nucleus. While, AC099850.3, AC090246.1, MIAT, AC079921.1 and lINC00524 were localized in the cytoplasm (Figs. 6 I, J).
Differential immune cell infiltration and function between the two groups
The association of the risk model with immune cell infiltration was explored. The immune cell infiltration analysis results for both groups are presented in a heatmap (Fig. 7A). Vertical coordinates represent immune cell infiltration results for both groups as predicted by different software. Furthermore, immune functions were compared between the two groups, and differences between the groups in immune-related functions, including T-cell costimulation, T-cell coinhibition, CCR, and HLA were identified (Fig. 7B). Thus, the results of the ssGSEA of immune infiltration suggested that immune status was significantly different between the two groups, suggesting a need to develop individualized immunotherapy for OSCC patients.
In addition to differences in immune function and immune cell infiltration, we also examined differences in m6A-associated genes and immune checkpoints between the two groups. A total of 48 immune checkpoints were analysed, and only 29 checkpoint genes were found to be expressed significantly differently between the groups, as shown in Fig. 7C. The expression of m6A-related genes, including ALKBH5, HNRNPC, and YTHDF1, exhibited significant upregulation in high-risk patients (p < 0.05), whereas YTHDC2 gene expression was significantly downregulated in high-risk patients (p < 0.01) (Fig. 7D).
The KEGG enrichment analysis was performed to assess differences in the pathways enriched between the two groups. Based on the findings, 10 active pathways were identified in the high-risk patients and as many as 24 active signalling pathways were identified in the low-risk patients (p < 0.05). Figure 8 shows the key enrichment results. More active pathways in high-risk patients were related to metabolism, such as spliceosome, pyrimidine metabolism, and purine metabolism. On the other hand, significant enrichment in the low-risk group was identified in immune-associated biological process terms: B-cell receptor pathway, T-cell receptor pathway, and FcεRI pathway.
L1000FWD analysis led to the identification of potential target drugs
We searched for potential target drugs in OSCC by uploading the up- and downregulated differentially expressed ferroptosis genes to the L1000FWD database. The top ten drug candidates were obtained, and the basic information of the drugs is shown in Table 2. Treatment with these drugs led to differences in gene enrichment, and thus, MEK inhibitors, oestrogen receptor agonists, RAF inhibitors, etc., were identified Therefore, these small-molecule drugs may be candidate drugs for OSCC treatment and be references for the development of new individualized small-molecule drugs. Among these small-molecule drugs, we selected the three most promising for visualization, and the 2D and 3D structures of KM-03949SC, RJC-00245SC and BRD-K82185908 are shown in Fig. 9.
Internal validation and real-time PCR
We also analysed differences in the expression of ferroptosis-related lncRNAs with respect to different clinical features (Figs. 10A-G). FIRRE, LINC01305, AC099850.3, AL512274.1, AC090246.1, MIAT, AC079921.2 and LINC00524 were differentially expressed in tumour and normal tissues (Fig. 10A). In addition, AL512274.1 and AC090246.1 was differentially expressed in N stage tumours (Fig. 10C); AL512274.1 and MIAT were differentially expressed in lymphovascular invasion (Fig. 10D); LINC01305, AL512274.1 and AC079921.2 were differentially expressed in different grades (Fig. 10F); and AC099850.3, AL512274.1 and MIAT expression was strongly correlated with OS events in OSCC patients (Fig. 10G). In addition, we also detected the expression levels of eight lncRNAs in four pairs of matched OSCC (T), adjacent normal tissues (N). As shown in Figs. 10H-O, the relative expression levels of FIRRE, LINC01305, AC099850.3, AC090246.1, MIAT, AC079921.2 and LINC00524 in OSCC tissues were higher than those in adjacent normal tissues, while the relative expression level of AL512274.1 was lower than those in adjacent normal tissues. Therefore, the expression levels of the eight lncRNAs were consistent with the results of our model analysis.
Patients with OSCC, a common head and neck cancer, have an overall poor prognosis. According to the latest NCCN dental guidelines, surgery, chemotherapy and radiotherapy are recommended for OSCC. Through individualized therapy, treatments are selected on the basis of different states of the disease . In recent years, the use of multidisciplinary therapies has enabled OSCC patients to obtain optimal treatment options with minimal risk of complications. A pathologist is responsible for a definitive diagnosis, a surgeon completely removes the lesion, and a radiologist performs precision radiotherapy. In addition to advising on postoperative repair and care, our goal is to defeat OSCC . In this process, an effective predictor of OSCC prognosis is critical for treatment decisions. However, the most commonly used TNM staging and relevant clinical features of patients lead to certain limitation in the analysis [32,33]. TNM staging is based solely on the number and size of positive lymph nodes and does not pay account for negative lymph node number. Therefore, an increasing number of studies has emphasized the importance of negative lymph node and total lymph node numbers in the prognosis of OSCC. These two indicators have been shown to be independent prognostic factors in several malignancies [34,35,36]. A retrospective case analysis of 120 patients with OSCC in Europe indicated that there were insignificant differences in prognosis between young patients and patients 40 years or older . A similar retrospective analysis from Taiwan did not reveal a difference in survival between males and females with OSCC . These results suggested that TNM stage and clinical features are not ideal factors for predicting OSCC prognosis. Therefore, it is necessary to find other prognostic indicators to accurately determine the prognosis and guide treatment for patients. In this research, a bioinformatics method was adopted to explore the effect of ferroptosis-related lncRNAs on OSCC prognosis. A prognostic risk model was established for validating model accuracy, which offered a new perspective on the effect of ferroptosis on OSCC pathogenesis and prognosis.
The prognostic model developed in this study included eight ferroptosis-related lncRNAs (FIRRE, LINC01305, AC099850.3, AL512274.1, AC090246.1, MIAT, AC079921.2, and LINC00524). FIRRE is related to breast cancer , gallbladder cancer  and diffuse large B-cell lymphoma . In another study, FIRRE inhibited proinflammatory factor production by decreasing the expression of HMGB1, thus relieving neuropathic pain in female rats . In oesophageal squamous carcinoma, LINC01305 was found to regulate HTR3A mRNA, thereby promoting cancer cell metastasis and proliferation . In addition, LINC01305 regulated the epithelial-mesenchymal transition (EMT) in cervical cancer  and lung cancer  through different pathways. MIAT is also an important lncRNA that was first found to be associated with myocardial infarction  and later it was found to be involved in the progression of tumours, including retinoblastoma , smooth muscle tumour , and nasopharyngeal carcinoma . AC099850.3, AL512274.1, AC090246.1, AC079921.2, and LINC00524 have rarely been studied in solid tumours, but different bioinformatics analysis methods suggest their potential prognostic value in different cancers (AC099850.3  in non-small-cell lung cancer, AL512274.1  in OSCC, and LINC00524  in clear cell renal cell carcinoma). Thus, the eight lncRNAs identified in our study are coexpressed with many ferroptosis genes, and they are previously identified lncRNAs. Unfortunately, they have not been extensively studied in OSCC, and their biological functions in OSCC have not been reported.
Ferroptosis is inextricably linked to cellular metabolism. When the dynamic balance of cellular metabolism in the body is disrupted, ferroptosis leads to many diseases caused by metabolic imbalance, such as heart disease and brain injury . By performing a functional enrichment analysis, we found that the pathways associated with high-risk cases were primarily cellular metabolism-related pathways such as purine metabolism, pyrimidine metabolism, and amino acid metabolism, indicating that ferroptosis is closely related to cellular metabolism in OSCC, which aligns with the results reported in the literature .
Immune cell infiltration results demonstrated that the prognostic model correlated with immune activity. Numerous investigations have been conducted to explore the correlation of ferroptosis with tumour immunity, but these findings were insufficient to determine the exact relationship between the immune system and ferroptosis. The pathways associated with low-risk cases were primarily T and B-cell receptor (TCR, BCR) pathways. TCR is the molecule on the T-cell surface that specifically recognizes antigens and determines the diversity of T cells [54,55]. Signal transduction initiated by T-cell antigen receptors is at the core of T-cell activation, which is necessary for acquired immunity . CD4+ T cells and CD8+ T cell activation initiates the TCR signalling pathway, allowing the body to respond with the appropriate immune response . BCR is a molecule on the surface of B cells that specifically recognizes antigens, and the activation signal generated by its binding to antigens is the first signal in B-cell activation by regulating B-cell gene transcription. These results suggest that the acquired immune system is more active and immune activity is stronger in low-risk OSCC patients than in high-risk OSCC patients. Additionally, differences in immune checkpoint, immune function, and m6A-associated genes between the two groups offer a theoretical foundation for the development of individualized immune-targeted therapies for OSCC patients. Small-molecule drugs predicted on the basis of ferroptosis-related DEGs have not yet been experimentally validated, and their specific role in OSCC still needs further exploration.
Although we constructed a highly accurate prognostic model for OSCC patients, it must be admitted that there were still some limitations in this study. This model was derived from the comprehensive genome analysis of OSCC cases in TCGA database, and lacked the ability to specifically recognize tumor cells. However, OSCC was known to be highly heterogeneous and to have a poor prognosis . This heterogeneity can occur not only among individuals, but also within the same individual or even within the same organization. Therefore, further study on heterogeneity is necessary. Our model needs to be combined with single-cell sequencing to better understand the heterogeneity between cells, and its applicability and accuracy also need to be further explored in clinical patients with OSCC. In conclusion, the prognostic model of ferroptosis-related lncRNAs was expected to be a novel biomarker for OSCC diagnosis and treatment decision making. However, further in vitro and in vivo studies are needed to clarify the role and mechanism of ferroptosis-related lncRNAs in the OSCC occurrence and development.
We constructed a new prognostic model of OSCC based on ferroptosis-related lncRNAs. The model is valuable for prognostic prediction and immune evaluation, laying a foundation for the study of ferroptosis-related lncRNAs in OSCC.
Availability of data and materials
This study followed the policies and guidelines for data access and publication specified by The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). The public databases involved in this research are as follows: FerrDb website (http://www.zhounan.org/ferrdb/), LNCipedia (https://lncipedia.org/), lncLocator database (http://www.csbio.sjtu.edu.cn/bioinf/lncLocator/).
Long non-coding RNAs
Oral squamous cell carcinoma
The Cancer Genome Atlas
Genome Data Sharing
Least absolute shrinkage and selection operator
Receiver Operating Characteristic
Decision Curve Analysis
Areas under the curve
Kyoto Encyclopedia of Genes and Genomes
T cell receptor
B cell receptor; EMT: epithelial-mesenchymal transition
Sung H, Ferlay J, Siegel RL, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
César R. Essentials of Oral Cancer. Int J Clin Exp Pathol. 2015;8(9):11884–94. https://doi.org/10.5281/zenodo.192487.
Machiels JP, Leemans CR, Golusinski W, et al. Squamous cell carcinoma of the oral cavity, larynx, oropharynx and hypopharynx: EHNS–ESMO–ESTRO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2020;31(11):1462–75. https://doi.org/10.1016/j.annonc.2020.07.011.
Beltramini GA, Belloni LM, Fusco N, et al. Comparing Prognostic Utility Between the 8th Edition of TNM Staging System and the Lymph Node Ratio for Oral Squamous Cell Carcinoma. Head Neck. 2021;43(10):2876–82. https://doi.org/10.1002/hed.26769.
Gamez ME, Kraus R, Hinni ML, et al. Treatment outcomes of squamous cell carcinoma of the oral cavity in young adults. Oral Oncol. 2018;87:43–8. https://doi.org/10.1016/j.oraloncology.2018.10.014.
Xie Y, Hou W, Song X, et al. Ferroptosis: process and function. Cell Death Differ. 2016;23:369–79. https://doi.org/10.1038/cdd.2015.158.
Tang D, Kang R, Berghe TV, et al. The molecular machinery of regulated cell death. Cell Res. 2019;29(5):347–64. https://doi.org/10.1038/s41422-019-0164-5.
Dixon SJ, Lemberg KM, Lamprecht MR, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–72. https://doi.org/10.1016/j.cell.2012.03.042.
Luo WJ, Wang J, Xu WH, et al. LncRNA RP11-89 facilitates tumorigenesis and ferroptosis resistance through PROM2-activated iron export by sponging miR-129-5p in bladder cancer. Cell Death Dis. 2021;12(1043):1–11. https://doi.org/10.1038/s41419-021-04296-1.
Li DX, Zhang ML, Chao HT. Significance of glutathione peroxidase 4 and intracellular iron level in ovarian cancer cells-"utilization" of ferroptosis mechanism. Inflamm Res. 2021;70(10–12):1177–89. https://doi.org/10.1007/s00011-021-01495-6.
Yadav P, Sharma P, Sundaram S, et al. SLC7A11/xCT is a target of miR-5096 and its restoration partially rescues miR-5096-mediated ferroptosis and anti-tumor effects in human breast cancer cells. Cancer Lett. 2021;2021(522):211–24. https://doi.org/10.1016/j.canlet.2021.09.033.
Wu J, Minikes AM, Gao M, et al. Intercellular interaction dictates cancer cell ferroptosis via NF2–YAP signalling. Nature. 2019;572(7769):402–6. https://doi.org/10.1038/s41586-019-1426-6.
Wenz C, Faust D, Linz B, et al. Cell–cell contacts protect against t-BuOOH-induced cellular damage and ferroptosis in vitro. Arch Toxicol. 2019;93(5):1265–79. https://doi.org/10.1007/s00204-019-02413-w.
Yang WH, Cornelia Ding CK, Sun T, et al. The hippo pathway effector TAZ regulates Ferroptosis in renal cell carcinoma. Cell Rep. 2019;28(10):2501–8. https://doi.org/10.1016/j.celrep.
Li Z, Jiang L, Chew SH, et al. Carbonic anhydrase 9 confers resistance to ferroptosis/apoptosis in malignant mesothelioma under hypoxia. Redox Biol. 2019;26: 101297. https://doi.org/10.1016/j.redox.2019.101297.
Zhang H, Deng T, Liu R, et al. CAF secreted miR-522 suppresses ferroptosis and promotes acquired chemo-resistance in gastric cancer. Mol Cancer. 2020;19(1):43–59. https://doi.org/10.1186/s12943-020-01168-8.
Lei G, Zhang Y, Koppula P, et al. The role of ferroptosis in ionizing radiation-induced cell death and tumor suppression. Cell Res. 2020;30(2):146–62. https://doi.org/10.1038/s41422-019-0263-3.
Yang J, Cao XH, Luan KF, Huang YD, et al. Circular RNA FNDC3B protects oral squamous cell carcinoma cells from Ferroptosis and contributes to the malignant progression by regulating miR-520d-5p/SLC7A11 Axis. Front Oncol. 2021;11:672724. https://doi.org/10.3389/fonc.2021.672724.
Fukuda M, Ogasawara Y, Hayashi H, et al. Down-regulation of glutathione peroxidase 4 in oral cancer inhibits tumor growth through SREBP1 Signaling. Anticancer Res. 2021;41(4):1785–92. https://doi.org/10.21873/anticanres.14944.
Kung JT, Colognori D, Lee JT. Long Noncoding RNAs: Past, Present, and Future. Genetics. 2013;193(3):651–69. https://doi.org/10.1534/genetics.112.146704.
Iyer MK, Niknafs YS, Malik R, et al. The landscape of long noncoding RNAs in the human transcriptome. Nat Genet. 2015;47(3):199–208. https://doi.org/10.1038/ng.3192.
Slack FJ, Chinnaiyan AM. The role of non-coding RNAs in oncology. Cell. 2019;179(5):1033–55. https://doi.org/10.1016/j.cell.2019.10.017.
Qiu YT, Li HB, Zhang Q, et al. Ferroptosis-related long noncoding RNAs as prognostic marker for colon adenocarcinoma. Appl Bionics Biomech. 2022;2022:5220368. https://doi.org/10.1155/2022/5220368.
Xu ZJ, Jiang SX, Ma J, et al. Comprehensive analysis of Ferroptosis-related LncRNAs in breast cancer patients reveals prognostic value and relationship with tumor immune microenvironment. Front Surg. 2021;8(742360):1–16. https://doi.org/10.3389/fsurg.2021.742360.
Biancur DE, Kapner KS, Yamamoto K, et al. Functional genomics identifies metabolic vulnerabilities in pancreatic cancer. Cell Metabol. 2021;33(1):199-210.e8. https://doi.org/10.1016/j.cmet.2020.10.018.
Zheng Q, Min S, Zhou Q. Identification of potential diagnostic and prognostic biomarkers for LUAD based on TCGA and GEO databases. Biosci Rep. 2021;41(6):BSR20204370. https://doi.org/10.1042/BSR20204370.
Xie JH, Li HB, Chen L, et al. A novel pyroptosis-related lncRNA signature for predicting the prognosis of skin cutaneous melanoma. Int J Gen Med. 2021;14:6517–27. https://doi.org/10.2147/IJGM.S335396.
Li T, Wang Y, Xiang XW, et al. Development and validation of a Ferroptosis-related lncRNAs prognosis model in oral squamous cell carcinoma. Front genet. 2022;13: 847940. https://doi.org/10.3389/fgene.2022.847940.
Lu R, Li ZY, Yi SC. Constructing a Ferroptosis-related long non-coding RNA signature to predict the prognostic of head and neck squamous cell carcinoma patients by bioinformatic analysis. Biochem Genet. 2022;60(5):1825–44. https://doi.org/10.1007/s10528-021-10176-2.
National Comprehensive Cancer Network Guidelines Head and Neck Version 1.2020. http://www.nccn.org/
Felice FD, Musio D, Terenzi V, et al. Treatment improvement and better patient care: which is the most important one in oral cavity cancer? Radiat Oncol. 2014;9(1):263–9. https://doi.org/10.1186/s13014-014-0263-x.
Jan JC, Hsu WH, Liu SA, et al. Prognostic factors in patients with buccal squamous cell carcinoma: 10-year experience. J Oral Maxillofac Surg. 2011;69(2):396–404. https://doi.org/10.1016/j.joms.2010.05.017.
Shah JP, Gil Z. Current concepts in management of oral cancer–surgery. Oral Oncol. 2009;45(4–5):394–440. https://doi.org/10.1016/j.oraloncology.2008.05.017.
Huang CM, Lin JX, Zheng CH, et al. Prognostic impact of dissected lymph node count on patients with node-negative gastric cancer. World J Gastroenterol. 2009;15(31):3926–30. https://doi.org/10.3748/wjg.15.3926.
Lisovsky M, Schutz SN, Drage MG, et al. Number of lymph nodes in primary nodal basin and a “Second Look” protocol as quality indicators for optimal nodal staging of colon cancer. Arch Pathol Lab Med. 2017;141(1):125–30. https://doi.org/10.5858/arpa.2015-0401-OA.
Tarantino I, Warschkow R, Hackert T, et al. Staging of pancreatic cancer based on the number of positive lymph nodes. Br J Surg. 2017;104(5):608–18. https://doi.org/10.1002/bjs.10472.
Blanchard P, Belkhir F, Temam S, et al. Outcomes and prognostic factors for squamous cell carcinoma of the oral tongue in young adults: a single-institution case-matched analysis. Eur Arch Oto-Rhino-L. 2016;274(3):1683–90. https://doi.org/10.1007/s00405-016-4419-1.
Lin NC, Hsu JT, Tsai KY. Difference between female and male patients with oral squamous cell carcinoma: a single-center retrospective study in Taiwan. Int J Env Res Pub He. 2020;17(11):3978. https://doi.org/10.3390/ijerph17113978.
An M, Zang XW, Wang JM, et al. Comprehensive analysis of differentially expressed long noncoding RNAs, miRNAs and mRNAs in breast cancer brain metastasis. Epigenomics. 2021;13(14):1113–28. https://doi.org/10.2217/epi-2021-0152.
Wang SQ, Wang Y, Wang SH, et al. Long non-coding RNA FIRRE acts as a miR-520a-3p sponge to promote gallbladder cancer progression via Mediating YOD1 expression. Front Genet. 2021;12(674653):1–13. https://doi.org/10.3389/fgene.2021.674653.
Shi X, Cui ZG, Liu XD, et al. LncRNA FIRRE is activated by MYC and promotes the development of diffuse large B-cell lymphoma via Wnt/β-catenin signaling pathway. Biochem Biophys Res Commun. 2021;510(4):594–600. https://doi.org/10.1016/j.bbrc.2019.01.105.
Wen YY, Fan XC, Bu HL, et al. Downregulation of lncRNA FIRRE relieved the neuropathic pain of female mice by suppressing HMGB1 expression. Mol Cell Biochem. 2021;476(2):841–52. https://doi.org/10.1007/s11010-020-03949-7.
Huang GW, Chen QQ, Ma CC, et al. Linc01305 promotes metastasis and proliferation of esophageal squamous cell carcinoma through interacting with IGF2BP2 and IGF2BP3 to stabilize HTR3A mRNA. Int J Biochem Cell Biol. 2021;136(106015):1–12. https://doi.org/10.1016/j.biocel.2021.106015.
Yan SP, Chu DX, Qiu HF, et al. LncRNA LINC01305 silencing inhibits cell epithelial-mesenchymal transition in cervical cancer by inhibiting TNXB-mediated PI3K/Akt signalling pathway. J Cell Mol Med. 2019;23(4):56–2666. https://doi.org/10.1111/jcmm.14161.
Yan F, Liu SW, Li XY, et al. Silencing LncRNA LINC01305 inhibits epithelial mesenchymal transition in lung cancer cells by regulating TNXB-mediated PI3K/Akt signaling pathway. J Biol Regul Homeost Agents. 2020;34(2):499–508. https://doi.org/10.23812/20-73-A-33.
Ishii N, Ozaki K, Sato H, et al. Identification of a novel non-coding RNA, MIAT, that confers risk of myocardial infarction. J Hum Genet. 2006;51(12):1087–99. https://doi.org/10.1007/s10038-006-0070-9.
Xu XB, Zhao YD, Duan G, et al. Downregulation of MIAT reduces the proliferation and migratory and invasive abilities of retinoblastoma cells by sponging miR-665 and regulating LASP1. Exp Ther Med. 2021;22(5):1342–52. https://doi.org/10.3892/etm.2021.10777.
Chuang TD, Quintanilla D, Boos D, et al. Long Noncoding RNA MIAT modulates the extracellular matrix deposition in Leiomyomas by Sponging MiR-29 Family. Endocrinology. 2021;162(11):bqab186. https://doi.org/10.1210/endocr/bqab186.
Zhu XW, Liu L, Wang Y, et al. LncRNA MIAT/HMGB1 axis is involved in cisplatin resistance via regulating IL6-mediated activation of the JAK2/STAT3 pathway in nasopharyngeal carcinoma. Front Oncol. 2021;11(651693):1–14. https://doi.org/10.3389/fonc.2021.651693.
Zhou JL, Zhang MY, Dong HH, et al. Comprehensive analysis of acetylation-related lncRNAs and Identified AC0998503 as prognostic biomarker in non-small cell lung cancer. J Oncol. 2021;2021(4405697):1–19. https://doi.org/10.1155/2021/4405697.
Jiang QK, Xue DF, Shi FZ, et al. Prognostic significance of an autophagy-related long non-coding RNA signature in patients with oral and oropharyngeal squamous cell carcinoma. Oncol Lett. 2021;21(1):29–38. https://doi.org/10.3892/ol.2020.12290.
Wang S, Chai KQ, Chen JB. A novel prognostic nomogram based on 5 long non-coding RNAs in clear cell renal cell carcinoma. Oncol Lett. 2019;18(6):6605–13. https://doi.org/10.3892/ol.2019.11009.
Stockwell BR, Friedmann AJP, Bayir H, et al. Ferroptosis: a regulated cell death nexus linking metabolism. Cell. 2017;171(1):273–85. https://doi.org/10.1016/j.cell.2017.09.021.
Samelson LE. Signal transduction mediated by the T cell antigen receptor: the role of adapter proteins. Annu Rev Immunol. 2002;20: 371394. https://doi.org/10.1146/annurev.
Call ME, Wucherpfennig KW. The T cell receptor: critical role of the membrane environment in receptor assembly and function. Annu Rev Immunol. 2005;23:101–25. https://doi.org/10.1146/annurev.
Garcia KC, Adams JJ, Feng D, et al. The molecular basis of TCR germline bias for MHC is surprisingly simple. Nat Immunol. 2009;10(2):143–7. https://doi.org/10.1038/ni.f.219.
Schamel WW, Alarcon B, Minguet S. The TCR is an allosterically regulated macromolecular machinery changing its conformation while working. Immunol Rev. 2019;291(1):8–25. https://doi.org/10.1111/imr.12788.
Cleary AS, Leonard TL, Gestl SA, et al. Tumour cell heterogeneity maintained by cooperating subclones in Wnt-driven mammary cancers. Nature. 2014;508:113–7. https://doi.org/10.1038/nature13187.
This work was supported by the National Nature Science Foundation of China (Grant numbers 81072214, 30371547) and the National Key R&D Program of China (Grant number 2016YFC1102603).
Ethics approval and consent to participate
The studies involving human participants were reviewed and approved by the Biomedical Ethics Committee of Peking University Stomatological Hospital. The ethical code is PKUSSIRB-202274058.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Qiu, L., Tao, A., Liu, F. et al. Potential prognostic value of a eight ferroptosis-related lncRNAs model and the correlative immune activity in oral squamous cell carcinoma. BMC Genom Data 23, 80 (2022). https://doi.org/10.1186/s12863-022-01097-z
- Oral squamous cell carcinoma
- Long non-coding RNAs
- Immune activity