Original Research ARTICLE
UBE2C Overexpression Aggravates Patient Outcome by Promoting Estrogen-Dependent/Independent Cell Proliferation in Early Hormone Receptor-Positive and HER2-Negative Breast Cancer
- 1Laboratory of Cancer Genomics and Molecular Pathology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea
- 2Laboratory of Molecular Pathology and Cancer Genomics, Department of Pharmacy, College of Pharmacy, Seoul National University, Seoul, South Korea
- 3Gencurix, Inc., Seoul, South Korea
- 4The Center for Companion Diagnostics, LOGONE Bio Convergence Research Foundation, Seoul, South Korea
- 5College of Pharmacy, Daegu Catholic University, Gyeongsan-si, South Korea
- 6Department of Pathology and Translational Genomics, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea
- 7Department of Health Sciences and Technology, SAIHST, Sungkyunkwan University, Seoul, South Korea
- 8The Center for Anti-cancer Companion Diagnostics, BioMAX/N-Bio, Seoul National University, Seoul, South Korea
- 9Department of Molecular Medicine and Biopharmaceutical Sciences, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, South Korea
We previously showed that UBE2C mRNA expression is significantly associated with poor prognosis only in patients with hormone receptor (HR)+/human epidermal growth factor receptor 2 (HER2)– breast cancer. In this study, we further reanalyzed the correlation between UBE2C mRNA expression and clinical outcomes in patients with HR+/HER2– breast cancer, and we investigated the molecular mechanism underlying the role of UBE2C modulation in disease progression in this subgroup of patients. Univariate and multivariate analyses showed that high UBE2C expression was associated with significantly shorter survival of breast cancer patients with pN0 and pN1 tumors but not pN2/N3 tumors (P < 0.05). In vitro functional experiments in HR+/HER2– breast cancer cells showed that UBE2C expression is a tumorigenic factor, and that estrogen upregulated UBE2C mRNA and protein by directly binding to the UBE2C promoter region. UBE2C knockdown inhibited cell proliferation by affecting cell cycle progression, and UBE2C overexpression was associated with estrogen-independent growth. UBE2C depletion markedly increased the cytotoxicity of tamoxifen by inducing apoptosis. The present findings suggest that UBE2C overexpression is correlated with relapse and promotes estrogen-dependent/independent proliferation in early HR+/HER2– breast cancer.
Breast cancer is classified according to estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) expression status into hormone receptor HR+/HER2– (ER+ or PR+/HER2–), HR+/HER2+ (ER+ or PR+/HER2+), HR–/HER2+ (ER–/PR–/HER2+), or triple-negative breast cancer (TNBC; ER–/PR–/HER2–); the HR+/HER2– subtype accounts for 25–40% of breast cancers (1, 2). Above 70% of patients with HR+/HER2– breast cancer remain disease-free without chemotherapy after surgery. However, more than 5% of patients with HR+/HER2– breast cancer have a higher of late recurrence beyond 5 years after primary adjuvant hormone therapy than those with other subtypes (3–5). Therefore, an exact estimation of the risk of recurrence could be useful for decision-making regarding the duration of adjuvant hormone therapy or treatment with adjuvant chemotherapy.
Because gene expression provides significant prognostic or predictive information for patients with HR+ breast cancer, several commercial assays such as Oncotype DX, MammaPrint, and EndoPredict, which are based on the expression of multiple genes in frozen or formalin-fixed, paraffin-embedded (FFPE) samples, have been developed (6–8). These assays predict the risk of distant recurrence after hormone therapy and help to identify patients who may benefit from adjuvant chemotherapy by discriminating between high- and low-risk patients with early HR+ breast cancer. However, these assays were validated in only Western countries, therefore the prognostic value of these signatures in Asian breast cancer patients remains unclear as well as Korean. In a previous study, we identified 384 candidate prognostic genes associated with distant metastasis in patients with lymph node (LN) negative early breast cancer (9). Of these 384 genes, 16 candidate prognostic genes were selected, and the association between the expression of these genes determined by quantitative real-time reverse transcription-PCR (qRT-PCR) and clinical outcomes was evaluated in 997 FFPE tissue samples from 819 patients with breast cancer (9, 10). Six genes were selected to develop a prognostic model for predicting the risk of distant recurrence or distant metastasis, and the performance of the model was compared with that of classical clinicopathological risk factors (10–13).
One of the six genes included in the model is ubiquitin-conjugating enzyme E2C (UBE2C), which plays an important role in the ubiquitin-proteasome system. The ubiquitin-proteasome system precisely regulates the cell cycle through proteasome-mediated protein degradation pathways in eukaryotes. Dysregulation of ubiquitin-dependent processes is involved in the occurrence and progression of cancer, and the ubiquitin-proteasome system has thus become a new druggable target in certain cancers such as bladder carcinoma, cervical cancer, colorectal carcinoma, esophageal squamous cell carcinoma, lung cancer, and pancreatic carcinoma (14–20). High UBE2C expression is associated with a high degree of malignancy, low differentiation, high metastatic tendency, and poor patient survival in a wide range of solid tumors including breast cancer (21–25). We previously reported that UBE2C is significantly associated with poor prognosis in patients with HR+/HER2– breast cancer (9, 10). Consistent with our report, Filipits et al. reported the use of UBE2C mRNA expression as a marker in the EndoPredict assay for predicting the risk of recurrence or distant metastasis in patients with HR+/HER2– breast cancer (8). However, the clinical and functional significance of UBE2C expression in HR+/HER2– breast cancer remains unknown.
In this study, we examined the correlation between UBE2C mRNA expression and clinical outcomes in patients with HR+/HER2– breast cancer. We also evaluated the expression status of UBE2C and investigated the molecular mechanism underlying the role of UBE2C regulation in HR+/HER2– breast cancer progression.
Materials and Methods
A total of 997 FFPE tissue specimens were obtained from patients with breast cancer who underwent curative resection of primary tumors with LN dissection at Samsung Medical Center (SMC, Korea) between 1994 and 2002. The protocol for the present study was approved by the SMC Institutional Review Board (IRB file No. 2008-12-035). Tumor size and LN involvement were evaluated according to the American Joint Committee on Cancer 7th TNM Staging System, and tumor histological grades were determined according to the Bloom–Richardson grading scheme. Paraffin-embedded tissue samples (mounted on slides) were analyzed to define tumor regions and select representative tumor areas for further analysis. Breast cancer specimens were classified into subtypes using an immunohistochemical assay with ER, PR, and HER2 as markers.
qRT-PCR Analysis of Patient Samples
RNA was isolated from patient-derived FFPE samples using a tissue preparation system (Siemens AG), and qRT-PCR was performed to measure the expression levels of UBE2C (Roche Applied Science). The results of qRT-PCR were expressed as cycle threshold (Ct) values. The Ct value for UBE2C was normalized to a relative expression value (ΔCt value) using three reference genes (CTBP1, CUL1, and UBQLN1), and it was calculated as follows:
Expression values were compared using the independent t-test and Chi-square test. Univariate and multivariate survival analyses were performed using the Cox proportional hazard model. Study endpoints were disease-free survival (DFS, time from primary surgery to any first event of recurrence), distant metastasis-free survival (DMFS, time from primary surgery until distant metastasis), and overall survival (OS, time from primary surgery until death from any cause). The genomic and clinical variables that were significant in the univariate analysis were included in the multivariate analysis to identify independent predictors of DFS, DMFS, and OS. Survival rates were estimated by the Kaplan–Meier (KM) method, and the log-rank test was used for assessing the difference in survival among patient groups. To classify patients into low- and high-expression groups for distant metastasis, the cut-off value was set as the value that maximized the sum of sensitivity and specificity. Statistical significance was considered at P < 0.05. All statistical analyses were performed using R 3.5.1 (http://r-project.org).
The human breast cell lines were obtained from the American Type Culture Collection and Korean Cell Line Bank. All cell lines were cultured according to the manufacturers' recommendations. Cell lines were validated by human cell line authentication (STR DNA profiling) using the AmpFLSTR™ Identifiler PCR Amplification Kit (Thermo Fisher Scientific).
Real-Time qRT-PCR in Cells
The expression levels of UBE2C mRNA were measured by real-time qRT-PCR. Total RNA was isolated using RiboEx (GeneAll) and the Hybrid-R kit (GeneAll) followed by the Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science), according to the manufacturers' instructions. qRT-PCR was performed on cDNA using a LightCycler 480 System (Roche Applied Science). The UBE2C primers used were as follows: 5′-TGCCGAGCTCTGGAAAAA-3′ (forward primer) and 5′-AAAAGACGACACAAGGACAGG-3′ (reverse primer). The amplified UBE2C cDNAs obtained using these primers consisted of five transcript isoforms among seven coding sequence (CDS) transcripts (https://www.ensembl.org). The HPRT primers were used as a control. Commercial Universal Probe Library (UPL) probes were purchased from Roche Applied Science.
Western Blot Analysis
Cells were lysed with RIPA buffer [20 mM Tris-HCl (pH 8.0), 150 mM NaCl, 10% glycerol, 1% NP40, and 2 mM EDTA]. Equal amounts of protein were subjected to 10% SDS-PAGE and transferred to a nitrocellulose membrane (Millipore). The membrane was incubated overnight at 4°C with primary anti-UBE2C (Abcam, ab187181), anti-cleaved PARP (Cell Signaling, #9541), anti-cleaved CASPASE 3 (Cell Signaling, #9661), anti-ERα (Santa Cruz Biotechnology, sc-8002), or anti-β-ACTIN (Santa Cruz Biotechnology, sc-47778) antibody, and then washed for 30 min with TBST. The membrane was then incubated for 1 h with horseradish peroxidase-conjugated secondary goat anti-rabbit or anti-mouse antibody (diluted 1:1,000 in 5% non-fat milk) before washing for 30 min with TBST. Signals were detected using an ECL solution (Promega).
Cloning and siRNA Knockdown
A UBE2C cDNA expression construct was generated by amplifying the region encoding the primary transcript from human UBE2C mRNA (NM_007019.4). The PCR products were cloned into the N-terminal pEF/HIS A vector (Addgene). The full sequence of wild-type UBE2C cDNA was confirmed by Sanger sequencing. UBE2C mRNA knockdown was performed using On-TARGET plus Human UBE2C specific siRNA (Dharmacon). Four target sequences of UBE2C were used for mixed siRNA.
Cell Viability, Proliferation, and Migration Assays
Cell viability was assessed using the live cell imaging tool (IncuCyte ZOOM, Essen BioScience). For proliferation assays using estrogen, equal numbers of cells were plated in 12-well plates and treated with 17-beta-estradiol (E2) added to the media at 10 nM every 48 h after siRNA treatment. The total number of cells was counted. The migration of MCF-7 and T47D cells was analyzed using Transwell chambers (Corning) and the wound healing assay in T47D cells. Treated and control cells were serum starved for 24 h and plated in equal numbers in Transwell chambers. After 24 h of incubation, cells that migrated to the bottom surface were detached and lysed using a cell detachment solution and 4× cell lysis buffer (Millipore). The cell lysate was incubated with CyQuant GR Dye 1 (Millipore), and fluorescence was measured at 485/535 nm. Cell migration was assessed using the wound healing assay in T47D cells. Equal numbers of treated and control cells were cultured to 90–95% confluency in 96-well plates. Cell monolayers were wounded using a sterilized wound maker (Cell player), and gap recovery was measured using live cell imaging (IncuCyte ZOOM, Essen BioScience).
Clonogenic Assay in Two-Dimensions (2D) and Three-Dimensions (3D), Apoptosis Assay, and Cell Cycle Distribution Analysis
For the clonogenic assay in 2D, cell numbers were assessed indirectly by staining with crystal violet because the total amount of cell-bound dye is proportional to cell number. Equal numbers of treated and control cells were seeded in six-well plates at a density of 5,000 and 10,000 cells per well in medium supplemented with 10% FBS and 1% antibiotics for 14 days. Viable cells were stained with 0.5% crystal violet dye in 25% methanol for 5 min. After staining, the fixed cells were washed with phosphate-buffered saline and photographed. The dye was eluted with 10% acetic acid, and absorbance was measured at 570 nm. For the clonogenic assay in 3D, cells were cultured in Matrigel (BD Science). The bottom layer of each well was coated with 45 μL of Matrigel and allowed to form a gel by incubation for 30 min at 37°C. Then, 20,000 cells resuspended in 200 μL of matrigel were loaded onto the bottom layer of each well. Complete medium was then overlaid onto the gel and replaced every 2 days. The images of cell colonies were obtained after culturing for 7 days. Cell apoptosis was determined using the FITC Annexin V Apoptosis Detection Kit I (BD Biosciences) following treatment with 1 μM staurosporine (STS, Sigma-Aldrich). Cell apoptosis was analyzed with a FACSAria system (BD Biosciences). Harvested cells were stained with propidium iodide (PI) and subjected to flow cytometric analysis (BD FACSAria).
MCF-7 cells were treated with siRNA against UBE2C and tamoxifen, and growth inhibition was measured using the WST assay. MCF-7 cells were seeded into six-well plates and treated with tamoxifen at 48 h after siRNA transfection. The pharmacological interaction between UBE2C-specific siRNA and tamoxifen was assessed by Chou–Talalay analysis (26, 27). This method assesses synergy and antagonism by quantifying the divergence of the combination effect from the expected additive effect of the two therapeutic agents. A combination index (CI) was estimated from the dose-effect data. A CI of <1, =1, and >1 indicates synergy, additive effects, and antagonism, respectively.
Chromatin Immunoprecipitation (ChIP)-PCR
Cells were plated in a 150 mm dish at 80–90% density, and chromatin immunoprecipitation (ChIP) was performed using an EZ ChIP kit (Millipore) as described in the instruction manual. In brief, cells were fixed with 10 ml formaldehyde (37%, Sigma) for 10 min to crosslink proteins to DNA. Cell suspensions were then sonicated (Bioruptor), and the chromatin was sheared to a DNA size of ~200–1,000 base pairs. The cross-linked protein/DNA complexes were subjected to immunoprecipitation with anti-ERα (Santa Cruz Biotechnology, sc-8002) and normal mouse IgG as the negative control. Protein/DNA complexes were purified using a spin column. Following purification of associated DNA, semi-quantitative PCR was performed to quantify the DNA sequence associated with ERα using specific UBE2C promoter primers. The primer sequences used for ChIP-PCR assay were as follows (Fwd: forward, Rev: reverse): GREB1 (positive control, Fwd: 5′-CAGCTGACTGTCTTCCACCA-3′, Rev: 5′-CCACCGTTTCGTGTCTTCTT-3′); GAPDH (negative control, Fwd: 5′-TACTAGCGGTTTTACGGGCG-3′, Rev: 5′-TCGAACAGGAGGAGCAGAGAGCGA-3′) (28); and UBE2C (Fwd: 5′-ACCAATCGGTTGTCAGAAGC-3′, Rev: 5′-GGCAGAGAGACAGGAACTCG-3′).
Implications of UBE2C mRNA Expression in the Survival of Patients With HR+/HER2– Breast Cancer
The clinical significance of UBE2C in breast cancer was investigated by analyzing UBE2C mRNA expression in invasive breast cancer tissues. In a previous study, we performed subgroup analysis according to molecular breast cancer subtype using FFPE samples from 819 breast cancer patients. The results of univariate analysis showed that UBE2C mRNA expression as a continuous variable is statistically significantly associated with DFS, DMFS, and OS in only patients with HR+/HER2– breast cancer regardless of LN negativity or positivity (10). Of the previous cohort, 410 patients with HR+/HER2– breast cancer were included in this study; the detailed clinical characteristics of patients were shown in Table 1. Briefly, the median patient age was 47.3 years (range, 25.2–80.5), and 49.5% (203/410) were LN negative. Of 410 patients, 342 (83.4%) received adjuvant chemotherapy, and histologic grade 2 was predominant in the patient group (48.5%, 199/410).
Patients were classified into low- or high-UBE2C groups according to a cut-off for normalized expression of UBE2C mRNA of 29.05; 60.5% (n = 248) were assigned to the high-expression group and the remaining patients were classified as the low-expression group. The association of classical clinicopathological factors and UBE2C gene expression levels with clinical outcomes was analyzed. Univariate analysis of DFS, DMFS, and OS indicated that positive LN involvement, large tumor size, high histologic grade, and high UBE2C expression were significantly associated with poor survival in all patients (P < 0.05, Supplementary Table 1). However, univariate, Kaplan–Meier, and multivariate analyses of subgroups according to LN status revealed that high UBE2C expression was associated with significantly shorter DFS, DMFS, and OS in breast cancer with pN0 and pN1 tumors but not in pN2/N3 tumors (P < 0.05, Supplementary Table 1, Figure 1, and Table 2). The results indicate that high UBE2C expression has an unfavorable impact on DFS, DMFS, and OS only in patients with HR+/HER2– pN0 and pN1 breast cancer.
Figure 1. UBE2C mRNA expression is strongly associated with poor survival in patients with HR+/HER2– breast cancer. The impact of UBE2C mRNA expression on DFS (left panels), DMFS (middle panels), and OS (right panels) of (A) all, (B) pN0, (C) pN1, or (D) pN2/N3 patients with HR+/HER2– breast cancer was analyzed using the Kaplan–Meier method.
Table 2. Multivariate analysis of the effect of UBE2C mRNA expression on disease-free survival (DFS), distant metastasis-free survival (DMFS), and overall survival (OS) in patients with HR+/HER2– breast cancer.
UBE2C Expression Status in HR+/HER2– Breast Cancer Tissues and Cell Lines
On the basis of the association between UBE2C expression and clinical outcomes in HR+/HER2– breast cancer, we examined the expression status of UBE2C in different breast cancer subtypes. The mRNA expression levels of UBE2C in four breast cancer subtypes were compared using 819 FFPE tissues from patients with different LN status. The median UBE2C mRNA expression was lower in HR+/HER2– breast cancers with pN0 and pN1 tumors than in other subtypes (Figure 2A). However, UBE2C mRNA expression in pN2/N3 tumors did not differ between HR+/HER2– and HR+/HER2+ or HR–/HER2+ subtypes excluding TNBC (Figure 2A). Among patients with HR+/HER2– breast cancer, UBE2C expression did not differ between pN0 and pN1 tumors, whereas it was higher in pN2/N3 tumors than in pN0 and pN1 tumors (Figure 2B). HR+/HER2– breast cancer cell lines expressed the UBE2C protein at detectable levels (Figure 2C).
Figure 2. UBE2C is expressed in HR+/HER2– breast cancer tissues and cell lines. (A) UBE2C mRNA expression in FFPE breast cancer tissue was evaluated in a dataset of 819 samples from our cohort according to LN status, pN0 (left panels), pN1 (middle panel), and pN2/N3 (right panel). (B) UBE2C mRNA expression in FFPE breast cancer tissue was evaluated in a dataset of 410 samples from patients with HR+/HER2– breast cancer according to LN status, pN0, pN1, and pN2/N3. The bar indicates the median value. (C) UBE2C protein expression was quantified by immunoblotting in breast cancer cell lines. *P < 0.05, **P < 0.01, Student's t-test.
Effects of UBE2C Depletion or Overexpression on the Tumorigenicity of HR+/HER2– Breast Cancer Cells
To examine whether UBE2C regulates the growth of HR+/HER2– breast cancer cells, we used T47D cells: a HR+/HER2– breast cancer cell line with high UBE2C expression level. T47D cells were transfected with UBE2C siRNA, and UBE2C knockdown was confirmed by immunoblotting (Figure 3A, left panel). Loss of UBE2C expression decreased cell proliferation (Figure 3A, middle and right panels) and migration (Figure 3C). Next, an HR+/HER2– cell line with moderate UBE2C expression level, MCF-7, was transfected with HA (Hemagglutinin, control) or UBE2C-expressing cDNA, and UBE2C overexpression was confirmed by immunoblotting (Figure 3B, upper left panel). Overexpression of UBE2C increased cell proliferation in 2D (Figure 3B, lower left and upper right panels) and 3D (Figure 3B, lower right panel) cell cultures and increased migration (Figure 3D). These findings indicate that UBE2C expression contributes to the tumorigenicity of HR+/HER2– breast cancer cells.
Figure 3. UBE2C contributes to the tumorigenicity of HR+/HER2– breast cancer cells. (A) The efficacy of siRNA-mediated UBE2C knockdown was evaluated by immunoblotting in T47D cells (left panel). NT, non-targeting; N, siNT; U, siUBE2C. The effects of UBE2C knockdown on cell proliferation were evaluated using an IncuCyte analyzer (middle panel) and clonogenic assay (right panel). (B) Ectopic UBE2C overexpression was confirmed by immunoblotting in MCF-7 cells (upper left panel). The effects of UBE2C overexpression on cell proliferation were evaluated using an IncuCyte analyzer (lower left panel) and clonogenic assay in 2D (upper right panel) and 3D Matrigel (100×, lower right panel). (C) The effects of UBE2C knockdown on cell migration were evaluated with the wound healing assay using the IncuCyte analyzer in T47D cells. (D) The effects of ectopic UBE2C overexpression on cell migration were evaluated with the Boyden chamber assay in MCF-7 cells. n = 3; Bars, SE; *P < 0.05, **P < 0.01, Student's t-test.
The Growth of MCF-7 Cells Is Regulated by UBE2C Expression in an Estrogen-Dependent/Independent Manner
Because tamoxifen downregulated UBE2C expression in MCF-7 cells (Figure 5C), we examined the effect of estrogen on UBE2C expression in HR+/HER2– breast cancer cells. MCF-7 (Figure 4A) and T47D (Supplementary Figure 1) cells were treated with E2, and UBE2C expression was monitored at the indicated time points. Exposure to E2 upregulated UBE2C at the mRNA and protein levels in MCF-7 and T47D cells.
Figure 4. UBE2C is a target of ERα and is required for estrogen-induced cell proliferation. (A) The induction of UBE2C mRNA (upper panel) and protein (lower panel) expression by estrogen in MCF-7 cells was quantified by qRT-PCR and immunoblotting, respectively. (B) The ERα-bound UBE2C promoter region was assayed by conventional ChIP-PCR. The GREB1 and GAPDH promoters served as ER α-bound positive and negative controls, respectively. (C) MCF-7 cells were transfected with UBE2C siRNAs for 48 h, and estrogen was added to the medium at 48 h intervals. Cell proliferation was evaluated using an IncuCyte analyzer (upper panel). MCF-7 cells were transfected with UBE2C siRNAs for 48 h, and estrogen was added to the medium for 24 h. Then, cells were stained with propidium iodide, and DNA content and intensity were analyzed by flow cytometry. (D) Hormone-independent cell proliferation of UBE2C-overexpressing MCF-7 cells was monitored by clonogenic assay under DMEM (upper and lower panels) or charcoal stripped DMEM (middle and lower panels) culture conditions. n = 3; Bars: SE. *P < 0.05, **P < 0.01, Student's t-test.
The role of estrogen receptor α (ERα) in the regulation of UBE2C mRNA expression in HR+/HER2– breast cancer was evaluated by analyzing the recruitment of ERα to the promoter region of UBE2C in MCF-7 cells by ChIP-PCR (Figure 4B). A search of the public database (https://chip-atlas.org) was performed to predict the putative binding region of ERα in the UBE2C promoter, and specific sequences were selected (Supplementary Figure 2). ChIP analysis revealed increased ERα binding to UBE2C and GREB1, an ERα target gene, but not to GAPDH, a negative ERα target gene, in the presence of estrogen.
Given our observations suggesting direct regulation of UBE2C mRNA expression by ERα binding, we assessed the effect of UBE2C on cell proliferation in response to estrogen. MCF-7 cells were transfected with UBE2C siRNA and treated with E2, and cell growth and cell cycle progression were assessed by live cell imaging and flow cytometry, respectively (Figure 4C). Cell proliferation was significantly lower in UBE2C knockdown cells exposed to estrogen than in those cultured in the absence of estrogen. Cell cycle analysis showed that UBE2C depletion caused cell cycle arrest at the S and G2/M phases of the cell cycle, consistent with previous studies by Shen et al. (29) and Palumbo et al. (20). UBE2C knockdown strongly induced arrest at S and G2/M phases in MCF-7 cells exposed to estrogen compared with UBE2C knockdown in cells not exposed to estrogen. Taken together, these results show that loss of UBE2C inhibits cell proliferation by affecting cell cycle progression in HR+/HER2– breast cancer cells treated with estrogen.
The tumorigenicity of UBE2C led us to hypothesize that sustained UBE2C overexpression promotes estrogen-independent tumor growth. The proliferation rate of UBE2C-overexpressing MCF-7 cells cultured in estrogen-deprived medium was markedly higher than that of cells cultured in DMEM for a long period (Figure 4D). These data suggest that UBE2C overexpression contributes to estrogen-independent growth in HR+/HER2– breast cancer cells.
Effects of UBE2C Depletion in Combination With Tamoxifen on Apoptosis of HR+/HER2– Breast Cancer Cells
To evaluate the effect of UBE2C depletion in combination with pharmacological inhibition with tamoxifen, MCF-7 (Figure 5A) and T47D (Supplementary Figure 3) cells were transfected with UBE2C siRNA and treated with tamoxifen, and cell growth/viability was measured. Tamoxifen inhibited the growth of HR+/HER2– breast cancer cells, with inhibitory concentration (IC50) values of 8.26 μM and 7.61 μM in MCF-7 and T47D cells, respectively. Next, MCF-7 and T47D cells were transfected with various concentrations of UBE2C siRNA, followed by tamoxifen treatment at a dose corresponding to the IC50 value. The Chou–Talalay method was used to assess the synergistic effect of tamoxifen and UBE2C siRNA. The CI was 0.3472 ± 0.0046 in MCF-7 cells and 0.8200 ± 0.0083 in T47D cells, demonstrating the synergistic effect of tamoxifen and UBE2C siRNA. MCF-7 cells were more sensitive to combination treatment with tamoxifen and UBE2C siRNA than were T47D cells. These data suggest that tamoxifen acts synergistically with UBE2C depletion to decrease cell viability in HR+/HER2– breast cancer cells expressing UBE2C.
Figure 5. UBE2C depletion in combination with tamoxifen increases apoptosis of HR+/HER2– breast cancer cells. (A) Chou–Talalay drug interaction analysis in MCF-7 cells treated with UBE2C siRNA and tamoxifen. (B) The apoptotic potential of treated cells was determined by flow cytometry and Annexin V assay. The cells were stained with propidium iodide and Annexin V antibody, and DNA content and intensity were analyzed by flow cytometry. Numbers indicate the percentage of apoptotic cells. (C) The viability of treated MCF-7 cells was evaluated using a 3D Matrigel assay (40×). (D) Lysates from treated MCF-7 cells were subjected to western blotting to monitor the protein levels of cleaved PARP and cleaved CASPASE 3. Fold changes were determined by comparing the protein levels with those of β-ACTIN using ImageJ. Tam, tamoxifen.
To assess the effect of tamoxifen and UBE2C depletion on cell death, MCF-7 cells were transfected with UBE2C siRNA and treated with tamoxifen, and cell death was analyzed by flow cytometry (Figure 5B), 3D cell culture (Figure 5C), and immunoblotting (Figure 5D). UBE2C siRNA alone slightly increased the population of apoptotic cells compared with that in cells treated with non-targeting control siRNA, and tamoxifen markedly increased the population of apoptotic cells compared with the effect of UBE2C siRNA. Combination treatment with UBE2C siRNA and tamoxifen significantly increased apoptosis. The size and number of cell colonies was smaller in cells treated with UBE2C siRNA and tamoxifen than in those exposed to UBE2C siRNA alone or tamoxifen alone. The expression of cleaved PARP and cleaved CASPASE 3 was higher in cells treated with UBE2C siRNA and tamoxifen than in those exposed to UBE2C siRNA alone or tamoxifen alone. These findings suggest that UBE2C depletion markedly increased the cytotoxicity of tamoxifen by inducing apoptosis in HR+/HER2– breast cancer cells.
Many studies report that high UBE2C expression is associated with poor survival in patients with breast cancer irrespective of subtype (22, 24, 25, 30, 31). However, we previously showed that high UBE2C mRNA expression is significantly associated with poor DFS, DMFS, and OS only in patients with HR+/HER2– breast cancer among the four subtypes (10). This finding indicated that the association between UBE2C overexpression and the survival of patients with HR+/HER2– breast cancer might be estrogen related; thus, in this study, we assessed the effect of estrogen on UBE2C expression using HR+/HER2– breast cancer cells.
A literature review was performed to determine the correlation between the mRNA and protein expression of UBE2C. The results showed a correlation value of UBE2C mRNA-to-protein of 0.7127 (32), which is considered acceptable. Although the expression of UBE2C was only measured at the mRNA level and not at the protein level, we further assessed the clinical significance of UBE2C overexpression in HR+/HER2– breast cancer. The results of univariate and multivariate analyses showed that high UBE2C mRNA expression was strongly associated with poor DFS, DMFS, and OS only in patients with HR+/HER2– pN0 and pN1 (early) breast cancer but not in those with pN2/N3 (late) breast cancer. However, UBE2C mRNA expression was lower in pN0 and pN1 tumors from HR+/HER2– breast cancer patients than in tumors from other subtype patients. These findings indicate that the effect of UBE2C overexpression on the survival of patients with early HR+/HER2– breast cancer may be estrogen related, and UBE2C expression can predict the risk of distant recurrence among patients with early HR+/HER2– breast cancer. In addition, UBE2C mRNA expression was not associated with clinical outcomes in HR+/HER2– pN2/N3 breast cancer and other more aggressive subtypes (HR+/HER2+, HR–/HER2+, and TNBC). This finding suggests that although UBE2C expression correlates with worse prognosis of early HR+/HER2– breast cancer in this study, since most highly malignant breast cancer cases also show high UBE2C expression, it implies that the expression does not hold a prognostic value.
Consistent with the results of two studies by Takeo et al. and Chou et al. in MCF-7 cells (23, 33), functional in vitro experiments based on UBE2C overexpression or depletion showed that aberrantly high levels of UBE2C promoted cell proliferation and migration. These previous studies were based on literature reviews on the importance of UBE2C tumorigenicity as an essential part of the ubiquitin-conjugating enzyme complex, and they used the MCF-7 cell line as a representative breast cancer cell line. However, this study analyzed the association between UBE2C expression and clinical outcomes only in HR+/HER2– breast cancer through a parametric survival analysis using public datasets and a patient cohort. The MCF-7 cell line was used as a representative HR+/HER2– breast cancer cell line. Although MCF-7 cells are widely used in research, cancer cell lines are different from patient-derived cancer cells. We therefore confirmed the results obtained in MCF-7 cells by using a 3D cell culture system.
Estrogen is a key regulator of female development and reproductive function, and it is a causal factor in breast cancer. Estrogen-driven gene expression is mediated by the activity of ERα in normal breast cells and in breast cancer (34). We found that estrogen exposure highly upregulated UBE2C mRNA and protein, and that ERα directly bound to the UBE2C promoter region. Loss of UBE2C expression markedly decreased cell proliferation by affecting cell cycle progression under estrogen exposure conditions compared with UBE2C knockdown cells without estrogen exposure in HR+/HER2– breast cancer cells. These findings suggest that UBE2C is a target of ERα and is required for estrogen-induced cell proliferation, thereby activating estrogen-mediated signaling networks. This may explain why high UBE2C expression was strongly associated with poor DMFS, DFS, and OS only in patients with HR+/HER2– breast cancer. In addition, we found that UBE2C mRNA expression levels were high in certain triple negative cell lines. Nicolau-Neto et al. reported that UBE2C is a transcriptional target of the cell cycle regulator FOXM1, which was identified as a specific marker for TNBC (35, 36).
Tamoxifen is currently used for the management of ER+ breast cancers. We found that tamoxifen increased the apoptotic population of MCF-7 cells. The cytotoxicity of tamoxifen is attributed to the induction of apoptosis in breast cancer cells (37, 38). Despite its antitumor efficacy in the clinic, 20–30% of patients with ER+ breast cancer are resistant to tamoxifen treatment (39). However, the mechanisms underlying resistance to tamoxifen treatment are complex and not fully understood (40–42). Dysregulation of apoptosis or cell cycle regulators is suggested as one of several mechanisms of resistance to endocrine therapy (41). In this study, we showed that UBE2C is involved in cell cycle regulation, and that UBE2C overexpression contributes to estrogen-independent growth in HR+/HER2– breast cancer cells. The present findings indicate that UBE2C overexpression may induce resistance to tamoxifen treatment by promoting estrogen-independent growth. UBE2C depletion markedly increased the cytotoxicity of tamoxifen by inducing apoptosis in HR+/HER2– breast cancer cells. This suggests that tamoxifen treatment combined with inhibitors targeting UBE2C could be an effective strategy to increase the inhibitory effect of tamoxifen and overcome resistance in HR+/HER2– breast cancer. The combination of palbociclib, a selective inhibitor of the cyclin-dependent kinases CDK4 and CDK6, and ovarian suppression with letrozole significantly prolonged progression-free survival in women with ER+/HER2– advanced breast cancer in the PALOMA-3 (NCT01942135) study (43). Therefore, we suggest that the combination of palbociclib with tamoxifen is a promising treatment strategy in patients with HR+/HER2– breast cancer overexpressing UBE2C.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
The studies involving human participants were reviewed and approved by the Samsung Medical Center (SMC) Institutional Review Board (IRB file No. 2008-12-035). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Y-JK and YS: conception and design and final approval of the article. Y-JK, GL, JH, KS, J-SC, Y-LC, and YS: analysis and interpretation. Y-JK, GL, and JH: data collection. Y-JK and JH: writing the article.
This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean Government (grant number 2016R1A6A3A11935830). We are very grateful for the financial support of the Research Institute of Pharmaceutical Sciences, Seoul National University College of Pharmacy.
Conflict of Interest
JH was employed by the company Gencurix, Inc.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01574/full#supplementary-material
Supplementary Figure 1. UBE2C is induced by estrogen in T47D cells. The induction of UBE2C mRNA (left panel) and protein (right panel) expression by estrogen in T47D cells was quantified by qRT-PCR and immunoblotting, respectively. Fold changes were determined by comparing the protein levels with those of β-ACTIN using ImageJ.
Supplementary Figure 2. The putative binding region of ERα on the UBE2C promoter was determined by a review of https://chip-atlas.org. (A) A binding region of ERα on the promoter from the SRX1531791 ChIP-sequencing dataset using Integrative Genomics Viewer. (B) Design of the putative UBE2C promoter region ChIP-PCR primer pair.
Supplementary Figure 3. Chou–Talalay drug interaction analysis of UBE2C siRNA and tamoxifen in T47D cells.
Supplementary Table 1. Univariate analysis of the effect of UBE2C mRNA expression on disease-free survival (DFS), distant metastasis-free survival (DMFS), and overall survival (OS) in patients with HR+/HER2– breast cancer.
ER, Estrogen receptor; PR, Progesterone receptor; HER2, Human epidermal growth factor receptor 2; FFPE, Formalin-fixed paraffin-embedded; DMFS, Distant metastasis-free survival; DFS, Disease-free survival; OS, Overall survival.
1. Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thurlimann B, Senn HJ. 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:1736–47. doi: 10.1093/annonc/mdr304
2. Kaufmann M, Pusztai L. Use of standard markers and incorporation of molecular markers into breast cancer therapy: consensus recommendations from an International Expert Panel. Cancer. (2011) 117:1575–82. doi: 10.1002/cncr.25660
5. Early Breast Cancer Trialists' Collaborative Group (EBCTCG). Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. (2005) 365:1687–717. doi: 10.1016/S0140-6736(05)66544-0
6. Van De Vijver MJ, He YD, Van't Veer LJ, Dai H, Hart AA, Voskuil DW, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. (2002) 347:1999–2009. doi: 10.1056/NEJMoa021967
7. Paik S, Tang G, Shak S, Kim C, Baker J, Kim W, et al. Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol. (2006) 24:3726–34. doi: 10.1200/JCO.2005.04.7985
8. Filipits M, Rudas M, Jakesz R, Dubsky P, Fitzal F, Singer CF, et al. A new molecular predictor of distant recurrence in ER-positive, HER2-negative breast cancer adds independent information to conventional clinical risk factors. Clin Cancer Res. (2011) 17:6012–20. doi: 10.1158/1078-0432.CCR-11-0926
9. Oh E, Choi YL, Park T, Lee S, Nam SJ, Shin YK. A prognostic model for lymph node-negative breast cancer patients based on the integration of proliferation and immunity. Breast Cancer Res Treat. (2012) 132:499–509. doi: 10.1007/s10549-011-1626-8
10. Han J, Choi YL, Kim H, Choi JY, Lee SK, Lee JE, et al. MMP11 and CD2 as novel prognostic factors in hormone receptor-negative, HER2-positive breast cancer. Breast Cancer Res Treat. (2017) 164:41–56. doi: 10.1007/s10549-017-4234-4
11. Gong G, Kwon MJ, Han J, Lee HJ, Lee SK, Lee JE, et al. A new molecular prognostic score for predicting the risk of distant metastasis in patients with HR+/HER2– early breast cancer. Sci Rep. (2017) 7:45554. doi: 10.1038/srep45554
12. Kwon MJ, Lee SB, Han J, Lee JE, Lee JW, Gong G, et al. BCT score predicts chemotherapy benefit in Asian patients with hormone receptor-positive, HER2-negative, lymph node-negative breast cancer. PLoS ONE. (2018) 13:e0207155. doi: 10.1371/journal.pone.0207155
13. Kwon MJ, Lee JE, Jeong J, Woo SU, Han J, Kang BI, et al. Comparison of geneswell BCT score with oncotype DX recurrence score for risk classification in asian women with hormone receptor-positive, HER2-negative early breast cancer. Front Oncol. (2019) 9:667. doi: 10.3389/fonc.2019.00667
16. Chen S, Chen Y, Hu C, Jing H, Cao Y, Liu X. Association of clinicopathological features with UbcH10 expression in colorectal cancer. J Cancer Res Clin Oncol. (2010) 136:419–26. doi: 10.1007/s00432-009-0672-7
17. Rajkumar T, Sabitha K, Vijayalakshmi N, Shirley S, Bose MV, Gopal G, et al. Identification and validation of genes involved in cervical tumourigenesis. BMC Cancer. (2011) 11:80. doi: 10.1186/1471-2407-11-80
18. Zhao ZK, Wu WG, Chen L, Dong P, Gu J, Mu JS, et al. Expression of UbcH10 in pancreatic ductal adenocarcinoma and its correlation with prognosis. Tumour Biol. (2013) 34:1473–7. doi: 10.1007/s13277-013-0671-9
19. Kim WT, Jeong P, Yan C, Kim YH, Lee IS, Kang HW, et al. UBE2C cell-free RNA in urine can discriminate between bladder cancer and hematuria. Oncotarget. (2016) 7:58193–202. doi: 10.18632/oncotarget.11277
20. Palumbo A Jr, Da Costa NM, De Martino M, Sepe R, Pellecchia S, De Sousa VP, et al. UBE2C is overexpressed in ESCC tissues and its abrogation attenuates the malignant phenotype of ESCC cell lines. Oncotarget. (2016) 7:65876–87. doi: 10.18632/oncotarget.11674
22. Psyrri A, Kalogeras KT, Kronenwett R, Wirtz RM, Batistatou A, Bournakis E, et al. Prognostic significance of UBE2C mRNA expression in high-risk early breast cancer. A Hellenic Cooperative Oncology Group (HeCOG) study. Ann Oncol. (2012) 23:1422–7. doi: 10.1093/annonc/mdr527
23. Chou CP, Huang NC, Jhuang SJ, Pan HB, Peng NJ, Cheng JT, et al. Ubiquitin-conjugating enzyme UBE2C is highly expressed in breast microcalcification lesions. PLoS ONE. (2014) 9:e93934. doi: 10.1371/journal.pone.0093934
24. Parris TZ, Kovacs A, Aziz L, Hajizadeh S, Nemes S, Semaan M, et al. Additive effect of the AZGP1, PIP, S100A8 and UBE2C molecular biomarkers improves outcome prediction in breast carcinoma. Int J Cancer. (2014) 134:1617–29. doi: 10.1002/ijc.28497
25. Mo CH, Gao L, Zhu XF, Wei KL, Zeng JJ, Chen G, et al. The clinicopathological significance of UBE2C in breast cancer: a study based on immunohistochemistry, microarray and RNA-sequencing data. Cancer Cell Int. (2017) 17:83. doi: 10.1186/s12935-017-0455-1
26. Chou TC, Talalay P. Quantitative analysis of dose-effect relationships: the combined effects of multiple drugs or enzyme inhibitors. Adv Enzyme Regul. (1984) 22:27–55. doi: 10.1016/0065-2571(84)90007-4
27. Chou TC, Motzer RJ, Tong Y, Bosl GJ. Computerized quantitation of synergism and antagonism of taxol, topotecan, and cisplatin against human teratocarcinoma cell growth: a rational approach to clinical protocol design. J Natl Cancer Inst. (1994) 86:1517–24. doi: 10.1093/jnci/86.20.1517
28. Caizzi L, Ferrero G, Cutrupi S, Cordero F, Ballare C, Miano V, et al. Genome-wide activity of unliganded estrogen receptor-alpha in breast cancer cells. Proc Natl Acad Sci USA. (2014) 111:4892–7. doi: 10.1073/pnas.1315445111
29. Shen Z, Jiang X, Zeng C, Zheng S, Luo B, Zeng Y, et al. High expression of ubiquitin-conjugating enzyme 2C (UBE2C) correlates with nasopharyngeal carcinoma progression. BMC Cancer. (2013) 13:192. doi: 10.1186/1471-2407-13-192
30. Loussouarn D, Campion L, Leclair F, Campone M, Charbonnel C, Ricolleau G, et al. Validation of UBE2C protein as a prognostic marker in node-positive breast cancer. Br J Cancer. (2009) 101:166–73. doi: 10.1038/sj.bjc.6605122
31. Garcia-Escudero R, Martinez-Cruz AB, Santos M, Lorz C, Segrelles C, Garaulet G, et al. Gene expression profiling of mouse p53-deficient epidermal carcinoma defines molecular determinants of human cancer malignancy. Mol Cancer. (2010) 9:193. doi: 10.1186/1476-4598-9-193
32. Mertins P, Mani DR, Ruggles KV, Gillette MA, Clauser KR, Wang P, et al. Proteogenomics connects somatic mutations to signalling in breast cancer. Nature. (2016) 534:55–62. doi: 10.1038/nature18003
33. Fujita T, Ikeda H, Kawasaki K, Taira N, Ogasawara Y, Nakagawara A, et al. Clinicopathological relevance of UbcH10 in breast cancer. Cancer Sci. (2009) 100:238–48. doi: 10.1111/j.1349-7006.2008.01026.x
34. Carroll JS, Liu XS, Brodsky AS, Li W, Meyer CA, Szary AJ, et al. Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell. (2005) 122:33–43. doi: 10.1016/j.cell.2005.05.008
35. Nicolau-Neto P, Palumbo A, De Martino M, Esposito F, De Almeida Simao T, Fusco A, et al. UBE2C is a transcriptional target of the cell cycle regulator FOXM1. Genes. (2018) 9:188. doi: 10.3390/genes9040188
37. Zhang GJ, Kimijima I, Onda M, Kanno M, Sato H, Watanabe T, et al. Tamoxifen-induced apoptosis in breast cancer cells relates to down-regulation of bcl-2, but not bax and bcl-X(L), without alteration of p53 protein levels. Clin Cancer Res. (1999) 5:2971–7.
39. Ali S, Rasool M, Chaoudhry H, N Pushparaj P, Jha P, Hafiz A, et al. Molecular mechanisms and mode of tamoxifen resistance in breast cancer. Bioinformation. (2016) 12:135–9. doi: 10.6026/97320630012135
40. Shou J, Massarweh S, Osborne CK, Wakeling AE, Ali S, Weiss H, et al. Mechanisms of tamoxifen resistance: increased estrogen receptor-HER2/neu cross-talk in ER/HER2-positive breast cancer. J Natl Cancer Inst. (2004) 96:926–35. doi: 10.1093/jnci/djh166
42. Puyang X, Furman C, Zheng GZ, Wu ZJ, Banka D, Aithal K, et al. Discovery of selective estrogen receptor covalent antagonists for the treatment of ERalphaWT and ERalphaMUT breast cancer. Cancer Discov. (2018) 8:1176–93. doi: 10.1158/2159-8290.CD-17-1229
43. Cristofanilli M, Turner NC, Bondarenko I, Ro J, Im SA, Masuda N, et al. Fulvestrant plus palbociclib versus fulvestrant plus placebo for treatment of hormone-receptor-positive, HER2-negative metastatic breast cancer that progressed on previous endocrine therapy (PALOMA-3): final analysis of the multicentre, double-blind, phase 3 randomised controlled trial. Lancet Oncol. (2016) 17:425–39. doi: 10.1016/S1470-2045(15)00613-0
Keywords: UBE2C, HR+/HER2– breast cancer, estrogen, cell cycle, tamoxifen, apoptosis
Citation: Kim Y-J, Lee G, Han J, Song K, Choi J-S, Choi Y-L and Shin YK (2020) UBE2C Overexpression Aggravates Patient Outcome by Promoting Estrogen-Dependent/Independent Cell Proliferation in Early Hormone Receptor-Positive and HER2-Negative Breast Cancer. Front. Oncol. 9:1574. doi: 10.3389/fonc.2019.01574
Received: 30 August 2019; Accepted: 30 December 2019;
Published: 23 January 2020.
Edited by:Carmine De Angelis, University of Naples Federico II, Italy
Reviewed by:Tomohiro Chiba, Kyorin University, Japan
Mark Opdam, The Netherlands Cancer Institute (NKI), Netherlands
Copyright © 2020 Kim, Lee, Han, Song, Choi, Choi and Shin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.