Front. Oncol. Frontiers in Oncology Front. Oncol. 2234-943X Frontiers Media S.A. 10.3389/fonc.2022.844642 Oncology Original Research Development and Validation of Ferroptosis- and Immune-Related lncRNAs Signatures for Breast Infiltrating Duct and Lobular Carcinoma Wei Tao 1 † Zhu Ning 2 † Jiang Weihua 3 Xing Xiao-Liang 2 * 1 Department of Surgical Oncology, Urumqi Friendship Hospital, Urumqi, China 2 School of Public Health and Laboratory Medicine, Hunan Provincial Key Laboratory for Synthetic Biology of Traditional Chinese Medicine, Hunan University of Medicine, Huaihua, China 3 Department of Breast Surgery, The Affiliated Tumor Hospital of Xinjiang Medical University, Urumqi, China

Edited by: Karen Elizabeth Nava-Castro, National Autonomous University of Mexico, Mexico

Reviewed by: Shengchun Liu, First Affiliated Hospital of Chongqing Medical University, China; Feng Wang, Affiliated Hospital of Nantong University, China

*Correspondence: Xiao-Liang Xing, xiaoliangxinghnm@126.com

†These authors have contributed equally to this work

This article was submitted to Breast Cancer, a section of the journal Frontiers in Oncology

 04 04 2022 2022 12 844642 28 12 2021 09 03 2022 Copyright © 2022 Wei, Zhu, Jiang and Xing 2022 Wei, Zhu, Jiang and Xing

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.

 Background

Heterogeneity of breast cancer (BRCA) is significantly correlated with its prognosis. Target therapy for ferroptosis and immunity is a new cancer treatment option discovered in recent years. In the present study, we aimed to identify ferroptosis- and immune-related long non-coding RNAs (lncRNAs) to accurately predict the prognosis and diagnosis of patients with breast infiltrating duct and lobular carcinoma by integrated analyses.

 Methods

The corresponding data for the patients with breast infiltrating duct and lobular carcinoma by integrated analyses were obtained from The Cancer Genome Atlas (TCGA). Analyses of univariate and multivariate Cox regressions were used to identify the suitable candidate biomarkers.

 Results

We found that seven ferroptosis- and immune-related differentially expressed lncRNAs (FI-DELs) (AC007686.3, AC078883.1, ADAMTS9-AS1, AL035661.1, CBR3-AS1, FTX, and TMEM105) were correlated with the overall survival of patients with breast infiltrating duct and lobular carcinoma. The areas under the receiver operating characteristic (AUCs) value of the prognosis model were all over 0.6 in training, validation, and entire groups. The sensitivity and specificity of the diagnosis model was 87.84% and 97.06%, respectively.

 Conclusions

Through a series of bioinformatics analyses, we found that the seven FI-DELs could serve as prognostic and diagnostic biomarkers for patients with breast infiltrating duct and lobular carcinoma. However, whether these seven biomarkers could be really applied to the clinic requires further investigations.

 breast infiltrating duct and lobular carcinoma ferroptosis immune lncRNAs prognosis diagnosis Introduction

Breast cancer (BRCA) is the most common type of malignant tumor. Globally, BRCA has very high incidence rates (1). It is established that there are almost 2.3 million new cases and 70 thousands deaths in 2020 (1). As a benefit from the development of BRCA diagnosis and treatment, the 5-year survival rate of BRCA patients has increased to 90% in recent years (2–4). However, BRCA is often diagnosed at later stage because it is asymptomatic early, which can lead to serious consequences (5). Previous studies have also demonstrated that there are significant differences in the prognosis of BRCA due to the heterogeneous of tumor cell (6, 7). Therefore, it is necessary and important to identify novel prognostic and diagnosis signatures for different types of BRCA to guide clinical practice and improve overall survival.

Previous studies indicated that BRCA easily metastasizes to other organs, such as the bones and lungs (5). For non-metastatic BRCA, the primary goal of treatment is to eradicate the tumor from the breast and regional lymph nodes and prevent metastatic recurrence (8). For metastatic BRCA, treatment is aimed at prolonging life and relieving symptoms. At present, metastatic BRCA remains incurable in almost all patients. Neoadjuvant/adjuvant therapy has been used primarily in metastatic breast cancer (8). Evidence indicated that the response to neoadjuvant treatment and the prognosis of BRCA are positively influenced by the tumor-infiltrating lymphocytes (9, 10). The development and progression of BRCA are closely related to the immune microenvironment, such as CD8+ and CD4+ T cells (11–13). Ferroptosis, an iron-dependent form of regulated cell death, is first proposed by Dixon in 2012 (14). Cumulative evidence indicted that the ferroptosis and immune system can regulate each other and participate in the development and progression of cancers (15–20).

Long non-coding RNA (lncRNA) is a class of non-coding RNA molecules with transcription length over 200 nt (21). Previous studies have found that lncRNAs are closely involved in the development and progression of several cancers (including BRCA) in various ways and may serve as specific biomarkers (22–25). In this study, we aimed to identify suitable ferroptosis- and immune-related lncRNAs as prognosis and diagnosis biomarkers for patients with breast infiltrating duct and lobular carcinoma by integrated analyses.

 Material and Methods Data Acquisition

The RNA-seq counts data and their corresponding clinical information for 1,215 samples (102 normal and 826 patients with breast infiltrating duct and lobular carcinoma) were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). The list of recognized ferroptosis- and immune-related genes were downloaded from the database of FerrDb and ImmPort, respectively. The annotated lncRNAs were downloaded from Gene Transfer Format (GTF) file (http://asia.ensembl.org).

 Candidate Biomarkers Identification

The DESeq2 package in R (3.6.1) was used to screen the differentially expressed genes with the following criterion adj. p < 0.05, |logFC| ≥ 0.5 and base mean ≥ 50. Co-expression analyses for the genes and lncRNAs were determined by Pearson analyses with the following criteria: p < 0.05 and R ≥ 0.3. After dividing the patients with breast infiltrating duct and lobular carcinoma into two groups depending on the median value, univariate and multivariate Cox regressions were used to identify the overall survival (OS)-related signatures.

 Prognosis and Diagnosis Model Construction

According to previous studies, we constructed the risk assessment model and the diagnosis model (26, 27). Risk value = (−1.08) × Exp (AC007686.3) + (−0.85) × Exp (AC078883.1) + (−1.03) × Exp (ADAMTS9-AS1) + (1.48) × Exp (AL035661.1) + (1.18) × Exp (CBR3-AS1) + (1.56) × Exp (FTX) + (0.86) × Exp (TMEM105). The Logit value = 1.023 + (−0.150) × Exp (AC007686.3) + (−0.072) × Exp (AC078883.1) + (−0.132) × Exp (ADAMTS9-AS1) + (0.051) × Exp (AL035661.1) + (0.087) × Exp (CBR3-AS1) + (−0.104) × Exp (FTX) + (0.053) × Exp (TMEM105). The Express (Exp) values of the candidate biomarkers were obtained from the normalized value of DESeq2 analyses.

 Principal Component Analyses and Enrichment Analyses

Principal component analyses (PCAs) were conducted to reduce dimension and visualize the patients with breast infiltrating duct and lobular carcinoma with different risk values in R (3.6.1).

DAVID bioinformatics Resource 6.8 was used to perform the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses with the default parameters (28).

 Statistical Analyses

A repeated measures ANOVA followed by unpaired two-tailed Student’s t-test was used as indicated. All results are expressed as mean ± SEM.

 Results Identification of Ferroptosis- and Immune-Related lncRNAs as Biomarkers

First, 16,901 lncRNAs were identified by analyzing the RNA-seq data from patients with breast infiltrating duct and lobular carcinoma, of which 741 lncRNAs were differentially expressed between the normal and the cancer patients. A total of 1,038 prognostic lncRNAs were identified by univariate Cox regression analyses in all patients with breast infiltrating duct and lobular carcinoma. To obtain ferroptosis- and immune-related lncRNAs, we downloaded 329 ferroptosis-related genes and 2,483 immune-related genes from the database of FerrDb and ImmPort, respectively. We then performed co-expression analyses for these ferroptosis- and immune-related genes with the lncRNAs and found that 8,997 ferroptosis- and immune-related lncRNAs correlated with 322 ferroptosis- and immune-related genes. Finally, we found that 47 ferroptosis- and immune-related differentially expressed lncRNAs (FI-DELs) were present among differentially expressed lncRNAs, prognostic lncRNAs, and ferroptosis- and immune-related lncRNAs ( Figure 1A ). The expression of these 47 FI-DELs are shown in Figure 1B . To construct a risk assessment model and further verify it, 826 patients with breast infiltrating duct and lobular carcinoma were randomly separated into training (n=413) or validation (n=413) group. Table 1 presents the detail clinical characteristics of these patients.

Identification of ferroptosis- and immune-related lncRNAs as biomarkers. (A) Venn diagram to identify the common lncRNAs of differentially expressed lncRNAs, ferroptosis- and immune-related lncRNAs, and prognostic lncRNAs. (B) Heatmap of 47 common lncRNAs identified by Venn analyses. (C) Multivariate Cox regression analyses illustrated that seven FI-DELs were correlated with the OS of patients with breast infiltrating duct and lobular carcinoma independently. (D–J) Overall survival curve of these seven FI-DELs. (K) The correlation of these seven FI-DELs with the FI-DGEs.

Clinical characteristic of patients in the training and validation groups.

Variables Training group (n=413) Validation group (n=413)
No. % No. %
Age
 ≤60 237 57.38% 240 58.11%
 >60 176 42.62% 173 41.89%
Stage
 I 86 20.82% 62 15.01%
 II 225 54.48% 250 60.53%
 III 83 20.10% 83 20.10%
 IV 12 2.91% 6 1.45%
 X 7 1.69% 12 2.91%
T
 T1 129 31.23% 104 25.18%
 T2 229 55.45% 263 63.68%
 T3 36 8.72% 31 7.51%
 T4 18 4.36% 14 3.39%
 TX 1 0.24% 1 0.24%
N
 N0 191 46.25% 188 45.52%
 N1 143 34.62% 146 35.35%
 N2 47 11.38% 54 13.08%
 N3 24 5.81% 17 4.12%
 NX 8 1.94% 8 1.94%
M
 M0 356 86.20% 357 86.44%
 M1 12 2.91% 8 1.94%
 MX 45 10.90% 48 11.62%
Gender
 Female 406 98.30% 409 99.03%
 Male 7 1.70% 4 0.97%

Through multivariate Cox regression analyses in the training group, we found that seven FI-DELs were correlated independently with the overall survival of patients with breast infiltrating duct and lobular carcinoma ( Figures 1C–J ). The interactions of these seven FI-DELs with the FI-DEGs are displayed in Figure 1K .

 Construction and Validation of the Prognostic Model

According to previous studies, we first constructed a risk assessment model using these seven FI-DELs in the training group. The patients were further divided into high- and low-risk groups by the Youden index as the optimal cutoff value ( Supplementary Figure S1 ). The risk value (up) and survival status (down) for each patients in the training group was displayed in Figure 2A . The expression of AC078883.1 and ADAMTS9-AS1 were significantly decreased, while the expression of the AL035661.1, CBR3-AS1, FTX, and TMEM105 were significantly increased in the patient with breast infiltrating duct and lobular carcinoma with high-risk value ( Figure 2B ). Patients with breast infiltrating duct and lobular carcinoma with high-risk value had a higher probability of death than these patients with the low-risk value ( Figure 2C ). As shown in Figure 2D , the area under ROC curve (AUC) reached 0.67. The time-dependent AUC reached 0.75 at 1 year, 0.71 at 3 years, 0.63 at 5 years, and 0.68 at 10 years.

Developments and validations of prognosis model. (A, E, I) Risk value (up) and survival status (down) for each patients in the training group (A), validation group (E), and entire group (I). (B, F, J) Expression of these seven candidate biomarkers in the training group (B), validation group (F), and entire group (J). (C, G, K) Overall survival curve of patients with breast infiltrating duct and lobular carcinoma with different risk value in the training group (C), validation group (G), and entire group (J). (D, H, L) ROC curve of risk assessment model in the training group (D), validation group (H), and entire group (L). *p < 0.05, ***p < 0.001. ns, no significance.

We then performed similar investigations for the patients in the validation and entire groups. Similar results are shown in Figures 2E–L . All AUC values of the prognosis model in the training, validation, and entire group were over 0.60.

 Independent Prognostic Value of Risk Assessment Model

To determine whether the risk assessment model was an independent prognostic factor for the patients with breast infiltrating duct and lobular carcinoma, univariate and multivariate Cox regression analyses were performed among the clinical characteristics and the prognostic model. The results of the univariate Cox regression showed that the risk assessment model was significantly associated with the OS in the training, validation, and entire groups ( Figure 3A ; Supplementary Figures S2A, D ). The risk assessment model remained an independent predictor of OS in multivariate Cox regression analyses ( Figure 3B ; Supplementary Figures S2B, E ). Additionally, the ROC curve analyses showed that the AUC value of the risk assessment model in the training group, validation group, and entire group were 0.67, 0.60, and 0.63 respectively ( Figures 2D, H, L ), which were higher than the AUC values of other prognostic factors ( Figure 3C ; Supplementary Figures S2C, F ).

Independent prognostic values of these four FI-DELs. (A) Univariate Cox regression analyses illustrated four clinical characteristics and risk assessment model were correlated with OS of patients. (B) Multivariate Cox regression analyses illustrated that three clinical characteristics and risk assessment model were correlated with the OS of patients independently. (C) ROC curve of different clinical characteristics and risk assessment model. (D–L) Comparison of the risk score in different clinical characteristics groups. (M) Correlation analyses of these seven FI-DELs with the clinical characteristics. *p < 0.05, **p < 0.01, ***p < 0.001.

To investigate the relationships of these seven FI-DELs with different clinical characteristics, we performed differential expression analyses and found that the risk values displayed significant difference among the patients with different age, gender, pathological TN, and pathological stage ( Figures 3D–L ). We further performed the correlation analyses and found that AC078883.1, CBR3-AS1, and AL035661.1 were significantly correlated with the risk value ( Figure 3M ).

 The Immune Cell Infiltration Landscape in Breast Cancer

To further explore the relationship between FI-DELs and the immunity, we evaluated the immunity status using ESTIMATE in R (3.6.1). The estimate and stromal scores were significantly decreased, while the immune score and tumor purity was significantly increased in patients with breast infiltrating duct and lobular carcinoma ( Supplementary Figures S3A–D ). The estimate, immune, and stromal scores were significantly decreased, while the tumor purity was significantly increased in patients with breast infiltrating duct and lobular carcinoma with different risk score ( Figures 4A–D ). The correlation of immune status with these seven FI-DELs is displayed in Figure 4I .

Correlation analyses of risk assessment model with the immunity. (A–D) Comparison of the estimate score (A), immune score (B), stromal score (C), and tumor purity (D) between the patients with different risk value. (E–H) Comparison of the infiltrating score between the patients with different risk value. (E) B cell. (F) Macrophage. (G) Myeloid dendritic cell. (H) T cell CD8+. (I, J) Correlation analyses of these seven FI-DELs with the ESTIMATE score (I) and infiltrating score (J). **p < 0.01, ***p < 0.001.

Subsequently, we used the TIMER algorithm to investigate the landscape of the immune cell infiltration. The infiltration score of the macrophage, myeloid dendritic cell, and T cell CD8+ were significantly decreased, while the infiltration score of B cell was significantly increased in the patients with breast infiltrating duct and lobular carcinoma ( Supplementary Figures S3E–H ). The infiltration score of these four immune cells and factors were significantly decreased in patients with breast infiltrating duct and lobular carcinoma with high-risk value ( Figures 4E–H ). The correlation of immune cells and factors with these seven FI-DELs is displayed in Figure 4J .

 Immune-Related Pathways Were Activated in the Low-Risk Group

To explore the distribution of the patients with different risk score, we performed the PCA analyses using these seven FI-DELs and found that the patients with breast infiltrating duct and lobular carcinoma with low-risk value could be well separated from these patients with high-risk value in the training, validation, and entire groups ( Figures 5A–C ).

Principal component analyses and enrichment analyses. (A–C) Principal component analyses (PCA) plots depicted distinct distribution of high- and low-risk groups based on these seven FI-DEGs in the training group (A), validation group (B), and entire group (C). Blue means low-risk group. Red means high-risk group. (D–F) Top 10 significantly enriched GO. (D) For BP. (E) For CC. (F) For MF. (G) Top 10 significantly enriched KEGG.

We subsequently performed differentially expressed analyses for these patients with different risk score. A total of 2,582 genes (660 upregulated DEGs and 1,922 downregulated DGEs) were differentially expressed between the low- and high-risk group ( Supplementary Figure SS4 ). An enrichment analyses was carried out for these 2,582 DEGs using DAVID. A total of 563 BP, 89 CC, and 143 MF were enriched as measured by the false discovery rate (FDR) value <0.05 ( Supplementary Table S1 ), and the top 10 of these GOs are displayed in Figures 5D–F , including several immune response GO term. KEGG analyses showed that 110 pathways were enriched as measured by the FDR value <0.05 ( Supplementary Table S2 ), the top 10 of these pathways are displayed in Figure 5G , including ferroptosis- and immune-related pathway PI3K-Akt signaling pathway.

 Construction of a Diagnostic Model

To know the role in the diagnosis of patients with breast infiltrating duct and lobular carcinoma, we performed a stepwise logistic regression analyses for these seven FI-DELs ( Figure 6A ). The logit score of patients was significantly higher than that of the normal ( Figure 6B ). Correlation analyses showed that the logit score was significantly associated with AC078883.1, FTX, ADAMTS9-AS1, and AC007686.3 ( Figure 6C ). The sensitivity and specificity of the diagnose model were 87.53% and 97.06%, respectively ( Table 2 ). We also plot the ROC curve of the diagnose model, and the AUC value was 0.9702 ( Figure 6C ).

Construction of diagnosis model. (A) β-value of these seven FI-DELs analyzed by stepwise logistic regression. (B) Diagnosis values between normal and patients with breast infiltrating duct and lobular carcinoma. (C) Correlation analyses of these seven FI-DELs with diagnosis values. (D) ROC curves of the diagnostic model. ***p < 0.001.

Sensitivity and specificity of diagnosis model.

Real Cancer Real Normal
Predicted Cancer 723 3
Predicted Normal 103 99
Total 826 102
Correct 723 99
Sensitivity 0.8753
Specificity 0.9706
Discussions

Because BRCA is often asymptomatic in early stages, most patients were diagnosed in later stages (5). In addition, previous studies have shown that the heterogeneity of BRCA is significantly correlated with its prognosis. Therefore, it is important to screen new prognostic and diagnostic biomarkers for different types of breast cancer to guide clinical practice (6, 7).

In the present study, we aimed to identify suitable ferroptosis- and immune-related lncRNAs as prognosis and diagnosis biomarkers for patients with breast infiltrating duct and lobular carcinoma by integrated analyses. Finally, we found that AC007686.3, AC078883.1, ADAMTS9-AS1, AL035661.1, CBR3-AS1, FTX, and TMEM105 were significantly associated with the OS of patients with breast infiltrating duct and lobular carcinoma. The risk assessment model and diagnosis model could predict the outcome of these patients. In addition, we found that the expression ADAMTS9-AS1 was decreased significantly, which demonstrated in previous studies as suppressor in BRCA cell invasion, proliferation, and the aggression (29, 30). In another research, Chen et al. also found that ADAMTS9-AS1, as prognostic biomarker, could promote cell proliferation and EMT in colorectal cancer (31). These results suggested that ADAMTS9-AS1 may play different roles for different types of cancers (29–31).

Fang et al. found that AL035661.1 was upregulated in renal cell carcinoma and positively associated with the poor prognosis (32). Similarly, we also found that the expression of AL035661.1 was decreased in patients with breast infiltrating duct and lobular carcinoma, and these patients with high expression of AL035661.1 exhibited worse OS. However, Lv et al. found that the expression of AL035661.1 was downregulated and could serve as prognosis biomarkers of hepatocellular carcinoma (33). These results indicated that AL035661.1 was closely related to the development and progression of several cancers (32, 33) and also increased the possibility of AL035661.1 as prognosis biomarker for BRCA patients with breast infiltrating duct and lobular carcinoma.

Previous studies have demonstrated that CBR3-AS1 could promote lung adenocarcinoma cell proliferation, migration, and invasion (34, 35). In BRCA, several studies have also demonstrated that the expression of CBR3-AS was upregulated in BRCA tissues and cells (36, 37). Those patients with high expression displayed poor OS (36, 37). Consistent with previous studies, we also found that CBR3-AS was increased and negatively significantly associated with the OS. Our present study reinforced the potential role of CBR3-AS1 as a therapeutic target and potential prognostic for breast cancer, especially for these BRCA patients with infiltrating duct and lobular carcinoma.

LncRNA FTX was first identified in Xist gene locus. Evidence indicated that FTX could promote the migration, proliferation, and invasion for several caners, including colorectal cancer, gastric cancer, and lung cancer (38–40). In this study, we found that the expression of FTX was downregulated. However, the patients with low expression of FTX exhibited better OS. Therefore, we hypothesized that the FTX expression was just correlated with the OS of BRCA patients. The aberrant expression of FTX could be a result of the progression of cancer.

For AC007686.3, AC078883.1, and TMEM105, although there were no reported relationships between them and cancer at present, our research suggested that they were closely related to BRCA cancer. Whether their abnormal expressions were the cause or the compensatory results remains to be further studied. In the future, we will also carry out functional studies in the further study.

 Conclusions

Target therapy for ferroptosis and immunity are new cancer treatment options discovered in recent years. In the present study, through a series of bioinformatics analyses, we found that seven FI-DELs could serve as prognostic and diagnostic biomarkers for patients with breast infiltrating duct and lobular carcinoma. However, whether these seven biomarkers could be really applied to the clinic requires further investigations.

 Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: https://portal.gdc.cancer.gov/.

 Author Contributions

TW and NZ conceived and designed the experiments. XX performed the analyses. WJ helped in the analysis of the data. XX wrote the paper. All authors contributed to the article and approved the submitted version.

 Funding

This project was financially supported by Foundation of Hunan University of Medicine (2020122004 and 20KJPY02).

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.844642/full#supplementary-material

 References Sung H Ferlay J Siegel RL Laversanne M Soerjomataram I Jemal A . Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71:209–49. doi: 10.3322/caac.21660 Desantis CE Ma J Gaudet MM Newman LA Miller KD Goding Sauer A . Breast Cancer Statistics 2019. CA Cancer J Clin (2019) 69:438–51. doi: 10.3322/caac.21583 Ponde NF Zardavas D Piccart M . Progress in Adjuvant Systemic Therapy for Breast Cancer. Nat Rev Clin Oncol (2019) 16:27–44. doi: 10.1038/s41571-018-0089-9 Ruddy KJ Ganz PA . Treatment of Nonmetastatic Breast Cancer. JAMA (2019) 321:1716–7. doi: 10.1001/jama.2019.3927 Wu ZH Tang Y Yu H Li HD . The Role of Ferroptosis in Breast Cancer Patients: A Comprehensive Analysis. Cell Death Discovery (2021) 7:93. doi: 10.1038/s41420-021-00473-5 Masood S . Prognostic/predictive Factors in Breast Cancer. Clin Lab Med (2005) 25:809–825, viii. doi: 10.1016/j.cll.2005.08.012 Yersal O Barutca S . Biological Subtypes of Breast Cancer: Prognostic and Therapeutic Implications. World J Clin Oncol (2014) 5:412–24. doi: 10.5306/wjco.v5.i3.412 Waks AG Winer EP . Breast Cancer Treatment: A Review. JAMA (2019) 321:288–300. doi: 10.1001/jama.2018.19323 Savas P Salgado R Denkert C Sotiriou C Darcy PK Smyth MJ . Clinical Relevance of Host Immunity in Breast Cancer: From TILs to the Clinic. Nat Rev Clin Oncol (2016) 13:228–41. doi: 10.1038/nrclinonc.2015.215 Dieci MV Radosevic-Robin N Fineberg S Van Den Eynden G Ternes N Penault-Llorca F . Update on Tumor-Infiltrating Lymphocytes (TILs) in Breast Cancer, Including Recommendations to Assess TILs in Residual Disease After Neoadjuvant Therapy and in Carcinoma in Situ: A Report of the International Immuno-Oncology Biomarker Working Group on Breast Cancer. Semin Cancer Biol (2018) 52:16–25. doi: 10.1016/j.semcancer.2017.10.003 Smyth MJ Dunn GP Schreiber RD . Cancer Immunosurveillance and Immunoediting: The Roles of Immunity in Suppressing Tumor Development and Shaping Tumor Immunogenicity. Adv Immunol (2006) 90:1–50. doi: 10.1016/S0065-2776(06)90001-7 Schreiber RD Old LJ Smyth MJ . Cancer Immunoediting: Integrating Immunity’s Roles in Cancer Suppression and Promotion. Science (2011) 331:1565–70. doi: 10.1126/science.1203486 Boudreau A Van’t Veer LJ Bissell MJ . An “Elite Hacker”: Breast Tumors Exploit the Normal Microenvironment Program to Instruct Their Progression and Biological Diversity. Cell Adh Migr (2012) 6:236–48. doi: 10.4161/cam.20880 Dixon SJ Lemberg KM Lamprecht MR Skouta R Zaitsev EM Gleason CE . Ferroptosis: An Iron-Dependent Form of Nonapoptotic Cell Death. Cell (2012) 149:1060–72. doi: 10.1016/j.cell.2012.03.042 Matsushita M Freigang S Schneider C Conrad M Bornkamm GW Kopf M . T Cell Lipid Peroxidation Induces Ferroptosis and Prevents Immunity to Infection. J Exp Med (2015) 212:555–68. doi: 10.1084/jem.20140857 Wang D Dubois RN . Immunosuppression Associated With Chronic Inflammation in the Tumor Microenvironment. Carcinogenesis (2015) 36:1085–93. doi: 10.1093/carcin/bgv123 Wang D Xie N Gao W Kang R Tang D . The Ferroptosis Inducer Erastin Promotes Proliferation and Differentiation in Human Peripheral Blood Mononuclear Cells. Biochem Biophys Res Commun (2018) 503:1689–95. doi: 10.1016/j.bbrc.2018.07.100 Stockwell BR Jiang X . A Physiological Function for Ferroptosis in Tumor Suppression by the Immune System. Cell Metab (2019) 30:14–5. doi: 10.1016/j.cmet.2019.06.012 Wang W Green M Choi JE Gijon M Kennedy PD Johnson JK . CD8(+) T Cells Regulate Tumour Ferroptosis During Cancer Immunotherapy. Nature (2019) 569:270–4. doi: 10.1038/s41586-019-1170-y Li Z Rong L . Cascade Reaction-Mediated Efficient Ferroptosis Synergizes With Immunomodulation for High-Performance Cancer Therapy. Biomater Sci (2020) 8:6272–85. doi: 10.1039/D0BM01168A Cai P Otten ABC Cheng B Ishii MA Zhang W Huang B . A Genome-Wide Long Noncoding RNA CRISPRi Screen Identifies PRANCR as a Novel Regulator of Epidermal Homeostasis. Genome Res (2020) 30:22–34. doi: 10.1101/gr.251561.119 Kim J Piao HL Kim BJ Yao F Han Z Wang Y . Long Noncoding RNA MALAT1 Suppresses Breast Cancer Metastasis. Nat Genet (2018) 50:1705–15. doi: 10.1038/s41588-018-0252-3 Zhao W Geng D Li S Chen Z Sun M . LncRNA HOTAIR Influences Cell Growth, Migration, Invasion, and Apoptosis via the miR-20a-5p/HMGA2 Axis in Breast Cancer. Cancer Med (2018) 7:842–55. doi: 10.1002/cam4.1353 Kong X Duan Y Sang Y Li Y Zhang H Liang Y . LncRNA-CDC6 Promotes Breast Cancer Progression and Function as ceRNA to Target CDC6 by Sponging microRNA-215. J Cell Physiol (2019) 234:9105–17. doi: 10.1002/jcp.27587 Chang KC Diermeier SD Yu AT Brine LD Russo S Bhatia S . MaTAR25 lncRNA Regulates the Tensin1 Gene to Impact Breast Cancer Progression. Nat Commun (2020) 11:6438. doi: 10.1038/s41467-020-20207-y Tang B Zhu J Li J Fan K Gao Y Cheng S . The Ferroptosis and Iron-Metabolism Signature Robustly Predicts Clinical Diagnosis, Prognosis and Immune Microenvironment for Hepatocellular Carcinoma. Cell Commun Signal (2020) 18:174. doi: 10.1186/s12964-020-00663-1 Xing XL Xing C Huang Z Yao ZY Liu YW . Immune-Related lncRNAs to Construct Novel Signatures and Predict the Prognosis of Rectal Cancer. Front Oncol (2021) 11:661846. doi: 10.3389/fonc.2021.661846 Huang Da W Sherman BT Lempicki RA . Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat Protoc (2009) 4:44–57. doi: 10.1038/nprot.2008.211 Fang S Zhao Y Hu X . LncRNA ADAMTS9-AS1 Restrains the Aggressive Traits of Breast Carcinoma Cells via Sponging miR-513a-5p. Cancer Manag Res (2020) 12:10693–703. doi: 10.2147/CMAR.S266575 Chen J Cheng L Zou W Wang R Wang X Chen Z . ADAMTS9-AS1 Constrains Breast Cancer Cell Invasion and Proliferation via Sequestering miR-301b-3p. Front Cell Dev Biol (2021) 9:719993. doi: 10.3389/fcell.2021.719993 Chen W Tu Q Yu L Xu Y Yu G Jia B . LncRNA ADAMTS9-AS1, as Prognostic Marker, Promotes Cell Proliferation and EMT in Colorectal Cancer. Hum Cell (2020) 33:1133–41. doi: 10.1007/s13577-020-00388-w Fang X Liu X Lu L Liu G . Identification of a Somatic Mutation-Derived Long Non-Coding RNA Signatures of Genomic Instability in Renal Cell Carcinoma. Front Oncol (2021) 11:728181. doi: 10.3389/fonc.2021.728181 Lv Y Wei W Huang Z Chen Z Fang Y Pan L . Long non-Coding RNA Expression Profile can Predict Early Recurrence in Hepatocellular Carcinoma After Curative Resection. Hepatol Res (2018) 48:1140–8. doi: 10.1111/hepr.13220 Guan Y Yang J Liu X Chu L . Long Noncoding RNA CBR3 Antisense RNA 1 Promotes the Aggressive Phenotypes of Nonsmallcell Lung Cancer by Sponging Microrna5093p and Competitively Upregulating HDAC9 Expression. Oncol Rep (2020) 44:1403–14. doi: 10.3892/or.2020.7719 Hou M Wu N Yao L . LncRNA CBR3-AS1 Potentiates Wnt/beta-Catenin Signaling to Regulate Lung Adenocarcinoma Cells Proliferation, Migration and Invasion. Cancer Cell Int (2021) 21:36. doi: 10.1186/s12935-020-01685-y Xu L Zhu H Gao F Tang Y Zhu Y Sun Z . Upregulation of the Long non-Coding RNA CBR3-AS1 Predicts Tumor Prognosis and Contributes to Breast Cancer Progression. Gene X (2019) 2:100014. doi: 10.1016/j.gene.2019.100014 Zhang M Wang Y Jiang L Song X Zheng A Gao H . LncRNA CBR3-AS1 Regulates of Breast Cancer Drug Sensitivity as a Competing Endogenous RNA Through the JNK1/MEK4-Mediated MAPK Signal Pathway. J Exp Clin Cancer Res (2021) 40:41. doi: 10.1186/s13046-021-01844-7 Li H Yao G Zhai J Hu D Fan Y . LncRNA FTX Promotes Proliferation and Invasion of Gastric Cancer via miR-144/ZFX Axis. Onco Targets Ther (2019) 12:11701–13. doi: 10.2147/OTT.S220998 Jin S He J Zhou Y Wu D Li J Gao W . LncRNA FTX Activates FOXA2 Expression to Inhibit non-Small-Cell Lung Cancer Proliferation and Metastasis. J Cell Mol Med (2020) 24:4839–49. doi: 10.1111/jcmm.15163 Chen GQ Liao ZM Liu J Li F Huang D Zhou YD . LncRNA FTX Promotes Colorectal Cancer Cells Migration and Invasion by miRNA-590-5p/RBPJ Axis. Biochem Genet (2021) 59:560–73. doi: 10.1007/s10528-020-10017-8