- Open Access
Molecular characterization of breast cancer: a potential novel immune-related lncRNAs signature
Journal of Translational Medicine volume 18, Article number: 416 (2020)
Accumulating evidence has demonstrated that immune-related lncRNAs (IRLs) are commonly aberrantly expressed in breast cancer (BC). Thus, we aimed to establish an IRL-based tool to improve prognosis prediction in BC patients.
We obtained IRL expression profiles in large BC cohorts (N = 911) from The Cancer Genome Atlas (TCGA) database. Then, in light of the correlation between each IRL and recurrence-free survival (RFS), we screened prognostic IRL signatures to construct a novel RFS nomogram via a Cox regression model. Subsequently, the performance of the IRL-based model was evaluated through discrimination, calibration ability, risk stratification ability and decision curve analysis (DCA).
A total of 52 IRLs were obtained from TCGA. Based on multivariate Cox regression analyses, four IRLs (A1BG-AS1, AC004477.3, AC004585.1 and AC004854.2) and two risk parameters (tumor subtype and TNM stage) were utilized as independent indicators to develop a novel prognostic model. In terms of predictive accuracy, the IRL-based model was distinctly superior to the TNM staging system (AUC: 0.728 VS 0.673, P = 0.010). DCA indicated that our nomogram had favorable clinical practicability. In addition, risk stratification analysis showed that the IRL-based tool efficiently divided BC patients into high- and low-risk groups (P < 0.001).
A novel IRL-based model was constructed to predict the risk of 5-year RFS in BC. Our model can improve the predictive power of the TNM staging system and identify high-risk patients with tumor recurrence to implement more appropriate treatment strategies.
Globally, breast cancer (BC) is the most frequently common carcinoma in women [1, 2]. In light of the statistics of the American Cancer Society, approximately 279,100 new BC cases and 42,690 cancer deaths are estimated to occur in the United States in 2020 . Although various therapeutic strategies are currently being applied to improve BC prognosis, including surgery, chemotherapy, radiotherapy, endocrine therapy and targeted therapy, many BC patients still have poor survival outcomes because of recurrence . Survival prognosis prediction is mainly based on the American Joint Committee on Cancer (AJCC) tumor-node-metastasis (TNM) staging system. However, it should be noted that BC patients at the same TNM stage may have distinct survival outcomes . This highlights that the TNM staging system is inadequate to achieve accurate prognosis evaluation and that it is unable to reflect the biological heterogeneity of BC. For example, the system is mainly based on clinicopathological parameters rather than on molecular signatures [6, 7]. Increasing evidence has revealed that molecular biomarkers have the potential to improve prognostic assessment and identify high-risk cancer patients [8,9,10,11,12,13,14]. Accordingly, there is an urgent need to screen effective molecular biomarkers for improving survival prognosis prediction and identify high-risk BC patients with tumor recurrence.
With advancements in transcriptomics, long noncoding RNAs (lncRNAs) have captured considerable attention in human cancers in the past decade [15,16,17,18,19]. LncRNAs are a class of mRNA-like transcripts with a length of > 200 nucleotides that lack protein-coding ability. Accumulating evidence indicates that lncRNAs play important roles in various biological processes, including transcriptional modifications, cell proliferation, differentiation, epigenetic modulation, and immune system-modulated pathways [20,21,22,23,24,25]. For example, lncRNA SNHG1 enhances the differentiation of Treg cells to provoke immune escape in BC . LncRNA AC025580.2 has a positive influence on immune escape in pancreatic cancer . In addition, a large number of lncRNAs have been confirmed to act as diagnostic and prognostic biomarkers or potential therapeutic targets in multiple tumors [11, 28,29,30,31,32]. However, the role of immune-related lncRNAs (IRLs) in the prognostic evaluation of BC remains unclear.
Hence, the objective of this study was to establish and validate a novel IRL model to predict the 5-year RFS of BC patients. This tool is able to improve the predictive accuracy of the TNM staging system and identify high-risk BC patients with tumor recurrence to facilitate optimal therapeutic schemes.
Patients and study design
Gene and lncRNA expression profiles as well as corresponding BC cases information were acquired from The Cancer Genome Atlas (TCGA) database. The inclusion criteria were as follows: (1) pathologically confirmed invasive BC, (2) complete follow-up data, including RFS status and survival time, and (3) survival time > 1 month. The following clinical parameters for each patient were collected: age, estrogen receptor (ER) status, progesterone receptor (PR) status, T stage, N stage, TNM stage, human epithelial growth factor receptor 2 (HER2) status, tumor subtype, RFS, and survival time. The Research Ethics Committee of Guangdong Provincial People’s Hospital (GDPH) approved our study. As the study used data from the database TCGA, written informed consent was waived.
Establishment of IRL risk score and the IRL-based model
First, immune-related genes (IRGs) were retrieved from the Molecular Signatures Database (MSDB) v4.0 (immune response M19817, and immune system process M13664) [33, 34]. According to the GENCODE project, lncRNA expression data were obtained from the gene expression profile . Next, IRLs were defined based on association analysis between the mRNA expression level and lncRNA expression data (|R|> 0.7, P < 0.05). Subsequently, univariate Cox proportional hazards regression analysis (CPHRA) was conducted to explore the association between the IRL expression level and patient RFS. Multivariate CPHRA was carried out to investigate independent prognostic IRLs. Regarding multivariate analysis, the collinearity diagnosis was confirmed via the variance inflation factor (VIF), and the final independent prognostic IRLs were selected to develop the IRL-based risk score, as follows: IRL risk score = sum of coefficients × expression level of IRLs. Finally, a novel IRL-based model incorporating the IRL signature and clinical variables was constructed to improve the prediction power of TNM staging in BC patients.
Assessment of the IRL-based nomogram
The area under the curve (AUC) of the time-dependent receiver operating characteristic (ROC) curve for 5-year RFS was calculated to evaluate the performance of the IRL-based model. The nomogram score for each patient was allocated on the basis of the IRL-based nomogram. Patients were divided into low- and high-risk nomogram score groups using the median score as the cutoff point. Subgroup analysis was applied to assess the risk stratification ability of the IRL-based model. A calibration curve was drawn to estimate the calibration ability of the IRL-based nomogram. Decision curve analysis (DCA) was used to weigh the clinical practicability of the IRL-based nomogram [36,37,38,39,40,41].
Gene enrichment analysis of the IRL signature
Gene enrichment analysis of the four IRLs was applied using Metascape, a free online method for gene annotation . Functional pathway enrichment was assessed on the basis of co-expressing genes of the four IRLs in the same module.
The χ2 test, Fisher's exact test or Mann–Whitney U test was employed to explore differences in indicators. The optimal cut-off values of IRL-based nomogram scores were calculated via X-tile software, version 3.6.1 (Yale University, New Haven, CT, USA). Kaplan–Meier survival analysis and the log-rank test were used to examine the 5-year RFS difference between low- and high-risk groups. All statistical analyses were utilized via R software (www.r-project.org, version 3.6.1) and Stata/MP, version 14.0 (StataCorp LP, College Station, TX). A P value < 0.05 was defined as significant.
In total, 911 BC samples were collected as the primary cohort in our study. As in our previous studies, each patient can be assigned a computer-generated allocation number (0 or 1) based on Stata software, and two cohorts were differentiated from the primary cohort according to the allocation number. Of the total, 456 patients (allocation number = 0) were considered as the validation cohort [4, 5]. The baseline characteristics of all patients in the two datasets are shown in Table 1. There were no significant differences between the variables examined in the two independent cohorts (all P > 0.05). The median age of the included patients was 58 years (IQR: 48–66) and 58 years (IQR: 48–65) in the two datasets.
Identification of the IRL signature
First, 332 IRGs were collected from Molecular Signatures Database (MSDB), and the expression profiles of the 332 IRGs were downloaded from TCGA. We defined a lncRNA with (|R|> 0.7 and P < 0.05) expression levels between the lncRNAs and genes as IRLs. Thus, on the basis of Pearson correlation analysis, 52 lncRNAs were confirmed as IRLs (|R|> 0.7 and P < 0.05). After univariate CPHRA, 4 IRLs were identified in the primary cohort (P < 0.05). Subsequent multivariate CPHRA further screened four IRLs (A1BG-AS1, AC004477.3, AC004585.1 and AC004854.2) as independent prognostic factors in the primary cohort. The VIFs ranged from 1.01 to 1.03, indicating no collinearity among the predictors. The detailed characteristics of the four IRLs in the training dataset are listed in Table 2.
Development of the IRL-based risk score and prognostic nomogram
An IRL-based risk score formula was built on the basis of IRL expression levels and their estimated regression coefficients in the multivariate CPHRA, where IRL-based risk score = −0.451 × expressionA1BG-AS1 + 0.294 × expressionAC004477.3–0.450 × expressionAC004585.1 + 0.627 × expressionAC004854.2. Based on univariate and multivariate CPHRA (Table 3), three prognostic variables (IRL-based risk score, tumor subtype, and TNM stage) were selected as independent indicators to build an IRL-based model. The novel IRL-based nomogram is illustrated in Fig. 1.
Evaluation of the IRL-based prognostic nomogram
The AUCs of the IRL-based nomogram were 0.728 (95% CI: 0.658–0.797) and 0.751 (95% CI: 0.656–0.846) in the primary and validation datasets, respectively, showing that this model had good predictive performance (Fig. 2a). Calibration plots revealed ideal agreement between the IRL-based nomogram-predicted likelihoods and the actual observations of 5-year RFS, indicating that this tool had high calibration ability (Fig. 2b). In addition, compared with the predictive accuracy between the IRL-based model and TMN stage, the time-dependent ROC curve suggested that the IRL-based nomogram outperformed the TNM stage for both the primary and validation cohorts (P < 0.05, Fig. 3a, b). Moreover, DCA indicated that the IRL-based tool added more net benefit than did the TNM stage for the primary and validation datasets; therefore, this nomogram showed superior clinical usefulness (Fig. 3c, d). On the basis of the IRL-based model, patients in the primary and validation cohorts were stratified into low- or high-risk groups using the median nomogram score as the cutoff point. The distribution characteristics of the IRL-based model score, RFS, and RFS status are shown in Fig. 4, indicating that patients with lower nomogram scores had better 5-year RFS than those with higher nomogram scores (P < 0.001). Subgroup analyses demonstrated that the IRL-based model had good risk stratification ability for T1 (P < 0.0001), T2 (P < 0.0001), T3 (P = 0.0067), N0 (P = 0.029), N1 (P < 0.0001), N2 (P = 0.013), HR + /HER2- (P = 0.00033), HR-/HER2 + (P = 0.038), and TNBC (P < 0.0001) (Fig. 5).
Functional enrichment analysis of the IRL signature
To further explore the potential functional roles of the four IRLs, significantly associated IRGs (Pearson coefficient > 0.7 and P < 0.05) were included in pathway enrichment via Metascape. The top 20 highly enriched pathways are displayed in Fig. 6. The IRGs of the four IRLs clustered most significantly in lymphocyte activation, cytokine-mediated signaling pathway, regulation of cytokine production, negative regulation of immune system process, regulation of immune effector process, alpha–beta T cell activation, leukocyte migration, antigen receptor-mediated signaling pathway, B cell activation, TCR pathway, response to bacterium, T cell costimulation, regulation of peptidyl-tyrosine phosphorylation, calcium-mediated signaling, myeloid leukocyte activation, positive regulation of defense response, lymphocyte migration, interleukin-10 signaling, positive regulation of cytosolic calcium ion concentration, and regulation of interferon-gamma production categories. Moreover, according to LncRRIsearch, a web tool used for comprehensive prediction of lncRNA-mRNA interactions , the predicted targets of A1BG-AS1 are FKBP10, NFYC, UBR4, CEACAM19, MT-ND4, PKD1P6, and RNA28S5. The predicted target of AC004477.3 is SEH1L.
With the increasing development of clinical management and comprehensive treatment of BC, prognosis of the survival of these patients has greatly improved. Nevertheless, many BC patients experience tumor recurrence . High-throughput biological technologies have provided a platform for exploring the molecular characteristics of different tumors . Many lncRNAs have been used to predict BC prognosis [44, 46, 47], though few IRL signatures have been identified as improving BC survival prediction. Consequently, it is necessary to single out prognostic IRLs integrating TNM stage to predict 5-year RFS of BC patients. In this study, with the combination of four IRLs, TNM stage and tumor subtype, we efficiently constructed and validated a novel IRL-based nomogram to accurately predict 5-year RFS in BC patients. In addition, time-dependent ROC analysis revealed that the IRL-based nomogram had more favorable prognostic accuracy than did TNM stage. Moreover, the DCA indicated that the IRL-based model had better clinical application than TNM stage; our IRL-based tool also exhibited good risk stratification ability to significantly categorize BC patients into high- and low-risk groups.
To date, several prognostic tools have been established to predict survival in BC patients [44, 46,47,48,49,50,51,52,53,54,55]. Yao et al. identified five lncRNAs as vital prognostic factors to evaluate BC prognosis . However, these models ignored the prognostic value of IRLs for BC patients. Recent studies have suggested that the immune response plays a critical role in cancer progression and recurrence, and accumulating evidence has revealed that IRLs harbor vital predictive value for survival prognosis. For example, Zhang et al. found that si-lncMALAT1 both suppressed osteosarcoma progression and resulted in self-destruction of cells deep within the tumor. This would enable clearance of osteosarcoma cells from the body by the immune system. Therefore, lncMALAT1 was individually correlated with immune system activity and overall survival . LncRNAs cannot encode proteins but do regulate gene expression at different levels, including through epigenetic, transcriptional or posttranscriptional regulation . In the past decade, this new type of gene regulator, namely, lncRNA, has been associated with tumor development, progression, and prognosis, and many lncRNAs have been found to act as prognostic indicators in various tumors. Moreover, lncRNAs can regulate the immune microenvironment of BC . IRLs can be applied to characterize the infiltration of immune cells in tumors. In our study, the IRL signature developed had a better prognostic value than the IRG signature. Thus, IRLs should be integrated into nomograms to improve prognosis prediction in BC. The survival prediction based on the TNM staging system is not yet satisfactory in clinical practice , though it is universally acknowledged that the TNM staging system is beneficial to achieve optimal treatment choice for BC patients. The TNM staging system is based on tumor size, lymph node status and metastasis . Nonetheless, malignant behavior in BC is mainly determined by the molecular characteristics. Thus, incorporating the TNM staging system and molecular biomarkers is indispensable for developing an accurate tool for prognostic assessment in BC patients .
Some limitations should be noted in the present study. First, TCGA does not include postoperative treatment information (chemotherapy and hormone therapy). Therefore, we are unable to screen low-risk patients to tailor adjuvant treatment. Second, the molecular mechanism of the IRLs should be further explored by additional experiments. Third, our established model should be verified by prospective, large-scale multicenter datasets.
In summary, we identified IRLs signatures significantly associated with the 5-year RFS of BC and efficiently constructed a novel prognostic model for 5-year RFS prediction by incorporating four IRLs, TNM stage and tumor subtype. Our model can improve the predictive power of the AJCC TNM staging system and confirm high-risk patients with tumor recurrence to receive appropriate treatment.
Availability of data and materials
All data was obtained from The Cancer Genome Atlas (TCGA, https://tcgadatanci.nih.gov/tcga/).
The Cancer Genome Atlas
Decision curve analysis
American Joint Committee on Cancer
Long non-coding RNAs
Human epithelial growth factor receptor 2
Guangdong Provincial People’s Hospital
Molecular Signatures Database
Least absolute shrinkage and selection operator
Cox proportional hazards regression analysis
Area under the curve
Receiver operating characteristic
Yap Y-S, Lu Y-S, Tamura K, Lee JE, Ko EY, Park YH, Cao AY, Lin C-H, Toi M, Wu J, Lee S-C. Insights Into Breast Cancer in the East vs the West. JAMA Oncology. 2019;12:63.
Harbeck N, Penault-Llorca F, Cortes J, Gnant M, Houssami N, Poortmans P, Ruddy K. Tsang J. Cardoso F: Breast cancer. Nature Reviews Disease Primers; 2019. p. 5.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7–30.
Lai J, Chen B, Zhang G, Wang Y, Mok H, Wen L, Pan Z, Su F, Liao N. Identification of a novel microRNA recurrence-related signature and risk stratification system in breast cancer. Aging (Albany NY). 2019;11:7525–36.
Lai J, Wang H, Pan Z, Su F. A novel six-microRNA-based model to improve prognosis prediction of breast cancer. Aging (Albany NY). 2019;11:649–62.
Lai J, Pan Z, Chen P, Ye G, Chen K, Su F. Development and validation of a nomogram incorporating axillary lymph node ratio to predict survival in node-positive breast cancer patients after neoadjuvant chemotherapy. Jpn J Clin Oncol. 2018;49:22–8.
Lai J, Wang H, Peng J, Chen P, Pan Z. Establishment and external validation of a prognostic model for predicting disease-free survival and risk stratification in breast cancer patients treated with neoadjuvant chemotherapy. Cancer Manag Res. 2018;10:2347–56.
Zhu X, Tian X, Sun T, Yu C, Cao Y, Yan T, Shen C, Lin Y, Fang JY, Hong J, Chen H. GeneExpressScore signature: a robust prognostic and predictive classifier in gastric cancer. Mol Oncol. 2018;33:459.
Wen J, Wang G, Xie X, Lin G, Yang H, Luo K, Liu Q, Ling Y, Xie X, Lin P, et al. Prognostic Value of a Four-miRNA signature in patients with lymph node positive locoregional esophageal squamous cell carcinoma undergoing complete surgical resection. Ann Surg. 2019;26:1.
Mitra AP, Lam LL, Ghadessi M, Erho N, Vergara IA, Alshalalfa M, Buerki C, Haddad Z, Sierocinski T, Triche TJ, et al. Discovery and validation of novel expression signature for postcystectomy recurrence in high-risk bladder cancer. J Natl Cancer Inst. 2014;22:106.
Zhu X, Tian X, Yu C, Shen C, Yan T, Hong J, Wang Z, Fang JY, Chen H. A long non-coding RNA signature to improve prognosis prediction of gastric cancer. Mol Cancer. 2016;15:60.
Qu L, Wang Z, Chen Q, Li Y, He H, Hsieh JJ, Xue S, Wu Z, Liu B, Tang H, et al. Prognostic value of a long non-coding RNA signature in localized clear cell renal cell carcinoma. Eur Urol. 2018;32:98.
Zhang J-X, Song W, Chen Z-H, Wei J-H, Liao Y-J, Lei J, Hu M, Chen G-Z, Liao B, Lu J, et al. Prognostic and predictive value of a microRNA signature in stage II colon cancer: a microRNA expression analysis. Lancet Oncol. 2013;14:1295–306.
Tang X-R, Li Y-Q, Liang S-B, Jiang W, Liu F, Ge W-X, Tang L-L, Mao Y-P, He Q-M, Yang X-J, et al. Development and validation of a gene expression-based signature to predict distant metastasis in locoregionally advanced nasopharyngeal carcinoma: a retrospective, multicentre, cohort study. The Lancet Oncology. 2018;65:523.
Zhou W, Gong J, Chen Y, Chen J, Zhuang Q, Cao J, Mei Z, Hu B. Long noncoding RNA LINC00899 suppresses breast cancer progression by inhibiting miR-425. Aging (Albany NY). 2019;11:425.
Zhang N, Zeng X, Sun C, Guo H, Wang T, Wei L, Zhang Y, Zhao J, Ma X. LncRNA LINC00963 Promotes Tumorigenesis and Radioresistance in Breast Cancer by Sponging miR-324-3p and Inducing ACK1 Expression. Mol Ther Nucleic Acids. 2019;18:871–81.
Zhang G, Xu Y, Zou C, Tang Y, Lu J, Gong Z, Ma G, Zhang W, Jiang P. Long noncoding RNA ARHGAP27P1 inhibits gastric cancer cell proliferation and cell cycle progression through epigenetically regulating p15 and p16. Aging (Albany NY). 2019;11:9090–110.
Yang Z, An Y, Wang N, Dong X, Kang H. LINC02595 promotes tumor progression in colorectal cancer by inhibiting miR-203b-3p activity and facilitating BCL2L1 expression. J Cell Physiol. 2020;28:65.
Tan BS, Yang MC, Singh S, Chou YC, Chen HY, Wang MY, Wang YC, Chen RH. LncRNA NORAD is repressed by the YAP pathway and suppresses lung and breast cancer metastasis by sequestering S100P. Oncogene. 2019;38:5612–26.
Shen SN, Li K, Liu Y, Yang CL, He CY, Wang HR. Silencing lncRNAs PVT1 Upregulates miR-145 and confers inhibitory effects on viability, invasion, and migration in EC. Mol Ther Nucleic Acids. 2019;19:668–82.
Ren X, Chen C, Luo Y, Liu M, Li Y, Zheng S, Ye H, Fu Z, Li M, Li Z, Chen R. lncRNA-PLACT1 sustains activation of NF-kappaB pathway through a positive feedback loop with IkappaBalpha/E2F1 axis in pancreatic cancer. Mol Cancer. 2020;19:35.
Ji W, Diao YL, Qiu YR, Ge J, Cao XC, Yu Y. LINC00665 promotes breast cancer progression through regulation of the miR-379-5p/LIN28B axis. Cell Death Dis. 2020;11:16.
Guan Z, Wang Y, Wang Y, Liu X, Wang Y, Zhang W, Chi X, Dong Y, Liu X, Shao S, Zhan Q. Long non-coding RNA LOC100133669 promotes cell proliferation in oesophageal squamous cell carcinoma. Cell Prolif. 2020;2020:e12750.
Dong HT, Liu Q, Zhao T, Yao F, Xu Y, Chen B, Wu Y, Zheng X, Jin F, Li J, Xing P. Long non-coding RNA LOXL1-AS1 drives breast cancer invasion and metastasis by antagonizing miR-708-5p expression and activity. Mol Ther Nucleic Acids. 2019;19:696–705.
Fu J, Dong G, Shi H, Zhang J, Ning Z, Bao X, Liu C, Hu J, Liu M, Xiong B. LncRNA MIR503HG inhibits cell migration and invasion via miR-103/OLFM4 axis in triple negative breast cancer. J Cell Mol Med. 2019;2019:38.
Pei X, Wang X, Li H. LncRNA SNHG1 regulates the differentiation of Treg cells and affects the immune escape of breast cancer via regulating miR-448/IDO. Int J Biol Macromol. 2018;118:24–30.
Chen J, Wang S, Jia S, Ding G, Jiang G, Cao L. Integrated analysis of long non-coding RNA and mRNA expression profile in pancreatic cancer derived exosomes treated dendritic cells by microarray analysis. J Cancer. 2018;9:21–31.
Zhou M, Zhao H, Wang Z, Cheng L, Yang L, Shi H, Yang H, Sun J. Identification and validation of potential prognostic lncRNA biomarkers for predicting survival in patients with multiple myeloma. J Exp Clin Cancer Res. 2015;34:102.
Zhou M, Zhang Z, Zhao H, Bao S, Cheng L, Sun J. An immune-Related Six-lncRNA signature to improve prognosis prediction of glioblastoma multiforme. Mol Neurobiol. 2018;55:3684–97.
Zhou H, Zhang H, Chen J, Cao J, Liu L, Guo C, Huang G, Zeng D. A seven-long noncoding RNA signature predicts relapse in patients with early-stage lung adenocarcinoma. J Cell Biochem. 2019;12:315.
Zhang C, Li Z, Hu J, Qi F, Li X, Luo J. Identification of five long noncoding RNAs signature and risk score for prognosis of bladder urothelial carcinoma. Journal of Cellular Biochemistry. 2019. https://doi.org/10.1002/jcb.29330.
Zhang B, Wang H, Guo Z, Zhang X. Prediction of head and neck squamous cell carcinoma survival based on the expression of 15 lncRNAs. J Cell Physiol. 2019;36:15.
Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, Schinzel AC, Sandy P, Meylan E, Scholl C, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462:108–12.
Wang W, Zhao Z, Yang F, Wang H, Wu F, Liang T, Yan X, Li J, Lan Q, Wang J, Zhao J. An immune-related lncRNA signature for patients with anaplastic gliomas. J Neurooncol. 2018;136:263–71.
Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, Levy R. Prediction of survival in diffuse large-B-cell lymphoma based on the expression of six genes. N Engl J Med. 2004;350:1828–37.
Vickers AJ, Van Calster B, Steyerberg EW. Net benefit approaches to the evaluation of prediction models, molecular markers, and diagnostic tests. BMJ. 2016;352:i6.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. 2006;26:565–74.
Kastrinos F, Ojha RP, Leenen C, Alvero C, Mercado RC, Balmana J, Valenzuela I, Balaguer F, Green R, Lindor NM, et al. Comparison of prediction models for lynch syndrome among individuals with colorectal cancer. J Natl Cancer Inst. 2016;22:108.
Huang YQ, Liang CH, He L, Tian J, Liang CS, Chen X, Ma ZL, Liu ZY. Development and validation of a radiomics nomogram for preoperative prediction of lymph node metastasis in colorectal cancer. J Clin Oncol. 2016;34:2157–64.
Han SS, Rivera GA, Tammemagi MC, Plevritis SK, Gomez SL, Cheng I, Wakelee HA. Risk stratification for second primary lung cancer. J Clin Oncol. 2017;35:2893–9.
Kerr KF, Brown MD, Zhu K, Janes H. Assessing the clinical impact of risk prediction models with decision curves: guidance for correct interpretation and appropriate use. J Clin Oncol. 2016;34:2534–40.
Jassal B, Matthews L, Viteri G, Gong C, Lorente P, Fabregat A, Sidiropoulos K, Cook J, Gillespie M, Haw R, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020;48:D498-d503.
Fukunaga T, Iwakiri J, Ono Y, Hamada M. LncRRIsearch: A web server for lncRNA-RNA interaction prediction integrated with tissue-specific expression and subcellular localization data. Front Genet. 2019;10:462.
Tang J, Ren J, Cui Q, Zhang D, Kong D, Liao X, Lu M, Gong Y, Wu G. A prognostic 10-lncRNA expression signature for predicting the risk of tumour recurrence in breast cancer patients. J Cell Mol Med. 2019;2:35.
Lai J, Chen B, Mok H, Zhang G, Ren C, Liao N. Comprehensive analysis of autophagy-related prognostic genes in breast cancer. J Cell Mol Med. 2020;24:9145–53.
Yao Y, Zhang T, Qi L, Zhou C, Wei J, Feng F, Liu R, Sun C. Integrated analysis of co-expression and ceRNA network identifies five lncRNAs as prognostic markers for breast cancer. J Cell Mol Med. 2019;20:98.
Meng J, Li P, Zhang Q, Yang Z, Fu S. A four-long non-coding RNA signature in predicting breast cancer survival. J Exp Clin Cancer Res. 2014;33:84.
Wang K, Liao C, Zhong Q, Dong H, Zhang T, Jin R. CeNETs analysis reveals the prognostic value of a signature integration from five lncRNAs in breast cancer. J Cell Biochem. 2019;56:33.
Tian T, Gong Z, Wang M, Hao R, Lin S, Liu K, Guan F, Xu P, Deng Y, Song D, et al. Identification of long non-coding RNA signatures in triple-negative breast cancer. Cancer Cell Int. 2018;18:103.
Tang J, Cui Q, Zhang D, Kong D, Liao X, Ren J, Gong Y, Wu G. A prognostic eight-lncRNA expression signature in predicting recurrence of ER-positive breast cancer receiving endocrine therapy. J Cell Physiol. 2019;235:4746.
Sun M, Wu D, Zhou K, Li H, Gong X, Wei Q, Du M, Lei P, Zha J, Zhu H, et al. An eight-lncRNA signature predicts survival of breast cancer patients: a comprehensive study based on weighted gene co-expression network analysis and competing endogenous RNA network. Breast Cancer Res Treat. 2019;175:59.
Li J, Wang W, Xia P, Wan L, Zhang L, Yu L, Wang L, Chen X, Xiao Y, Xu C. Identification of a five-lncRNA signature for predicting the risk of tumor recurrence in patients with breast cancer. Int J Cancer. 2018;143:2150–60.
Li J, Gao C, Liu C, Zhou C, Ma X, Li H, Li J, Wang X, Qi L, Yao Y, et al. Four lncRNAs associated with breast cancer prognosis identified by coexpression network analysis. J Cell Physiol 2019.
He Y, Li X, Meng Y, Fu S, Cui Y, Shi Y, Du H. A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma. J Cell Biochem 2019.
Fan CN, Ma L, Liu N. Comprehensive analysis of novel three-long noncoding RNA signatures as a diagnostic and prognostic biomarkers of human triple-negative breast cancer. J Cell Biochem 2018.
Zhang ZC, Tang C, Dong Y, Zhang J, Yuan T, Li XL. Targeting LncRNA-MALAT1 suppresses the progression of osteosarcoma by altering the expression and localization of beta-catenin. J Cancer. 2018;9:71–80.
Cai P, Otten A, Cheng B, Ishii M, Zhang W, Huang B, Qu K, Sun B. PRANCRA genome-wide long noncoding RNA CRISPRi screen identifies as a novel regulator of epidermal homeostasis. Genome Res. 2020;30:22–34.
Li Z, Li Y, Wang X, Yang Q. Identification of a Six-Immune-Related Long Non-coding RNA Signature for Predicting Survival and Immune Infiltrating Status in Breast Cancer. Front Genet. 2020;11:680.
Pennathur A, Luketich JD. Resection for esophageal cancer: strategies for optimal management. Ann Thorac Surg. 2008;85:S751-756.
Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J CLIN ONCOL. 2008;26:1364–70.
Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16:e173–80.
This study was supported by the Doctor Launch Fund of Guangdong Provincial People's hospital (grant number 2020bq11), the National Natural Science Foundation of China (Grants 82002928; 81902828; 81602645; 81071851; 81001189), National Natural Science Foundation of Guangdong Province (Grant 2016A030313768; 2018A030313292), Medical Scientific Research Foundation of Guangdong Province (B2019039, Bo Chen), Research Funds from Guangzhou Municipal Science and Technology Project (201707010418; 201804010430).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lai, J., Chen, B., Zhang, G. et al. Molecular characterization of breast cancer: a potential novel immune-related lncRNAs signature. J Transl Med 18, 416 (2020). https://doi.org/10.1186/s12967-020-02578-4
- Breast cancer