AUTHOR=Chen Jiafeng , Li Xinrong , Yan Shuixin , Li Jiadi , Zhou Yuxin , Wu Minhua , Ding Jinhua , Yang Jiahui , Yuan Yijie , Zhu Ye , Wu Weizhu TITLE=An autophagy-related long non-coding RNA prognostic model and related immune research for female breast cancer JOURNAL=Frontiers in Oncology VOLUME=12 YEAR=2022 URL=https://www.frontiersin.org/journals/oncology/articles/10.3389/fonc.2022.929240 DOI=10.3389/fonc.2022.929240 ISSN=2234-943X ABSTRACT=Introduction

Breast cancer (BRCA) is the most common malignancy among women worldwide. It was widely accepted that autophagy and the tumor immune microenvironment play an important role in the biological process of BRCA. Long non-coding RNAs (lncRNAs), as vital regulatory molecules, are involved in the occurrence and development of BRCA. The aim of this study was to assess the prognosis of BRCA by constructing an autophagy-related lncRNA (ARlncRNA) prognostic model and to provide individualized guidance for the treatment of BRCA.

Methods

The clinical data and transcriptome data of patients with BRCA were acquired from the Cancer Genome Atlas database (TCGA), and autophagy-related genes were obtained from the human autophagy database (HADb). ARlncRNAs were identified by conducting co‑expression analysis. Univariate and multivariate Cox regression analysis were performed to construct an ARlncRNA prognostic model. The prognostic model was evaluated by Kaplan–Meier survival analysis, plotting risk curve, Independent prognostic analysis, clinical correlation analysis and plotting ROC curves. Finally, the tumor immune microenvironment of the prognostic model was studied.

Results

10 ARlncRNAs(AC090912.1, LINC01871, AL358472.3, AL122010.1, SEMA3B-AS1, BAIAP2-DT, MAPT-AS1, DNAH10OS, AC015819.1, AC090198.1) were included in the model. Kaplan–Meier survival analysis of the prognostic model showed that the overall survival(OS) of the low-risk group was significantly better than that of the high-risk group (p< 0.001). Multivariate Cox regression analyses suggested that the prognostic model was an independent prognostic factor for BRCA (HR = 1.788, CI = 1.534–2.084, p < 0.001). ROCs of 1-, 3- and 5-year survival revealed that the AUC values of the prognostic model were all > 0.7, with values of 0.779, 0.746, and 0.731, respectively. In addition, Gene Set Enrichment Analysis (GSEA) suggested that several tumor-related pathways were enriched in the high-risk group, while several immune‑related pathways were enriched in the low-risk group. Patients in the low-risk group had higher immune scores and their immune cells and immune pathways were more active. Patients in the low-risk group had higher PD-1 and CTLA-4 levels and received more benefits from immune checkpoint inhibitors (ICIs) therapy.

Discussion

The ARlncRNA prognostic model showed good performance in predicting the prognosis of patients with BRCA and is of great significance to guide the individualized treatment of these patients.