Abstract
Background: A surge in newly diagnosed breast cancer has overwhelmed the public health system worldwide. Joint effort had beed made to discover the genetic mechanism of these disease globally. Accumulated research has revealed autophagy may act as a vital part in the pathogenesis of breast cancer.
Objective: Aim to construct a prognostic model based on autophagy-related lncRNAs and investigate their potential mechanisms in breast cancer.
Methods: The transcriptome data and clinical information of patients with breast cancer were obtained from The Cancer Genome Atlas (TCGA) database. Autophagy-related genes were obtained from the Human Autophagy Database (HADb). Long non-coding RNAs (lncRNAs) related to autophagy were acquired through the Pearson correlation analysis. Univariate Cox regression analysis as well as the least absolute shrinkage and selection operator (LASSO) regression analysis were used to identify autophagy-related lncRNAs with prognostic value. We constructed a risk scoring model to assess the prognostic significance of the autophagy-related lncRNAs signatures. The nomogram was then established based on the risk score and clinical indicators. Through the calibration curve, the concordance index (C-index) and receiver operating characteristic (ROC) curve analysis were evaluated to obtain the model's predictive performance. Subgroup analysis was performed to evaluate the differential ability of the model. Subsequently, gene set enrichment analysis was conducted to investigate the potential functions of these lncRNAs.
Results: We attained 1,164 breast cancer samples from the TCGA database and 231 autophagy-related genes from the HAD database. Through correlation analysis, 179 autophagy-related lncRNAs were finally identified. Univariate Cox regression analysis and LASSO regression analysis further screened 18 prognosis-associated lncRNAs. The risk scoring model was constructed to divide patients into high-risk and low-risk groups. It was found that the low-risk group had better overall survival (OS) than those of the high-risk group. Then, the nomogram model including age, tumor stage, TNM stage and risk score was established. The evaluation index (C-index: 0.78, 3-year OS AUC: 0.813 and 5-year OS AUC: 0.785) showed that the nomogram had excellent predictive power. Subgroup analysis showed there were difference in OS between high-risk and low-risk patients in different subgroups (stage I-II, ER positive, Her-2 negative and non-TNBC subgroups; all P < 0.05). According to the results of gene set enrichment analysis, these lncRNAs were involved in the regulation of multicellular organismal macromolecule metabolic process in multicellular organisms, nucleotide excision repair, oxidative phosphorylation, and TGF-β signaling pathway.
Conclusions: We identified 18 autophagy-related lncRNAs with prognostic value in breast cancer, which may regulate tumor growth and progression in multiple ways.
Introduction
Breast cancer is the most common malignant tumor and the second leading cause of cancer death in females. In 2018, there were ~2.08 million new cases and 630,000 deaths worldwide (Bray et al., ). Accurately predicting the prognosis of breast cancer is essential for improving the prognosis and providing appropriate treatment for patients. With the advancement of high-throughput technologies, various prognostic models have been developed according to the sequencing data. Moreover, prognostic models based on the autophagy-related genes, immune-related genes, or immune-related lncRNAs have been studied before (Li B. et al., ; Lin et al., ; Shen et al., ). But none of them combined clinical data with sequencing data to build a predictive model for breast cancer patients. Our research aims to establish a prognostic model for patients with breast cancer utilizing the autophagy-related lncRNAs and clinical data, which may serve as an important supplement to prognosis prediction.
Long non-coding RNA (lncRNA) is a type of RNA with more than 200 nucleotides, which functions by regulating gene expression and interacting with specific proteins (Esteller, ; Kopp and Mendell, ; Choudhari et al., ). Substantial evidence showed that lncRNAs were involved in a series of biological processes, including proliferation and differentiation of tumor cells, chromosome remodeling, transcription, as well as post-transcriptional modification (Wang and Chang, ; Lau, ; Tomar et al., ; Wu et al., ). The role of lncRNA in breast cancer has been revealed gradually (Wang et al., ; Richard and Eichhorn, ). Wang et al. () found that lncRNA-encoded polypeptide ASRPS could inhibit breast cancer angiogenesis. Dong et al. () demonstrated that lncRNA TINCR promoted chemotherapy resistance and epithelial-mesenchymal transition of breast cancer by targeting microRNA-125b. Previous studies suggested that lncRNA has a potential impact on the prognosis and targeted therapy of breast cancer (Jiang et al., ; Yang et al., ; Du et al., ).
Autophagy is a cellular process in which damaged organelles and unnecessary macromolecules were degraded and recycled under the control of autophagy-related genes (Li X. et al., ), which played a unique role in different stages of tumor pathogenesis. Previous studies disclosed that decreased expression of autophagy-related gene Beclin-1 could inhibit autophagy of tumor cells and further accelerate breast cancer progression (Lisiak et al., ; Hamurcu et al., ). Wei et al. () indicated that Magnoflorine could induce autophagy in breast cancer cells through the AKT/mTOR signaling pathway, and thus increasing their sensitivity to doxorubicin. Lin et al. () proved that the prognostic model of breast cancer based on autophagy-related genes has satisfactory predictive power. Taken together, cell autophagy plays an important role in the progression, treatment, and prognosis of breast cancer.
Autophagy-related lncRNAs may have a potential effect on the prognosis of breast cancer patients. It has been reported that lncRNA H19 induces the tamoxifen-resistance in breast cancer through the activation of the H19/SAHH/DNMT3B axis, thereby facilitating cell autophagy (Wang et al., ). Both in vitro and in vivo experiments indicated that lncRNA NAMPT-AS regulated tumor cell autophagy and promoted cancer progression through the mTOR signaling pathway (Zhang et al., ). Luan et al. () screened 10 autophagy-related lncRNAs in glioma and proved the prognostic value of these lncRNAs. This paper aims to construct a prognostic model based on autophagy-related lncRNAs and investigate their potential mechanisms in breast cancer.
Materials and Methods
Sample Resources and Processing
The technical roadmap for the development of an autophagy-related lncRNA breast cancer prognosis model is shown in Figure 1. We downloaded the transcriptome data including mRNA and lncRNA as well as corresponding clinical information of breast cancer from the Cancer Genome Atlas database (TCGA, https://cancergenome.nih.gov/), and autophagy-related genes from the human autophagy database (HADb, http://www.autophagy.lu/). We used | log2FC | > 0.5 and false discovery rate (FDR) < 0.05 as the thresholds to identify differentially expressed autophagy-related genes and lncRNAs between cancer samples and normal samples. Pearson correlation analysis was utilized to calculate the correlation coefficient (R2) between lncRNAs and autophagy-related genes. Also, lncRNAs with | R2 | > 0.4 and p < 0.05 were defined as autophagy-related lncRNAs.
Figure 1
Construction of Risk Scoring Model
Univariate Cox regression analysis was performed on autophagy-related lncRNAs, and prognosis-related lncRNAs were recognized with criteria of p < 0.05. We performed the LASSO regression analysis to identify the crucial lncRNAs that were closely associated with OS. The selected lncRNAs were included in the multivariate Cox regression model to generate their coefficients (coefi). A risk scoring model containing both coefi and lncRNA expression levels (expi) was constructed to get the risk score for all samples. Patients were divided into high-risk group and low-risk group with the median risk score as the cut-off point. What's more, survival curve was executed to compare the prognosis of the two groups. Univariate and multivariate Cox regression analysis were exploited to evaluate the prognostic significance of risk score and clinical factors including age, tumor stage, and TNM stage. The area under the curve (AUC) was generated by ROC curve analysis to assess the predictive effectiveness of the model. Then we further evaluated the model by stratifying the patients into different subgroups by clinical stage, ER status, HER-2 status and TNBC.
Establishment and Evaluation of Prognostic Model
We established a nomogram based on the multivariate Cox regression results of age, tumor stage, TMN stage, and risk score. The 3- and 5-year OS of each patient was predicted based on the nomogram. Meanwhile, the concordance index (C-index), calibration curve, and ROC curve were generated to evaluate the predicted efficacy of the model.
Gene Set Enrichment Analysis
Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied to differentially expressed genes between high-risk and low-risk groups. We also analyzed the gene set to deduce their functions and determined whether the gene set differed significantly between the two groups. The research was to investigate whether the differentially expressed genes between the two groups were enriched during autophagy.
Results
Differentially Expressed Autophagy-Related Genes and lncRNAs
We extracted breast cancer samples from the TCGA database, including 1,053 tumor tissues and 111 normal tissues, and 231 autophagy-related genes were generated from the HADb database. Differentially expressed autophagy-related genes and lncRNAs were identified between breast cancer and para-cancer tissues with the help of the “limma” package in the R software. Seventy-five differentially expressed autophagy-related genes were found, of which 32 genes were down-regulated and 43 genes were up-regulated (Figure 2A). Besides, we recognized 3,770 differentially expressed lncRNAs with 2,068 lncRNAs down-regulated and 1,702 lncRNAs up-regulated (Figure 2B). Then, we attained 179 autophagy-related lncRNAs through Pearson correlation analysis with criteria of | R2 | > 0.4 and p < 0.05 (Supplementary Table 1).
Figure 2
Construction of Risk Scoring Model
Patients with follow-up time exceeding 30 days were included in univariate Cox regression analysis to explore the prognostic significance of autophagy-related lncRNAs. Twenty prognosis-associated lncRNAs were acquired with p < 0.05 (Figure 3A). We further screened 18 core lncRNAs tightly related to prognosis by LASSO regression analysis (Figures 3B,C). The multivariate Cox regression was employed to construct the risk scoring model as follows: Risk score . Patients were divided into high-risk group and low-risk group with the median risk score as the cut-off point. Survival curve result showed that the OS in the low-risk group was significantly higher than that in the high-risk group (p < 0.0001, Figure 4A). The risk curve and scatterplot showed that the low-risk group had lower risk coefficient and mortality rate compared to the high-risk group (Figures 4B,C). The heat map visualized the expression levels of the 18-lncRNAs between the high-risk and low-risk groups (Figure 4D). Univariate and multivariate Cox regression analysis results showed that 18-lncRNAs signature was an independent prognostic factor of patients with breast cancer (Figures 5A,B). The AUCs of 3- and 5-year OS generated by ROC curve analysis were 0.724 and 0.685, respectively (Figure 5C). In subgroup analysis, the risk scoring model showed satisfying differential ability in stage I-II, ER positive, Her-2 negative, and non-TNBC subgroups (P < 0.05). Due to small sample size of stage III-IV, ER negative, Her-2 positive and TNBC subgroups, the survival curve showed potential identification trends, but there were no statistical difference (Figure 6).
Figure 3
Figure 4
Figure 5
Figure 6
Establishment and Evaluation of the Prognostic Model
The nomogram was established according to the results of multivariate Cox regression including age, tumor stage, TMN stage, and risk score (Figure 7A). The calibration curve was shown in Figure 7. The C-index of the nomogram was 0.78. And the AUCs of the 3- and 5-year OS were 0.813 and 0.785, respectively (Figure 7B). Our results showed that the prognostic model presented excellent accuracy on prediction (Figures 7C,D).
Figure 7
Gene Set Enrichment Analysis
GO enrichment analysis and KEGG pathway analysis were conducted on differentially expressed genes between high-risk and low-risk groups. GO enrichment analysis indicated that the genes were enriched in angiogenesis, base excision repair, extracellular structure organization, assembly of mitochondrial respiratory chain complex, multicellular biological macromolecule metabolic process, negative regulation of cell-substrate adhesion, nucleotide excision repair, and positive regulation of protein acetylation (Figure 8A). KEGG pathway analysis denoted that these genes were involved in the modulation of cell adhesion molecules cams, ECM receptor interaction, focal adhesion, nucleotide excision repair, oxidative phosphorylation, ribosome, spliceosome, and TGF-β signaling pathway. All these contents may assist researchers to further explore the mechanisms of autophagy-related lncRNAs influencing breast cancer pathogenesis (Figure 8B).
Figure 8
Discussion
The incidence of newly diagnosed breast cancer is increasing, which severely affect women's life. Zhou et al. (2017) demonstrated that genetic variation of autophagy-related genes was associated with the prognosis of breast cancer patients. Autophagy can maintain the homeostasis and the survival of breast cancer cells by removing dysfunctional or unnecessary substances (Han et al., ; Ostendorf and Tavazoie, ). In vitro experiments showed that inhibiting autophagy could hinder the survival and metastasis of dormant breast cancer cells (Vera-Ramirez et al., ). Meanwhile, previous researchers have suggested that autophagy was tightly related to breast cancer prognosis (Gu et al., ; Lin et al., ; Zhong et al., ). Hence, finding reliable autophagy-related prognostic indicators may drive for early risk stratification and proper therapeutic decision for breast cancer.
The process of autophagy can be regulated by lncRNAs in a variety of ways. Particularly, lncRNAs could adsorb autophagy-related miRNA as a molecular sponge to regulate the expression of autophagy-related genes. It was found that LncRNA H19 could induce autophagy and suppress epithelial-mesenchymal transformation in breast cancer by targeting the let-7 miRNA family (Xiong et al., ). LncRNA HOTAIR may affect the autophagy of breast cancer cells by sponging miR-34a (Pawlowska et al., ). Luen et al. showed that the risk scoring model constructed by autophagy-related lncRNAs could effectively predict the prognosis of glioma patients (Luan et al., ). However, the prognostic role of autophagy-related lncRNAs in breast cancer has not been studied before.
In this study, we obtained breast cancer RNA-sequence data and clinical information from the TCGA database and autophagy-related genes from the HADb database. Univariate Cox analysis and LASSO regression analysis identified 18 autophagy-related lncRNAs with prognostic significance. A risk scoring model was constructed based on these lncRNAs to divide patients into high-risk and low-risk groups. Survival analysis revealed the low-risk sub-group had much more favorable prognosis. Sequentially, we established a nomogram based on the risk scores and clinical variables including age, tumor stage, and TMN stage to predict the 3- and 5-year OS of patients. The results of the C-index, calibration curve, and ROC curve validated the predictive performance of the model. The result of subgroup analysis implicated the model may have better differential ability in hormone positive and Her-2 negative breast cancer. Further studies are required to validate the result above.
Among the 18 vital lncRNAs, 9 lncRNAs (LIPE.AS1, AC061992.1, TFAP2A.AS1, LINC00667, AL512625.2, PDCD4.AS1, AC002398.1, AL451085 2, and AC110619.1) were protective factors for the prognosis of breast cancer, while other lncRNAs (LINC02169, LINC01235, LINC01614, AC109826.1, ST8SIA6.AS1, AP003119.3, MIR4435.2HG, AC147067.2, and AC098484.1) behaved the opposite (Supplementary Table 2). Previous studies have shown that AC061992.1 and LINC00667 are protective prognostic factors for breast cancer patients (Fan et al., ; Zhu et al., 2020). Zhou et al. (2019) reported that lncRNA TFAP2A.AS1 could inhibit the proliferation and invasion of breast cancer cells by adsorbing miR-933. Jadaliha et al. () found that lncRNA PDCD4.AS1 was capable of regulating the post-transcriptional process of oncogenes or tumor suppressor genes, thereby accelerating the progress of breast cancer. The above studies suggested that the screened lncRNAs played an important role in the progress of breast cancer.
The correlation frequency between autophagy-related genes and 18 prognosis-associated lncRNAs was shown in Supplementary Figure 1. The BAX expression in breast cancer cell line MCF-7 was upregulated when the apoptosis increasing significantly, caused by the transfection of HCCR-1 siRNA (Meng et al., ). Through real-time PCR, Shen et al. () revealed that ATG16L2, an autophagy-related gene, was closely related to cisplatin-induced apoptosis of breast cancer cell. Numerous studies have revealed the signaling pathway of autophagy and programmed cell death (Nagakannan et al., ; Chang et al., ). However, few studies could perfectly explain the role of BAX and ATG16L2 in breast cancer treatment resistance. ITGB1 was associated with EMT of breast cancer and the ITGB1/FAK/Src/AKT/β-catenin/MMP-9 signaling axis was an important signaling pathway (Du et al., ). ERBB2 was a star gene. Twenty-five percent of breast cancer over-expressed ERBB2, which conferred a more aggressive phenotype and targeted therapy was needed (D'Amato et al., ). Carolina et al. found the chromodomain helicase CHD4 regulated ERBB2 signaling pathway and autophagy, which represented a mechanism of resistance against Trastuzumab, a therapeutic anti-ERBB2 antibody (D'Alesio et al., ). MLST8 that controlled the rates of cell growth and proliferation was regulated by the mTOR pathway (McMahon et al., ). However, the function of MLST8 and SPNS1 in breast cancer was still unknown. In brief, the intrinsic mechanism of these genes in recurrence and metastasis of breast cancer needs to be clarified in our next step of research.
To further reveal the functions of 18-lncRNAs signature in breast cancer, we executed GO enrichment analysis and KEGG pathway analysis on the differentially expressed genes between high-risk and low-risk groups. We found that these lncRNAs were involved in regulating angiogenesis, repairing nucleotide excision, oxidative phosphorylation, TGF-β signaling pathway, and assembly of the mitochondrial respiratory chain. It is worth noting that the deletion of autophagy-related gene RB1CC1 could reduce the mitochondrial mass and oxidative respiration capacity of breast tumor cells, as well as inhibit the phosphorylation of mTOR substrate (Yeo et al., ). Breast cancer stem cells were reported to regulate the occurrence of autophagy through the TGF-β pathway (Yeo and Guan, ). Autophagy inducers could suppress breast cancer angiogenesis and thus inhibit the progression of breast tumors (Li et al., ). In summary, our findings indicated that the 18-lncRNAs signature was essential for autophagy modulation and tumor pathogenesis. However, the specific mechanism of 18 autophagy-related lncRNAs in breast cancer still needs further elucidation in vitro experiments.
There are still some limitations to the study. As it is a retrospective study derived from public data, it lacks some information such as further subdivision of ER subtypes and TNBC subtypes. Though the signature provides fine predictive effect on the prognosis of hormone receptor positive subgroup, it can not further verify its predictive performance in ER alpha and ER beta subtypes. As it is a retrospective study derived from public data, it lacks some information such as recurrence time and treatment records. Also, a relatively small size cohort of TNBC patients limits the predictive performance of the model in some subgroups analysis. Further in vivo or in vitro assays and prospective clinical trials are needed to verify our conclusions.
Conclusion
In summary, the present study constructed an 18-autophagy-related lncRNAs signature to predict the prognostic value for breast cancer patients. The prognostic model based on these lncRNAs has proved to be effective in predicting the OS. Moreover, gene set enrichment analysis revealed that the 18-lncRNAs signature might regulate the progress of breast cancer through a series of biological pathways.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
XL conceived the study, performed the data analysis, prepared the figures and tables, authored and reviewed drafts of the paper, and approved the final draft. JC collected the data from database, prepared the figures, reviewed drafts of paper and approved the final draft. JZ and WLo conceived the study, prepared the figures, reviewed drafts of the paper, and approved the final draft. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Guangdong Medical Research Foundation (Grant No. A2020622) and Elite Young Scholars Program of Jiangmen Central Hospital (Grant No. J201905).
Acknowledgments
We thank the The Cancer Genome Atlas (TCGA) Database for sharing the large amount of data.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.569318/full#supplementary-material
References
1
BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. 10.3322/caac.21492
2
ChangY. F.ZhangD.HuW. M.LiuD. X.LiL. (2020). Semaglutide-mediated protection against Abeta correlated with enhancement of autophagy and inhibition of apotosis. J. Clin. Neurosci. 81, 234–239. 10.1016/j.jocn.2020.09.054
3
ChoudhariR.SedanoM. J.HarrisonA. L.SubramaniR.LinK. Y.RamosE. I.et al. (2020). Long noncoding RNAs in cancer: from discovery to therapeutic targets. Adv. Clin. Chem.95, 105–147. 10.1016/bs.acc.2019.08.003
4
D'AlesioC.BelleseG.GaglianiM. C.LechiaraA.DameriM.GrasselliE.et al. (2019). The chromodomain helicase CHD4 regulates ERBB2 signaling pathway and autophagy in ERBB2(+) breast cancer cells. Biol Open8:4–8. 10.1242/bio.038323
5
D'AmatoV.RaimondoL.FormisanoL.GiulianoM.De PlacidoS.RosaR.et al. (2015). Mechanisms of lapatinib resistance in HER2-driven breast cancer. Cancer Treat. Rev. 41, 877–883. 10.1016/j.ctrv.2015.08.001
6
DongH.HuJ.ZouK.YeM.ChenY.WuC.et al. (2019). Activation of LncRNA TINCR by H3K27 acetylation promotes Trastuzumab resistance and epithelial-mesenchymal transition by targeting MicroRNA-125b in breast Cancer. Mol. Cancer18:3. 10.1186/s12943-018-0931-9
7
DuZ.GaoW.SunJ.LiY.SunY.ChenT.et al. (2019). Identification of long noncoding RNAmediated transcriptional dysregulation triplets reveals global patterns and prognostic biomarkers for ER+/PR+, HER2 and triple negative breast cancer. Int. J. Mol. Med. 44, 1015–1025. 10.3892/ijmm.2019.4261
8
EstellerM. (2011). Non-coding RNAs in human disease. Nat. Rev. Genet.12, 861–874. 10.1038/nrg3074
9
FanC. N.MaL.LiuN. (2018). Systematic analysis of lncRNA-miRNA-mRNA competing endogenous RNA network identifies four-lncRNA signature as a prognostic biomarker for breast cancer. J. Transl. Med. 16:264. 10.1186/s12967-018-1640-2
10
GuY.LiP.PengF.ZhangM.ZhangY.LiangH.et al. (2016). Autophagy-related prognostic signature for breast cancer. Mol. Carcinog.55, 292–299. 10.1002/mc.22278
11
HamurcuZ.DelibasiN.NalbantogluU.SenerE. F.NurdinovN.TasciB.et al. (2019). FOXM1 plays a role in autophagy by transcriptionally regulating Beclin-1 and LC3. genes in human triple-negative breast cancer cells. J. Mol. Med. 97, 491–508. 10.1007/s00109-019-01750-8
12
HanY.FanS.QinT.YangJ.SunY.LuY.et al. (2018). Role of autophagy in breast cancer and breast cancer stem cells (Review). Int. J. Oncol.52, 1057–1070. 10.3892/ijo.2018.4270
13
JadalihaM.GholamalamdariO.TangW.ZhangY.PetracoviciA.HaoQ.et al. (2018). A natural antisense lncRNA controls breast cancer progression by promoting tumor suppressor gene mRNA stability. PLoS Genet. 14:e1007802. 10.1371/journal.pgen.1007802
14
JiangY. Z.LiuY. R.XuX. E.JinX.HuX.YuK. D.et al. (2016). Transcriptome analysis of triple-negative breast cancer reveals an integrated mRNA-lncRNA signature with predictive and prognostic value. Cancer Res.76, 2105–2114. 10.1158/0008-5472.CAN-15-3284
15
KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs. Cell172, 393–407. 10.1016/j.cell.2018.01.011
16
LauE. (2014). Non-coding RNA: zooming in on lncRNA functions. Nat. Rev. Genet.15, 574–575. 10.1038/nrg3795
17
LiB.GengR.WuQ.YangQ.SunS.ZhuS.et al. (2020). Alterations in immune-related genes as potential marker of prognosis in breast cancer. Front. Oncol.10:333. 10.3389/fonc.2020.00333
18
LiP.GuoY.BledsoeG.YangZ.ChaoL.ChaoJ. (2016). Kallistatin induces breast cancer cell apoptosis and autophagy by modulating Wnt signaling and microRNA synthesis. Exp. Cell Res. 340, 305–314. 10.1016/j.yexcr.2016.01.004
19
LiX.HeS.MaB. (2020). Autophagy and autophagy-related proteins in cancer. Mol. Cancer19:12. 10.1186/s12943-020-1138-4
20
LinQ. G.LiuW.MoY. Z.HanJ.GuoZ. X.ZhengW.et al. (2020). Development of prognostic index based on autophagy-related genes analysis in breast cancer. Aging12, 1366–1376. 10.18632/aging.102687
21
LisiakN.TotonE.RybczynskaM. (2018). Autophagy as a potential therapeutic target in breast cancer treatment. Curr. Cancer Drug Targets18, 629–639. 10.2174/1568009617666171114143330
22
LuanF.ChenW.ChenM.YanJ.ChenH.YuH.et al. (2019). An autophagy-related long non-coding RNA signature for glioma. FEBS Open Bio.9, 653–667. 10.1002/2211-5463.12601
23
McMahonL. P.YueW.SantenR. J.LawrenceJ. J. (2005). Farnesylthiosalicylic acid inhibits mammalian target of rapamycin (mTOR) activity both in cells and in vitro by promoting dissociation of the mTOR-raptor complex. Mol. Endocrinol. 19, 175–183. 10.1210/me.2004-0305
24
MengK.YuanG.BaoH.WangL.MaR.YuB.et al. (2019). Interaction of HCCR-1 and Bax in breast cancer. J. BUON24, 1027–1037
25
NagakannanP.IqbalM. A.YeungA.ThliverisJ. A.RastegarM.GhavamiS.et al. (2016). Perturbation of redox balance after thioredoxin reductase deficiency interrupts autophagy-lysosomal degradation pathway and enhances cell death in nutritionally stressed SH-SY5Y cells. Free Radic. Biol. Med. 101, 53–70. 10.1016/j.freeradbiomed.2016.09.026
26
OstendorfB. N.TavazoieS. F. (2020). Autophagy suppresses breast cancer metastasis. Dev. Cell52, 542–544. 10.1016/j.devcel.2020.02.005
27
PawlowskaE.SzczepanskaJ.BlasiakJ. (2017). The long noncoding RNA HOTAIR in breast cancer: does autophagy play a role?Int. J. Mol. Sci. 18:11–18. 10.3390/ijms18112317
28
RichardJ.EichhornP. (2018). Deciphering the roles of lncRNAs in breast development and disease. Oncotarget9, 20179–20212. 10.18632/oncotarget.24591
29
ShenM.DuanW. M.WuM. Y.WangW. J.LiuL.XuM. D.et al. (2015). Participation of autophagy in the cytotoxicity against breast cancer cells by cisplatin. Oncol. Rep. 34, 359–367. 10.3892/or.2015.4005
30
ShenY.PengX.ShenC. (2020). Identification and validation of immune-related lncRNA prognostic signature for breast cancer. Genomics112, 2640–2646. 10.1016/j.ygeno.2020.02.015
31
TomarD.YadavA. S.KumarD.BhadauriyaG.KunduG. C. (2020). Non-coding RNAs as potential therapeutic targets in breast cancer. Biochim. Biophys. Acta Gene Regul. Mech.1863:194378. 10.1016/j.bbagrm.2019.04.005
32
Vera-RamirezL.VodnalaS. K.NiniR.HunterK. W.GreenJ. E. (2018). Autophagy promotes the survival of dormant breast cancer cells and metastatic tumour recurrence. Nat. Commun. 9:1944. 10.1038/s41467-018-04070-6
33
WangJ.XieS.YangJ.XiongH.JiaY.ZhouY.et al. (2019). The long noncoding RNA H19 promotes tamoxifen resistance in breast cancer via autophagy. J. Hematol. Oncol. 12:81. 10.1186/s13045-019-0747-0
34
WangJ.YeC.XiongH.ShenY.LuY.ZhouJ.et al. (2017). Dysregulation of long non-coding RNA in breast cancer: an overview of mechanism and clinical implication. Oncotarget8, 5508–5522. 10.18632/oncotarget.12537
35
WangK. C.ChangH. Y. (2011). Molecular mechanisms of long noncoding RNAs. Mol. Cell43, 904–914. 10.1016/j.molcel.2011.08.018
36
WangY.WuS.ZhuX.ZhangL.DengJ.LiF.et al. (2020). LncRNA-encoded polypeptide ASRPS inhibits triple-negative breast cancer angiogenesis. J. Exp. Med. 217:3–9. 10.1084/jem.20190950
37
WeiT.XiaojunX.PeilongC. (2020). Magnoflorine improves sensitivity to doxorubicin (DOX) of breast cancer cells via inducing apoptosis and autophagy through AKT/mTOR and p38 signaling pathways. Biomed. Pharmacother.121:109139. 10.1016/j.biopha.2019.109139
38
WuP.MoY.PengM.TangT.ZhongY.DengX.et al. (2020). Emerging role of tumor-related functional peptides encoded by lncRNA and circRNA. Mol. Cancer19:22. 10.1186/s12943-020-1147-3
39
XiongH.ShenJ.ChenZ.YangJ.XieB.JiaY.et al. (2020). H19/let7/Lin28 ceRNA network mediates autophagy inhibiting epithelialmesenchymal transition in breast cancer. Int. J. Oncol. 56, 794–806. 10.3892/ijo.2020.4967
40
YangR.XingL.WangM.ChiH.ZhangL.ChenJ. (2018). Comprehensive analysis of differentially expressed profiles of lncRNAs/mRNAs and miRNAs with associated ceRNA Networks in triple-negative breast cancer. Cell. Physiol. Biochem.50, 473–488. 10.1159/000494162
41
YeoS. K.GuanJ. L. (2016). Hierarchical heterogeneity in mammary tumors and its regulation by autophagy. Autophagy12, 1960–1961. 10.1080/15548627.2016.1213928
42
YeoS. K.PaulR.HaasM.WangC.GuanJ. L. (2018). Improved efficacy of mitochondrial disrupting agents upon inhibition of autophagy in a mouse model of BRCA1-deficient breast cancer. Autophagy14, 1214–1225. 10.1080/15548627.2018.1460010
43
ZhangH.ZhangN.LiuY.SuP.LiangY.LiY.et al. (2019). Epigenetic regulation of NAMPT by NAMPT-AS drives metastatic progression in triple-negative breast cancer. Cancer Res.79, 3347–3359. 10.1158/0008-5472.CAN-18-3418
44
ZhongS.ChenH.YangS.FengJ.ZhouS. (2020). Identification and validation of prognostic signature for breast cancer based on genes potentially involved in autophagy. PeerJ8:e9621. 10.7717/peerj.9621
45
ZhouB.GuoH.TangJ. (2019). Long Non-Coding RNA TFAP2A-AS1 Inhibits Cell Proliferation and Invasion in Breast Cancer via miR-933/SMAD2. Med. Sci. Monit. 25, 1242–1253. 10.12659/MSM.912421
46
ZhouJ.HangD.JiangY.ChenJ.HanJ.ZhouW.et al. (2017). Evaluation of genetic variants in autophagy pathway genes as prognostic biomarkers for breast cancer. Gene627, 549–555. 10.1016/j.gene.2017.06.053
47
ZhuM.LvQ.HuangH.SunC.PangD.WuJ. (2020). Identification of a four-long non-coding RNA signature in predicting breast cancer survival. Oncol. Lett.19, 221–228. 10.3892/ol.2019.11063
Summary
Keywords
breast cancer, autophagy, long non-coding RNA, nomogram, gene set enrichment analysis
Citation
Li X, Chen J, Yu Q, Huang H, Liu Z, Wang C, He Y, Zhang X, Li W, Li C, Zhao J and Long W (2021) A Signature of Autophagy-Related Long Non-coding RNA to Predict the Prognosis of Breast Cancer. Front. Genet. 12:569318. doi: 10.3389/fgene.2021.569318
Received
03 June 2020
Accepted
16 February 2021
Published
16 March 2021
Volume
12 - 2021
Edited by
Seyed Javad Mowla, Tarbiat Modares University, Iran
Reviewed by
Pieter Johan Adam Eichhorn, National University of Singapore, Singapore; Patricia Schoenlein, Augusta University, United States
Updates
Copyright
© 2021 Li, Chen, Yu, Huang, Liu, Wang, He, Zhang, Li, Li, Zhao and Long.
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: Wansheng Long jmlws2@163.comJinglin Zhao 403612238@qq.com
This article was submitted to RNA, a section of the journal Frontiers in Genetics
†These authors have contributed equally to this work
Disclaimer
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.