Skip to main content

A four-long non-coding RNA signature in predicting breast cancer survival



Many long non-coding RNAs(lncRNAs) have been found to be a good marker for several tumors. Using lncRNA-mining approach, we aimed to identify lncRNA expression signature that can predict breast cancer patient survival.


We performed LncRNA expression profiling in 887 breast cancer patients from Gene Expression Omnibus (GEO) datasets. The association between lncRNA signature and clinical survival was analyzed using the training set(n = 327, from GSE 20685). The validation for the association was performed in another three independent testing sets(252 from GSE21653, 204 from GSE12276, and 104 from GSE42568).


A set of four lncRNA genes (U79277, AK024118, BC040204, AK000974) have been identified by the random survival forest algorithm. Using a risk score based on the expression signature of these lncRNAs, we separated the patients into low-risk and high-risk groups with significantly different survival times in the training set. This signature was validated in the other three cohorts. Further study revealed that the four-lncRNA expression signature was independent of age and subtype. Gene Set Enrichment Analysis (GSEA) suggested that gene sets were involved in several cancer metastasis related pathways.


These findings indicate that lncRNAs may be implicated in breast cancer pathogenesis. The four-lncRNA signature may have clinical implications in the selection of high-risk patients for adjuvant therapy.


More than 50% of transcripts have no protein coding potential through the analysis of mammalian transcriptomes, a subset of these noncoding transcripts are termed long non-coding RNAs (lncRNAs) that range from 200 nucleotides to multiple kilobases in length [1]. These long, polyadenylated RNAs do not code for proteins, but function directly as RNAs. Many lncRNAs have already been associated with various disease processes, and cancer features prominently among these. In addition to the classic protein coding mRNAs, recent studies have revealed the contribution of lncRNAs as protooncogenes, tumor suppressor genes, and drivers of metastatic transformation at the transcriptional, post-transcriptional, and epigenetic levels [2]-[5]. Accumulating evidence indicates that lncRNAs are linked to a diverse range of functions in cellular development and their misregulation has also been implicated in various types of cancers [6]-[8]. In most cases, these transcripts are aberrantly expressed in cancers, which may indicate their potential as possible biomarkers and can be predictive of clinical outcome.

Breast cancer is a heterogeneous disease composed of multiple molecular alterations. Molecular differences between histologically similar tumors make clinical outcomes difficult to predict and treatment imperfectly adapted [9]. Breast cancers of varying histological subtypes and risk stratification are traditionally diagnosed based on their histopathological features, including tumor size, grade and lymph node status. Over the past decade, the “intrinsic” molecular subtypes of breast cancer: luminal A and B, basal, ERBB2 and normal-like, exhibit different histo-clinical features and treatment sensitivity [10],[11]. Given the heterogeneity of breast cancer and the multitude of variables influencing clinical evolution, the multi-gene signatures provide further prognostic and predictive information. One of the examples is a 21-gene classifer (Oncotype DX) ,which classifies breast tumors into low-,intermediate- and high-risk groups as to the advisability of adjuvant chemotherapy for patients in high-risk group [12],[13]. The utility of such gene signature might have clinical potential to predict patient outcome and aid in treatment choice [14].

In breast cancer, several lncRNA transcripts were involved in the biology of tumorigenesis. Furthermore, certain lncRNAs exhibit distinct expression patterns between primary tumors and metastases. A 2.2 kb lncRNA, HOTAIR has been shown to be an independent predictor of breast cancer survival. Elevated HOTAIR expression levels correlate with breast cancer, and are linked to poor prognosis and metastasis [3]. This lncRNA may induce metastases by remodeling the epigenetic machinery to repress metastasis suppressor genes (e.g., HOXD10 ). Another lncRNA, MALAT-1 (metastasis associated lung adenocarcinoma transcript 1) is overexpressed in many different cancer types, including lung, breast, colon, prostate, pancreatic, and hepatocellular carcinomas [15]-[17]. This highly conserved 8kb lncRNA is upregulated in invasive breast carcinomas and correlates with tumor grade [18]. GAS5 (growth arrest-specific 5) was found to be downregulated in breast cancer tissues, and overexpression of this lncRNA in the MCF-7 breast cancer cell line furthered growth arrest and apoptosis [19]. LSINCT5, the stress-regulated lncRNA, is overexpressed in breast and ovarian cancer cell lines and tumor tissues. In addition, LSINCT5 has been proved to play a role in cellular proliferation and also in the development of breast and ovarian cancers [20]. Transcriptional profiling has revealed highly aberrant lncRNA expression in human cancers [21]. However, the prognostic significance of lncRNAs in breast cancer has not been investigated.

Recently, the methodology of repurposing microarray data for probing lncRNA expression is well-established [22]-[24]. For instance, Du et al. used a large dataset of microarrays to build a resource of clinically relevant lncRNAs for the development of lncRNA biomarkers and the identification of lncRNA therapeutic targets [25]. Zhang et al. correlates lncRNA expression profiles with malignancy grade and histological differentiation in human gliomas by re-annotating Affymetrix HG-U133 Plus 2.0 array [26],[27]. Furthermore, several studies do have discovered new biomarkers to predict survival by re-annotation of previous microarray data. A six-lncRNA signature has been identified to predict survival of patients with glioblastoma multiforme [27], while a three-lncRNA signature has been shown to be associated with the prognosis of patients with oesophageal squamous cell carcinoma [28].

In this study, we aimed at profiling the lncRNA expression signatures by analyzing a cohort of previously published breast cancer gene expression profiles from the Gene Expression Omnibus (GEO), as well as another three independent data sets as testing sets. We identified a four-lncRNA signature associated with survival, and then established a risk score formula using the expressions of these four lncRNAs. The prognostic value of the signature was further confirmed in the testing cohorts. Our findings suggest that lncRNA signatures can be predictive of clinical outcome and they may be useful as biomarkers.

Materials and methods

GEO breast cancer gene expression data

Breast cancer gene expression data and corresponding clinical data used in this study were obtained from the publicly available GEO databases. To analyze the correlation of lncRNA expression signatures with survival endpoints for breast cancer as a whole (disease-free survival, metastasis-free survival and overall survival), we selected those data sets that included more than 100 patients with their survival status information. We followed the strategy of using the largest data set (GSE20685) as training set. This training set from GSE20685 [29] was first used to identify the gene expression signature. Another three independent data sets from GSE12276 [30], GSE21653 [31],[32], GSE42568 [33] were included in this study as testing sets. After filtering out samples without clinical survival information, there were a total of 887 samples, including 327 from GSE20685, 252 from GSE21653, 204 from GSE12276, and 104 from GSE42568, respectively. Figure 1 depicts the diagram of the study.

Figure 1
figure 1

Diagram of the study. The order of analyses to develop the risk score model and validate the efficiency of the gene signature to predict prognostic outcomes.

Microarray data processing and LncRNA profile mining

The raw CEL files were downloaded from GEO database and background adjusted using Robust Multichip Average. GATExplorer [22] was used to process microarrays on a local computer for gene expressions of lncRNAs. This GATExplorer provides a series of R packages, designed to be used with BioConductor tools, that allow to apply in a simple way the probe mapping data included in GATExplorer. A type of files called ncRNA Mapper were also obtained from GATExplorer, which include the probes that do not map to any coding region but that were mapped to a database for non-coding RNA of human and mouse (derived from RNAdb [34]). A customized R scripts was used to perform a microarray expression calculation according to the re-mapping data(file ncrnamapperhgu133plus2cdf_3.0). Each LncRNA should include at least a minimum of 3 probes mapping in the corresponding ncRNAs entity. All of the four lncRNAs were verified online in the ncRNA Expression Database ( [35], which provides gene expression information for thousands of long ncRNAs in human and mouse (Additional file 1: Table S1). We created a risk-score formula according to the expressions of these four lncRNAs for survival prediction. Patients having higher risk scores are expected to have poorer survival outcomes. The risk scores are calculated as follows: Risk score = (-0.35717× expression value of AK024118) + (0.518242 × expression value of U79277) + (-0.48664 × expression value of BC040204) + (-0.48122 × expression value of AK000974). In addition, the coding potential analysis of the lncRNAs was carried out by CNCI to classify protein-coding or noncoding transcripts [36].

Gene set enrichment analysis (GSEA)

GSEA was performed by the JAVA program ( using MSigDB C2 CP: Canonical pathways gene set collection(1320 gene sets available). Gene sets with a false discovery rate(FDR) value <0.05 after performing 1,000 permutations were considered to be significantly enriched [37]. Cytoscape and Enrichment Map were used for visualization of the GSEA results.

Statistical analysis

The association between the lncRNA gene expression and patient’s survival was assessed by univariable Cox regression analysis along with a permutation test using Biometric Research Branch-Array tools package [38] in the training set. With a parametric test(p ≤ 0.001), we identified a set of 30 lncRNA expressions strongly correlated with survival. Considering that a smaller number of genes in the model would make the model more practical, we then performed the random survival forests-variable hunting (RSFVH) algorithm [39]. We followed Kawaguchi [39],[40] for the parameters in the algorithm. In brief, the number of Monte Carlo iterations (nrep) was set as nrep = 100 and value controlling the step size used in the forward process (nstep) was set as nstep = 5. A set of 4 lncRNAs genes have been identified in which expressions were strongly and consistently related to patient survival (Figure 2).

Figure 2
figure 2

Error rate for the data as a function of trees (left) and out-of-bag importance values for predictors (right).

Using these four genes selected fitted in a multivariable Cox regression model, we constructed a formula that would predict survival in the training set. Each patient was then assigned a risk score that is a linear combination of the expression levels of the significant lncRNAs weighted by their respective Cox regression coefficients [41],[42]. According to this risk score, patients in the training set were divided into low-risk and high-risk groups using the median risk score as the cut-off. The Kaplan-Meier method was used to estimate survival time for the other three testing groups. Differences in survival times between the low-risk and high-risk groups in each set were then compared using the two-sided log rank test. Furthermore, we also used Cox multivariate analysis to test whether the risk score was independent of patient age and subtype with the available data. Oncotype DX score was implemented in an R package called genefu, available from the Comprehensive R Archive Network. We used receiver operating characteristic (ROC) curves to compare the sensitivity and specificity of the survival prediction of the lncRNA risk score. Area under the curve (AUC) values were calculated from the ROC curves [43]. All the data were analyzed by R program ( The significance was defined as p values being less than 0.05.


Identification of lncRNA genes from the training set

As summarized in the workflow (Figure 1), all analyses were performed in the training data set (GSE20685) first and then validated in the test data set (GSE21653, GSE12276, GSE42568). The training set(n = 327) was analyzed for the detection of prognostic lncRNA genes. By subjecting the lncRNA expression data derived from the training set to univariable Cox proportional hazards regression analysis using the BRB-Array Tools, we identified a set of 30 lncRNAs that were strongly correlated with patients’ overall survival (P ≤0.001). On the basis of the random survival forests model (see Materials and methods), four genes were selected as the predictors. Table 1 shows a list of these four genes with their obtained variable importance values. As depicted in Figure 2, from the plot we can see that AK024118 has dramatically larger importance value than other predictors. Of these, a positive coefficient of U79277 indicated that its higher level of expression was associated with shorter survival. The negative coefficients of the other genes (AK024118, BC040204, AK000974) indicated that their higher levels of expression were associated with longer survival. All of the four lncRNAs have been verified in the ncRNA Expression Database ( and these four transcripts were classified as ncRNAs in this website [35]. As coding potential analysis is commonly used to classify whether a transcript is of coding potential or not [25], we also used another tool, CNCI, developed by Sun to test those four transcripts [36]. This tool also suggests that all the four transcripts are non-coding transcripts with no coding potential.

Table 1 Four LncRNAs significantly associated with the overall survival in the training-set patients (n = 327)

The association of four-lncRNA signature and patient’s survival in the training set

With the risk score formula (see Materials and methods), we calculated the four-lncRNA expression signature risk score for each patient in the training set. The patients were then ranked according to their risk scores. Using the median risk score as cut-off in the training set, the patients were divided into low-risk(n = 164) and high-risk (n = 163) groups. Patients in the high-risk group had significantly shorter overall survival than those in the low-risk group (log-rank test P < 0.0001) (Figure 3A). Overall survival in the training set was 96.95% at 3 years, 91.89% at 6 years, 88.08% at 9 years and 84.86% at 12 years in the low risk group, versus 83.99%, 70.27%, 62.02% and 51.56% in the high risk group respectively. The correlation of the four-lncRNA risk score with overall survival was significant when it was analyzed as a continuous variable in the univariable Cox regression model.

Figure 3
figure 3

Kaplan-Meier estimates of the survival of Gene Expression Omnibus (GEO) patients using the four-lncRNA signature. The Kaplan-Meier plots were used to visualize the survival probabilities for the low-risk versus high-risk group of GEO patients determined on the basis of the median risk score from the training set patients. (A) Kaplan-Meier curves for GSE20685 training-set patients (n = 327); (B) Kaplan-Meier curves for GSE21653 testing-set patients (n = 252); (C) Kaplan-Meier curves for GSE12276 testing-set patients (n = 204); (D) Kaplan-Meier curves for GSE42568 testing-set patients (n = 104). The tick marks on the Kaplan-Meier curves represent the censored subjects. The differences between the two curves were determined by the two-sided log-rank test. The number of patients at risk was listed below the survival curves.

Validation of the four-lncRNA signature for survival prediction in the testing sets

In order to confirm our findings, we calculated the risk score for the testing sets including GSE21653(n = 252), GSE12276 (n = 204) and GSE42568(n = 104). By using the same cut-off value as the training set, the patients from each testing set were separately classified into low-risk and high-risk groups and subjected to survival comparison. As overall survival information was unavailable in GSE21653 and GSE12276, disease-free survival(DFS) and metastasis-free survival(MFS) was evaluated, respectively. Similar to the findings obtained from the training set, patients in the high-risk group had shorter survival time than patients in the low-risk group (Figures 3B, 3C, 3D). In consistence with the results described above, patient survival in the low-risk group was better than that in the high-risk group throughout the follow-up. In the univariable Cox regression model, the similar correlation of the risk score with overall survival was noted with the high-risk group having a shorter overall survival than the low-risk group. The distribution of patient risk scores (Z-score transformed), survival status and lncRNA values were analyzed independently for the training set (Figure 4). We found that patients with high-risk scores tended to have higher expression of U79277 and lower expression of the remaining genes(AK024118, BC040204, AK000974). Detail survival information of individual lncRNA in each data set and the gene signature in the context of different tumor subtypes was shown in Additional file 2: Figure S1 and Additional file 3: Figure S2, respectively.

Figure 4
figure 4

LncRNA risk score analysis of GEO patients. The distribution of four-lncRNA risk score, patients’ survival status and lncRNA expression signature were analyzed in the training set patients (n = 327). (A) LncRNA risk score distribution; (B) patients’ overall survival status and time; (C) heatmap of the lncRNA expression profiles. Rows represent lncRNAs, and columns represent patients. The black dotted line represents the median lncRNA risk score cutoff dividing patients into low-risk and high-risk groups.

Multivariate regression analysis shows that the four-lncRNA expression signature is independent of age and subtype

We carried out Cox multivariate analysis to ascertain whether the four-lncRNA expression signature was an independent predictor of breast cancer patient’s survival. Four-lncRNA risk score, age (available in GSE20685, GSE21653 and GSE42568) and subtype(only available in GSE21653) were defined as covariates. The effect of risk score, age and subtype on breast cancer patient survival time was further evaluated by multivariate Cox proportional hazard model. The results showed that risk score is an independent predictor of patient survival when adjusted by age or subtype in every cohort (Table 2).

Table 2 Univariable and multivariable Cox regression analyses in each data set

Evaluation of the risk score performance by receiver operating characteristic (ROC) curve analysis

As GSE21653 was the only data set with disease-free survival information(DFS), we performed receiver operating characteristic (ROC) analysis to compare the sensitivity and specificity of survival prediction between our model and Oncotype DX [13]. The area under receiver operating characteristic (AUROC) was determined and compared between these two gene signatures. As seen in Figure 5, ROC curves indicated that AUROC of four-lncRNA gene signature and Oncotype DX was 0.603 and 0.675, respectively. No significant difference(p=0.0837) was observed between the Oncotype group and the four-lncRNA gene signature group in terms of disease-free survival (DFS) (Figure 3).

Figure 5
figure 5

Receiver operating characteristic (ROC) analysis of sensitivity and specificity by four-lncRNA risk score, oncotype in predicting disease-free survival (DFS). In GSE21653, the score performance was assessed by calculating the area under the ROC (AUROC) of four-lncRNA risk score versus Oncotype, which was 0.603 and 0.675, respectively. (p = 0.0837).

Identification of four-lncRNA signature associated biological pathways and processes

Gene Set Enrichment Analysis(GSEA) was carried out to identify associated biological processes and signaling pathway [37]. We compared the gene expression profile of breast cancer patients with high-risk and low-risk group classified by four-lncRNA gene signature in the training set(GSE 20685). The gene sets with significantly different expression (FDR < 0.01, p <0.005) were picked up for Gene set enrichment analysis (GSEA). Several cancer related pathways such as epithelial mesenchymal transition(EMT) [44], cell cycle signaling and DNA replication [45] were enriched in the high risk group, which implies that the signature might be involved in the metastasis related pathways (Figure 6). The associated biological pathway with each lncRNA was shown in Additional file 4: Figure S3.

Figure 6
figure 6

Gene set enrichment analysis delineates biological pathways and processes associated with risk score in the training set(GSE 20685). GSEA validated enhanced activity of (A) epithelial mesenchymal transition pathway (B) cell cycle signaling pathway (C) DNA replication pathway in high risk score group.


The discovery of multiple functional regulatory lncRNAs has lead to genome-wide searches in multiple species as well as for transcripts that are aberrantly expressed in various types of cancers. Similar to protein-coding genes and miRNAs, lncRNAs play key roles in tumorigenesis. They are involved in a number of fundamental processes associated with cancer including cell cycle regulation, apoptosis, the DNA damage response, and metastasis [3],[46]. The expression of highly conserved lncRNAs is also altered in breast cancers [17]. Our recent study achieved the lncRNA profiling by mining the existing microarray gene expression data as reported [47],[48]. Except for several recent researches on the roles of lncRNAs in breast cancer, the prognostic value of lncRNA signatures have not been investigated. To our knowledge, this is the first report of a lncRNA expression signature predicting breast cancer patient survival.

In this study, we have identified a four-lncRNA expression signature that is associated with survival of breast cancer patients. We further revealed that the four-lncRNA signature is an independent predictor of breast cancer patient survival.

As for the characteristics of the four genes, the overexpression of U79277 was found to be correlated with shorter survival while three of these four lncRNAs identified (AK024118, BC040204, AK000974) were downregulated in the high-risk group compared to low-risk group. The functional study in cancer of these genes has not been reported so far. Nevertheless, our present study demonstrated the associations between the expressions of these genes and survival time. Interestingly, the locations of those putative lncRNAs overlap with many transcripts including some well-known oncogene and tumor suppress genes. AK024118 is located within the intron of BCL2 which is a known driver of lymphoma. U79277 is transcribed from the minus strand on human chromosome 8 and overlaps with YWHAZ 3’UTR. AK000974 which overlaps with many transcripts including CCNJ mRNA. We found that it is very common for ncRNAs. The lncRNAs were categorized as intergenic or genic. The genic lncRNAs were further classified as being exonic, intronic or overlapping and sense or antisense according to their relation with neighboring protein coding genes [49],[50]. Although some of lncRNAs may overlap with neighboring protein coding genes, most of them have their own function. Some lncRNAs regulate the transcription of nearby genes in cis, while others act in trans [28]. A concrete example is HOTAIR, a well-studied lncRNA, within the HOXC cluster was shown to help silence HOXD cluster genes in trans [51]. It may be worthwhile to further investigate these lncRNAs for the purpose of better understanding of their roles in determining breast cancer prognosis.

The median risk score was used as a cutoff point for two reasons. First, a previous lncRNA risk score formula used the median as a cutoff point for classifying patients into two groups [27]. Second, the most common approach for dichotomizing continuous variables was to take the sample median due to the absence of a prior cutpoint [52].

By performing multivariable Cox regression analysis that included age and subtype(when available) as covariables, we analyzed whether the prognostic value of the four-lncRNA signature was independent of age and subtype. The age at diagnosis exercises a complex influence on breast cancer prognosis. Young age at diagnosis influences negatively the prognosis [53]-[55], whereas breast cancer in elderly women is associated with an inferior prognosis when compared to that of middle-aged women [54]. Observational data in breast cancer patients is suggestive of an increased risk of disease specific mortality with increasing age [56],[57]. Several observations suggest that the percentage of deaths attributed to breast cancer decreased with age [58],[59]. These inconsistency in findings could explain the results that age was not significant prognostic factor when assessed in the univariable Cox regression analysis in our study. Nonetheless, we could conclude that the risk score obtained by the four-lncRNA signature was independent of age in the present study.

Breast cancer is clinically heteregeneous due to molecular differences between histologically similar tumors. Luminal, Her2 enriched, basal-like (Triple-negative) subgroups were identified and were shown to have different long-term survivals [10],[11],[60],[61]. There were few reports about the correlation between lncRNAs and molecular subtype of breast cancer. A newly identified lncRNA, LOC554202, has been found to express abundantly in the non-invasive breast cancer cell lines like luminal subtype, but the expression is lost in more aggressive triple-negative breast cancer cell lines of basal subtype [62]. It was therefore of interest to determine if our four-lncRNA signature was associated with this strong prognostic factor. As the data on molecular subtype was only available in GSE21653, we performed the analysis of multivariable Cox regression including risk score and subtype in this testing group. Because of the small sample size in some subgroups, we did not observe significant difference in either univariable or multivariable Cox regression analysis.

Further ROC analysis demonstrated that four-lncRNA gene signature was comparable with Oncotype DX (p = 0.0837). Although Oncotype DX is the most accepted in clinical practice for decision making as to the advisability of adjuvant chemotherapy for breast cancer patients [12], the test is not financially feasible for every patient in developing countries. As shown in this study, a small number of genes (4 genes) could be sufficient to predict the prognostic, using simply reverse transcription polymerase chain reaction (RT-PCR). Clinically, risk score may provide clues on biological behaviors as well as prognostic characteristics of tumors. Patients belonging to high-risk group may need more effective adjuvant therapy in addition to the standard treatment protocol. In addition to the current prognostic model, the four-lncRNA signature may develop easy-to-use prognostic model in order to facilitate further stratification of patients.

Moreover, Gene set enrichment analysis (GSEA) was performed aiming at analyzing coordinate expression changes at a pathway level. The associated molecular pathways, namely, epithelial mesenchymal transition(EMT) [44], cell cycle signaling and DNA replication revealed the four-lncRNA signature might be involved with cancer metastasis. Hence, these findings are likely to be implicated in the development of new targeted anti-cancer therapies. In breast cancer, it has been shown that knock-down of lncRNA HOTAIR with specific siRNAs may limit the metastatic potential of breast cancer cells [3]. The therapeutic potential of targeting regulatory lncRNAs in order to increase the expression of specific genes has also recently emerged [1],[63]. The four prognostic lncRNAs may have therapeutic potential as novel molecular targets.

Several limitations to this study need to be acknowledged. First, in our study, only a fraction of human lncRNA (5635 out of 15000+) were included in the analysis. So, the prognostic lncRNA genes identified here may not represent all the lncRNA candidates that are potentially correlated with breast cancer overall survival. Secondly, we lack information on the mechanisms behind the prognostic values of these four lncRNAs in breast cancer, and experimental studies on these lncRNAs might provide important information to further our understanding of their functional roles. Finally, although we recapitulated our findings in three published datasets to the extent possible based on data availability, the signature has not yet been tested prospectively in a clinical trial. Despite these drawbacks, however, the significant and consistent correlation of our four-lncRNA signature with overall survival in several independent data sets indicates that it is a potentially powerful prognostic marker for breast cancer.


In summary, we have identified a set of four-lncRNA signature, which predicts the overall survival in three independent cohorts. Further analysis revealed that the prognostic value was independent of age and subtype. Clinically, the identification of poor or good prognosis cases may help select the appropriate treatment. The identification of the prognostic lncRNAs indicates the potential roles of lncRNAs in breast cancer pathogenesis. The four-lncRNA signature may have clinical implications as molecular diagnosis markers and therapeutic targets.

Additional files



Long non-coding RNAs


Gene expression omnibus


Gene set enrichment analysis


False discovery rate


Random survival forests-variable hunting


Receiver operating characteristic


Area under the curve


Disease-free survival


Gene set enrichment analysis


Epithelial mesenchymal transition


Reverse transcription polymerase chain reaction


  1. Wahlestedt C: Targeting long non-coding RNA to therapeutically upregulate gene expression. Nat Rev Drug Discov. 2013, 12: 433-446. 10.1038/nrd4018.

    Article  CAS  PubMed  Google Scholar 

  2. Beltran M, Puig I, Pena C, Garcia JM, Alvarez AB, Pena R, Bonilla F, de Herreros AG: A natural antisense transcript regulates Zeb2/Sip1 gene expression during Snail1-induced epithelial-mesenchymal transition. Genes Dev. 2008, 22: 756-769. 10.1101/gad.455708.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, Wang Y, Brzoska P, Kong B, Li R, West RB, van de Vijver MJ, Sukumar S, Chang HY: Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010, 464: 1071-1076. 10.1038/nature08975.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E, Chang HY: Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007, 129: 1311-1323. 10.1016/j.cell.2007.05.022.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Viny AD, Clemente MJ, Jasek M, Askar M, Ishwaran H, Nowacki A, Zhang A, Maciejewski JP: MICA polymorphism identified by whole genome array associated with NKG2D-mediated cytotoxicity in T-cell large granular lymphocyte leukemia. Haematologica. 2010, 95: 1713-1721. 10.3324/haematol.2010.021865.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Huarte M, Rinn JL: Large non-coding RNAs: missing links in cancer?. Hum Mol Genet. 2010, 19: R152-R161. 10.1093/hmg/ddq353.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Wapinski O, Chang HY: Long noncoding RNAs and human disease. Trends Cell Biol. 2011, 21: 354-361. 10.1016/j.tcb.2011.04.001.

    Article  CAS  PubMed  Google Scholar 

  8. Li L, Sun R, Liang Y, Pan X, Li Z, Bai P, Zeng X, Zhang D, Zhang L, Gao L: Association between polymorphisms in long non-coding RNA PRNCR1 in 8q24 and risk of colorectal cancer. J Exper Clin Can Res. 2013, 32: 104-10.1186/1756-9966-32-104.

    Article  CAS  Google Scholar 

  9. Cancer Genome Atlas N: Comprehensive molecular portraits of human breast tumours. Nature 2012, 490:61–70. 10.1038/nature11412.

  10. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lønning PE, Børresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.

    Article  CAS  PubMed  Google Scholar 

  11. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Lønning PE, Børresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001, 98: 10869-10874. 10.1073/pnas.191367098.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Paik S: Development and clinical utility of a 21-gene recurrence score prognostic assay in patients with early breast cancer treated with tamoxifen. Oncologist. 2007, 12: 631-635. 10.1634/theoncologist.12-6-631.

    Article  CAS  PubMed  Google Scholar 

  13. Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, Baehner FL, Walker MG, Watson D, Park T, Hiller W, Fisher ER, Wickerham DL, Bryant J, Wolmark N: A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004, 351: 2817-2826. 10.1056/NEJMoa041588.

    Article  CAS  PubMed  Google Scholar 

  14. Yin ZQ, Liu JJ, Xu YC, Yu J, Ding GH, Yang F, Tang L, Liu BH, Ma Y, Xia YW, Lin XL, Wang HX: A 41-gene signature derived from breast cancer stem cells as a predictor of survival. J Exper Clin Can Res. 2014, 33: 49-10.1186/1756-9966-33-49.

    Article  CAS  Google Scholar 

  15. Ji P, Diederichs S, Wang W, Boing S, Metzger R, Schneider PM, Tidow N, Brandt B, Buerger H, Bulk E, Thomas M, Berdel WE, Serve H, Müller-Tidow C: MALAT-1, a novel noncoding RNA, and thymosin beta4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 2003, 22: 8031-8041. 10.1038/sj.onc.1206928.

    Article  PubMed  Google Scholar 

  16. Schmidt LH, Spieker T, Koschmieder S, Schaffers S, Humberg J, Jungen D, Bulk E, Hascher A, Wittmer D, Marra A, Hillejan L, Wiebe K, Berdel WE, Wiewrodt R, Muller-Tidow C: The long noncoding MALAT-1 RNA indicates a poor prognosis in non-small cell lung cancer and induces migration and tumor growth. J Thorac Oncol. 2011, 6: 1984-1992. 10.1097/JTO.0b013e3182307eac.

    Article  PubMed  Google Scholar 

  17. Perez DS, Hoage TR, Pritchett JR, Ducharme-Smith AL, Halling ML, Ganapathiraju SC, Streng PS, Smith DI: Long, abundantly expressed non-coding transcripts are altered in cancer. Hum Mol Genet. 2008, 17: 642-655. 10.1093/hmg/ddm336.

    Article  CAS  PubMed  Google Scholar 

  18. Iacoangeli A, Lin Y, Morley EJ, Muslimov IA, Bianchi R, Reilly J, Weedon J, Diallo R, Bocker W, Tiedge H: BC200 RNA in invasive and preinvasive breast cancer. Carcinogenesis. 2004, 25: 2125-2133. 10.1093/carcin/bgh228.

    Article  CAS  PubMed  Google Scholar 

  19. Mourtada-Maarabouni M, Pickard MR, Hedge VL, Farzaneh F, Williams GT: GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene. 2009, 28: 195-208. 10.1038/onc.2008.373.

    Article  CAS  PubMed  Google Scholar 

  20. Silva JM, Boczek NJ, Berres MW, Ma X, Smith DI: LSINCT5 is over expressed in breast and ovarian cancer and affects cellular proliferation. RNA Biol. 2011, 8: 496-505. 10.4161/rna.8.3.14800.

    Article  CAS  PubMed  Google Scholar 

  21. Gibb EA, Vucic EA, Enfield KS, Stewart GL, Lonergan KM, Kennett JY, Becker-Santos DD, MacAulay CE, Lam S, Brown CJ, Lam WL: Human cancer long non-coding RNA transcriptomes. PLoS One. 2011, 6: e25915-10.1371/journal.pone.0025915.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Risueno A, Fontanillo C, Dinger ME, de Las Rivas J: GATExplorer: genomic and transcriptomic explorer; mapping expression probes to gene loci, transcripts, exons and ncRNAs. BMC Bioinformatics. 2010, 11: 221-10.1186/1471-2105-11-221.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Liao Q, Xiao H, Bu D, Xie C, Miao R, Luo H, Zhao G, Yu K, Zhao H, Skogerbo G, Chen R, Wu Z, Liu C, Zhao Y: ncFANs: a web server for functional annotation of long non-coding RNAs. Nucleic Acids Res. 2011, 39: W118-W124. 10.1093/nar/gkr432.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Gellert P, Ponomareva Y, Braun T, Uchida S: Noncoder: a web interface for exon array-based detection of long non-coding RNAs. Nucleic Acids Res. 2013, 41: e20-10.1093/nar/gks877.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Du Z, Fei T, Verhaak RG, Su Z, Zhang Y, Brown M, Chen Y, Liu XS: Integrative genomic analyses reveal clinically relevant long noncoding RNAs in human cancer. Nat Struct Mol Biol. 2013, 20: 908-913. 10.1038/nsmb.2591.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Zhang X, Sun S, Pu JK, Tsang AC, Lee D, Man VO, Lui WM, Wong ST, Leung GK: Long non-coding RNA expression profiles predict clinical phenotypes in glioma. Neurobiol Dis. 2012, 48: 1-8. 10.1016/j.nbd.2012.06.004.

    Article  PubMed  Google Scholar 

  27. Zhang XQ, Sun S, Lam KF, Kiang KM, Pu JK, Ho AS, Lui WM, Fung CF, Wong TS, Leung GK: A long non-coding RNA signature in glioblastoma multiforme predicts survival. Neurobiol Dis. 2013, 58: 123-131. 10.1016/j.nbd.2013.05.011.

    Article  CAS  PubMed  Google Scholar 

  28. Li J, Chen Z, Tian L, Zhou C, He MY, Gao Y, Wang S, Zhou F, Shi S, Feng X, Sun N, Liu Z, Skogerboe G, Dong J, Yao R, Zhao Y, Sun J, Zhang B, Yu Y, Shi X, Luo M, Shao K, Li N, Qiu B, Tan F, Chen R, He J: LncRNA profile study reveals a three-lncRNA signature associated with the survival of patients with oesophageal squamous cell carcinoma.Gut 2014, doi:10.1136/gutjnl-2013-305806. [Epub ahead of print].,

  29. Kao KJ, Chang KM, Hsu HC, Huang AT: Correlation of microarray-based breast cancer molecular subtypes and clinical outcomes: implications for treatment optimization. BMC Cancer. 2011, 11: 143-10.1186/1471-2407-11-143.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Bos PD, Zhang XH, Nadal C, Shu W, Gomis RR, Nguyen DX, Minn AJ, van de Vijver MJ, Gerald WL, Foekens JA, Massague J: Genes that mediate breast cancer metastasis to the brain. Nature. 2009, 459: 1005-1009. 10.1038/nature08021.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Sabatier R, Finetti P, Cervera N, Lambaudie E, Esterni B, Mamessier E, Tallet A, Chabannon C, Extra JM, Jacquemier J, Viens P, Birnbaum D, Bertucci F, Viens P, Birnbaum D, Bertucci F: A gene expression signature identifies two prognostic subgroups of basal breast cancer. Breast Cancer Res Treat. 2011, 126: 407-420. 10.1007/s10549-010-0897-9.

    Article  CAS  PubMed  Google Scholar 

  32. Sabatier R, Finetti P, Adelaide J, Guille A, Borg JP, Chaffanet M, Lane L, Birnbaum D, Bertucci F: Down-regulation of ECRG4, a candidate tumor suppressor gene, in human breast cancer. PLoS One. 2011, 6: e27656-10.1371/journal.pone.0027656.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Clarke C, Madden SF, Doolan P, Aherne ST, Joyce H, O’Driscoll L, Gallagher WM, Hennessy BT, Moriarty M, Crown J, Kennedy S, Clynes M: Correlating transcriptional networks to breast cancer survival: a large-scale coexpression analysis. Carcinogenesis. 2013, 34 (10): 2300-2308. 10.1093/carcin/bgt208.

    Article  CAS  PubMed  Google Scholar 

  34. Pang KC, Stephen S, Engstrom PG, Tajul-Arifin K, Chen W, Wahlestedt C, Lenhard B, Hayashizaki Y, Mattick JS: RNAdb-a comprehensive mammalian noncoding RNA database. Nucleic Acids Res. 2005, 33: D125-D130. 10.1093/nar/gki089.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Dinger ME, Pang KC, Mercer TR, Crowe ML, Grimmond SM, Mattick JS: NRED: a database of long noncoding RNA expression. Nucleic Acids Res. 2009, 37: D122-D126. 10.1093/nar/gkn617.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, Liu Y, Chen R, Zhao Y: Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41: e166-10.1093/nar/gkt646.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y: Analysis of gene expression data using BRB-ArrayTools. Cancer Inform. 2007, 3: 11-17.

    PubMed Central  PubMed  Google Scholar 

  39. Ishwaran H, Kogalur UB: Consistency of random survival forests. Stat Probabil Lett. 2010, 80: 1056-1064. 10.1016/j.spl.2010.02.020.

    Article  Google Scholar 

  40. Kawaguchi A, Iwadate Y, Komohara Y, Sano M, Kajiwara K, Yajima N, Tsuchiya N, Homma J, Aoki H, Kobayashi T, Sakai Y, Hondoh H, Fujii Y, Kakuma T, Yamanaka R: Gene expression signature-based prognostic risk score in patients with primary central nervous system lymphoma. Clin Cancer Res. 2012, 18: 5672-5681. 10.1158/1078-0432.CCR-12-0596.

    Article  CAS  PubMed  Google Scholar 

  41. Alizadeh AA, Gentles AJ, Alencar AJ, Liu CL, Kohrt HE, Houot R, Goldstein MJ, Zhao S, Natkunam Y, Advani RH, Gascoyne RD, Briones J, Tibshirani RJ, Myklebust JH, Plevritis SK, Lossos IS, Levy R: Prediction of survival in diffuse large B-cell lymphoma based on the expression of 2 genes reflecting tumor and microenvironment. Blood. 2011, 118: 1350-1358. 10.1182/blood-2011-03-345272.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Bralten LB, French PJ: Genetic alterations in glioma. Cancers. 2011, 3: 1129-1140. 10.3390/cancers3011129.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Kang J, D’Andrea AD, Kozono D: A DNA repair pathway-focused score for prediction of outcomes in ovarian cancer treated with platinum-based chemotherapy. J Natl Cancer Inst. 2012, 104: 670-681. 10.1093/jnci/djs177.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Kalluri R, Weinberg RA: The basics of epithelial-mesenchymal transition. J Clin Invest. 2009, 119: 1420-1428. 10.1172/JCI39104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Vogelstein B, Kinzler KW: Cancer genes and the pathways they control. Nat Med. 2004, 10: 789-799. 10.1038/nm1087.

    Article  CAS  PubMed  Google Scholar 

  46. Cui Z, Ren S, Lu J, Wang F, Xu W, Sun Y, Wei M, Chen J, Gao X, Xu C, Mao JH, Sun Y: The prostate cancer-up-regulated long noncoding RNA PlncRNA-1 modulates apoptosis and proliferation through reciprocal regulation of androgen receptor. Urol Oncol. 2013, 31 (7): 1117-1123. 10.1016/j.urolonc.2011.11.030.

    Article  CAS  PubMed  Google Scholar 

  47. Johnson R: Long non-coding RNAs in Huntington’s disease neurodegeneration. Neurobiol Dis. 2012, 46: 245-254. 10.1016/j.nbd.2011.12.006.

    Article  CAS  PubMed  Google Scholar 

  48. Michelhaugh SK, Lipovich L, Blythe J, Jia H, Kapatos G, Bannon MJ: Mining Affymetrix microarray data for long non-coding RNAs: altered expression in the nucleus accumbens of heroin abusers. J Neurochem. 2011, 116: 459-466. 10.1111/j.1471-4159.2010.07126.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, Lagarde J, Veeravalli L, Ruan X, Ruan Y, Lassmann T, Carninci P, Brown JB, Lipovich L, Gonzalez JM, Thomas M, Davis CA, Shiekhattar R, Gingeras TR, Hubbard TJ, Notredame C, Harrow J, Guigó R: The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012, 22: 1775-1789. 10.1101/gr.132159.111.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Fachel AA, Tahira AC, Vilella-Arias SA, Maracaja-Coutinho V, Gimba ER, Vignal GM, Campos FS, Reis EM, Verjovski-Almeida S: Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations. Mol Cancer. 2013, 12: 140-10.1186/1476-4598-12-140.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Schorderet P, Duboule D: Structural and functional differences in the long non-coding RNA hotair in mouse and human. PLoS Genet. 2011, 7: e1002071-10.1371/journal.pgen.1002071.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Altman DG, Royston P: The cost of dichotomising continuous variables. BMJ. 2006, 332: 1080-10.1136/bmj.332.7549.1080.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Winchester DP, Osteen RT, Menck HR: The National Cancer Data Base report on breast carcinoma characteristics and outcome in relation to age. Cancer. 1996, 78: 1838-1843. 10.1002/(SICI)1097-0142(19961015)78:8<1838::AID-CNCR27>3.0.CO;2-Y.

    Article  CAS  PubMed  Google Scholar 

  54. Adami HO, Malker B, Holmberg L, Persson I, Stone B: The relation between survival and age at diagnosis in breast cancer. N Engl J Med. 1986, 315: 559-563. 10.1056/NEJM198608283150906.

    Article  CAS  PubMed  Google Scholar 

  55. Fredholm H, Eaker S, Frisell J, Holmberg L, Fredriksson I, Lindman H: Breast cancer in young women: poor survival despite intensive treatment. PLoS One. 2009, 4: e7695-10.1371/journal.pone.0007695.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Bastiaannet E, Liefers GJ, de Craen AJ, Kuppen PJ, van de Water W, Portielje JE, van der Geest LG, Janssen-Heijnen ML, Dekkers OM, van de Velde CJ, Westendorp RG: Breast cancer in elderly compared to younger patients in the Netherlands: stage at diagnosis, treatment and survival in 127,805 unselected patients. Breast Cancer Res Treat. 2010, 124: 801-807. 10.1007/s10549-010-0898-8.

    Article  CAS  PubMed  Google Scholar 

  57. Chapman JA, Meng D, Shepherd L, Parulekar W, Ingle JN, Muss HB, Palmer M, Yu C, Goss PE: Competing causes of death from a randomized trial of extended adjuvant endocrine therapy for breast cancer. J Natl Cancer Inst. 2008, 100: 252-260. 10.1093/jnci/djn014.

    Article  PubMed Central  PubMed  Google Scholar 

  58. Schairer C, Mink PJ, Carroll L, Devesa SS: Probabilities of death from breast cancer and other causes among female breast cancer patients. J Natl Cancer Inst. 2004, 96: 1311-1321. 10.1093/jnci/djh253.

    Article  PubMed  Google Scholar 

  59. Du XL, Fox EE, Lai D: Competing causes of death for women with breast cancer and change over time from 1975 to 2003. Am J Clin Oncol. 2008, 31: 105-116. 10.1097/COC.0b013e318142c865.

    Article  PubMed Central  PubMed  Google Scholar 

  60. Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, Deng S, Johnsen H, Pesich R, Geisler S, Demeter J, Perou CM, Lønning PE, Brown PO, Børresen-Dale AL, Botstein D: Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003, 100: 8418-8423. 10.1073/pnas.0932692100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Sotiriou C, Neo SY, McShane LM, Korn EL, Long PM, Jazaeri A, Martiat P, Fox SB, Harris AL, Liu ET: Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc Natl Acad Sci U S A. 2003, 100: 10393-10398. 10.1073/pnas.1732912100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Augoff K, McCue B, Plow EF, Sossey-Alaoui K: miR-31 and its host gene lncRNA LOC554202 are regulated by promoter hypermethylation in triple-negative breast cancer. Mol Cancer. 2012, 11: 5-10.1186/1476-4598-11-5.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Schwartz JC, Younger ST, Nguyen NB, Hardy DB, Monia BP, Corey DR, Janowski BA: Antisense transcripts are targets for activating small RNAs. Nat Struct Mol Biol. 2008, 15: 842-848. 10.1038/nsmb.1444.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


This article is funded by the research grant (YG2012ZD02, 2JC1407400) from Shanghai Jiaotong University and the Science and Technology Commission of Shanghai, China and the grant (81272506, 61227017) from National Natural Science Foundation of China.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Shen Fu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JM, PL and QZ identified all of the public datasets, JM, PL,QZ and ZY carried out all of the biostatistical and informatics analyses, helped develop the method, JM formulated the study conclusions and draft the manuscript; SF initiated and coordinated the project, guided the study design, supervised all data curation and analysis, finalized all study conclusions and manuscript writing. All coauthors reviewed and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Meng, J., Li, P., Zhang, Q. et al. A four-long non-coding RNA signature in predicting breast cancer survival. J Exp Clin Cancer Res 33, 84 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: