Systematic analysis of lncRNA–miRNA–mRNA competing endogenous RNA network identifies four-lncRNA signature as a prognostic biomarker for breast cancer
Journal of Translational Medicine volume 16, Article number: 264 (2018)
Increasing evidence has underscored the role of long non-coding RNAs (lncRNAs) acting as competing endogenous RNAs (ceRNAs) in the development and progression of tumors. Nevertheless, lncRNA biomarkers in lncRNA-related ceRNA network that can predict the prognosis of breast cancer (BC) are still lacking. The aim of our study was to identify potential lncRNA signatures capable of predicting overall survival (OS) of BC patients.
The RNA sequencing data and clinical characteristics of BC patients were obtained from the Cancer Genome Atlas database, and differentially expressed lncRNA (DElncRNAs), DEmRNAs, and DEmiRNAs were then identified between BC and normal breast tissue samples. Subsequently, the lncRNA–miRNA–mRNA ceRNA network of BC was established, and the gene oncology enrichment analyses for the DEmRNAs interacting with lncRNAs in the ceRNA network was implemented. Using univariate and multivariate Cox regression analyses, a four-lncRNA signature was developed and used for predicting the survival in BC patients. We applied receiver operating characteristic analysis to assess the performance of our model.
A total of 1061 DElncRNAs, 2150 DEmRNAs, and 82 DEmiRNAs were identified between BC and normal breast tissue samples. A lncRNA–miRNA–mRNA ceRNA network of BC was established, which comprised of 8 DEmiRNAs, 48 DElncRNAs, and 10 DEmRNAs. Further gene oncology enrichment analyses revealed that the DEmRNAs interacting with lncRNAs in the ceRNA network participated in cell leading edge, protease binding, alpha-catenin binding, gamma-catenin binding, and adenylate cyclase binding. A univariate regression analysis of the DElncRNAs revealed 7 lncRNAs (ADAMTS9-AS1, AC061992.1, LINC00536, HOTAIR, AL391421.1, TLR8-AS1 and LINC00491) that were associated with OS of BC patients. A multivariate Cox regression analysis demonstrated that 4 of those lncRNAs (ADAMTS9-AS1, LINC00536, AL391421.1 and LINC00491) had significant prognostic value, and their cumulative risk score indicated that this 4-lncRNA signature independently predicted OS in BC patients. Furthermore, the area under the curve of the 4-lncRNA signature associated with 3-year survival was 0.696.
The current study provides novel insights into the lncRNA-related ceRNA network in BC and the 4 lncRNA biomarkers may be independent prognostic signatures in predicting the survival of BC patients.
Breast cancer (BC) is a heterogeneous and malignant neoplasm derived from breast tissue, and accounts for about 16% of all cancers and 22.9% of invasive cancers in women . The most common cause of the BC-related mortality is metastasis [2, 3]. Currently, BC diagnosis and prognosis is evaluated on the basis of disease stage, histological grade, and the expression level of hormone receptors . However, clinical and pathological symptoms have limited predictive value in detecting early BC, and the clinical outcomes are highly variable on account of its heterogeneity. In addition, the underlying molecular mechanisms of BC still remain unclear. Therefore, it is vital to identify potential molecular diagnostic markers and/or therapeutic targets to combat BC, especially the invasive form.
Long non-coding RNAs (lncRNAs) is a class of ncRNA over 200 nucleotides long , and are reportedly involved in a number of cellular processes, for example, transcriptional and post-transcriptional regulation [5, 6]. Due to their strong tissue specificity, lncRNAs are potentially effective early diagnostic biomarkers of various cancers . Identification of a BC specific lncRNA biomarker may therefore be of clinical significance for the diagnosis and prognosis of BC. Several lncRNAs have been reported to be associated with BC initiation and progression [8, 9], and although some have been found to predict clinical outcomes for BC, the results are inconsistent due to limited tissue samples. Furthermore, studies without large sample size are also not able to determine with statistical power whether these lncRNAs are associated with survival or other clinical factors. The Cancer Genome Atlas (TCGA) is an open-access and large-scale database which can provide multidimensional molecular profiles for a large number tumor samples. To increase the statistical reliability of our studies, we identified BC specific lncRNAs using data obtained from TCGA database.
The competing endogenous RNA (ceRNA) hypothesis presented by Salmena et al.  was proposed as a novel regulatory mechanism between ncRNA and coding messenger RNA. LncRNAs contain miRNA-response elements (MREs) which function as ceRNAs, and play a key role in various pathological processes like tumorigenesis . Zhang et al.  have proved the biological role of lncRNA related-ceRNAs in glioblastomas. A recent study has demonstrated that lncRNA NUTF2P3-001 acts as a ceRNA to communicate with KRAS by competitively binding to hsa-mir-3923, and the up-regulation of NUTF2P3-001 reverses the suppressive effect of hsa-mir-3923 on KRAS, leading to the proliferation and invasion of pancreatic cancer . In addition, the aberrant expression of 7-lncRNA signature (called LncRisk-7) led to differential gene expression via a dysregulated lncRNA-associated ceRNA network, contributing to pancreatic ductal adenocarcinoma progression . Collectively, these findings show that dysregulation of important lncRNAs in the ceRNA network also disrupt the miRNA-mediated lncRNA/mRNA ceRNA interactions and therefore contribute to cancer initiation and progression [15, 16]. Nevertheless, very little information is available on BC ceRNAs.
In our work, RNA sequencing data of 1109 BC samples and 113 adjacent non-tumor breast tissues samples were retrieved from the TCGA database. To the best of our knowledge, this is the first study to use large scale sequencing database (TCGA) and ceRNA network to identify BC-specific lncRNAs. This new approach of predicting cancer specific lncRNA and ceRNA networks can elucidate the lncRNA-mediated ceRNA regulatory mechanisms in the development and prognosis of BC, and identify novel lncRNAs as potential diagnostic biomarkers or therapeutic targets.
Patients and samples from the TCGA database
RNA sequencing (RNA-Seq) data associated with BC were retrieved from the TCGA database (https://portal.gdc.cancer.gov/, version 10.1, release time: February 15, 2018). A total of 1222 individuals with BC were included in the current study. The exclusion criteria were (1) histological diagnosis negating BC, (2) presence of a malignancy other than BC, and (3) lack of complete clinical data. The gene expression profiles of 1109 BC and 113 adjacent normal breast tissues, and miRNA data of 1103 BC and 104 adjacent normal breast tissues were downloaded. In addition to the RNA expression data, clinical information of BC patients were also downloaded from the TCGA database. No approval from the ethics committee was needed because all the information was required from the TCGA database. The clinical characteristics for BC patients are listed in Table 1.
RNA sequence data processing and differential expression analysis
The raw RNA sequencing (lncRNA, miRNA, and mRNA) reads were post-processed and normalized using the trimmed mean of M-values (TMM) method.
EdgeR package in R (version 3.4.1) was used to identify the differentially expressed mRNAs (DEmRNAs), lncRNAs (DElncRNAs) and miRNAs (DEmiRNAs) between the BC and adjacent-normal breast tissues , and the cut-off criteria were set as P < 0.01 and |logFC| > 2. Volcano plots were visualized using the ggplot2packages in R . The heat map was plotted using the pheatmap function of pheatmap package version 1.0.8 .
Establishment of the ceRNA network
The lncRNA–miRNA–mRNA ceRNA network was constructed based on the hypothesis that lncRNAs directly interact with and regulate the activity of mRNAs by acting as miRNA sponges . Based on this hypothesis, we established the lncRNA–miRNA–mRNA ceRNA network in three steps: (1) BC-specific RNAs (lncRNA, mRNA, and miRNA) with P < 0.01, and |logFC| > 2 were reserved, (2) the potential miRNAs targeted by DElncRNAs and the lncRNA–miRNA interactions were predicted by the miRcode online tool (http://www.mircode.org), and (3) the MiRDB (http://www.mirdb.org/), miRTarBase (http://mirtarbase.mbc.nctu.edu.tw//), and Targetscan (http://www.targetscan.org//) programs were used to predict the target mRNAs of miRNAs. Finally, the miRNAs that were negatively regulated by the lncRNAs and mRNAs were selected to build the ceRNA network. Cytoscape (version 3.5.1) was used to visualize the lncRNA–miRNA–mRNA ceRNA network.
Functional enrichment analysis
Gene oncology (GO) is widely used as functional enrichment analysis for a large number of genes . The putative biological roles of DElncRNAs corresponds to that of their associated mRNAs. GO function analyses were therefore conducted for the DEmRNAs in the ceRNA network using R clusterProfiler package . Fisher’s test was used to identify the significant GO terms, and GO categories with P < 0.05 were considered statistically significant.
Construction of the BC-specific prognostic signatures
Kaplan–Meier and log-rank test was used to determine the association between the DEmRNAs, DElncRNAs and DEmiRNAs in the ceRNA network and the overall survival (OS) of BC patients, and statistical significance was set at P < 0.05. Univariate Cox proportional hazards regression method was implemented to analyze the relationship between the DElncRNAs and OS when a significant level was set at 0.05, in order to determine those with a prognostic value in BC. This was followed by multivariate Cox hazards regression model to determine the independent prognostic factors for BC, and the prognostic risk score for predicting OS was as follows:
(where “exp” denotes the expression level of DElncRNAs, and “β” is the regression coefficient obtained from the multivariate Cox regression model) . Using the median risk score as the threshold, the BC patients were stratified into the high- and low-risk groups. The “survival ROC” package in R was used to construct the time-dependent receiver operating characteristic (ROC) curves within 3 years as the defining point, and to measure the risk prediction rate of specific lncRNAs between the two groups. In addition, the univariate and multivariate analyses were used to evaluate the effects of other clinical variables of BC patients on OS risk scores. The R software (version 3.4.1) was used for all statistical analyses.
Identification of DEmRNAs, DElncRNAs, and DEmiRNAs
We identified the DEmRNAs, DElncRNAs, and DEmiRNAs in BC and adjacent-normal breast tissues using the TCGA database, with P < 0.01 and |logFC| > 2 as the thresholds. A total of 2150 DEmRNAs (1368 up- and 782 down-regulated), 1061 DElncRNAs (839 up-, and 222 down-regulated), and 82 DEmiRNAs (62 up- and 20 down-regulated) were identified between BC and normal samples. Volcano plots displaying the distribution of the DElncRNAs, DEmiRNAs, and DEmRNAs were generated, as shown in Fig. 1a. The heat map showed clear separation and consistency in the expression profiles of the BC and normal samples (Fig. 1b).
MiRNA predicted target analysis and ceRNA network establishment
The differentially expressed RNAs identified above were selected, and the lncRNAs and mRNAs targeted by miRNAs were extracted to establish the lncRNA–miRNA–mRNA ceRNA network. The relationships among 1061 DElncRNAs and 82 DEmiRNAs were first evaluated. Since lncRNAs might interact with the miRNAs through MREs, the miRcode tool was then used to detect the potential MREs; 18 BC-specific miRNAs that putatively target 70 BC-specific lncRNAs were then identified (Additional file 1: Table S1). The MiRDB, miRTarBase and Targetscan programs were then used to determine the relationship between the 82 DEmiRNAs and 2150 DEmRNAs, and predict the mRNA targets of miRNAs. The results indicated that 8 BC specific miRNAs targeted 10 BC-specific mRNAs (Additional file 1: Table S2).
On the basis of the above data, the lncRNA–miRNA–mRNA ceRNA network was established and plotted using Cytoscape 3.5.1. Overall, 8 miRNAs (6 up-, and 2 down-regulated, Table 2), 48 lncRNAs (32 up-, and 16 down-regulated, Table 3), and 10 mRNAs (2 up-, and 8 down-regulated, Table 4) were involved in the proposed ceRNA network (Fig. 2). Based on the expression levels of DEmRNAs, DElncRNAs and DEmiRNAs, two ceRNA networks including under-expressed (Fig. 2a) and over-expressed (Fig. 2b) networks were constructed.
Delineation of GO analysis
In order to better understand the role of the DElncRNAs in BC, we analyzed the mRNAs of the ceRNA network and identified the lncRNA regulated GO terms. The results of GO analysis are shown in Additional file 2: Figure S1. Our data showed that the mRNAs associated to cellular component (CC), was cell leading edge. Ge et al.  have demonstrated that trypsin secreted from MDA MB-231 BC cells activates the protease-activated receptor-2 and the activated protease-activated receptor-2 can promote cell migration based on ERK1/2-dependent pathway, involving the formation of a scaffolding complex at the cell leading edge. Meanwhile, the mRNAs related to molecular function (MF) were most relevant to protease binding, alpha-catenin binding, gamma-catenin binding, and adenylate cyclase binding. Proteases provide the cancer a characteristic of being able to invade into other tissues, and protease-activated receptor-1 is involved in the migration and invasion of breast cancer cells . The loss of alpha-catenin has been implicated to be related to the metastasis and poor survival in BC . The expression of gamma-catenin was also reported to be associated with the metastasis in human BC . The activity of the adenylate cyclase is positively linked to the inhibition of cell proliferation, as well as induction of apoptosis in human BC MCF-7 cells . Demonstrated herein, these GO terms are associated with BC pathogenesis and prognosis.
Correlations between BC specific signatures and OS
Kaplan–Meier and log-rank test were used to determine the relationship between the DEmRNAs, DElncRNAs and DEmiRNAs in the ceRNA network and the OS of BC patients with a cut-off threshold of P < 0.05. Totally, 4 DElncRNAs (ADAMTS9-AS1, AL356479.1, CHL1-AS2, and LINC00536; Fig. 3a), 3 DEmiRNAs (hsa-miR-204, hsa-miR-210, and hsa-miR-429; Fig. 3b), and 1 DEmRNA (KPNA2; Fig. 3c) were found to be related to OS.
Establishment of the 4-lncRNAs prognostic model
Univariate regression analysis was used to identify the lncRNAs associated with the OS of BC patients. With the significance level cutoff threshold set at P < 0.05, a group of lncRNA signatures including ADAMTS9-AS1, AC061992.1, LINC00536, HOTAIR, AL391421.1, TLR8-AS1, and LINC00491 lncRNAs was detected to have significant prognostic value (Additional file 1: Table S3). Significantly, we found that ADAMTS9-AS1 was simultaneously identified to be connected with OS in Kaplan–Meier (log-rank test) and univariate Cox regression analysis.
All the above lncRNAs were then fitted into the multivariate Cox regression model, which indicated that only four lncRNAs—ADAMTS9-AS1, LINC00536, AL391421.1 and LINC00491—had a significant prognostic value in BC (Additional file 3: Figure S2C), and these four lncRNAs were used to develop an lncRNA prognostic model. A risk score analysis of the four lncRNAs was performed for each patient, and based on the risk scores, the patients were divided into the “low risk” and “high risk” groups (Additional file 3: Figure S2A). The mortality rate of the high risk patients was significantly higher compared to the low risk patients (12.99% vs 5.94%; P < 0.05; Additional file 3: Figure S2B).
In addition, the high risk group was correlated with worse prognosis compared to the low risk group (Fig. 4a). The 3 year survival correlation of the 4-lncRNA signature was analyzed by ROC and AUC was computed to assess the discriminatory capacity of lncRNA signature (Fig. 4b). The AUC of the 4-lncRNA signature was 0.696 indicating its utility as a prognostic model for predicting the survival status of BC.
The 4-lncRNA expression was then analyzed in the tumor and normal tissues, and in the high- and low-risk patient groups (Fig. 5). ADAMTS9-AS1, and AL391421.1 were expressed at high levels in patients with low-risk scores, whereas LINC00536 and LINC00491 were up-regulated in the high-risk patients. Furthermore, LINC00491, AL391421.1, and LINC00536 were expressed at high levels, and ADAMTS9-AS1 was expressed at low levels in the BC patients.
Prognostic value of the four-lncRNA signature in BC
Univariate and multivariate regression models were used to assess the prognostic power of the 4-lncRNA signature. Univariate analysis indicated that age, pathological stage, N stage, M stage, ER, PR, Her2, and risk scores were significantly correlated with OS of BRCA patients (P < 0.05). Similarly, Kaplan–Meier analysis demonstrated that clinical factors (age, pathological stage, N stage, M stage, ER, PR, and Her2) were significantly correlated to OS, which was consistent with the univariate analysis. Kaplan–Meier curves of the clinical characteristics are shown in Fig. 6. Multivariate analysis indicated that only age, pathological stage, Her2, and risk scores were independent prognostic factors of OS (P= 0.003, 0.024, 0.034, and 0.017 respectively; Table 1).
BC is a common malignant gynecological cancer, and is one of the main causes for the cancer-related deaths in women . The lack of specific diagnostic and prognostic biomarkers may contribute to the current low survival rate among BC patients. To improve clinical outcomes therefore, it is essential to explore the exact regulatory mechanisms of BC initiation and progression, and to identify the potential BC-related prognostic signatures that predict those outcomes. Growing experimental evidence indicates that lncRNAs play important roles in many biological processes, and ceRNA activity is closely related to the development of cancers [30, 31].
In recent years, some studies have investigated the ceRNAs in BC. For instance, Chen et al.  analyzed the BC ceRNA network on the basis of common miRNAs as well as co-expression, but did not consider miRNA expression. Another study also established a BC specific ceRNA network to investigate its underlying molecular mechanisms based on the PCC of miRNA–mRNA pairs . However, both studies focused on the roles of mRNAs rather than that of lncRNAs in the BC ceRNA networks. In 2018, Zhou et al.  constructed four BC-related ceRNA networks by combining the miRNA targets and the expression data of lncRNAs, miRNAs and mRNA, but they did not take into account the relationship between survival and lncRNAs, nor construct the prognostic signature. In our study, in addition to constructing the ceRNA networks by combining lncRNA, miRNA, and mRNA expression data, we also investigated the association of lncRNA and OS in BC patients. Furthermore, based on the theory of ceRNA network, we established the 4-lncRNAs prognostic signature. With the goal of identifying lncRNAs significantly associated with OS, we established an lncRNA–miRNA–mRNA ceRNA network using the information obtained from the TCGA database. Univariate regression analysis on the DElncRNAs of the ceRNA network identified 7 lncRNAs—ADAMTS9-AS1, AC061992.1, LINC00536, HOTAIR, AL391421.1, TLR8-AS1 and LINC00491—that were associated with OS. Multivariate analysis showed significant prognostic value of 4 of those lncRNAs (ADAMTS9-AS1, LINC00536, AL391421.1 and LINC00491) in the OS of BC patients. A cumulative risk score of the 4 lncRNAs was calculated, which indicated that this 4-lncRNA signature independently predicted OS in BC patients. To the best of our knowledge, this is the first report integrating a ceRNA network with TCGA data to build an lncRNA-related risk score, and evaluate the OS of BC patients. Our study will help improve the understanding of lncRNA-mediated ceRNA regulatory mechanisms in BC and identify novel lncRNAs as therapeutic targets.
In the current study, among this 4-lncRNA signature, ADAMTS9-AS1 was demonstrated to play important roles in the progression and prognosis of cancer. ADAMTS9-AS1 is an antisense lncRNA, and growing evidence has implicated that a large amount of antisense lncRNAs play crucial roles in the cancer [35, 36]. Li et al.  reported that ADAMTS9-AS1 could predict the survival status of patients with esophageal squamous cell carcinoma. Another study has also indicated a prognostic role of ADAMTS9-AS1 in patients with colon adenocarcinoma . In addition, ADAMTS9-AS1 has been demonstrated to be a risk lncRNA in ovarian cancer, which is involved in the progression of ovarian cancer . In our study, we noticed that ADAMTS9-AS1 with low-expression could compete with up-regulated miRNAs (hsa-mir-182, and hsa-mir-21), to regulate the expression of the target genes such as CHL1, SPRY2, and TCEAL7 involved in the ceRNA network. Previous studies have shown high expression of hsa-mir-182 in MCF-7 breast cancer cells [40, 41]. In addition, the high-expression of hsa-mir-21 was reported to be correlated to the metastasis and poor prognosis of BC patients . The remaining three lncRNAs of the ceRNA network (LINC00536/AL391421.1/LINC00491) were up-regulated and competed with the decreased hsa-mir-204 and hsa-mir-144 levels. Down-regulation of has-mir-204 has been suggested to enhance cell proliferation and invasion in gastric cancer , and low-expression of has-mir-204 is related to the poor prognosis of acute myeloid leukemia patients . In addition, decreased expression of has-mir-144 is strongly correlated with the progression of colorectal cancer . No study so far has reported any association of LINC00536, AL391421.1 or LINC00491 with cancer. This is the first study to show aberrant expression of ADAMTS9-AS1, LINC00536, AL391421.1 and LINC00491 in BC, and indicates a potential prognostic role of this 4-lncRNA signature in BC. In addition, the bioinformatics based investigation of lncRNAs will be helpful in future experimental studies.
Although the findings of our study have important clinical implications, the limitations must also be noted. First, a longer follow-up duration is required to verify our results, and second, the findings based on the TCGA database will need to be verified using other experimental methods. In addition, the biological roles of ADAMTS9-AS1, LINC00536, AL391421.1, and LINC00491 in BC also need to be further investigated.
Taken together, we have identified a 4-lncRNA signature as a potential prognostic predictor for BC patients by analyzing the genome-wide lncRNA expression data from the TCGA database based on a ceRNA network. The current findings provide novel insights into the lncRNA-related ceRNA network in BC and identify potential diagnostic and prognostic biomarkers. Further functional studies are needed to elucidate the molecular mechanisms underlying lncRNA function in BC.
long non-coding RNAs
competing endogenous RNAs
the Cancer Genome Atlas
differentially expressed lncRNA
receiver operating characteristic
Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013;63(1):11.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646.
Gupta GP, Massagué J. Cancer metastasis: building a framework. Cell. 2006;127(4):679.
Vuong D, et al. Molecular classification of breast cancer. Virchows Archiv. 2014;465(1):1–14.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009;10(3):155–9.
Fatica A, Bozzoni I. Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014;15(1):7–21.
Hung T, Chang HY. Long noncoding RNA in genome regulation: prospects and mechanisms. RNA Biol. 2010;7(5):582.
Su X, et al. Comprehensive analysis of long non-coding RNAs in human breast cancer clinical subtypes. Oncotarget. 2014;5(20):9864.
Jin M, et al. A four-long non-coding RNA signature in predicting breast cancer survival. J Exp Clin Cancer Res. 2014;33(1):84.
Salmena L, et al. ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Song X, et al. Analysing the relationship between lncRNA and protein-coding gene and the role of lncRNA as ceRNA in pulmonary fibrosis. J Cell Mol Med. 2014;18(6):991–1003.
Zhang K, et al. Identification and functional characterization of lncRNAs acting as ceRNA involved in the malignant progression of glioblastoma multiforme. Oncol Rep. 2016;36(5):2911–25.
Xiang L, et al. Hypoxia-induced lncRNA-NUTF2P3-001 contributes to tumorigenesis of pancreatic cancer by derepressing the miR-3923/KRAS pathway. Oncotarget. 2016;7(5):6000–14.
Zhou M, et al. Construction and analysis of dysregulated lncRNA-associated ceRNA network identified novel lncRNA biomarkers for early diagnosis of human pancreatic cancer. Oncotarget. 2016;7(35):56383–94.
Karreth FA, Pandolfi PP. ceRNA cross-talk in cancer: when ce-bling rivalries go awry. Cancer Discov. 2013;3(10):1113–21.
Tan JY, et al. Extensive microRNA-mediated crosstalk between lncRNAs and mRNAs in mouse embryonic stem cells. Genome Res. 2015;25(5):655.
Robinson MD, McCarthy DJ, Smyth GK. EdgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26(1):139–40.
Wickham H. ggplot2: elegant graphics for data analysis. Houston: Springer Publishing Company, Incorporated; 2009. p. 180–5.
Kolde R. pheatmap: pretty heatmaps; 2015.
Guo LL, et al. Competing endogenous RNA networks and gastric cancer. World J Gastroenterol. 2015;21(41):11680.
Ashburner M, et al. Gene ontology: tool for the unification of biology. Gene ontology consortium. Nat Genet. 2000;25(1):25–9.
Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics-a J Integr Biol. 2012;16(5):284–7.
Zeng JH, et al. Comprehensive investigation of a novel differentially expressed lncRNA expression profile signature to assess the survival of patients with colorectal adenocarcinoma. Oncotarget. 2017;8(10):16811.
Ge L, et al. Constitutive protease-activated receptor-2-mediated migration of MDA MB-231 breast cancer cells requires both beta-arrestin-1 and -2. J Biol Chem. 2004;279(53):55419–24.
Kamath L, et al. Signaling from protease-activated receptor-1 inhibits migration and invasion of breast cancer cells. Can Res. 2001;61(15):5933.
Ashida K, et al. Expression of e-cadherin, α-catenin, β-catenin, and CD44 (standard and variant isoforms) in human cholangiocarcinoma: an immunohistochemical study. Hepatology. 2010;27(4):974.
Bukholm IK, et al. E-cadherin and α-, β-, and γ-catenin protein expression in relation to metastasis in human breast carcinoma. J Pathol. 2015;185(3):262–6.
Jr LW, Fu R, Banko M. Growth factor-induced transcription via the serum response element is inhibited by cyclic adenosine 3′,5′-monophosphate in MCF-7 breast cancer cells. Endocrinology. 1997;138(6):2219.
Siegel R, et al. Cancer statistics, 2014. CA Cancer J Clin. 2014;64(1):9–29.
Tang J, et al. A novel biomarker Linc00974 interacting with KRT19 promotes proliferation and metastasis in hepatocellular carcinoma. Cell Death Dis. 2014;5(12):e1549.
Rutnam ZJ, et al. The pseudogene TUSC2P promotes TUSC2 function by binding multiple microRNAs. Nat Commun. 2011;5(1):2914.
Chen J, et al. Competing endogenous RNA network analysis identifies critical genes among the different breast cancer subtypes. Oncotarget. 2017;8(6):10171–84.
Zhou X, Liu J, Wang W. Construction and investigation of breast-cancer-specific ceRNA network based on the mRNA and miRNA expression data. IET Syst Biol. 2014;8(3):96–103.
Zhou S, et al. Systematical analysis of lncRNA–mRNA competing endogenous RNA network in breast cancer subtypes. Breast Cancer Res Treat. 2018;169(2):267–75.
Li T, et al. Upregulation of long noncoding RNA ZEB1-AS1 promotes tumor metastasis and predicts poor prognosis in hepatocellular carcinoma. Oncogene. 2016;35(12):1575.
Yuan SX, et al. Antisense long non-coding RNA PCNA-AS1 promotes tumor growth by regulating proliferating cell nuclear antigen in hepatocellular carcinoma. Cancer Lett. 2014;349(1):87–94.
Li Z, et al. Comprehensive analysis of differential co-expression patterns reveal transcriptional dysregulation mechanism and identify novel prognostic lncRNAs in esophageal squamous cell carcinoma. Oncotargets Ther. 2017;10:3095–105.
Xing Y, et al. Comprehensive analysis of differential expression profiles of mRNAs and lncRNAs and identification of a 14-lncRNA prognostic signature for patients with colon adenocarcinoma. Oncol Rep. 2018;39(5):2365–75.
Wang H, et al. LncRNAs expression profiling in normal ovary, benign ovarian cyst and malignant epithelial ovarian cancer. Sci Rep. 2016;6:38983.
Guttilla IK, White BA. Coordinate regulation of FOXO1 by miR-27a, miR-96, and miR-182 in breast cancer cells. J Biol Chem. 2009;284(35):23204–16.
Li P, et al. MiR-183/-96/-182 cluster is up-regulated in most breast cancers and increases cell proliferation and migration. Breast Cancer Res. 2014;16(6):473.
Yan LX, et al. MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA. 2008;14(11):2348.
Zhou X, et al. Decreased miR-204 in H. pylori-associated gastric cancer promotes cancer cell proliferation and invasion by targeting SOX4. PLoS ONE. 2014;9(7):e101457.
Butrym A, et al. Low expression of microRNA-204 (miR-204) is associated with poor clinical outcome of acute myeloid leukemia (AML) patients. J Exp Clin Cancer Res. 2015;34(1):1–5.
Iwaya T, et al. Downregulation of miR-144 is associated with colorectal cancer progression via activation of mTOR signaling pathway. Carcinogenesis. 2012;33(12):2391–7.
CNF and LM conceived and designed the experiments. CNF, LM and NL analyzed data. LM and NL wrote this manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
miRNAs targeting lncRNAs and mRNAs ofBC, as well as the prognostic value of the lncRNAs obtained from the univariate Cox’s analysis.
Gene Ontology (GO) analysis. GO results for aberrantly expressed mRNAs with significant Enrichment score covering domains of cellular components (CCs) and molecular functions (MFs). The bar plot devotes the enrichment scores of the significant GO terms.
Four-lncRNA signature ((ADAMTS9-AS1, AL391421.1, LINC00491, and LINC00536)) predicted OS in BC cohort. A. Risk-score distribution. Red demonstrating higher expression while blue representing lower expression. Risk scores for all BC patients were created in ascending order and blue is marked as low risk or red is labeled as high risk. B. Patients’ survival status with blue devoting dead, and red standing for alive. C. Heat map of the four-lncRNA expression profiles in BC patients.
About this article
Cite this article
Fan, CN., Ma, L. & Liu, N. 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 (2018). https://doi.org/10.1186/s12967-018-1640-2
- Breast cancer
- Long non-coding RNA
- Competing endogenous RNA network
- Overall survival