- 1Department of Blood Transfusion, The Third Xiangya Hospital of Central South University, Changsha, China
- 2College of Engineering and Computer Science, Australian National University, Canberra, ACT, Australia
Lung adenocarcinoma (LUAD) is a highly heterogeneous malignancy, which makes prognosis prediction of LUAD very challenging. Ferroptosis is an iron-dependent cell death mechanism that is important in the survival of tumor cells. Long non-coding RNAs (lncRNAs) are considered to be key regulators of LUAD development and are involved in ferroptosis of tumor cells, and ferroptosis-related lncRNAs have gradually emerged as new targets for LUAD treatment and prognosis. It is essential to determine the prognostic value of ferroptosis-related lncRNAs in LUAD. In this study, we obtained RNA sequencing (RNA-seq) data and corresponding clinical information of LUAD patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and ferroptosis-related lncRNAs by co-expression analysis. The best predictors associated with LUAD prognosis, including C5orf64, LINC01800, LINC00968, LINC01352, PGM5-AS1, LINC02097, DEPDC1-AS1, WWC2-AS2, SATB2-AS1, LINC00628, LINC01537, LMO7DN, were identified by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis, and the LUAD risk prediction model was successfully constructed. Kaplan-Meier analysis, receiver operating characteristic (ROC) time curve analysis and univariate and multivariate Cox regression analysis and further demonstrated that the model has excellent robustness and predictive ability. Further, based on the risk prediction model, functional enrichment analysis revealed that 12 prognostic indicators involved a variety of cellular functions and signaling pathways, and the immune status was different in the high-risk and low-risk groups. In conclusion, a risk model of 12 ferroptosis related lncRNAs has important prognostic value for LUAD and may be ferroptosis-related therapeutic targets in the clinic.
Introduction
Lung cancer is the malignancy with the highest mortality rate worldwide, with approximately 25% of cancer patients dying from lung cancer (1). Non-small-cell lung cancer (NSCLC) is the most common type, accounting for approximately 85% of all lung cancers (2). NSCLC is divided into three main subtypes: squamous cell carcinoma, adenocarcinoma and large-cell carcinoma (3). Lung adenocarcinoma (LUAD) is the most prevalent histotype of NSCLC (4). Despite the rapid development of molecular biology and precision medicine, the 5-year overall survival (OS) rate for patients with LUAD is only 15% (2). Furthermore, because lung cancer is a highly heterogeneous tumor with complex molecular mechanisms, many molecular targeted drugs are ineffective in some patients, which poses a great challenge for the treatment of lung cancer (5, 6). Therefore, it is of great importance to identify specific tumor factors and provide new biomarkers for accurate diagnosis, individualized treatment, and prognosis prediction of LUAD.
Long non-coding RNAs (lncRNAs) are RNAs with a molecular weight of more than 200 nucleotides. Although lncRNAs are not involved in protein translation, they play a key role in gene regulation. In different biological and physiopathological contexts, lncRNAs can affect gene expression by regulating chromatin function, regulating the assembly and function of membrane-free nucleosomes, as well as changing the stability and translation of cytoplasmic mRNAs and interfering with signaling pathways (7). In recent years, with the development of high-throughput sequencing, people have discovered a large number of non-coding genes play an important role in the development and progression of tumors (8). Studies have demonstrated that lncRNAs are involved in cell growth, invasion and metastasis (9, 10), as well as tumor angiogenesis (11). Many lncRNAs are specifically expressed in different tumors and are involved in complex biological processes by interacting with DNA, RNA and proteins (12). The key role of lncRNAs in tumors indicates their promise as new targets for precise cancer therapy.
Ferroptosis is a unique cell death mechanism distinguished from apoptosis. Ferroptosis is characterized by iron-dependent peroxidation and is associated with oxidative polyunsaturated fatty acids (PUFAs), reactive oxygen species (ROS) and Lipoperoxide (13). In recent years, studies have demonstrated that ferroptosis is an important regulatory mechanism for tumor growth and is important for chemoradiotherapy and immunotherapy of tumors (14). Therefore, combinations with agents targeting iron death signaling could improve the anti-tumor efficacy of these therapeutic approaches. Peng Chen et al. demonstrated that Erianin, a natural product isolated from Dendrobium Chrysotoxum Lindl, can induce ferroptosis in lung cancer cells through Ca2 +/CaM signaling pathway and exert anti-tumor effects (15). In addition, it was demonstrated that lncRNA LINC00336 is highly expressed in lung cancer, and acts as a competitive endogenous RNA to function as an oncogene. Notably, the carcinogenesis of lncRNA LINC00336 is associated with ferroptosis (16). At present, the mechanism of ferroptosis-related lncRNAs in lung cancer progression has not been clarified, and their important significance in lung cancer treatment and prognosis needs to be further elucidated.
The tumor immune microenvironment (TIME) is associated with iron metabolism and homeostasis in vivo, while ferroptosis plays a key role in tumor immunity (17). Ferroptosis can expose tumor antigens thereby improving the immunogenicity of TIME and the effect of immunotherapy (18). Studies have shown that inhibition of system (a cystine/glutamate antiporter) expression by interferon-γ (IFNγ) and CT8 + T cells in TIME can induce ferroptosis in cancer cells (19, 20). The effect of ferroptosis on tumors involves a variety of immune cells and immune factors. Therefore, it is important to explore biomarkers associated with tumor immunity and ferroptosis for immunotherapy of lung cancer.
Aberrant expression of lncRNAs and abnormal alterations in ferroptosis are both common phenomena in tumor cells and are associated with tumor progression (14, 21). However, their direct interconnection and role in lung cancer are unknown and require further investigation. The purpose of this study is to explore the significance and molecular mechanism of ferroptosis-related lncRNAs in LUAD for the first time by bioinformatics methods and find biomarkers with potential utilization value so as to improve the clinical efficacy and improve the prognosis of LUAD patients. In this study, we used TCGA and GEO databases to search for ferroptosis-related lncRNAs that are differentially expressed in LUAD, and identified 12 lncRNAs that are closely related to LUAD prognosis and successfully constructed a LUAD risk prediction model. In addition, the acting mechanism of ferroptosis-related lncRNAs in tumor progression was further mined by functional analysis and immune infiltration analysis to provide new insights into the prognosis and immunotherapy of LUAD.
Materials and Methods
Patients and Datasets
The RNA sequencing (RNA-seq) data and corresponding clinical information of patients with LUAD were downloaded from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database. The data from TCGA database was used as a training cohort to establish the multi-lncRNA prognosis model, and the data from GEO (GSE30219, GES31210, GSE31546) as a validation cohort was used to test the predictive power of the risk score.
Identification of Ferroptosis-Related LncRNAs
Ferroptosis-related genes obtained by summarizing previous literature (22–24) and were provided in Table S1. First, the lncRNA and mRNA expression profiles were extracted separately from TCGA database. Then, the ferroptosis-related lncRNA profiles were excerpted from the ferroptosis-related genes based on the co-expression analysis (25). The ferroptosis-related lncRNAs was obtained through correlation analysis between the expression levels of the lncRNAs and ferroptosis-related genes, using Pearson correlation coefficients (coefficients > 0.40, P < 0.05). In addition, we explored regulatory mechanisms of ferroptosis-related lncRNA. We extracted ferroptosis-related lncRNA significantly correlated with co-expression and ferroptosis genes to construct the regulatory network and used Cytoscape software to display the ferroptosis genes and ferroptosis-related lncRNA regulatory network.
Construction of the Risk Score
First, differential expression analysis of ferroptosis-related lncRNAs was performed by using a “limma” R package, and differentially expressed lncRNAs were selected with an absolute log2-fold change (FC) > 1 and an adjusted P value < 0.01 (26). Next, overall survival was set as the clinical endpoint in our study. Then, univariate Cox regression was applied to identify ferroptosis-related lncRNAs related to overall survival for model construction. LncRNAs with a P value less than 0.01 were considered as significant prognostic signatures. A hazard ratio (HR) was used to determine whether lncRNAs are risk factors (HR>1) or protective factors (HR<1).
The set of differentially expressed lncRNAs and the set of prognosis-related ones were intersected, and differentially expressed lncRNAs with prognosis value were selected. Based on these genes, we used the Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis by using the “glmnet” R package to identify the best prognostic signature (27, 28). Specifically, the LASSO Cox regression was used to fit the overall survival of LUAD patients based on selected lncRNAs. During regression, the key lncRNAs were assigned non-zero coefficients and were selected to develop a risk score. Then, our risk score was constructed by using LASSO regression coefficients and expression levels of the key lncRNAs (29, 30). LUAD patients from the TCGA database, as a training cohort, were divided into high-risk group and low-risk group for further study based on the median score as the cutoff. Last, the patients in the GEO cohorts were also divided into two risk groups based on a cohort-specific cut-off.
Validation of the Risk Score
Using the “survival” R package and “survivalROC” R package, the area under the time dependent receiver operating characteristic (AUROC) were assessed to evaluate the predictive power of this prognostic model (26). Besides, the survival analysis of both high-risk group and low-risk group were performed by “Kaplan–Meier” and the “survival” package log-rank test (30). Based on the expression of prognostic ferroptosis-related lncRNA, principal component analysis (PCA) was performed with the “prcomp” function of the “stats” R package (31). Besides, t-distributed stochastic neighbor embedding (t-SNE) was implemented to explore the distribution of different groups using the “Rtsne” R package (23). The same analyses were performed on the validation cohort (the data from GEO database) for validation.
Independent Prognostic Analysis
Univariate and multivariate Cox regression analyses was performed to evaluate whether the prognostic model was independent of other traditional clinical characteristics (including age, gender, smoking, race and TNM stage) in predicting OS of patients with LUAD. The independent prognostic analysis was done on the TCGA cohort and on the GEO cohorts.
Functional Enrichment Analysis
After identifying the lncRNA signature correlating with LUAD prognosis, the gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathways were used to assess the biological roles of the prognostic candidates by the “clusterProfiler” R package (31). To further investigate the potential immunomodulatory mechanism of lncRNA in the regulation of tumor-infiltrating immune cells, the single-sample gene set enrichment analysis (ssGSEA) was performed to assess the infiltration abundance between the high-risk group and low-risk group of the TCGA cohort.
Immunogenomic Landscape Analyses
Computational methods were used to evaluate the immune infiltration and functions, including ESTIMATE (32), TIMER (33), MCP-counter (34), CIBERSORTx (35), and single-sample gene set enrichment analysis (ssGSEA) in an attempt to comprehensively analyze the immune differences between the two groups of the TCGA cohort (35). The immune infiltration and functions were compared between the high- and the low-risk groups by using the two-sample Wilcoxon test.
Results
Differentially Expressed mRNAs and LncRNAs Associated With Prognosis
Expression data and clinical information of four cohorts were downloaded from TCGA and GEO. The baseline characteristics were summarized in Table 1. The first cohort consisted of 535 LUAD and 59 normal samples from TCGA, and was used for differential expression analysis. 17122 lncRNAs were screened, and 282 lncRNAs were differentially expressed between LUAD and normal groups. The expression profiles of differentially expressed lncRNAs (DElncRNAs) were visualized in the form of volcano map and heatmap (Figures 1A and S1, respectively). As shown in Figures 1A and S1, 171 lncRNAs were upregulated, and 111 lncRNAs were downregulated in LUAD samples. 119 DElncRNAs were associated with the overall survival of LUAD patients using univariate Cox regression. The lncRNAs with the 10 highest and the 10 lowest Hazard Ratios (HRs) were reported in Figure 1B.
Figure 1 Differentially expressed ferroptosis-related lncRNA between LUAD and normal samples (log2 fold change>2, adjusted p-value < 0.01). Compared with normal samples, 49 ferroptosis-related lncRNAs were upregulated and 43 ferroptosis-related lncRNA downregulated in LUAD samples. (A) is a volcano map for differentially expressed ferroptosis-related lncRNAs. Red dots represent significantly upregulated expressed genes and green dots represent significantly downregulated expressed genes. (B) indicates the HR (95% Cl) and p-value of selected lncRNAs by univariate Cox proportional hazards. Blue dots represent protective factors, and red dots represent risk factors.
Construction of a Prognostic Model in the TCGA Cohort
12 ferroptosis-related lncRNAs were identified in the LASSO Cox regression model. According to this risk model, we used a special formula (Table S2) to calculate the risk score of each sample. Samples were classified into high- and low-risk groups by comparing their risk scores to the median score in the TCGA cohort and GEO cohort. As shown in Figures 2A, B, the areas under the time-dependent ROC of the TCGA cohort are 0.756, 0.713, and 0.694 for 1-, 3-, and 5-years survival, respectively. The areas under the time-dependent ROC of the GEO cohort are 0.728, 0.754, and 0.734 for 1-, 3-, and 5-years survival. The Kaplan Meier analysis with log-rank tests revealed a significant difference in the overall survival (p<0.01) between the two groups (shown in Figures 2C, D) in two cohorts. The patients in the high-risk group had a higher probability of death earlier than patients in the low-risk group. The univariate and multivariate Cox regression analyses were performed to evaluate whether clinical parameters and the risk score are independent prognostic factors of OS (Figures 2E–H). As shown in Figures 2G, H the stage and the risk score were proven to be independent factors for LUAD in TCGA cohort (HR>1, p<0.05). The age, gender, stage, and risk score were proven to be independent factors for LUAD in the GEO cohort (HR>1, p<0.05). We also indicated that LUAD patients in different risk groups were distributed in two directions by PCA and t-SNE analyses (Figure 3).
Figure 2 Prognostic analysis of the risk score in the TCGA and three GEO cohorts. (A, B) Time-dependent receiver operating characteristic curves assess the prognostic performance of the risk score in TCGA and three GEO cohorts respectively. (C, D) Kaplan-Meier curves display the overall survival of patients in the high- and low-risk group in TCGA and three GEO cohorts respectively. (E, F) The univariate Cox regression analysis of the associations between the risk scores and clinical parameters and the overall survival (OS) of patients in TCGA and three GEO cohort. (G, H) The multivariate Cox regression analysis of the associations between the risk scores and clinical parameters and the OS of patients in TCGA and three GEO cohorts.
Figure 3 (A, C) Principal component analysis (PCA) plot in the TCGA and three GEO cohorts. (B, D) t-distributed stochastic neighbor embedding (tSNE) analysis in the TCGA and three GEO cohorts.
Functional Enrichment Analyses
GO functional enrichment analyses, and KEGG pathway enrichment analyses were performed on the DElncRNAs between the high- and low-risk groups. The GO analysis results indicated that the 12 prognostic lncRNAs mainly focused on DNA replication, DNA-dependent DNA replication, ribosome biogenesis, and mitochondrial gene expression (Figures 4A, B). The KEGG analysis results show that the 12 prognostic lncRNAs are mainly enriched in DNA replication pathway, B cell receptor signaling pathways, hematopoietic cell lineage pathway, and cell cycle pathway (Figures 4C, D).
Figure 4 Further comparison between high- and low-risk groups of the TCGA cohort. (A) and (B) show the bar plot and cluster plot of significant GO functional items, respectively. (C) and (D) show the bar plot and cluster plot of significant KEGG pathways, respectively.
Immunity Analyses
To further explore the relationship between the LUAD prognosis and immune status, we quantified the infiltrating scores of immune cell- and immunity-related functions in two groups with ESTIMATE, TIMER, MCP counter, CIBERSORTx, and single-sample gene set enrichment analysis(ssGSEA) algorithms. There was a remarkable difference between the two groups, as shown in Figure 5. The ESTIMATE score was significantly higher in the low-risk group than in the high-risk group, indicating a higher overall immune level and immunogenicity of the tumor microenvironment (TME) in the low-risk group (p<0.001). Correlation analysis of immune cell subpopulations based on ssGSEA of TCGA-LUAD data revealed that the scores of immune cell, including aDCs, B cells, DCs, iDCs, Macrophages, NK-cells, T helper cells and Treg were significantly different between the low- and high-risk groups (p<0.001). In addition, there were significant differences in immune function scores including APC co-stimulation, inflammation promoting, check point, T cell co−inhibition T cell co−stimulation and Type II IFN Response between the low-risk and high-risk groups (p<0.01).
Figure 5 Comparison of the immunity analysis between different risk groups. (A–F) indicates the ESTIMATE, TIMER, MCP counter, CIBERSORTx, and single-sample gene set enrichment analysis (ssGSEA) algorithms to compare the cellular components or cell immune responses between two groups. CCR, cytokine-cytokine receptor. Adjusted P values were showed as: ns, not significant; *P < 0.05; **P < 0.01; ***P < 0.001.
Discussion
LUAD is characterized by a high degree of heterogeneity, and studies have demonstrated that a variety of genetic mutations and epigenetic changes drive the development and progression of LUAD, including gene fusion, major chromosomal events, single nucleotide changes (SNVs) and copy number alterations (CNAs) (36, 37). In recent years, with the in-depth understanding of the changes in the molecular mechanism of LUAD, many molecular targeted drugs for the precise treatment of LUAD have been developed, such as tyrosine kinase inhibitors (TLIs) of epidermal growth factor receptor (EGFR) and PDL1/PDL inhibitors (38, 39), which bring hope to LUAD patients. However, these mutations are not ubiquitous in all LUAD patients, and existing targeted drugs are only effective in some patients, and even some patients develop resistance and relapse during treatment. Therefore, further investigation of possible molecular mechanisms in LUAD is needed to provide new targets and biomarkers for individualized treatment and prognosis of LUAD. Recently, studies have demonstrated that the therapeutic effect of tumors can be improved by regulating the cell death process (40). Ferroptosis is an iron-dependent cell death program that has been shown to be involved in the development of tumors and the response to anti-tumor treatment (14, 19, 41). In addition, many lncRNAs are known to participate in malignant tumor progression and tumor resistance, and become new biomarkers and therapeutic targets for cancer diagnosis and treatment (42–44). Many lncRNAs can affect the ferroptosis of cancer cell and tumor immunity (45–47), being promising to become tumor markers with clinical value.
LncRNAs have been found to play a key role in ferroptosis of cancer cells. For example, lncRNAs engage the development and progression of NSCLC by regulating ferroptosis process (48). SLC7A11, the target gene of lncRNA, is a key gene of ferroptosis that can be downregulated by XAV939 (an NSCLC inhibitor), which inhibits the development of NSCLC through a ferroptosis mediated pathway (49). Jie Wu et al. found that lncRNA APCDC1L-AS can induce LUAD resistance to icotinib (a kind of TLIs) by inhibiting autophagic degradation of EGFR (44). The lncRNA P53RRA in the cytoplasm affects some metabolic gene transcription and promotes ferroptosis by activating P53 (45). P53RRA enhances erastin-induced growth inhibition while also increasing the concentration of intracellular iron and lipid ROS, which correlates with ferroptosis in NSCLC. LncRNALINC00336 is a key suppressor of ferroptosis and exerts antitumor effects in lung cancer by reducing intracellular iron and lipid ROS levels through interaction with ELAVL1 (50). These evidences suggest that ferroptosis-related lncRNAs are key regulators of ferroptosis and may be effective targets and specific markers for LUAD treatment and prognosis.
Here, we explored the association between ferroptosis-related lncRNAs and LUAD prognosis. The LUAD risk prediction model was successfully constructed with 12 ferroptosis-related lncRNAs, including five protective factors (C5orf64, LINC01800, LINC00968, LINC01352, PGM5-AS1) and seven risk molecules (LINC02097, DEPDC1-AS1, WWC2-AS2, SATB2-AS1, LINC00628, LINC01537, LMO7DN).C5orf64, LINC00968, WWC2-AS2, LINC00628, and LMO7DN have been shown to be associated with progression or prognosis of LUAD. We found that C5orf64 was the most significant protective factor, and Zhaofei Pang et al. found that C5orf64 was also an immune-related lncRNA that positively correlated with M2 macrophage, monocyte, eosinophil, and neutrophil infiltration in TIME of LUAD, but negatively correlated with Tregs and plasma cell infiltration; in addition, PD-1, PD-L1, and CTLA-4 were highly expressed in patients with high levels of C5orf64 (51). LINC00968 was prove to be associated with tumor progression and drug resistance (52–55). LINC00968 is under expressed in LUAD and inhibits tumor cell growth through the linc00968/miR-9-5p/CPEB3 (56) and miR-21-5p/SMAD7 (57) axes.LINC00628 is a tumor suppressor gene in multiple tumors, including gastric cancer (58), breast cancer (59), osteosarcoma (60) and hepatocellular carcinoma(HCC) (61). However, it was found that LINC00628 silencing inhibited proliferation, migration, and invasion of LUAD cells (62),which is consistent with our findings that LINC00628 acts as a risk factor for LUAD. Yonathan Brhane et al. identified that LMO7DN was associated with LUAD prognosis at the genome-wide level (63), suggesting that LMO7DN may play an important role in LUAD progression. Xiuqing Shen et al. demonstrated that WWC2-AS2 can be used as a predictor of LUAD, but its role in LUAD is unknown (64). Although there is little evidence of direct association of LINC01352, PGM5-AS1, and SATB2-AS1 with LUAD, many studies have shown that they play a key role in different tumors. Pinbo Huang et al. found that HBV/HBx (HBV X protein) complex can reduce the expression level of LINC01352 which in turn inhibits LINC01352-miR-135b-APC axis and promotes the development of HCC.LINC01352 is an independent prognostic molecule for HCC (65). PGM5-AS1 has different effects in different tumors. For example, PGM5-AS1 upregulates GDF10 gene expression by competitively binding to miR-587, inhibiting prostate cancer cell growth and accelerating apoptosis (66). In contrast, PGM5-AS1 can weaken the inhibitory effect of fibrillin-1 which mediated by miR-140-5P, and promote epithelial-mesenchymal transition, invasion and metastasis of osteosarcoma cells (67). Low SATB2-AS1 expression is associated with poor prognosis in colorectal cancer (CRC), and it can inhibit CRC metastasis and regulate TH1-type chemokine expression in TIME (68). However, we have not yet found studies on the significance of LINC02097, DEPDC1-AS1, and LINC01537 in LUAD or other tumors. Our findings demonstrate for the first time that these three lncRNAs are associated with the prognosis of LUAD, and their mechanism of action in LUAD needs to be explored.
In addition, we constructed the LUAD risk model with good predictive power and robustness, and it can be used as an independent prognostic factor for LUAD. This predictive model can effectively classify patients into low-risk group and high-risk group. Whether in the training cohort or validation cohort, we found a better prognosis in the low-risk group, which illustrates that 12 ferroptosis-related lncRNAs may be good prognostic factors for LUAD.
Due to different genetic mutations in different patients, targeted therapy is more likely to be a precise individualized treatment. We need to understand the molecular mechanisms that may occur in LUAD and find new therapeutic targets. Therefore, we further revealed the possible signaling pathways involved in ferroptosis-related lncRNAs by functional enrichment analysis, which are associated with a variety of life activities of cells, including cellular immunity, cell cycle, DNA replication, mitochondrial function and so on, indicating that ferroptosis-related lncRNAs may be involved in the development and progression of LUAD. The effect of immunotherapy depends on the immunogenicity of TME, so understanding TIME is the key to evaluating immunotherapy (69). By immune infiltration analysis, comparing to high-risk group, we found that the abundance of immune cells (dendritic cells, B cells, CD8 + cells, macrophages and neutrophils) in low-risk group was higher. Besides, a variety of immune system process pathways are involved, suggesting that low-risk patients have higher immunogenicity in TEM and better response to immunotherapy. The results of our study may help guide immunotherapy for LUAD.
However, our study has some shortcomings. First, based on traditional statistical methods, we constructed and evaluated 12 ferroptosis-related lncRNA risk prediction models. Although many studies have demonstrated the feasibility of these methods, more high-level techniques and methods need to be developed to improve the accuracy and robustness of prognostic models. In addition, this model has only been validated in GEO datasets, and external and practical validation is needed in the future to evaluate whether it can be applicable to clinical patients. Finally, we only preliminarily explored the signaling pathways involved in ferroptosis-related lncRNAs, however, the specific mechanism of ferroptosis-related lncRNAs in LUAD and their interconnection with TIME and ferroptosis are not yet fully understood, and more studies are needed to validate our findings.
Conclusion
In conclusion, we found 12 ferroptosis-related lncRNAs associated with LUAD prognosis. A risk score constructed based on these 12 biomarkers, can independently predict the prognosis of LUAD patients. Furthermore, by functional enrichment analysis, we explored the role of these 12 biomarkers in the immunity and ferroptosis of LUAD, providing new insights for further understanding of the molecular mechanisms in the development and progression of LUAD. These biomarkers have many potential therapeutic and prognostic implications for LUAD patient management.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: The Cancer Genome Atlas database (https://portal.gdc.cancer.gov/) and Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/).
Author Contributions
RG and Q-YZ conceived and designed the study. LL and L-PL provided equal contributions to research design, data analysis and article writing. Q-QZ revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (Grant No. 81573091) and the Fundamental Research Funds for the Central Universities of Central South University under Grant (No. 2020zzts892; 2021zzts1093).
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.
Acknowledgments
The authors gratefully acknowledge contributions from the GEO and TCGA network.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.675545/full#supplementary-material
References
1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer Statistics, 2021. CA Cancer J Clin (2021) 71:7–33. doi: 10.3322/caac.21654
2. Spella M, Stathopoulos GT. Immune Resistance in Lung Adenocarcinoma. Cancers (Basel) (2021) 13(3). doi: 10.3390/cancers13030384
3. Zappa C, Mousa SA. Non-Small Cell Lung Cancer: Current Treatment and Future Advances. Transl Lung Cancer Res (2016) 5:288–300. doi: 10.21037/tlcr.2016.06.07
4. Couraud S, Zalcman G, Milleron B, Morin F, Souquet PJ. Lung Cancer in Never Smokers–A Review. Eur J Cancer (2012) 48:1299–311. doi: 10.1016/j.ejca.2012.03.007
5. Calvayrac O, Pradines A, Pons E, Mazieres J, Guibert N. Molecular Biomarkers for Lung Adenocarcinoma. Eur Respir J (2017) 49(4). doi: 10.1183/13993003.01734-2016
6. Devarakonda S, Morgensztern D, Govindan R. Genomic Alterations in Lung Adenocarcinoma. Lancet Oncol (2015) 16:e342–51. doi: 10.1016/S1470-2045(15)00077-7
7. Statello L, Guo CJ, Chen LL, Huarte M. Gene Regulation by Long Non-Coding RNAs and Its Biological Functions. Nat Rev Mol Cell Biol (2021) 22:96–118. doi: 10.1038/s41580-020-00315-9
8. Gao N, Li Y, Li J, Gao Z, Yang Z, Li Y, et al. Long Non-Coding RNAs: The Regulatory Mechanisms, Research Strategies, and Future Directions in Cancers. Front Oncol (2020) 10:598817. doi: 10.3389/fonc.2020.598817
9. Hung CL, Wang LY, Yu YL, Chen HW, Srivastava S, Petrovics G, et al. A Long Noncoding RNA Connects c-Myc to Tumor Metabolism. Proc Natl Acad Sci U S A (2014) 111:18697–702. doi: 10.1073/pnas.1415669112
10. Kim T, Cui R, Jeon YJ, Lee JH, Lee JH, Sim H, et al. Long-Range Interaction and Correlation Between MYC Enhancer and Oncogenic Long Noncoding RNA CARLo-5. Proc Natl Acad Sci U S A (2014) 111:4173–8. doi: 10.1073/pnas.1400350111
11. Prensner JR, Zhao S, Erho N, Schipper M, Iyer MK, Dhanasekaran SM, et al. RNA Biomarkers Associated With Metastatic Progression in Prostate Cancer: A Multi-Institutional High-Throughput Analysis of SChLAP1. Lancet Oncol (2014) 15:1469–80. doi: 10.1016/S1470-2045(14)71113-1
12. Schmitt AM, Chang HY. Long Noncoding RNAs in Cancer Pathways. Cancer Cell (2016) 29:452–63. doi: 10.1016/j.ccell.2016.03.010
13. Perez MA, Magtanong L, Dixon SJ, Watts JL. Dietary Lipids Induce Ferroptosis in Caenorhabditiselegans and Human Cancer Cells. Dev Cell (2020) 54:447–54 e4. doi: 10.1016/j.devcel.2020.06.019
14. Chen X, Kang R, Kroemer G, Tang D. Broadening Horizons: The Role of Ferroptosis in Cancer. Nat Rev Clin Oncol (2021) 18(5):280–96. doi: 10.1038/s41571-020-00462-0
15. Chen P, Wu Q, Feng J, Yan L, Sun Y, Liu S, et al. Erianin, A Novel Dibenzyl Compound in Dendrobium Extract, Inhibits Lung Cancer Cell Growth and Migration Via Calcium/Calmodulin-Dependent Ferroptosis. Signal Transduct Target Ther (2020) 5:51. doi: 10.1038/s41392-020-0149-3
16. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, et al. Correction to: Long Noncoding RNA LINC00336 Inhibits Ferroptosis in Lung Cancer by Functioning as a Competing Endogenous RNA. Cell Death Differ (2020) 27:1447. doi: 10.1038/s41418-019-0394-6
17. Wang Y, Yu L, Ding J, Chen Y. Iron Metabolism in Cancer. Int J Mol Sci (2018) 20(1). doi: 10.3390/ijms20010095
18. Zhang F, Li F, Lu GH, Nie W, Zhang L, Lv Y, et al. Engineering Magnetosomes for Ferroptosis/Immunomodulation Synergism in Cancer. ACS Nano (2019) 13:5662–73. doi: 10.1021/acsnano.9b00892
19. Jiang X, Stockwell BR, Conrad M. Ferroptosis: Mechanisms, Biology and Role in Disease. Nat Rev Mol Cell Biol (2021) 22(4):266–82. doi: 10.1038/s41580-020-00324-8
20. Wang W, Green M, Choi JE, Gijon M, Kennedy PD, Johnson JK, et al. CD8(+) T Cells Regulate Tumour Ferroptosis During Cancer Immunotherapy. Nature (2019) 569:270–4. doi: 10.1038/s41586-019-1170-y
21. Gui CP, Cao JZ, Tan L, Huang Y, Tang YM, Li PJ, et al. A Panel of Eight Autophagy-Related Long Non-Coding RNAs Is a Good Predictive Parameter for Clear Cell Renal Cell Carcinoma. Genomics (2021) 113:740–54. doi: 10.1016/j.ygeno.2021.01.016
22. Li H, Li L, Xue C, Huang R, Hu A, An X, et al. A Novel Ferroptosis-Related Gene Signature Predicts Overall Survival of Breast Cancer Patients. Biology (Basel) (2021) 10(2). doi: 10.3390/biology10020151
23. Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, et al. A Novel Ferroptosis-Related Gene Signature for Overall Survival Prediction in Patients With Hepatocellular Carcinoma. Int J Biol Sci (2020) 16:2430–41. doi: 10.7150/ijbs.45050
24. Stockwell BR, Friedmann Angeli JP, Bayir H, Bush AI, Conrad M, Dixon SJ, et al. Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell (2017) 171:273–85. doi: 10.1016/j.cell.2017.09.021
25. Yuan M, Wang Y, Sun Q, Liu S, Xian S, Dai F, et al. Identification of a Nine Immune-Related LncRNA Signature as a Novel Diagnostic Biomarker for Hepatocellular Carcinoma. BioMed Res Int (2021) 2021:9798231. doi: 10.1155/2021/9798231
26. Li M, Liang M, Lan T, Wu X, Xie W, Wang T, et al. Four Immune-Related Long Non-coding RNAs for Prognosis Prediction in Patients With Hepatocellular Carcinoma. Front Mol Biosci (2020) 7:566491. doi: 10.3389/fmolb.2020.566491
27. Gupta S, Lee REC, Faeder JR. Parallel Tempering With Lasso for Model Reduction in Systems Biology. PloS Comput Biol (2020) 16:e1007669. doi: 10.1371/journal.pcbi.1007669
28. Zhang M, Zhu K, Pu H, Wang Z, Zhao H, Zhang J, et al. An Immune-Related Signature Predicts Survival in Patients With Lung Adenocarcinoma. Front Oncol (2019) 9:1314. doi: 10.3389/fonc.2019.01314
29. Lian P, Wang Q, Zhao Y, Chen C, Sun X, Li H, et al. An Eight-Long Non-Coding RNA Signature as a Candidate Prognostic Biomarker for Bladder Cancer. Aging (Albany NY) (2019) 11:6930–40. doi: 10.18632/aging.102225
30. Wei C, Liang Q, Li X, Li H, Liu Y, Huang X, et al. Bioinformatics Profiling Utilized a Nine Immune-Related Long Noncoding RNA Signature as a Prognostic Target for Pancreatic Cancer. J Cell Biochem (2019) 120:14916–27. doi: 10.1002/jcb.28754
31. Du X, Zhang Y. Integrated Analysis of Immunity- and Ferroptosis-Related Biomarker Signatures to Improve the Prognosis Prediction of Hepatocellular Carcinoma. Front Genet (2020) 11:614888. doi: 10.3389/fgene.2020.614888
32. Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
33. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for Analysis of Tumor-Infiltrating Immune Cells. Nucleic Acids Res (2020) 48:W509–W14. doi: 10.1093/nar/gkaa407
34. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol (2016) 17:218. doi: 10.1186/s13059-016-1070-5
35. Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, et al. Determining Cell Type Abundance and Expression From Bulk Tissues With Digital Cytometry. Nat Biotechnol (2019) 37:773–82. doi: 10.1038/s41587-019-0114-2
36. Campbell JD, Alexandrov A, Kim J, Wala J, Berger AH, Pedamallu CS, et al. Distinct Patterns of Somatic Genome Alterations in Lung Adenocarcinomas and Squamous Cell Carcinomas. Nat Genet (2016) 48:607–16. doi: 10.1038/ng.3564
37. Cancer Genome Atlas Research N. Comprehensive Molecular Profiling of Lung Adenocarcinoma. Nature (2014) 511:543–50. doi: 10.1038/nature13385
38. Imyanitov EN, Iyevleva AG, Levchenko EV. Molecular Testing and Targeted Therapy for Non-Small Cell Lung Cancer: Current Status and Perspectives. Crit Rev Oncol Hematol (2021) 157:103194. doi: 10.1016/j.critrevonc.2020.103194
39. Ortega-Franco A, Calvo V, Franco F, Provencio M, Califano R. Integrating Immune Checkpoint Inhibitors and Targeted Therapies in the Treatment of Early Stage Non-Small Cell Lung Cancer: A Narrative Review. Transl Lung Cancer Res (2020) 9:2656–73. doi: 10.21037/tlcr-20-546
40. Koren E, Fuchs Y. Modes of Regulated Cell Death in Cancer. Cancer Discov (2021) 11:245–65. doi: 10.1158/2159-8290.CD-20-0789
41. Yan HF, Zou T, Tuo QZ, Xu S, Li H, Belaidi AA, et al. Ferroptosis: Mechanisms and Links With Diseases. Signal Transduct Target Ther (2021) 6:49. doi: 10.1038/s41392-020-00428-9
42. Huang FT, Chen WY, Gu ZQ, Zhuang YY, Li CQ, Wang LY, et al. The Novel Long Intergenic Noncoding RNA UCC Promotes Colorectal Cancer Progression by Sponging miR-143. Cell Death Dis (2017) 8:e2778. doi: 10.1038/cddis.2017.191
43. Sun Y, Hu B, Wang Q, Ye M, Qiu Q, Zhou Y, et al. Long Non-Coding RNA HOTTIP Promotes BCL-2 Expression and Induces Chemoresistance in Small Cell Lung Cancer by Sponging miR-216a. Cell Death Dis (2018) 9:85. doi: 10.1038/s41419-017-0113-5
44. Wu J, Zheng C, Wang Y, Yang Z, Li C, Fang W, et al. LncRNA APCDD1L-AS1 Induces Icotinib Resistance by Inhibition of EGFR Autophagic Degradation Via the miR-1322/miR-1972/miR-324-3p-SIRT5 Axis in Lung Adenocarcinoma. Biomark Res (2021) 9:9. doi: 10.1186/s40364-021-00262-3
45. Mao C, Wang X, Liu Y, Wang M, Yan B, Jiang Y, et al. A G3BP1-Interacting LncRNA Promotes Ferroptosis and Apoptosis in Cancer Via Nuclear Sequestration of P53. Cancer Res (2018) 78:3484–96. doi: 10.1158/0008-5472.CAN-17-3454
46. Qi W, Li Z, Xia L, Dai J, Zhang Q, Wu C, et al. LncRNA GABPB1-AS1 and GABPB1 Regulate Oxidative Stress During Erastin-Induced Ferroptosis in HepG2 Hepatocellular Carcinoma Cells. Sci Rep (2019) 9:16185. doi: 10.1038/s41598-019-52837-8
47. Wang Z, Chen X, Liu N, Shi Y, Liu Y, Ouyang L, et al. A Nuclear Long Non-Coding RNA LINC00618 Accelerates Ferroptosis in a Manner Dependent Upon Apoptosis. Mol Ther (2021) 29:263–74. doi: 10.1016/j.ymthe.2020.09.024
48. Wu Y, Zhang S, Gong X, Tam S, Xiao D, Liu S, et al. The Epigenetic Regulators and Metabolic Changes in Ferroptosis-Associated Cancer Progression. Mol Cancer (2020) 19:39. doi: 10.1186/s12943-020-01157-x
49. Yu H, Han Z, Xu Z, An C, Xu L, Xin H. RNA Sequencing Uncovers the Key Long Non-Coding RNAs and Potential Molecular Mechanism Contributing to XAV939-mediated Inhibition of Non-Small Cell Lung Cancer. Oncol Lett (2019) 17:4994–5004. doi: 10.3892/ol.2019.10191
50. Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, et al. Long Noncoding RNA LINC00336 Inhibits Ferroptosis in Lung Cancer by Functioning as a Competing Endogenous RNA. Cell Death Differ (2019) 26:2329–43. doi: 10.1038/s41418-019-0304-y
51. Pang Z, Chen X, Wang Y, Wang Y, Yan T, Wan J, et al. Long Non-Coding RNA C5orf64 Is a Potential Indicator for Tumor Microenvironment and Mutation Pattern Remodeling in Lung Adenocarcinoma. Genomics (2021) 113:291–304. doi: 10.1016/j.ygeno.2020.12.010
52. Liu G, Yuan D, Sun P, Liu W, Wu PF, Liu H, et al. LINC00968 Functions as an Oncogene in Osteosarcoma by Activating the PI3K/AKT/mTOR Signaling. J Cell Physiol (2018) 233:8639–47. doi: 10.1002/jcp.26624
53. Sun X, Huang T, Zhang C, Zhang S, Wang Y, Zhang Q, et al. Long Non-Coding RNA LINC00968 Reduces Cell Proliferation and Migration and Angiogenesis in Breast Cancer Through Up-Regulation of PROX1 by Reducing hsa-miR-423-5p. Cell Cycle (2019) 18:1908–24. doi: 10.1080/15384101.2019.1632641
54. Xiu DH, Liu GF, Yu SN, Li LY, Zhao GQ, Liu L, et al. Long Non-Coding RNA LINC00968 Attenuates Drug Resistance of Breast Cancer Cells Through Inhibiting the Wnt2/Beta-Catenin Signaling Pathway by Regulating WNT2. J Exp Clin Cancer Res (2019) 38:94. doi: 10.1186/s13046-019-1100-8
55. Yao N, Sun JQ, Yu L, Ma L, Guo BQ. LINC00968 Accelerates the Progression of Epithelial Ovarian Cancer Via Mediating the Cell Cycle Progression. Eur Rev Med Pharmacol Sci (2019) 23:4642–9. doi: 10.26355/eurrev_201906_18043
56. Tang H, Han X, Feng Y, Hao Y. linc00968 Inhibits the Tumorigenesis and Metastasis of Lung Adenocarcinoma Via Serving as a ceRNA Against miR-9-5p and Increasing CPEB3. Aging (Albany NY) (2020) 12:22582–98. doi: 10.18632/aging.103833
57. Zhu Y, Bo H, Chen Z, Li J, He D, Xiao M, et al. LINC00968 can Inhibit the Progression of Lung Adenocarcinoma Through the miR-21-5p/SMAD7 Signal Axis. Aging (Albany NY) (2020) 12:21904–22. doi: 10.18632/aging.104011
58. Zhang ZZ, Zhao G, Zhuang C, Shen YY, Zhao WY, Xu J, et al. Long Non-Coding RNA LINC00628 Functions as a Gastric Cancer Suppressor Via Long-Range Modulating the Expression of Cell Cycle Related Genes. Sci Rep (2016) 6:27435. doi: 10.1038/srep27435
59. Chen DQ, Zheng XD, Cao Y, He XD, Nian WQ, Zeng XH, et al. Long Non-Coding RNA LINC00628 Suppresses the Growth and Metastasis and Promotes Cell Apoptosis in Breast Cancer. Eur Rev Med Pharmacol Sci (2017) 21:275–83.
60. He R, Wu JX, Zhang Y, Che H, Yang L. LncRNA LINC00628 Overexpression Inhibits the Growth and Invasion Through Regulating PI3K/Akt Signaling Pathway in Osteosarcoma. Eur Rev Med Pharmacol Sci (2018) 22:5857–66. doi: 10.26355/eurrev_201809_15915
61. Chen Q, Wang D, Li Y, Yan S, Dang H, Yue H, et al. LINC00628 Suppresses Migration and Invasion of Hepatocellular Carcinoma by its Conserved Region Interacting With the Promoter of VEGFA. J Cell Physiol (2019) 234(9):15751–62. doi: 10.1002/jcp.28233
62. Xu SF, Zheng Y, Zhang L, Wang P, Niu CM, Wu T, et al. Long Non-coding RNA LINC00628 Interacts Epigenetically With the LAMA3 Promoter and Contributes to Lung Adenocarcinoma. Mol Ther Nucleic Acids (2019) 18:166–82. doi: 10.1016/j.omtn.2019.08.005
63. Brhane Y, Yang P, Christiani DC, Liu G, McLaughlin JR, Brennan P, et al. Genetic Determinants of Lung Cancer Prognosis in Never Smokers: A Pooled Analysis in the International Lung Cancer Consortium. Cancer Epidemiol Biomarkers Prev (2020) 29:1983–92. doi: 10.1158/1055-9965.EPI-20-0248
64. Shen XQ, Wu QM, Yang CH, Yan QD, Cao PJ, Chen FL. Four Low Expression LncRNAs Are Associated With Prognosis of Human Lung Adenocarcinoma. Clin Lab (2020) 66(10). doi: 10.7754/Clin.Lab.2020.200211
65. Huang P, Xu Q, Yan Y, Lu Y, Hu Z, Ou B, et al. HBx/Eralpha Complex-Mediated LINC01352 Downregulation Promotes HBV-Related Hepatocellular Carcinoma Via the miR-135b-APC Axis. Oncogene (2020) 39:3774–89. doi: 10.1038/s41388-020-1254-z
66. Du L, Gao Y. PGM5-AS1 Impairs miR-587-Mediated GDF10 Inhibition and Abrogates Progression of Prostate Cancer. J Transl Med (2021) 19:12. doi: 10.1186/s12967-020-02572-w
67. Liu W, Liu P, Gao H, Wang X, Yan M. Long Non-Coding RNA PGM5-AS1 Promotes Epithelial-Mesenchymal Transition, Invasion and Metastasis of Osteosarcoma Cells by Impairing miR-140-5p-Mediated FBN1 Inhibition. Mol Oncol (2020) 14:2660–77. doi: 10.1002/1878-0261.12711
68. Xu M, Xu X, Pan B, Chen X, Lin K, Zeng K, et al. LncRNA SATB2-AS1 Inhibits Tumor Metastasis and Affects the Tumor Immune Cell Microenvironment in Colorectal Cancer by Regulating SATB2. Mol Cancer (2019) 18:135. doi: 10.1186/s12943-019-1063-6
Keywords: lung adenocarcinoma, ferroptosis, lncRNAs, prognosis, risk score
Citation: Lu L, Liu L-P, Zhao Q-Q, Gui R and Zhao Q-Y (2021) Identification of a Ferroptosis-Related LncRNA Signature as a Novel Prognosis Model for Lung Adenocarcinoma. Front. Oncol. 11:675545. doi: 10.3389/fonc.2021.675545
Received: 03 March 2021; Accepted: 08 June 2021;
Published: 23 June 2021.
Edited by:
Wilfried Ernst Erich Eberhardt, Essen University Hospital, GermanyReviewed by:
Daniel Christoph, Kliniken Essen-Mitte, GermanyMarcel Wiesweg, Universitätsklinikum Essen, Germany
Copyright © 2021 Lu, Liu, Zhao, Gui and Zhao. 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: Rong Gui, Z3Vpcm9uZ0Bjc3UuZWR1LmNu; Qin-Yu Zhao, cWlueXUuemhhb0BhbnUuZWR1LmF1
†These authors have contributed equally to this work and share first authorship