- 1Department of Surgical Oncology, Urumqi Friendship Hospital, Urumqi, China
- 2School of Public Health and Laboratory Medicine, Hunan Provincial Key Laboratory for Synthetic Biology of Traditional Chinese Medicine, Hunan University of Medicine, Huaihua, China
- 3Department of Breast Surgery, The Affiliated Tumor Hospital of Xinjiang Medical University, Urumqi, China
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.
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.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g001.jpg)
Figure 1 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.
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.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g002.jpg)
Figure 2 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).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g003.jpg)
Figure 3 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.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g004.jpg)
Figure 4 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).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g005.jpg)
Figure 5 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).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/844642/fonc-12-844642-HTML/image_m/fonc-12-844642-g006.jpg)
Figure 6 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.
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
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, 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:209–49. doi: 10.3322/caac.21660
2. Desantis CE, Ma J, Gaudet MM, Newman LA, Miller KD, Goding Sauer A, et al. Breast Cancer Statistics 2019. CA Cancer J Clin (2019) 69:438–51. doi: 10.3322/caac.21583
3. 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
4. Ruddy KJ, Ganz PA. Treatment of Nonmetastatic Breast Cancer. JAMA (2019) 321:1716–7. doi: 10.1001/jama.2019.3927
5. 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
6. Masood S. Prognostic/predictive Factors in Breast Cancer. Clin Lab Med (2005) 25:809–825, viii. doi: 10.1016/j.cll.2005.08.012
7. 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
8. Waks AG, Winer EP. Breast Cancer Treatment: A Review. JAMA (2019) 321:288–300. doi: 10.1001/jama.2018.19323
9. Savas P, Salgado R, Denkert C, Sotiriou C, Darcy PK, Smyth MJ, et al. 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
10. Dieci MV, Radosevic-Robin N, Fineberg S, Van Den Eynden G, Ternes N, Penault-Llorca F, et al. 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
11. 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
12. 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
13. 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
14. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: An Iron-Dependent Form of Nonapoptotic Cell Death. Cell (2012) 149:1060–72. doi: 10.1016/j.cell.2012.03.042
15. 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
16. Wang D, Dubois RN. Immunosuppression Associated With Chronic Inflammation in the Tumor Microenvironment. Carcinogenesis (2015) 36:1085–93. doi: 10.1093/carcin/bgv123
17. 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
18. 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
19. 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
20. 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
21. Cai P, Otten ABC, Cheng B, Ishii MA, Zhang W, Huang B, et al. 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
22. Kim J, Piao HL, Kim BJ, Yao F, Han Z, Wang Y, et al. Long Noncoding RNA MALAT1 Suppresses Breast Cancer Metastasis. Nat Genet (2018) 50:1705–15. doi: 10.1038/s41588-018-0252-3
23. 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
24. Kong X, Duan Y, Sang Y, Li Y, Zhang H, Liang Y, et al. 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
25. Chang KC, Diermeier SD, Yu AT, Brine LD, Russo S, Bhatia S, et al. MaTAR25 lncRNA Regulates the Tensin1 Gene to Impact Breast Cancer Progression. Nat Commun (2020) 11:6438. doi: 10.1038/s41467-020-20207-y
26. Tang B, Zhu J, Li J, Fan K, Gao Y, Cheng S, et al. 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
27. 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
28. 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
29. 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
30. 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
31. Chen W, Tu Q, Yu L, Xu Y, Yu G, Jia B, et al. 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
32. 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
33. Lv Y, Wei W, Huang Z, Chen Z, Fang Y, Pan L, et al. 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
34. 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
35. 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
36. Xu L, Zhu H, Gao F, Tang Y, Zhu Y, Sun Z, et al. 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
37. Zhang M, Wang Y, Jiang L, Song X, Zheng A, Gao H, et al. 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
38. 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
39. 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
Keywords: breast infiltrating duct and lobular carcinoma, ferroptosis, immune, lncRNAs, prognosis, diagnosis
Citation: Wei T, Zhu N, Jiang W and Xing X-L (2022) Development and Validation of Ferroptosis- and Immune-Related lncRNAs Signatures for Breast Infiltrating Duct and Lobular Carcinoma. Front. Oncol. 12:844642. doi: 10.3389/fonc.2022.844642
Received: 28 December 2021; Accepted: 09 March 2022;
Published: 04 April 2022.
Edited by:
Karen Elizabeth Nava-Castro, National Autonomous University of Mexico, MexicoReviewed by:
Shengchun Liu, First Affiliated Hospital of Chongqing Medical University, ChinaFeng Wang, Affiliated Hospital of Nantong University, China
Copyright © 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.
*Correspondence: Xiao-Liang Xing, xiaoliangxinghnm@126.com
†These authors have contributed equally to this work