- Research article
- Open Access
- Open Peer Review
Prognostic and therapeutic significance of ribonucleotide reductase small subunit M2 in estrogen-negative breast cancers
BMC Cancervolume 14, Article number: 664 (2014)
Ribonucleotide reductase (RR) is an essential enzyme involved in DNA synthesis. We hypothesized that RR subunit M2 (RRM2) might be a novel prognostic and predictive biomarker for estrogen receptor (ER)-negative breast cancers.
Individual and pooled survival analyses were conducted on six independent large-scale breast cancer microarray data sets; and findings were validated on a human breast tissue set (ZJU set).
Gene set enrichment analysis revealed that RRM2-high breast cancers were significantly enriched for expression of gene sets that increased in proliferation, invasiveness, undifferentiation, embryonic stem/progenitor-like phenotypes, and poor patient survival (p < 0.01). Independent and pooled analyses verified that increased RRM2 mRNA levels were associated with poor patient outcome in a dose-dependent manner. The prognostic power of RRM2 mRNA was comparable to multiple gene signatures, and it was superior to TNM stage. In ER-negative breast cancers, RRM2 showed more prognostic power than that in ER-positive breast cancers. Further analysis indicated that RRM2 was a more accurate prognostic biomarker for ER-negative breast cancers than the pathoclinical indicators and uPA. A new RR inhibitor, COH29, could significantly enhance the chemosensitivity to doxorubicin in ER-negative MDA-MB-231 cells, but not in ER-positive MCF-7 cells.
RRM2 appears to be a promising prognostic biomarker and therapeutic target for ER-negative breast cancer patients.
Over a million new cases of breast cancer are diagnosed and ~400,000 deaths occur annually . Breast cancer is a heterogeneous disease that has variable gene expression and different outcomes that cannot be predicted by pathologic grade or clinical stages [2, 3]. A comprehensive gene expression signature has identified four major molecular subtypes: luminal A, luminal B, HER2-enriched and basal-like breast cancers, which are largely comprised of the triple negative breast cancer (TNBC) subtype. Each subtype has a distinct clinical behavior and response to therapy . Among the four, the TN and HER2-enriched subtypes are considered to be ER-negative, and account for 15-17% and 15-20% of invasive breast carcinomas respectively . Recently, therapies that specifically target the HER2 receptor have significantly increased the survival of patients with HER2-enriched breast cancers, and PARP (poly-ADP ribose polymerase) inhibition holds some promise as a targeted therapy for TNBC [6, 7]. Gene signatures that include the 70-gene signature , 21-gene recurrence score (commercially developed as Oncotype Dx) , PI3K signature , core serum response signature (CSR) , and the grade signature , have been developed to predict the survival of breast cancer patients, and have been used to predict the outcome for ER-positive patients. A 7-gene immune response module (IRM)  and a HER2-derived prognostic predictor (HDPP, 158 genes)  were developed to identify ER-negative breast cancers that were associated with poor prognosis. Multiple-gene-based signatures potentially enhance the accuracy of prediction. However, their disadvantages include higher cost and lack of specific targets for chemotherapeutic agent selection.
Ribonucleotide reductase (RR) is a rate-limiting enzyme essential for DNA replication. RR catalyzes the formation of 2′-deoxyribonucleoside diphosphates from corresponding ribonucleotide diphosphates . Human RR is a heterodimeric tetramer that consists of two large RRM1 subunits and two small subunits, RRM2 and/or RRM2B . RRM2 overexpression induces cell proliferation and invasion [17, 18] and is positively correlated with higher grade and muscularis propria invasion in bladder  and gastric cancers . Previously, we showed that high RRM2 expression is a prognostic indicator for poor survival for several kinds of cancers, including lung and colorectal cancers [21–23]. RRM2 enhances the invasive capacity of pancreatic adenocarcinoma cells in a NF-kappaB-dependent manner , and RRM2 overexpression can increase angiogenesis by down-regulating thrombspondin-1 and increasing production of VEGF . Knocking down RRM2 also sensitizes cancer cells to the cytotoxic effects of the nucleoside analogs . The RRM2 gene (mRNA) has been included in some gene signatures, such as the 3D-culture signature  and Five-gene Molecular Grade Index (MGI) [28, 29], and the prognostic significance of RRM2 mRNA has been evaluated [30, 31]. Independent prognostic significance of RRM2 protein needs to be further investigated, and the prognostic performance of RRM2 (mRNA and protein) also needs to be compared with other existing breast cancer biomarkers or gene signatures. Especially, the prognostic value of RRM2 in ER-negative breast cancers is worth further evaluation.
In this study, we evaluated the prognostic value of RRM2 by analyzing seven independent breast cancer data sets, both individually and in a pooled manner. RRM2 expression was significantly correlated with poor survival in a dose-dependent manner, particularly for ER-negative breast cancers. Moreover, the prognostic value of RRM2 was comparable to the 70-gene signature and 21-gene recurrence score for breast cancer overall, and superior to these biomarkers and pathoclinical indicators for ER-negative breast cancers.
Microarray data sets
A total of 10 published microarray data sets including: Ivshina (GSE4922), Chin (E-TABM-158), Wang (GSE2034), Pawitan (GSE1456), Desmedt (GSE7390), Expo (GSE2109), Huang , Bild (GSE3143), Sortiriou (GSE2990) and NKI  with clinical annotations were downloaded from the combined microarray dataset BRAVO (Biomarker recognition and validation on-line). The Expo set was excluded because it lacks survival data. Based on BLAST results, the data sets of Huang, Bild and Sortiriou were excluded because of using low similarity probes (36922_at) (Additional file 1: Figure S1A). The NKI set was selected for RRM2 prognostic evaluation because the probe (Agilent Technologies) for RRM2 also is 100% similarity to sequence of > gi|260064012|ref|NM_001165931.1|, and it contains information of most gene signatures’ classification. Finally, six independent microarray data sets include Ivshina , Chin , Wang , Pawitan , Desmedt  and NKI  were chosen for this study. The summarized microarray data sets are shown in Additional file 2: Table S1. The significant correlation between signals yielded from two RRM2 probes (201890_at and 209773_s_at) in the Ivshina set was shown in Additional file 1: Figure S1B (R2 = 0.736, p < 0.001). Neither 201890_at nor 209773_s_at correlated with the expression of RRM2B or RRM1 (Additional file 1: Figure S1C to 1E), indicating that the signal was specific for RRM2 expression.
The protocol for the use of human tissues was reviewed and approved by the Institutional Review Board (IRB) of the Second Affiliated Hospital of Zhejiang University (ZJU) (Zhejiang, China). Prior to the study, all patients gave their written informed consent to allow us to use leftover tissue samples for scientific research. All eligible participants had received modified radical mastectomy and the primary tumor samples were obtained from surgical specimens. The exclusion criteria were: 1) no informed consent obtained, 2) multiple cancers, 3) lack of histological diagnosis, and 4) no follow-up information. After applying the selection criteria, a total of 175 breast cancer patients who were diagnosed from 2002 to 2006 were enrolled in the ZJU set (Additional file 2: Table S2). All patients recruited were Chinese females with a median age of 50.0 years (range 20–84 years). The pathoclinical and demographic data was collected by reviewing the hospital records. Of the 175 patients, 46 patients with local advanced breast cancer( Stage IIIA, IIIB and IIIC) received neo-adjuvant chemotherapy before surgery, 140 patients with positive lympho node or aggressive potential(Including: Age <35, Tumor size equal or more than 2 cm, grade 2–3, lymphovascular invasion , and ER/PR negative) were treated with adjuvant chemotherapy (primarily an anthracycline-based regimen) and 9 patients with large tumor size (equal or lager than 5 cm) received post-surgical radiotherapy. After surgery, all patients in the ZJU set were followed up twice a year until September 2010. Median follow-up time was 73.6 months (range: 8.9-104.9) for overall survival (OS) and 71.2 months (range: 2.5-104.9) for progression-free survival (PFS).
This is a population-based, retrospective outcome study (Additional file 1: Figure S2). The PFS was defined as the time from initial surgery until tumor recurrence, including local relapse and metastasis. The OS period was calculated as the time from initial surgery to the date the patient was last seen. Only deaths from metastasis and local relapse of breast cancer were considered as the end of the survival period. According to calculations conducted using nQuery Advisor 6.01 software, 170 cases would be sufficient for 80% study power with a two-sided α of 0.05, so the sample size in the ZJU set is sufficient to reach meaningful conclusions.
To normalize the mRNA expression levels among the above data sets, we re-stratified the scores of RRM2 and other related genes in to four grades (Q1, Q2, Q3 and Q4) based on the percentile for each independently downloaded data set. For Cox analysis, less than the value of the median was regarded as RRM2-low, and greater or equal to the median was RRM2-high.
Gene set enrichment analysis (GSEA)
The detailed GSEA protocol can be obtained from the Broad Institute Gene Set Enrichment Analysis website (http://www.broad.mit.edu/gsea) and related references . The GSEA software v2.0.13 was run using a JAVA 7.0 platform. The dataset (.gct) and phenotype label (.cls) files were created and loaded into GSEA software, and gene sets were downloaded from the Board Institute website. The number of permutation was set to 1000, and the phenotype label was RRM2-high versus RRM2-low. The ranked-list metric was generated by calculating the signal-to-noise ratio, which is based on the difference of means scaled according to the standard deviation.
Quantitative immunohistochemistry (IHC) assays
The specimen of primary breast cancers specimens in the ZJU set was assembled and built into a Multiple-tissue array (MTA) as previously described . Antibodies against RRM2 and RRM2B were previously generated, selected and tested in our lab . Antibodies against ER (Clone: SP1), PR (Clone: PgR636), HER2 (Clone: A0485) and Ki67 (Clone: MIB-1) were purchased from Dako company, and the p53 antibody (Clone: DO-7) was purchased from Vector Lab. Proteins were stained in tissue sections from the MTAs using a deparaffinization and staining protocol previously described . The cut-off values for ER, PR, HER2 or Ki67 positivity were based on previous reports [2, 39]. Based on the above IHC staining, all participants were classified into four intrinsic subtypes Luminal A, Luminal B, HER2 positive and Basal-like TNBC . For breast cancers, the CD44+/CD24-/low was an indicator for the tumor stem/progenitor cells . Double staining of CD44/CD24 was performed using the Autostainer (Dako) and EnVision G|2 Doublestain System for Rabbit/Mouse (DAB+/Permanent Red) (Dako) according to the manufacturer’s instructions. CD44 was detected with diaminobenzidine (DAB) and CD24 was detected with Permanent Red.
RRM2 was located mainly in cytoplasm, but some nuclear staining was also observed. The immunoreactivity in the cytoplasm was scored according to percentage and intensity of the staining. The staining intensity scored as 0: negative, 1: weak, 2: moderate, and 3: strong. Because of heterogeneity of cancer cells, the IHC staining in each cell was varying. In our score system, only more than 10% cells of whole cells with positive staining were taken into consideration. Highest staining staining would be scored final score. Finally, the expression level of the proteins was divided into four subgroups: negative (-), weak positive (+), positive (++), and strong positive (+++). In Cox analysis, (-) and (+) were defined as RRM2-low, whereas (++) and (+++) were defined as RRM2-high.
To avoid observer bias, all the slides were evaluated independently by two different observers with training in pathology in a double-blind manner. Discrepancies were jointly reviewed by the two readers, and missing samples were left blank.
Small interfering RNA, invasion and cytotoxicity assay
RRM2 and scrambled control siRNA were purchased from Santa Cruz Biotechnology Inc. Cells were seeded at 2 × 105 cells per well in 6-well plates filled with 2 mL antibiotic-free normal growth medium supplemented with FBS and then incubated at 37°C in a 5% CO2 incubator for 24 hours. A total of 4 μL of 10 nmol/mL RRM2 or control siRNA was transfected into MDA-MB-231 or ZR-75-1 cells by using the transfection RNAiMAX reagent (Invitrogen, Carlsbad, CA). Cells were incubated in the transfection medium for 6 hours and then placed in normal cell culture medium. The inhibition of RRM2 and RRM2B was measured by using Western blot analysis.
For the invasion assay , 2.5 × 104 cells were seeded on the Matrigel™ (BD company) insert of the 24-well chamber. After incubation for 72 hours in 5% CO2 at 37°C, the top of the Matrigel™ inserts were wiped with a cotton-tipped swab to remove cells that had not migrated through the membrane. The cells on the lower surface of the membrane were stained with 0.5% Coomassie blue (dissolved in 50% ethanol) and counted. Each experiment was performed three times.
Cytotoxicity was assessed using the MTS assay, according to the manufacturer’s instructions (CellTiter 96 AQueous Assay reagent; Promega) on 10 replicates of 2,500 cells per well in a 96-well plate treated with test drugs for 72 hours.
Data management and statistical methods
The database was downloaded, converted, constructed and managed by using MS-Excel. The JMP 10.0 software (SAS Institution, Cary, NC) was used for general statistical analysis. Group comparisons for continuous data were performed using t-tests for independent means or one-way analyses of variance. Categorical variables were compared using χ2 analysis, Fisher’s exact test or the binomial test of proportions. Kaplan-Meier analysis and Cox hazard proportional hazard models were used to analyze OS and PFS. In PFS analysis, patients with distant metastatic breast cancer not completely resected were excluded. Multivariate Cox analysis was applied to adjust for covariate effects, and stratification analysis was used to reduce the confounding effect on estimation of Hazard Ratio (HR). Missing data were coded and excluded from the analysis. Meanwhile the Comprehensive Meta-analysis V2 software was used for meta-analysis. A Fixed-effect model was used to adjusted the HR of RRM2.
Oncofetal characteristics of RRM2in breast cancers
The GSEA was conducted to determine the correlation between RRM2 expression and oncogenic pathways in the Desmedt set. Samples were re-stratified as RRM2-high (equal or greater than median, n = 80) or RRM2-low (less than median, n = 79), according to the expression levels of RRM2 mRNA. Genes that were significantly enriched in RRM2-high breast cancer subgroups controlled cancer undifferentiation, Sarrio’s epithelial-mesenchymal transition, breast cancer progenitor-like phenotypes, proliferation gene signature, Chang’s core serum response gene signature, and the Naderi breast cancer prognosis (up) gene signature (Figure 1A through 1 F). These findings suggested that RRM2 might be associated with poor differentiation and proliferation, and promote the invasiveness of breast cancer.
The association of RRM2 protein levels with clinicopathological features of breast cancer was analyzed in the ZJU set. The Figure 2A shows the IHC staining of RRM2-, +, ++ and +++ samples. Increased RRM2 protein levels were associated with larger tumor size, positive lymph nodes and relapse/metastasis (p < 0.05), but not with patient age or histological type (Additional file 2: Table S2). CD44+/CD24-/low is considered to be the biomarker for the tumor stem/progenitor cell phenotype in breast cancers . Representative CD44+/CD24+ and CD44+/CD24-/low cells are displayed in Figure 2B. The high expression of RRM2 protein was significantly correlated with CD44+/CD24-/low, a tumor stem/progenitor cell indicator (Trends p = 0.044), and Ki-67, a proliferation indicator (p = 0.018) (Additional file 2: Table S2, and Figure 2C and D). In addition, RRM2 levels were significantly correlated with poorly differentiated breast cancers (p < 0.001) (Figure 2E). Moreover, the percentage of RRM2-high in Luminal A, Luminal B/HER2-, Luminal B/HER2+, TNBC and HER2-positive molecular subtypes were 46.0%, 55.9%, 70.0%, 65.4% and 66.7%, respectively (Additional file 2: Table S2). It was indicating that RRM2 protein levels were higher in the molecular subtypes with poor outcomes. The relation between RRM2 mRNA levels and molecular subtypes of breast cancer also could be seen in the Pawitan set (Figure 2F) (p < 0.05).
The prognostic significance of RRM2for breast cancers
Uni- and multivariate Cox analyses of the data were used to evaluate the prognostic value of RRM2 for breast cancer patient outcome. The analysis of OS and PFS revealed that high levels of RRM2 were significantly associated with poor survival in the Desmedt, Pawitan, Wang, Ivshina and NKI sets, except in Chin set. The adjusted HR values for OS were 2.21(95% CI 1.16-4.24), 0.76(95% CI 0.38-1.49), 2.48(95% CI 1.03-6.36) and 2.31(95% CI 1.38-4.02) among the Desmedt, Chin, Pawitan and NKI sets respectively (Table 1), and the adjusted HR values for PFS were 2.16(95% CI 1.32-3.54), 0.80(95% CI 0.40-1.62), 3.28(95% CI 1.43-7.89), 2.20(95% CI 1.44-3.39), 2.21(95% CI 1.37-3.64) and 1.96(95% CI 1.31-2.98) among Desmedt, Chin, Pawitan, Wang, Ivshina and NKI sets respectively. The HR values from the pooled analysis were 2.19 (95%CI 1.28-3.49) and 2.34 (95% CI 1.62-3.41) for OS and PFS respectively. Meta-analysis (Fixed effect model) also yielded similar result (Table 1). It was revealed that RRM2 mRNA levels were significantly associated with poor outcome in breast cancer patients.
RRM2 was negatively correlated with the survival of breast cancer patients in a dose-dependent manner. Kaplan-Meier analysis revealed that RRM2 mRNA levels (either 209773_s_at or 201890_at) significantly and negatively impacted the PFS of breast cancer patients in the Pawitan and Ivshina sets (Figure 3A and B). As RRM2 mRNA levels increased, the outcome for the breast cancer patient became worse, and the relative risk of death increased in a dose-dependent manner (Figure 3C). Similar results were also observed for all others except the Chin set (Additional file 1: Figure S3A to J).
The prognostic performance of RRM2 mRNA levels was also compared with some of the prognostic gene signatures and the TNM stage in the NKI dataset. The Cox analyses revealed that the increase in HR was steadily and significantly correlated with increases in RRM2 mRNA levels (Figure 3C). The prognostic performance of RRM2 was similar to the 70-genes, wound response genes and 21 gene recurrence score. Also, both the RRM2 mRNA levels and the gene signatures were more efficient prognostic indicators than the TNM stage.
The prognostic significance of RRM2 protein was also validated in the ZJU set. Multivariate Cox proportional hazard analysis confirmed that high levels of RRM2 protein were significantly and negatively associated with OS of breast cancer patients. After adjusting for RRM2B, ER, PR, pre-operation chemotherapy, adjuvant chemotherapy, radiotherapy and other factors, the HR of RRM2 was 9.84 (95% CI 2.00-80.03, P = 0.003) (Figure 3D and E).
The prognostic significance of RRM2in ER-negative breast cancers
ER-negative breast cancers (including the HER2-positive and TNBC subtypes) have poorer prognosis . A Kaplan-Meier analysis of the Pawitan set revealed that RRM2 mRNA levels, similar to uPA, which is a validated prognostic biomarker for ER-negative breast cancer , were significantly associated with poor OS and PFS in a dose-dependent manner in ER-negative, but not ER-positive, breast cancers (Additional file 1: Figure S4). Pooled analysis results also suggested that RRM2 mRNA levels were significantly correlated with poor OS and PFS (Figure 4A and B). Further multivariate Cox analysis indicated that the HR value for OS increased steadily as the levels of RRM2 increased in the ER-negative, but not in ER-positive, subgroup (Figure 4C). For PFS, RRM2 was prognostic for poor survival in both the ER-negative and ER-positive breast cancers, but was more significant for the ER-negative breast cancers (Figure 4D). Similar results were also seen from the ZJU set. A multivariate Cox analysis indicated that the protein levels of RRM2 were significantly associated with poor OS in the ER-negative subgroup (n = 64, HR = 22.4, 95% CI 2.2-285.0) (Figure 4E), but not in the ER-positive subgroup (n = 96) (Figure 4F).
The prognostic value of RRM2 for the HER2-positive and TNBC subtypes was also evaluated. In the ZJU set, a total of 12 HER2-positive breast cancer and 28 TNBCs were identified. Among these patients, none of the RRM2-low cases (n = 15) died of breast cancer during the following up period (Log rank P = 0.025). However, due to the small number of cases there was no statistical significance in the HER2-positive or TNBC subsets when separated.
Prognostic performance of RRM2in ER-negative breast cancers
The efficiency of RRM2 in predicting the outcome for patients with ER-negative breast cancer was compared with other prognostic markers for this subtype including ki-67, HER2, tumor invasiveness, lymph node involvement, distant organ metastasis, Elson histological stage and uPA. The uPA marker has been reported as a powerful prognostic predictor of breast cancer  and its performance was comparable to the Immune Response Module (IRM, 7 genes) and HER2- Derived Prognostic Predictor (HDPP, 158 genes) signatures . In the pooled dataset, the HR of RRM2 increased steadily in a dose-dependent manner both for OS and PFS (Figure 5A and B). It was comparable to distant organ metastasis for the analysis of both OS and PFS. Furthermore, RRM2 was more accurate than ki-67, HER2, lymph node involvement, tumor grade and uPA for classifying a subgroup of ER-negative breast cancer patients that were at higher risk. Similar results were obtained using the ZJU set (Figure 5C). However, because there are no data about distant metastasis in the ZJU patient set, we could not make this comparison.
Inhibition of RRM2reduced invasion ability and enhanced the drug sensitivity to doxorubicin in ER-negative breast cancer cells
The above findings suggested that RRM2 inhibitors might be a promising anti-cancer agent for reducing invasion ability and enhancing the chemosensitivity in ER-negative breast cancers. Here, the RRM2 siRNA was used to inhibit the expression of RRM2 in MDA-MB-231(ER-negative) and ZR-75-1(ER-positive) breast cancer cell lines. It was shown that the protein expression levels of RRM2 were significantly reduced by RRM2 siRNA in both cell lines (Figure 6A). As with inhibition of RRM2, the invasion abilities of cells were significantly reduced about 50% in both MDA-MB-231 and ZR-75-1 cells (p < 0.05) (Figure 6B and C). This phenomenon also could be seen in MCF-7, another ER-positive breast cancer cell line (Data not shown). This is compatible with our previous findings [23, 50].
A novel RR inhibitor, COH 29, was reported to inhibit RR activity by blocking RRM2 binding to RRM1, thus inactivating the RR holoenzyme . Here, we have preliminarily demonstrated that inhibition of RRM2 by COH29 could significantly reduce the growth of MCF-7 and MDA-MB-231 breast cancer cells. Also, we found that COH29 also could significantly enhance the chemosensitivity of breast cancer cells to doxorubicin in ER-negative MDA-MB-231 cells, but not in ER-positive MCF-7 cells (Figure 6D). This synergetic effect also could not be seen in ZR-75-1 cells (Data not shown). Therefore, the RRM2 might be served as a therapeutic biomarker for application of RR inhibitor to treatment of ER-negative breast cancers.
A recent investigation demonstrated that the reproducibility of current preclinical cancer researches was only 11% . A similar result was also seen in another report . Therefore, third party validations play a very important role in proving the reliability of research results. Open access public microarray databases provided an opportunity for researchers to validate their findings in a non-biased platform. Here, we demonstrated that RRM2 impacted the survivability of breast cancer patients in 6 of 7 independently microarray data sets (except Chin’s set), which indicates the high reproducibility of this result. The reason why RRM2 did not impact poor survival in Chin’s set might be of the small sample size (only 118 cases), which did not have enough power to detect significant differences. We showed that increased RRM2 expression was accompanied by increased expression of other widely used biomarkers and gene sets of proliferation, undifferentiation, and stem/progenitor phenotypes. The role of RRM2 in malignancy has been demonstrated in several solid cancers including bladder , pancreatic , gastric  and colorectal cancers , and has been attributed to increased proliferation and invasion . Previously, gene expression profiling has uncovered molecular signatures that influence the choice for breast cancer treatment. Four distinct subgroups are characterized by different gene expression patterns and show diverse epidemiologic, histopathology and clinical features, respond to different therapeutic strategies. Multiple studies have shown the outcomes are much poorer for patients with ER-negative breast cancer . This study demonstrated that either protein or mRNA of RRM2 was highly expressed in the molecular subtypes of TNBC, HER-2 positive and Luminal B, which were subtypes with relative poorer outcome. The protein expression of RRM2 was also positively and significantly associated with CD44+/CD24-/low, a tumor stem/progenitor cell biomarker. In addition, RRM2 levels were significantly correlated with poorly differentiated breast cancers (p < 0.001). Overall, RRM2 was significantly related to aggressiveness and associated with poor outcome in breast cancers, especially in ER-negative breast cancer. These findings were validated in many independent studies including participants with different genetic and socio-economic background, which indicates the reliability of this study. In 3D–signature study, the RRM2 could not impact poor survivability in ER negative breast cancer . This is because there were insufficient ER negative samples in their study.
Currently, prognostic/predictive factors including primary tumor size, lymph node stage, histological grade and hormone receptor status are used routinely in clinical practice to choose systemic therapies . However, breast cancers are significantly heterogeneous even between patients whose tumors have similar clinical and pathological characteristics. Therefore, distinguishing patients who will be more prone to relapse and will need more intensive therapy is an urgent problem for oncologists. Based on high throughput gene expression, a series of multiple gene signatures have been developed for predicting the survival of breast cancer patients. RRM2 also has been included in some gene signatures [27–30]. Overall, our data also indicate that the 70-gene, 21-gene recurrence score and wound-response gene signatures are more efficient than TNM stage for predicting patient outcome. Although multiple-gene-based signatures significantly enhance the accuracy of prediction, they also dramatically increase the diagnostic cost. Interestingly, the prognostic performance of RRM2 was comparable to all of the above gene signatures for breast cancers overall in the NKI set. Furthermore, RRM2 was more accurate for predicting the outcome of patients with ER-negative breast cancer than uPA, Ki-67, HER2 and Elson grade. Therefore, determining the RRM2 mRNA or protein levels could be more accurate and cost-effective for predicting the outcome of patients with ER-negative breast cancers.
Out of all breast cancers diagnosed, approximately 15% to 20% are HER2-enriched and another 15%-17% are TNBC . These two types typically experience a poorer prognosis than the ER-positive subtypes and biomarkers predicting their prognosis and therapeutic response are urgently needed. Several studies have revealed that RRM2 expression is also related to resistance to gamma radiation and chemotherapeutic agents [56, 57], which is presumably due to the essential role of RRM2 in DNA repair . Many studies have shown that using small interfering RNA (siRNA) to inhibit RRM2 overexpression significantly reverses cancer cell resistance to chemotherapy agents and gamma radiation [56, 57]. In a previous study, we demonstrated that using hydroxyurea to inhibit RRM2 significantly enhanced the chemosensitivity of KB cells (a head and neck cancer cell line) to gemcitabine . Cell culture studies also indicate that expression of siRNAs targeting RRM2 significantly reduced the invasion ability in MDA-MB-231 and ZR-75-1 breast cancer cells, which was consistent with our previous studies on other cancers [23, 50]. Similar to the way that HER2 is used to identify subtypes that are likely to respond to trastuzumab, RRM2 may be a predictive marker that can stratify the more fatal types of ER-negative breast cancers to identify breast cancer subtypes that will respond to RRM2 inhibitors. Alternatively, sorting patients who have ER-negative breast cancers according to their levels of RRM2 expression may prevent over-treatment of the low expressers. The COH29 is a newly developed RR inhibitor that inhibits RR activity by blocking RRM2 binding to RRM1, thus inactivating the RR holoenzyme . Here, our pilot study indicated that inhibition of RRM2 by COH29 could significantly reduce the growth of breast cancer cells. Also, it could significantly enhance the sensitivity of MDA-MB-231 (ER-negative) cell to doxorubicin, a common therapeutic agent for breast cancers. Nevertheless, we will perform further clinical trials using novel RRM2 inhibitors (such as RRM2 siRNA nanoparticles and COH29) to address this hypothesis. According to our data, further investigating and understanding the effect of RRM2 expression on breast cancer malignancy could lead to a novel approach for treating ER-negative breast cancers.
Our study revealed RRM2 was not only an indicator for breast cancers’ proliferation and invasiveness, but also a biomarker for undifferentiation and stem/progenitor-like phenotypes. The prognostic performance of RRM2 was superior to the current pathology stages and biomarkers in ER-negative breast cancers. Targeting on RRM2 might be a promising therapeutic method to enhance the chemosensitivity for ER-negative breast cancer patients.
Human ribonucleotide reductase large subunit M1
Human ribonucleotide reductase small subunit M2
a p53 dependent human ribonucleotide reductase small subunit R2 (p53R2)
Triple negative breast cancer
Formalin-fixed, paraffin-embedded tissue block
Poly ADP ribose polymerase
short interfering RNA
Gene set enrichment analysis
Progression- free survival
Proportional hazard ratio
- 95% CI:
95% confidence interval
Immune response rmodule
HER2- derived prognostic predictor.
Siegel R, Naishadham D, Jemal A: Cancer statistics, 2013. CA Cancer J Clin. 2013, 63 (1): 11-30. 10.3322/caac.21166.
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 (19): 10869-10874. 10.1073/pnas.191367098.
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 (6797): 747-752. 10.1038/35021093.
Sorlie T, Wang Y, Xiao C, Johnsen H, Naume B, Samaha RR, Borresen-Dale AL: Distinct molecular mechanisms underlying clinically relevant subtypes of breast cancer: gene expression analyses across three different platforms. BMC Genomics. 2006, 7: 127-10.1186/1471-2164-7-127.
Romond EH, Perez EA, Bryant J, Suman VJ, Geyer CE, Davidson NE, Tan-Chiu E, Martino S, Paik S, Kaufman PA, Swain SM, Pisansky TM, Fehrenbacher L, Kutteh LA, Vogel VG, Visscher DW, Yothers G, Jenkins RB, Brown AM, Dakhil SR, Mamounas EP, Lingle WL, Klein PM, Ingle JN, Wolmark N: Trastuzumab plus adjuvant chemotherapy for operable HER2-positive breast cancer. N Engl J Med. 2005, 353 (16): 1673-1684. 10.1056/NEJMoa052122.
Andre F, Pusztai L: Molecular classification of breast cancer: implications for selection of adjuvant chemotherapy. Nat Clin Pract. 2006, 3 (11): 621-632. 10.1038/ncponc0636.
Nanda R: “Targeting” triple-negative breast cancer: the lessons learned from BRCA1-associated breast cancers. Semin Oncol. 2011, 38 (2): 254-262. 10.1053/j.seminoncol.2011.01.007.
van ’t Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.
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 (27): 2817-2826. 10.1056/NEJMoa041588.
Saal LH, Johansson P, Holm K, Gruvberger-Saal SK, She QB, Maurer M, Koujak S, Ferrando AA, Malmstrom P, Memeo L, Isola J, Bendahl PO, Rosen N, Hibshoosh H, Ringnér M, Borg A, Parsons R: Poor prognosis in carcinoma is associated with a gene expression signature of aberrant PTEN tumor suppressor pathway activity. Proc Natl Acad Sci U S A. 2007, 104 (18): 7564-7569. 10.1073/pnas.0702507104.
Chang HY, Nuyten DS, Sneddon JB, Hastie T, Tibshirani R, Sorlie T, Dai H, He YD, van’t Veer LJ, Bartelink H, van de Rijn M, Brown PO, van de Vijver MJ: Robustness, scalability, and integration of a wound-response gene expression signature in predicting breast cancer survival. Proc Natl Acad Sci U S A. 2005, 102 (10): 3738-3743. 10.1073/pnas.0409462102.
Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, Wong JE, Liu ET, Bergh J, Kuznetsov VA, Miller LD: Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006, 66 (21): 10292-10301. 10.1158/0008-5472.CAN-05-4414.
Teschendorff AE, Miremadi A, Pinder SE, Ellis IO, Caldas C: An immune response gene expression module identifies a good prognosis subtype in estrogen receptor negative breast cancer. Genome Biol. 2007, 8 (8): R157-10.1186/gb-2007-8-8-r157.
Staaf J, Ringner M, Vallon-Christersson J, Jonsson G, Bendahl PO, Holm K, Arason A, Gunnarsson H, Hegardt C, Agnarsson BA, Luts L, Grabau D, Fernö M, Malmström PO, Johannsson OT, Loman N, Barkardottir RB, Borg A: Identification of subtypes in human epidermal growth factor receptor 2–positive breast cancer reveals a gene signature prognostic of outcome. J Clin Oncol. 2010, 28 (11): 1813-1820. 10.1200/JCO.2009.22.8775.
Nordlund P, Reichard P: Ribonucleotide reductases. Annu Rev Biochem. 2006, 75: 681-706. 10.1146/annurev.biochem.75.103004.142443.
Jordan A, Reichard P: Ribonucleotide reductases. Annu Rev Biochem. 1998, 67: 71-98. 10.1146/annurev.biochem.67.1.71.
Heidel JD, Liu JY, Yen Y, Zhou B, Heale BS, Rossi JJ, Bartlett DW, Davis ME: Potent siRNA inhibitors of ribonucleotide reductase subunit RRM2 reduce cell proliferation in vitro and in vivo. Clin Cancer Res. 2007, 13 (7): 2207-2215. 10.1158/1078-0432.CCR-06-2218.
Zhou BS, Tsai P, Ker R, Tsai J, Ho R, Yu J, Shih J, Yen Y: Overexpression of transfected human ribonucleotide reductase M2 subunit in human cancer cells enhances their invasive potential. Clin Exp Metastasis. 1998, 16 (1): 43-49.
Morikawa T, Maeda D, Kume H, Homma Y, Fukayama M: Ribonucleotide reductase M2 subunit is a novel diagnostic marker and a potential therapeutic target in bladder cancer. Histopathology. 2010, 57 (6): 885-892. 10.1111/j.1365-2559.2010.03725.x.
Morikawa T, Hino R, Uozaki H, Maeda D, Ushiku T, Shinozaki A, Sakatani T, Fukayama M: Expression of ribonucleotide reductase M2 subunit in gastric cancer and effects of RRM2 inhibition in vitro. Hum Pathol. 2010, 41 (12): 1742-1748. 10.1016/j.humpath.2010.06.001.
Hsu NY, Wu JY, Liu X, Yen Y, Chen CY, Chou MC, Lin CH, Lee H, Cheng YW: Expression status of ribonucleotide reductase small subunits hRRM2/p53R2 as prognostic biomarkers in stage I and II non-small cell lung cancer. Anticancer Res. 2011, 31 (10): 3475-3481.
Liu X, Lai L, Wang X, Xue L, Leora S, Wu J, Hu S, Zhang K, Kuo ML, Zhou L, Zhang H, Wang Y, Wang Y, Zhou B, Nelson RA, Zheng S, Zhang S, Chu P, Yen Y: Ribonucleotide reductase small subunit M2B prognoses better survival in colorectal cancer. Cancer Res. 2011, 71 (9): 3202-3213. 10.1158/0008-5472.CAN-11-0054.
Liu X, Zhang H, Lai L, Wang X, Loera S, Xue L, He H, Zhang K, Hu S, Huang Y, Nelson RA, Zhou B, Zhou L, Chu P, Zhang S, Zheng S, Yen Y: Ribonucleotide reductase small subunit M2 serves as a prognostic biomarker and predicts poor survival of colorectal cancers. Clin Sci (Lond). 2013, 124 (9): 567-578. 10.1042/CS20120240.
Duxbury MS, Whang EE: RRM2 induces NF-kappaB-dependent MMP-9 activation and enhances cellular invasiveness. Biochem Biophys Res Commun. 2007, 354 (1): 190-196. 10.1016/j.bbrc.2006.12.177.
Zhang K, Hu S, Wu J, Chen L, Lu J, Wang X, Liu X, Zhou B, Yen Y: Overexpression of RRM2 decreases thrombspondin-1 and increases VEGF production in human cancer cells in vitro and in vivo: implication of RRM2 in angiogenesis. Mol Cancer. 2009, 8: 11-10.1186/1476-4598-8-11.
Reid G, Wallant NC, Patel R, Antonic A, Saxon-Aliifaalogo F, Cao H, Webster G, Watson JD: Potent subunit-specific effects on cell growth and drug sensitivity from optimised siRNA-mediated silencing of ribonucleotide reductase. J RNAi Gene Silencing. 2009, 5 (1): 321-330.
Martin KJ, Patrick DR, Bissell MJ, Fournier MV: Prognostic breast cancer signature identified from 3D culture model accurately predicts clinical outcome across independent datasets. PLoS One. 2008, 3 (8): e2994-10.1371/journal.pone.0002994.
Habel LA, Sakoda LC, Achacoso N, Ma XJ, Erlander MG, Sgroi DC, Fehrenbacher L, Greenberg D, Quesenberry CP: HOXB13:IL17BR and molecular grade index and risk of breast cancer death among patients with lymph node-negative invasive disease. Breast Cancer Res. 2013, 15 (2): R24-10.1186/bcr3402.
Ma XJ, Salunga R, Dahiya S, Wang W, Carney E, Durbecq V, Harris A, Goss P, Sotiriou C, Erlander M, Sgroi D: A five-gene molecular grade index and HOXB13:IL17BR are complementary prognostic factors in early stage breast cancer. Clin Cancer Res. 2008, 14 (9): 2601-2608. 10.1158/1078-0432.CCR-07-5026.
Li Z, Peng L, Han S, Huang Z, Shi F, Cai Z, Li X, Zhang P, Zhu H, Jin W: [Screening molecular markers in early breast cancer of the same pathological types but with different prognoses using Agilent gene chip]. Nan Fang Yi Ke Da Xue Xue Bao. 2013, 33 (10): 1483-1488.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, Desmedt C, Larsimont D, Cardoso F, Peterse H, Nuyten D, Buyse M, Van de Vijver MJ, Bergh J, Piccart M, Delorenzi M: Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006, 98 (4): 262-272. 10.1093/jnci/djj052.
Huang E, Cheng SH, Dressman H, Pittman J, Tsou MH, Horng CF, Bild A, Iversen ES, Liao M, Chen CM, West M, Nevins JR, Huang AT: Gene expression predictors of breast cancer outcomes. Lancet. 2003, 361 (9369): 1590-1596. 10.1016/S0140-6736(03)13308-9.
van de Vijver MJ, He YD, van’t Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, Parrish M, Atsma D, Witteveen A, Glas A, Delahaye L, van der Velde T, Bartelink H, Rodenhuis S, Rutgers ET, Friend SH, Bernards R: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347 (25): 1999-2009. 10.1056/NEJMoa021967.
Chin SF, Wang Y, Thorne NP, Teschendorff AE, Pinder SE, Vias M, Naderi A, Roberts I, Barbosa-Morais NL, Garcia MJ, Iyer NG, Kranjac T, Robertson JF, Aparicio S, Tavaré S, Ellis I, Brenton JD, Caldas C: Using array-comparative genomic hybridization to define molecular portraits of primary breast cancers. Oncogene. 2007, 26 (13): 1959-1970. 10.1038/sj.onc.1209985.
Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365 (9460): 671-679. 10.1016/S0140-6736(05)17947-1.
Smeds J, Miller LD, Bjohle J, Hall P, Klaar S, Liu ET, Pawitan Y, Ploner A, Bergh J: Gene profile and response to treatment. Ann Oncol. 2005, 16 Suppl 2: ii195-ii202.
Desmedt C, Piette F, Loi S, Wang Y, Lallemand F, Haibe-Kains B, Viale G, Delorenzi M, Zhang Y, d’Assignies MS, Bergh J, Lidereau R, Ellis P, Harris AL, Klijn JG, Foekens JA, Cardoso F, Piccart MJ, Buyse M, Sotiriou C, TRANSBIG Consortium: Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series. Clin Cancer Res. 2007, 13 (11): 3207-3214. 10.1158/1078-0432.CCR-06-2765.
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 (43): 15545-15550. 10.1073/pnas.0506580102.
Millar EK, Graham PH, McNeil CM, Browne L, O’Toole SA, Boulghourjian A, Kearsley JH, Papadatos G, Delaney G, Fox C, Nasser E, Capp A, Sutherland RL: Prediction of outcome of early ER + breast cancer is improved using a biomarker panel, which includes Ki-67 and p53. Br J Cancer. 2011, 105 (2): 272-280. 10.1038/bjc.2011.228.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thurlimann B, Senn HJ, Panel M: Strategies for subtypes--dealing with the diversity of breast cancer: highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol. 2011, 22 (8): 1736-1747. 10.1093/annonc/mdr304.
Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF: Prospective identification of tumorigenic breast cancer cells. Proc Natl Acad Sci U S A. 2003, 100 (7): 3983-3988. 10.1073/pnas.0530291100.
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression. Proc Natl Acad Sci U S A. 2004, 101 (25): 9309-9314. 10.1073/pnas.0401994101.
Sarrio D, Rodriguez-Pinilla SM, Hardisson D, Cano A, Moreno-Bueno G, Palacios J: Epithelial-mesenchymal transition in breast cancer relates to the basal-like phenotype. Cancer Res. 2008, 68 (4): 989-997. 10.1158/0008-5472.CAN-07-2017.
Zhang Y, Crump M, Berge SA: Purging of contaminating breast cancer cells from hematopoietic progenitor cell preparations using activation enhanced cell death. Breast Cancer Res Treat. 2002, 72 (3): 265-278. 10.1023/A:1014965726663.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA: An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008, 40 (5): 499-507. 10.1038/ng.127.
Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, Chi JT, van de Rijn M, Botstein D, Brown PO: Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004, 2 (2): E7-10.1371/journal.pbio.0020007.
Naderi A, Teschendorff AE, Barbosa-Morais NL, Pinder SE, Green AR, Powe DG, Robertson JF, Aparicio S, Ellis IO, Brenton JD, Caldas C: A gene-expression signature to predict survival in breast cancer across independent data sets. Oncogene. 2007, 26 (10): 1507-1516. 10.1038/sj.onc.1209920.
Carey LA, Perou CM, Livasy CA, Dressler LG, Cowan D, Conway K, Karaca G, Troester MA, Tse CK, Edmiston S, Deming SL, Geradts J, Cheang MC, Nielsen TO, Moorman PG, Earp HS, Millikan RC: Race, breast cancer subtypes, and survival in the Carolina Breast Cancer Study. JAMA. 2006, 295 (21): 2492-2502. 10.1001/jama.295.21.2492.
Urban P, Vuaroqueaux V, Labuhn M, Delorenzi M, Wirapati P, Wight E, Senn HJ, Benz C, Eppenberger U, Eppenberger-Castori S: Increased expression of urokinase-type plasminogen activator mRNA determines adverse prognosis in ErbB2-positive primary breast cancer. J Clin Oncol. 2006, 24 (26): 4245-4253. 10.1200/JCO.2005.05.1912.
Huang Y, Liu X, Wang YH, Yeh SD, Chen CL, Nelson RA, Chu P, Wilson T, Yen Y: The prognostic value of ribonucleotide reductase small subunit M2 in predicting recurrence for prostate cancers. Urol Oncol. 2014, 32 (1): 51-e59-19
Zhou B, Su L, Hu S, Hu W, Yip ML, Wu J, Gaur S, Smith DL, Yuan YC, Synold TW, Horne D, Yen Y: A small-molecule blocking ribonucleotide reductase holoenzyme formation inhibits cancer cell growth and overcomes drug resistance. Cancer Res. 2013, 73 (21): 6484-6493. 10.1158/0008-5472.CAN-13-1094.
Begley CG, Ellis LM: Drug development: raise standards for preclinical cancer research. Nature. 2012, 483 (7391): 531-533. 10.1038/483531a.
Prinz F, Schlange T, Asadullah K: Believe it or not: how much can we rely on published data on potential drug targets?. Nat Rev Drug Discov. 2011, 10 (9): 712-10.1038/nrd3439-c1.
Duxbury MS, Ito H, Zinner MJ, Ashley SW, Whang EE: RNA interference targeting the M2 subunit of ribonucleotide reductase enhances pancreatic adenocarcinoma chemosensitivity to gemcitabine. Oncogene. 2004, 23 (8): 1539-1548. 10.1038/sj.onc.1207272.
Anderson TJ, Alexander FE, Lamb J, Smith A, Forrest AP: Pathology characteristics that optimize outcome prediction of a breast screening trial. Br J Cancer. 2000, 83 (4): 487-492. 10.1054/bjoc.2000.1286.
Aimiuwu J, Wang H, Chen P, Xie Z, Wang J, Liu S, Klisovic R, Mims A, Blum W, Marcucci G, Chan KK: RNA-dependent inhibition of ribonucleotide reductase is a major pathway for 5-azacytidine activity in acute myeloid leukemia. Blood. 2012, 119 (22): 5229-5238. 10.1182/blood-2011-11-382226.
Zhang YW, Jones TL, Martin SE, Caplen NJ, Pommier Y: Implication of checkpoint kinase-dependent up-regulation of ribonucleotide reductase R2 in DNA damage response. J Biol Chem. 2009, 284 (27): 18085-18095. 10.1074/jbc.M109.003020.
Niida H, Katsuno Y, Sengoku M, Shimada M, Yukawa M, Ikura M, Ikura T, Kohno K, Shima H, Suzuki H, Tashiro S, Nakanishi M: Essential role of Tip60-dependent recruitment of ribonucleotide reductase at DNA damage sites in DNA repair during G1 phase. Genes Dev. 2010, 24 (4): 333-338. 10.1101/gad.1863810.
Zhou B, Mi S, Mo X, Shih J, Tsai J, Hu E, Hsu M, Kay K, Yen Y: Time and sequence dependence of hydroxyurea in combination with gemcitabine in human KB cells. Anticancer Res. 2002, 22 (3): 1369-1377.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/14/664/prepub
This work was supported by Grants from the Major State Basic Research Development Program of China, 973 Program (2004CB518707); the Zhejiang Provincial Natural Science Foundation of China (R2090353); US National Institutes of Health Grant (R01 CA127541), and the China Scholarship Council. We thank Drs. Xutao Deng, Zheng Liu, Rui Ba, Huarong Chen, Dan Li and Jiaping Peng for help with clinical data collection, tissue arrays construction and microarray data downloading. We also thank Margaret Morgan, Nancy Linford and Mansze Kong for English editing.
The authors declare that they have no competing interests.
HZ carried out breast cancer patients’ database creation, multiple tissue array construction, data collection and manuscript writing. XL was charged for study design, data analysis and manuscript writing. CDW conducted public breast cancer data downloading and normalization. YH scored IHC staining. SL conducted IHC staining. LX carried out cell cytotoxicity assay and provided tissue samples. SZ coordinated the project. PC is a pathologist who helped with IHC staining and scoring. SZ provided breast cancer patients clinical, pathology and following-up data. YY was in charge of study design and manuscript writing. All authors read and approved the final manuscript.
Hang Zhang, Xiyong Liu contributed equally to this work.