Abstract
Estrogen receptor (ER) β1 and its splice variants are expressed both in ovary and ovarian cancer. We studied the role of ERβ1 and two of its splice variants in regulation of gene expression, cellular proliferation, apoptosis, and migration of an ovarian cancer cell line. In this study, we transfected SK-OV-3 ovarian cancer cells with vectors coding for ERβ1 or its splice variants ERβ-δ125 and ERβ-δ1256, and tested their response to estrogen and tamoxifen in comparison with the untransfected cells. Heterologous expression of ERβ1, but not of the exon-deleted ERβ variants resulted in notably slower cell growth of SK-OV-3 ovarian cancer cells, an effect accompanied by more than tenfold increase of cyclin-dependent kinase inhibitor p21(WAF1) transcript levels and a significant reduction of cyclin A2 mRNA levels. SK-OV-3 cells stably overexpressing ERβ1 ligand independently also exhibited an increased apoptosis rate and a significantly decreased motility, an effect accompanied by upregulation of fibulin 1c. Our data demonstrate that ERβ1, but not the exon-deleted isoforms tested exerts multiple antitumoral effects on SK-OV-3 ovarian cancer cells even in the absence of estradiol or functional ERα.
Introduction
Although 40–60% of ovarian cancers express estrogen receptor (ER) α (Greenlee et al. 2000, Havrilesky et al. 2001), only a minor proportion of patients (ranging from 7 to 18%) respond clinically to treatment with selective ER modulator tamoxifen (Hatch et al. 1991, Scambia et al. 1995). However, the role of estrogens has been recently highlighted by the results of three large prospective studies showing that estradiol uptake in postmenopausal stage increased the risk of ovarian cancer incidence and mortality in women who used long-term estrogen replacement therapy (Rodriguez et al. 2001, Lacey et al. 2002, Anderson et al. 2003). Estrogen effects are mediated by two ER types, named ERα and ERβ (Gustafsson 1999, Pettersson & Gustafsson 2001, Osborne & Schiff 2005). Although particularly the molecular mechanisms of ERβ function in ovary and ovarian cancer are still poorly elucidated, it is becoming increasingly clear that both receptor types are responsible for different biological functions, as indicated by their specific expression patterns and different effects of their gene knockout (Merchenthaler & Shugrue 1999, Couse et al. 2000). Besides their different physiological functions, recent studies have suggested that ERβ, in contrast to ERα, might act as a tumor suppressor in breast or prostate cancer cells (Lazennec et al. 2001, Cheng et al. 2004), whereas other studies did not come to such conclusions (Burns et al. 2003). Given that ERβ is able to counteract ERα signaling in some settings, loss of ERβ is thought to enhance ERα-mediated proliferation of hormone-dependent cancer cells (Lindberg et al. 2003). Furthermore, recent studies suggested that ERβ signaling might affect cellular apoptosis (Cheng et al. 2004). An interesting feature of both ERs is the variety of their mRNA isoforms resulting from differential splicing (Price et al. 2000, 2001, Speirs et al. 2000, Poola et al. 2002a,b, Herynk & Fuqua 2004). The so far identified ERβ splice variants are characterized by alternative 3′-exons (ERβ2, ERβ3, ERβ4, ERβ5) or by deletion of single or multiple exons (e.g., ERβΔ2, ERβΔ5/6). Some of these mRNA isoforms were demonstrated to code for ERβ proteins, which are characterized by impaired estrogen or DNA binding or altered cofactor interaction (Sierens et al. 2004, Zhao et al. 2005). The emerging picture of multiple ERβ mRNA isoforms, and thus also the multitude of differentially built proteins, strongly suggests their synthesis to be considered as another level of complexity of estrogen signaling.
A loss of ERβ expression or a decrease in ERβ/ERα ratio in epithelial ovarian cancer as compared with normal tissues has been reported consistently by several groups (Pujol et al. 1998, Rutherford et al. 2000). However, the role of ERβ and particularly its splice variants in ovarian carcinogenesis is not fully understood and the suggested role of ERβ as a tumor suppressor raised from observations on breast and prostate cancer cells (Merchenthaler & Shugrue 1999, Lazennec et al. 2001) has to be tested with regard to ovarian cancer. In this study, we engineered SK-OV-3 ovarian cancer cells heterologously expressing ERβ1 or the exon-skipped ERβ splice variants ERβ-δ125 and ERβ-δ1256 recently identified by our group (Treeck et al. 2007). The predicted proteins coded by these novel ERβ isoforms lack the activation function 1 (AF-1) domain and have large deletions both in the ligand-binding domain (LBD) and the DNA-binding domain (DBD), and thus are expected to exhibit a drastically changed function profile in comparison with ERβ1. In this study, we examined to what extent expression of ERβ1 and two of its splice variants are able to modulate cellular proliferation, apoptosis, motility, and gene expression of ERα-negative, estrogen unresponsive SK-OV-3 ovarian cancer cells.
Materials and Methods
Materials
Phenol red-free Dulbecco’s Modified Eaglo’s Medium (DMEM) culture medium was obtained from Invitrogen; Fetal Calf Serum (FCS) was purchased from PAA Laboratories GmbH (Pasching, Austria). 17-β Estradiol (E2), 4-OH tamoxifen (4-OH TAM), ICI 182 780, staurosporine, and serum replacement 2 (SR2) were obtained from Sigma, SK-OV-3 and OVCAR-3 ovarian cancer cells were obtained from American Type Culture Collection (Manassas, VA, USA). M-MLV-P reverse transcriptase, Cell Titer Blue kit, Caspase-Glo 3/7 kit and ImProm-II Reverse Transkriptase were purchased from Promega. RNeasy Mini Kit, RNase-Free DNase Set, and Quantitect SYBR Green PCR Kit were obtained from Qiagen. PCR primers were synthesized at Metabion (Planegg-Martinsried, Germany). Transfectin reagent was obtained from Bio-Rad. Platinum Pfx Polymerase and OptiMEM medium were purchased from Invitrogen. Annexin V-FLUOS Staining Kit was obtained from Roche. Rapid-Scan gene expression panel was obtained from Origene (Rockville, MD, USA). siRNAs were obtained from Ambion (Austin, TX, USA).
Plasmids
Vector pTARGET (Promega) allows cloning in Escherichia coli and additionally carries the human cytomegalovirus immediate-early enhancer/promoter region to promote constitutive expression of cloned DNA inserts in mammalian cells. This vector also contains the neomycin phosphotransferase gene, a selectable marker for mammalian cells. pTARGET derivatives containing ORFs of ERβ1, ERβ-δ125, or ERβ?δ1256 were used for heterologous expression in SK-OV-3 cells. Vector pEGFP-N2 (Clontech) codes for the green fluorescent protein (GFP) for visualization of transfection efficacy using a fluorescence microscope. Vector pTAL-SEAP (Clontech) constitutively codes for the secreted alkaline phosphatase (SEAP) protein and served as positive control for the SEAP assay and the pTAL-estrogen response element (ERE)-SEAP is a reporter gene vector containing EREs in the promotor of the SEAP gene. Both vectors were used for the reporter gene assays performed in this study. Vector pSV-β-GAL (Promega) constitutively codes for the β-galactosidase enzyme and was used as internal control for transfection efficacy in the reporter gene assays.
Cell culture, transfections, and siRNA
SK-OV-3 cells were maintained in phenol red-free DMEM/F12 medium supplemented with 10% FCS. Cells were cultured with 5% CO2 at 37 °C in a humidified incubator. For transfection, 4 × 105 SK-OV-3 cells per well of a 6-well dish were seeded in DMEM/F12 10% FCS. The next day, 2 ml fresh culture medium was added to the cells and the transfection solution was prepared by mixing 5 μl Transfectin reagent (Bio-Rad) and 1 μg plasmid DNA or 30 nM siRNA in OptiMEM reduced serum medium (Invitrogen) and added to the cultured cells. For generation of stable clones, G418 selection (300 μg/ml) was started 48 h after transfection. For analysis of mRNA levels in siRNA-treated cells or for subsequent proliferation or apoptosis assays of siRNA-treated cells, transfected cells were harvested 24–48 h later. The siRNA sequence for knockdown of ERβ1 was 5′-CCUUACCUGUAAACAGAGAtt-3′, the sequence of the scrambled negative-control siRNA was 5′-CCAGAUU-CAGACCAAAUGUtt-3′ (Ambion).
RT and PCR
Total RNA was isolated by means of the RNeasy kit (Qiagen) according to the manufacturer’s instructions. From 1 μg total RNA, cDNA was synthesized using 100 U M-MLV-P reverse transcriptase (Promega), 2.5 mM dNTP mixture, and 50 pM random primers (Invitrogen). For detection of ERβ splice variants by standard RT-PCR, 2 μl cDNA was amplified in a reaction mix of 1 U platinum polymerase (Invitrogen), 20 pmol of each primer, 1 × PCR-buffer, 1.5 mM MgCl2, and 2.5 mM of each dNTP. The cDNA was amplified in 35 cycles (1 cycle = 1 min at 94 °C melting, 2 min at 56 °C annealing, 3 min at 72 °C extension). All PCR primers were designed intron-spanning, sequences are indicated in Table 1, position of ERβ primers is illustrated in Fig. 1.
For real-time PCR detection of ERβ isoforms or estrogen target genes, 2 μl cDNA were amplified using the Quantitect SYBR Green PCR Kit (Qiagen) and the LightCyler PCR device (Roche Diagnostics). The PCR program was 95 °C for 15 min, followed by 35 PCR cycles (95 °C for 10 s, 56 °C for 30 s, 72 °C for 30 s) and a final extension for 5 min at 72 °C, followed by a standard melting curve analysis. In all RT-PCR experiments, a 190 bp β-actin fragment was amplified as reference gene using intron-spanning primers actin-2573 and actin-2876. After performing dilution experiments with sample cDNA over a 100-fold range confirming the PCR efficiencies of all primer pairs to be approximately equal (Ståhlberg et al. 2003), data were analyzed using the comparative ΔΔCT method (Livak & Schmittgen 2001) calculating the difference between the threshold cycle (CT) values of the target and reference gene of each sample and then comparing the resulting ΔCT values between different samples. In these experiments, mRNA not subjected to RT was used as a negative control to distinguish cDNA and vector or genomic DNA amplification.
Antibodies and western blot analysis
SK-OV-3 cells were lysed in RIPA buffer 1% (v/v) Igepal CA-630, 0.5% (w/v) sodium deoxycholate, 0.1%(w/v) SDS in PBS containing aprotonin and sodium orthovanadate. Aliquots containing 15 μg proteins were resolved by 10% (w/v) SDS–PAGE, followed by electrotransfer to a PVDF hybond (Amersham) membrane. Immunodetection was carried out using ERβ antibody (GR-39, Oncogene) or β-actin antibody (8226, ABCAM, Cambridge, UK) diluted 1:5000 in PBS containing 5% skim milk (w/v) followed by horseradish peroxidase-conjugated secondary antibody, which was detected using a chemiluminescence (ECL) system (Amersham).
Cell viability assay
SK-OV-3 wild-type (WT) cells and SK-OV-3 clones cultured in DMEM containing 10% FCS or 1 × SR2 were seeded in 96-well plates in triplicates (1000 cells/well), and serum-free cultured cells were treated with 1 nM E2 alone or in combination with 4-OH TAM (0.5 or 5 μM). After 72, 96, 120, and 144 h, relative numbers of viable cells were measured in comparison with the untreated control and the solvent control using the fluorimetrical, resazurin-based Cell Titer Blue assay (Promega) according to the manufacturer’s instructions at 560Ex/590Em nm in a Victor3 multilabel counter (Perkin–Elmer, Waltham, MD, USA). Cell growth was expressed as percentage of the untreated medium control. Statistical analysis of the data was performed by one-way ANOVA using Prism 2.0 Software (Graph pad, San Diego, CA, USA), with statistical significance accepted at P < 0.05.
Apoptosis assays
SK-OV-3 WT cells and SK-OV-3 clones cultured in DMEM supplemented with 1 × SR2 (Sigma) were seeded in 96-well plates (5 × 103 cells/well) and treated with 1 nM E2, 10 μM 4-OH TAM or apoptosis inductor staurosporine (0.1 μM) as a positive control. After 6 h treatment, cellular apoptosis was determined by measurement of caspase 3 and 7 activity by means of the luminometric Caspase-Glo 3/7 assay (Promega) according to the manufacturer’s protocol using a Victor3 multilabel counter (Perkin–Elmer). Additionally, apoptosis was measured by means of the Annexin V-FLUOS Staining Kit (Roche). Cells were treated with Annexin V and propidium iodide (PI) according to the manufacturer’s protocol, and apoptotic cells exhibiting positive green Annexin V fluorescence but no red PI staining were counted. Cellular apoptosis was expressed as percentage of the untreated medium control or as percentage of the SK-OV-3 WT cells. Statistical analysis of the data was performed by one-way ANOVA using Prism 2.0 Software (Graph pad), with statistical significance accepted at P < 0.05.
Wound-healing assay
SK-OV-3 cells were plated in 6-well dishes (3 × 105 cells/well) in DMEM/F12 containing 1 × SR2 (Sigma). The next morning, cells were treated with 1 nM E2 or ethanol as negative control. After 24 h of treatment, wound-induced migration was triggered by scraping the cells with a blue tip, and the scratch was pictured immediately (day 0). The cells were pictured again 48 h later. The percentage of wound filling was calculated by computer-aided measuring of the remaining gap space on the pictures using the software Adobe Photoshop Elements 2.0.
Reporter gene assays
SK-OV-3 and OVCAR-3 WT cells were seeded in 6-well plates in DMEM/F12 supplemented with 5% FCS (4 × 105 cell per well), 5 h later serum concentration was reduced to 1% and 0.5 × serum-free SR2 medium was added. The next day, the prior to transfection medium was changed to 1 × SR2. Transfections were carried out mixing 10 μl Transfectin reagent (Bio-Rad) in a total volume of 250 μl OptiMEM medium with 5 μg pEGFP-N2 vector (Clontech) foreasy visualization of transfection efficacy using a fluorescence microscope, 5 μg pTAL-SEAP vector (Clontech) as positive control for the SEAP assay, or 10 μg reporter gene vector pTAL-ERE-SEAP (Clontech). Generally, 5 μg pSV-β-GAL vector (Promega) was added to the transfection solution serving as internal control for transfection efficacy. 24 h after adding the 250 μl transfection solution to the medium, cells were stimulated with 100 nM E2 alone or in combination with 1 μM 4-OH TAM in fresh DMEM/F12 containing 1 × SR2. The next day, the medium was removed and 20 μl of it was subjected to the Phospha-Light Assay (Applied Biosystem) for luminometric quantification of secreted SEAP protein in the culture supernatant according to the instructions of the manufacturer. Cells were lysed using the β-Glo Assay (Promega) and subjected to this assay for luminometric determination of transfected β-galactosidase enzyme as internal control for the transfection efficacy. Both luminometric SEAP and β-GAL quantification were carried out using a VICTOR3 multilabel plate reader (Perkin–Elmer). To normalize the data, SEAP values are expressed in relation to the measured β-GAL values.
Results
Expression of ERβ1, ERβ-δ125, and ERβ-δ1256 in human ovary and ovarian cancer
First, we examined whether ERβ1 and the exon-skipped ERβ splice variants (Fig. 1) cloned from human breast cancer cells (Treeck et al. 2007) were also expressed in ovary and in ovarian cancer. For this purpose, a cDNA pool from epithelial ovarian cancer tissue specimen and from normal ovary (Rapid-Scan Kit, Origene) was screened for ERβ1, ERβδ125, and ERβδ1256 transcripts by means of RT-PCR. To confirm specificity of amplification of the exon-skipped variants, a set of isoform-specific PCR primers was used annealing at the junction of the 5′-UTR and exon 3 (primer δ12) and the junction of exon 4 and 6 (primer δ5) or exon 4 and 7 (primer δ56) respectively and identity of the resulting amplicons was confirmed by sequencing. ERβ1, the specific 438 bp ERβ-δ125 amplicon and the 450 bp ERβ-δ1256 cDNA fragment were detected both in normal ovary and in ovarian cancer tissue (Fig. 2a). In all samples, detection of β-actin generally was used as positive control (not shown).
Heterologous expression of ERβ1, ERβ-δ125, and ERβ-δ1256 in SK-OV-3 cells and siRNA-triggered knockdown of ERβ1
To confirm estrogen-unresponsiveness of SK-OV-3 cells, we assessed ERE activation in this cell line and in OVCAR-3 cells by means of reporter gene assays. Stimulation by estradiol resulted in ERE activation in ERα-positive OVCAR-3 cells, but not in SK-OV-3 ovarian cancer cells (Fig. 2b).
Given that real time PCR analysis of ERβ1, ERβ-δ125 and ERβ-δ1256 mRNA levels revealed a very weak expression of these receptor isoforms in SK-OV-3 cells (Fig. 2c), which was 20- to 50-fold lower than in the pooled ovarian cancer samples (data not shown), we used this cell line to elucidate the function of ERβ1 and the exon-skipped ERβ splice variants in ovarian cancer cells by means of heterologous gene expression. SK-OV-3 cells were transfected with pTARGET mammalian expression vectors (Promega) containing the coding region of ERβ1, ERβ-δ125, or ERβ-δ1256 or the original pTARGET vector as negative control. After verification of their expression in transient transfection assays on mRNA level by means of RT-PCR (data not shown), SK-OV-3 clones stably expressing the transfected pTARGET derivatives were generated by G418 selection (300 μg/ml). About 6 weeks after transfection, 3–6 SK-OV-3 clones per derivative were isolated using cloning disks and propagated. In these clones, mRNA levels of ERβ1, ERβ-δ125, or ERβ-δ1256 respectively, was quantified in relation to β-actin expression by means of real time RT-PCR, avoiding false-positive signals from vector DNA by comparison to a sample, which was not reversely transcribed. Heterologous expression of ERβ isoforms in SK-OV-3 cells was additionally verified by sequencing of the amplified cDNA. SK-OV-3 clones mock transfected with the original pTARGET vector as negative control were identified by detection of mRNA transcribed from the neomycin resistance gene of this vector by means of RT-PCR (primers pTAR1 and pTAR2). Two clones, H (higher expression) and L (lower expression) from SK-OV-3/β1, SK-OV-3/δ125, SK-OV-3/δ1256 cells exhibiting ERβ-isoform mRNA levels similar to the ovarian cancer samples, were chosen for further characterization. Additionally, we chose an RNAi approach to confirm specificity of the observed ERβ1 effects on proliferation and apoptosis. We used siRNA specific for ERβ1 to knockdown ERβ1 expression in SK-OV-3/ERβ1-H cells. Transfection of these cells with 30 nM ERβ1 siRNA resulted in a significant reduction of ERβ1 mRNA levels down to 10% of the respective level in control cells (transfected with scrambled siRNA; Fig. 2c). For the first time we succeeded in detection of the δ125 and δ1256 isoforms on protein level in transfected SK-OV-3 clones by means of western blot analysis. The use of ERβ antibody GR39/Ab-2 (Oncogene) succeeded in detection of weak bands of expected size in transfected SK-OV-3 cells, but not in WT or vector-transfected cells, confirming heterologous expression of the exon-skipped isoforms on protein level (Fig. 2d).
Proliferation of SK-OV-3 cells heterologously expressing ERβ1, ERβ-δ125 or ERβ-δ1256
Given that ERs are known to regulate cellular proliferation by different molecular mechanisms, we examined the effect of heterologous ERβ isoform expression on cellular proliferation of SK-OV-3 cells. For this purpose, both vector-transfected and ERβ-transfected SK-OV-3 cells were cultured in serum-free SR2 medium and treated with E2 (1 nM) alone or in combination with 4-OH TAM (1 μM) for up to 6 days. In serum-free culture medium, SK-OV-3 cells containing higher ERβ1 transcript levels exhibited a significantly reduced proliferation if compared with vector-transfected cells even in the absence of E2 (P < 0.01 versus vector control). This effect was not observed in cells expressing lower ERβ1 mRNA levels or the exon-skipped variants. Addition of 1 nM E2, a concentration, which was observed to exert the strongest effect on cell growth of ERβ1-H cells (Fig. 3b), further slowed cell growth of SK-OV-3 cells stably expressing higher levels of ERβ1 mRNA, but did not affect proliferation of SK-OV-3 cells expressing δ125 or δ1256 transcript isoforms or lower levels of ERβ1. Addition of 4-OH TAM to the E2-containing culture medium was not able to affect E2-triggered growth inhibition in SK-OV-3 cells overexpressing ERβ1 (Fig. 3a). In contrast, addition of pure antiestrogen ICI 182 780 (100 nM) significantly inhibited the effect of E2 on growth of these cells (Fig. 3c). To confirm specificity of the observed antiproliferative effect of ERβ1 expression, we used an RNAi approach to knockdown ERβ1 expression in SK-OV-3/ERβ1-H cells. Growth analysis of siRNA-treated cells exhibiting significantly reduced ERβ1 levels (Fig. 2c) revealed that ERβ1 knockdown reverted the growth inhibitory effect of ERβ1 overexpression (Fig. 3d).
Apoptosis of SK-OV-3 cells heterologously expressing ERβ1, ERβ-δ125, or ERβ-δ1256
The decreased cell growth observed in SK-OV-3 cells transfected with ERβ isoforms could result not only from cell-cycle blockage, but also from increased apoptosis. To examine the effect of ERβ isoforms on apoptosis of SK-OV-3 cells, basal caspase 3/7 activity was analyzed. We observed significantly increased caspase 3/7 activation in SK-OV-3 cells expressing higher levels of ERβ1 mRNA grown in serum-free medium even in the absence of E2, but not in cells expressing lower levels of ERβ1 or the exon-skipped isoforms. Addition of E2 (Fig. 4a) or 4-OH TAM (not shown) did not affect the increased apoptosis of SK-OV-3 cells overexpressing ERβ1. In two of four experiments, we additionally performed experiments comparing the apoptotic cell membrane phosphatidylserine translocation in the different SK-OV-3 clones by means of double staining with Annexin V and PI confirming the results from the caspase assays (data not shown).
To confirm specificity of the observed apoptotic effect of ERβ1 on SK-OV-3 ovarian cancer cells, again we used an RNAi approach. SK-OV-3/ERβ1-H cells treated with ERβ1 siRNA 24 h prior to apoptosis detection and tested for ERβ1 knockdown to about 15% (Fig. 2c) exhibited a significantly reduced apoptosis when compared with mock-transfected cells (Fig. 4b).
Motility of SK-OV-3 cells heterologously expressing ERβ1, ERβ-δ125, or ERβ-δ1256
Because estrogen signaling is known not only to affect cell growth and apoptosis, but also cellular migration, it was important to determine whether ERβ1 or the exon-skipped ERβ isoforms could also affect motility of SK-OV-3 ovarian cells. For this purpose, we performed wound healing assays (Fig. 5a). In serum-free medium, SK-OV-3 WT and vector transfected control cells had filled about 75% of the wound after 4 days. In contrast, SK-OV-3 cells expressing higher levels of ERβ1 exhibited significantly slower migration ability as they filled only about 5% of the gap. SK-OV-3 cells expressing ERβδ125, high levels of ERβδ1256 or low levels of ERβ1 also exhibited a significantly diminished motility as they filled about 40–50% of the wound (Fig. 5b). Migration of SK-OV-3 cells expressing lower levels of ERβδ1256 did not differ from WT or control cells. The motility of all SK-OV-3 clones was not affected by treatment with 1 nM E2.
Expression of estrogen-responsive genes in SK-OV-3 cells heterologously expressing ERβ1, ERβ-δ125, or ERβ-δ1256
Given that ERs are ligand-inducible transcription factors directly regulating gene transcription, we studied the effect of ERβ1 and the exon-skipped ERβ splice variants on expression of 15 estrogen-responsive genes in SK-OV-3 cells (progesterone receptor (PR), cyclin D1, CDK2, autotaxin, PS2, ERα, FAS ligand, HER2, cathepsin D, EGFR, IGFBP-4, WISP-2, p21(WAF1), cyclin A2, and fibulin 1c). For this purpose, we analyzed expression of these genes in WT, control, and ERβ-transfected SK-OV-3 cells cultured in serum-free medium (± 1 nM E2 for 24 h) on mRNA level by means of real time RT-PCR.
Three of the analyzed genes, p21(WAF1), cyclin A2, and fibulin-1c exhibited altered mRNA levels in ERβ1-transfected SK-OV-3 ovarian cancer cells (Fig. 6). In SK-OV-3 cells overexpressing ERβ1, the up to fourfold elevated basal p21(WAF1) transcript levels were further increased after addition of E2 or 4-OH TAM. In SK-OV-3 cells expressing the δ1256 isoform or higher levels of the δ125 isoform, p21(WAF1) levels were also slightly increased, but were not affected by treatment with E2 or tamoxifen. In contrast, cyclin A2 mRNA level was decreased in SK-OV-3/ERβ1-H cells in the absence of E2, and mRNA levels of this gene were further reduced after addition of E2 or 4-OH TAM. Fibulin-1c transcript levels were fourfold elevated in SK-OV-3 cells expressing higher levels of ERβ1, and were further increased after treatment with E2 or tamoxifen.
Irrespective of estrogen or tamoxifen treatment, we did not observe any significant differences between the different clones and WTor vector-transfected SK-OV-3 cells regarding the expression of PR, cyclin D1, CDK2, autotaxin, PS2, ERα, FAS ligand, HER2, cathepsin D, EGFR, IGFBP-4, or WISP-2.
Discussion
The aim of this study was to determine the extent to which ERβ1 and two of its exon-skipped isoforms modulate basic features of SK-OV-3 ovarian cancer cells like proliferation, motility and apoptosis, and to study changes in gene expression as potential underlying molecular mechanism. Recently, we have identified the two novel exon-skipped ERβ transcript variants ERβ-δ125 and ERβ-δ1256 in human breast cancer cells (Treeck et al. 2007). In this study, we detected ERβ1, ERβ-δ125, and ERβ-δ1256 transcripts both in human ovary and ovarian cancer tissue, but we did not measure a notable expression of these ERβ types in SK-OV-3 ovarian cancer cells. In contrast to other ovarian adenocarcinoma lines like OVCAR-3 or BG-1 expressing functional ERα and other steroid hormone receptors like PR and AR, or normal ovarian epithelial cells, this cell line derived from an epithelial ovarian tumor is estrogen-unresponsive and HER2-overexpressing, and represents a relatively aggressive and fast growing ovarian cancer cell type (Lau et al. 1999). Our data demonstrating ERE activation by E2 in OVCAR-3, but not in SK-OV-3 cells confirm estrogen unresponsiveness of this cell line on molecular level.
In this study, we stably introduced cDNA coding for ERβ1, ERβ-δ125, and ERβ-δ1256 into SK-OV-3 cells to study the function of ERβ in this ovarian cancer model. For further characterization, we have chosen not the stably transfected clones exhibiting the highest expression levels, but the ones with lower overexpression levels comparable with the respective expression we measured in OVCAR-3 cells.
Several in vitro studies show evidence that ERβ may negatively regulate cellular proliferation, promote apoptosis and thus may have a protective role in normal breast and prostate. The same studies also suggested that antitumoral effects of ERβ are not necessarily dependent on the presence of ERα (Brandenberger et al. 1998, Lazennec et al. 2001, Cheng et al. 2004). Supporting these reports, in this study, we demonstrate that ERβ exerts antitumoral effects on SK-OV-3 ovarian cancer cells not expressing functional ERα (Lau et al. 1999). A loss of ERβ expression or increased ERα/ERβ ratio in epithelial ovarian cancer as compared with normal tissues has been reported consistently by several groups (Brandenberger et al. 1998, Pujol et al. 1998, Rutherford et al. 2000). A loss of ERβ expression could thus constitute a crucial step in ovarian carcinogenesis and hormone unresponsiveness. However, the role of ERβ in ERα-positive or -negative ovarian cancer cells is not completely understood. In this regard, the specific role of ERβ splice isoforms also remains unclear, though many ERβ splice variants are expressed in ovary and ovarian cancer (Poola et al. 2002a,b). Three different ERβ variant mRNAs that have deletions in exon 5 or 6 or exons 5/6 have been identified in human breast, uterus, and ovarian tissues (Lu et al. 1998, Vladusic et al. 1998, Speirs et al. 2000). A recent study examined the function of one of these exon-skipped variants, ERβ-δ5, suggesting that this isoform might act as a dominant negative receptor on ERα and ERβ pathways (Helguero et al. 2005). In another study, an ERβ isoform lacking the exons 2, 5, and 6 was identified and it was stated that deletion of these exons would cause a frame shift mutation resulting in premature termination of translation (Poola et al. 2002a,b). The exon-skipped variants ERβ-δ125 and ERβ-δ1256 we examined here use a different translation initiation codon in the beginning of exon 3 allowing translation in the same reading frame as ERβ1. The proteins coded by these variants are predicted not to contain the AF-1 domain mediating the ligand-independent transcriptional activity of ERβ and are predicted to have deletions both in the DBD and LBD. Thus, it is expected that both the ligand-dependent and ligand-independent activity of the deduced proteins are significantly diminished.
SK-OV-3 ovarian cancer cells are reported to be estrogen-unresponsive because they express a truncated, dysfunctional ERα, and were described to express very low levels of ERβ (Jones et al. 1994, Lau et al. 1999). Our findings demonstrating a reduced proliferation of SK-OV-3 cells stably expressing higher levels of ERβ1 even in absence of E2 are in agreement with previous studies showing ligand-independent antiproliferative effects of this receptor on tumor cells of different origin (Lazennec et al. 2001, Cheng et al. 2004). Growth of SK-OV-3 cells overexpressing ERβ1 was further reduced by E2 also demonstrating a ligand-dependent action of this receptor in our cellular system. The results of our RNAi approach demonstrating a reversion of ERβ1-triggered growth inhibition clearly confirm the antiproliferative action of this receptor in SK-OV-3 ovarian cancer cells. The absence of any significant effect of E2 or tamoxifen on SK-OV-3 cells stably expressing the exon-skipped variants could be explained by their LBD deletions, which are expected to impair ligand binding.
It was also of great interest to analyze the effects of ERβ1 and the two exon-skipped isoforms on apoptosis, because cell growth results from the balance of both cell cycle events and apoptosis regulation. We decided to examine the intrinsic apoptotic pathway, because ERβ previously was reported to promote apoptosis in a caspase 3-dependent manner in breast and prostate cancer cells (Cheng et al. 2004, Mak et al. 2006). And in fact, introduction of ERβ1 in ovarian cancer cells led to an increased basal apoptosis rate as shown both by caspase 3 and 7 activation and Annexin V staining. Apoptosis rate was not significantly elevated in SK-OV-3 cells expressing the exon-skipped isoforms, suggesting that the AF-1 or LBD domain are important for the apoptotic effect of ERβ. Interestingly, in SK-OV-3 cells overexpressing ERβ1 we did not observe any effect of E2 or tamoxifen on apoptosis suggesting that ERβ induction of apoptosis in this ovarian cancer cell line is a ligand-independent effect. The data from our siRNA approach demonstrating apoptosis reduction after treatment of SK-OV-3/ERβ1-H cells with ERβ1 siRNA clearly supports specificity of apoptotic ERβ1 action in this ovarian cancer model.
We also investigated the potential modulation of motility by ERβ1 and both exon-skipped isoforms, as another key event occurring during tumor development. We indeed observed that ERβ1 drastically inhibits motility of the ovarian cancer cell line as observed previously in ERα-negative breast and prostate cancer models exogenously expressing ERβ (Lazennec et al. 2001, Cheng et al. 2004). Again, these effects were estrogen-independent and motility of ERβ-δ125- and ERβ-δ1256-transfected SK-OV-3 cells was altered only marginally suggesting that the exons deleted in both variants are important to confer ERβ inhibition of motility.
Little is known about the gene regulatory function of ERβ in ovarian tissue. To analyze the molecular mechanisms underlying the observed alterations in proliferation, apoptosis, and motility of SK-OV-3 ovarian cancer cells expressing higher levels of ERβ1, we examined expression of a set of 15 estrogen-responsive genes on mRNA level. Even in the absence of E2, transcript levels of cell cycle inhibitor p21(WAF1) were strongly elevated in ERβ1-expressing cells and to a smaller extent also in SK-OV-3 cells expressing the exon-skipped isoforms, suggesting that ERβ is able to increase p21(WAF1) levels in a ligand-independent manner. Only in SK-OV-3 cells expressing full-length ERβ1, p21(WAF1) transcript levels were further elevated after treatment with E2 alone or in combination with 4-OH TAM, supporting previous studies suggesting an involvement of p21(WAF1) in cellular estrogen response (Planas-Silva & Weinberg 1997, Thomas et al. 1998). Given that our results demonstrate that p21(WAF1) mRNA level is increased in Sk-OV-3 cells expressing higher levels of ERβ1 and, which exhibit a reduced cell growth, it is tempting to speculate that p21(WAF1) might be a key mediator of the antiproliferative effect of ERβ1 in this ovarian cancer model.
Cyclin A2 is a cell cycle regulator, which is known to be estrogen responsive (Vendrell et al. 2004). The observed reduction of cyclin A2 mRNA levels in SK-OV-3 cells expressing ERβ1 could be another molecular mechanism underlying the growth inhibitory action of this receptor on SK-OV-3 cells. Downregulation of cyclin A2 in SK-OV-3 cells overexpressing ERβ1 after treatment with E2 clearly corresponds with the observed growth inhibitory effects.
The third gene exhibiting altered transcript levels in ERβ1-overexpressing SK-OV-3 cells was fibulin-1c, an extracellular matrix protein, which is overexpressed in epithelial ovarian and breast cancers and is involved in the regulation of cellular motility (Hayashido et al. 1998). Previous studies demonstrated that in ERα-positive ovarian and breast cancer cell lines, fibulin-1c mRNA levels are markedly increased by estrogens (Moll et al. 2002, Bardin et al. 2005). Our data demonstrate that heterologous expression of ERβ1 in ovarian cancer cells also is able to strongly increase fibulin-1c transcript levels even in absence of functional ERα or E2, suggesting that the regulation of fibulin-1c by ER pathways is more complex than assumed. Indeed, recent studies reported that fibulin-1c is a gene, which is estrogen-responsive not through classical ERα binding to ERE, but by E2-triggered activation of specificity protein 1 binding sites (Moll et al. 2002). Thus, the observed upregulation of fibulin 1c expression in SK-OV-3/ERβ1-H cells could be at least one molecular mechanism underlying the decreased motility of these cells.
In this study, we analyzed the functions of ERβ1 and two exon-skipped ERβ splice isoforms by means of heterologous gene expression in estrogen-unresponsive SK-OV-3 ovarian cancer cells. Particularly overexpression of ERβ1 exerted strong antitumoral effects on SK-OV-3 cells in terms of inhibition of growth and motility and induction of apoptosis, accompanied by specific changes in gene expression. Our results clearly suggest that the tumor suppressor function of ERβ1 in ovarian cancer cells is not necessarily dependent on ERα expression.
Primer sequences used for RT-PCR amplification
Oligonucleotides | Sequences 5′ –3′ | |
---|---|---|
Target | ||
ERβ-δ125 | δ12 | GGTGTGTTTATCTGCAAGGACA |
δ5 | CTCATCCCGGGAATCTTCTT | |
ERβ-δ1256 | δ12 | GGTGTGTTTATCTGCAAGGACA |
δ56 | CCAGAGGGTACATACCGGGAA | |
β-Actin | Actin-2573 | CTGTGGCATCCACGAAACTA |
Actin-2876 | CGCTCAGGAGGAGCAATG | |
neor | pTAR1 | ATGACTGGGCACAACAGACA |
pTAR2 | CTCGTCCTGCAGTTCATTCA | |
Fibulin-1c | Fib1 | CGAGTGCCCTGAGAACTACC |
Fib2 | GATGTTGGTGGGGAAAGAGA | |
Cyclin A2 | CYCA-1 | CTGCTGCTATGCTGTTAGCC |
CYCA-2 | TGTTGGAGCAGCTAAGTCAAAA | |
p21(WAF1) | WAF1 | GCATGACAGATTTCTACCACTCC |
WAF2 | AAGATGTAGAGCGGGCCTTT |
(O Treeck and G Pfeiler contributed equally to this study)
We thank Angelika Vollmer, Bettina Ederhofer and Helena Houlihan for their expert technical assistance. The authors declare that there is no conflict of interest that would prejudice the impartiality of this scientific work.
References
Anderson GL, Judd HL, Kaunitz AM, Barad DH, Beresford SA, Pettinger M, Liu J, McNeeley SG & Lopez AM 2003 Effects of estrogen plus progestin on gynecologic cancers and associated diagnostic procedures: the women’s health initiative randomized trial. JAMA 290 1739–1748.
Bardin A, Moll F, Margueron R, Delfour C, Chu ML, Maudelonde T, Cavailles V & Pujol P 2005 Transcriptional and posttranscriptional regulation of fibulin-1 by estrogens leads to differential induction of messenger ribonucleic acid variants in ovarian and breast cancer cells. Endocrinology 146 760–768.
Brandenberger AW, Tee MK & Jaffe RB 1998 Estrogen receptor alpha (ER-alpha) and beta (ER-beta) mRNAs in normal ovary, ovarian serous cystadenocarcinoma and ovarian cancer cell lines: down-regulation of ER-beta in neoplastic tissues. Journal of Clinical Endocrinology and Metabolism 83 1025–1028.
Burns KH, Agno JE, Chen L, Haupt B, Ogbonna SC, Korach KS & Matzuk MM 2003 Sexually dimorphic roles of steroid hormone receptor signaling in gonadal tumorigenesis. Molecular Endocrinology 17 2039–2052.
Cheng J, Lee EJ, Madison LD & Lazennec G 2004 Expression of estrogen receptor beta in prostate carcinoma cells inhibits invasion and proliferation and triggers apoptosis. FEBS Letters 566 169–172.
Couse JF, Curtis Hewitt S & Korach KS 2000 Receptor null mice reveal contrasting roles for estrogen receptor alpha and beta in reproductive tissues. Journal of Steroid Biochemistry and Molecular Biology 74 287–296.
Greenlee RT, Murray T, Bolden S & Wingo PA 2000 Cancer statistics. CA: A Cancer Journal for Clinicians 50 7–33.
Gustafsson JA 1999 Estrogen receptor beta–a new dimension in estrogen mechanism of action. Journal of Endocrinology 163 379–383.
Hatch KD, Beecham JB, Blessing JA & Creasman WT 1991 Responsiveness of patients with advanced ovarian carcinoma to tamoxifen. A Gynecologic Oncology Group Study of second-line therapy in 105 patients. Cancer 68 269–271.
Havrilesky LJ, McMahon CP, Lobenhofer EK, Whitaker R, Marks JR & Berchuck A 2001 Relationship between expression of coactivators and corepressors of hormone receptors and resistance of ovarian cancers to growth regulation by steroid hormones. Journal of the Society for Gynecologic Investigation 8 104–113.
Hayashido Y, Lucas A, Rougeot C, Godyna S, Argraves WS & Rochefort H 1998 Estradiol and fibulin-1 inhibit motility of human ovarian- and breast-cancer cells induced by fibronectin. International Journal of Cancer 75 654–658.
Helguero LA, Faulds MH, Gustafsson JA & Haldosen LA 2005 Estrogen receptors alfa (ERalpha) and beta (ERbeta) differentially regulate proliferation and apoptosis of the normal murine mammary epithelial cell line HC11. Oncogene 24 6605–6616.
Herynk MH & Fuqua SA 2004 Estrogen receptor mutations in human disease. Endocrine Reviews 25 869–898.
Jones J, Lagasse LD & Karlan BY 1994 Steroid hormonal independence of HER-2/neu mRNA expression in four human ovarian carcinoma cell lines. Gynecologic Oncology 53 421–426.
Lacey JV Jr, Mink PJ, Lubin JH, Sherman ME, Troisi R, Hartge P, Schatzkin A & Schairer C 2002 Menopausal hormone replacement therapy and risk of ovarian cancer. JAMA 288 334–341.
Lau KM, Mok SC & Ho SM 1999 Expression of human estrogen receptor-alpha and -beta, progesterone receptor, and androgen receptor mRNA in normal and malignant ovarian epithelial cells. PNAS 96 5722–5727.
Lazennec G, Bresson D, Lucas A, Chauveau C & Vignon F 2001 ER beta inhibits proliferation and invasion of breast cancer cells. Endocrinology 142 4120–4130.
Lindberg MK, Moverare S, Skrtic S, Gao H, Dahlman-Wright K, Gustafsson JA & Ohlsson C 2003 Estrogen receptor (ER)-beta reduces ERalpha-regulated gene transcription, supporting a ‘ying yang’ relationship between ERalpha and ERbeta in mice. Molecular Endocrinology 17 203–208.
Livak KJ & Schmittgen TD 2001 Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔ CT method. Methods 25 402–408.
Lu B, Leygue E, Dotzlaw H, Murphy LJ, Murphy LC & Watson PH 1998 Estrogen receptor-beta mRNA variants in human and murine tissues. Molecular and Cellular Endocrinology 138 199–203.
Mak P, Leung YK, Tang WY, Harwood C & Ho SM 2006 Apigenin suppresses cancer cell growth through ERbeta. Neoplasia 8 896–904.
Merchenthaler I & Shugrue PJ 1999 Estrogen receptor-beta: a novel mediator of estrogen action in brain and reproductive tissues. Morphological considerations. Journal of Endocrinological Investigation 22 10–12.
Moll F, Katsaros D, Lazennec G, Hellio N, Roger P, Giacalone PL, Chalbos D, Maudelonde T, Rochefort H & Pujol P 2002 Estrogen induction and overexpression of fibulin-1C mRNA in ovarian cancer cells. Oncogene 21 1097–1107.
Osborne CK & Schiff R 2005 Estrogen-receptor biology: continuing progress and therapeutic implications. Journal of Clinical Oncology 23 1616–1622.
Pettersson K & Gustafsson JA 2001 Role of estrogen receptor beta in estrogen action. Annual Review of Physiology 63 165–192.
Planas-Silva MD & Weinberg RA 1997 Estrogen-dependent cyclin E-cdk2 activation through p21 redistribution. Molecular and Cellular Biology 17 4059–4069.
Poola I, Abraham J & Liu A 2002a Estrogen receptor beta splice variant mRNAs are differentially altered during breast carcinogenesis. Journal of Steroid Biochemistry and Molecular Biology 82 169–179.
Poola I, Abraham J & Baldwin K 2002b Identification of ten exon deleted ERbeta mRNAs in human ovary, breast, uterus and bone tissues: alternate splicing pattern of estrogen receptor beta mRNA is distinct from that of estrogen receptor alpha. FEBS Letters 516 133–138.
Price RH Jr, Lorenzon N & Handa RJ 2000 Differential expression of estrogen receptor beta splice variants in rat brain: identification and characterization of a novel variant missing exon 4. Brain Research. Molecular Brain Research 80 260–268.
Price RH Jr, Butler CA, Webb P, Uht R, Kushner P & Handa RJ 2001 A splice variant of estrogen receptor beta missing exon 3 displays altered subnuclear localization and capacity for transcriptional activation. Endocrinology 142 2039–2049.
Pujol P, Rey JM, Nirde P, Roger P, Gastaldi M, Laffargue F, Rochefort H & Maudelonde T 1998 Differential expression of estrogen receptor-alpha and -beta messenger RNAs as a potential marker of ovarian carcinogenesis. Cancer Research 58 5367–5373.
Rodriguez C, Patel AV, Calle EE, Jacob EJ & Thun MJ 2001 Estrogen replacement therapy and ovarian cancer mortality in a large prospective study of US women. JAMA 285 1460–1465.
Rutherford T, Brown WD, Sapi E, Aschkenazi S, Munoz A & Mor G 2000 Absence of estrogen receptor-beta expression in metastatic ovarian cancer. Obstetrics and Gynecology 96 417–421.
Scambia G, Benedetti-Panici P, Ferrandina G, Distefano M, Salerno G, Romanini ME, Fagotti A & Mancuso S 1995 Epidermal growth factor, oestrogen and progesterone receptor expression in primary ovarian cancer: correlation with clinical outcome and response to chemotherapy. British Journal of Cancer 72 361–366.
Sierens JE, Scobie GA, Wilson J & Saunders PT 2004 Cloning of oestrogen receptor beta from old and new world primates: identification of splice variants and functional analysis. Journal of Molecular Endocrinology 32 703–718.
Speirs V, Adams IP, Walton DS & Atkin SL 2000 Identification of wild-type and exon 5 deletion variants of estrogen receptor beta in normal human mammary gland. Journal of Clinical Endocrinology and Metabolism 85 1601–1605.
Ståhlberg A, Åman P, Ridell B, Mostad P & Kubista M 2003 A quantitative real-time PCR method for detection of B-Lymphocyte monoclonality by comparison of kappa and lambda immunoglobulin light chain expression. Clinical Chemistry 49 51–59.
Thomas TJ, Faaland CA, Adhikarakunnathu S, Watkins LF & Thomas T 1998 Induction of p21 (CIP1/WAF1/SID1) by estradiol in a breast epithelial cell line transfected with the recombinant estrogen receptor gene: a possible mechanism for a negative regulatory role of estradiol. Breast Cancer Research and Treatment 47 181–193.
Treeck O, Pfeiler G, Federhofer B, Houlihan H, Vollmer A, Horn F & Ortmann O 2007 Novel estrogen receptor beta transcript variants identified in human breast cancer cells affect cell growth and apoptosis of COS-1 cells. Molecular and Cellular Endocrinology 264 50–60.
Vendrell JA, Magnino F, Danis E, Duchesne MJ, Pinloche S, Pons M, Birnbaum D, Nguyen C, Theillet C & Cohen PA 2004 Estrogen regulation in human breast cancer cells of new downstream gene targets involved in estrogen metabolism, cell proliferation and cell transformation. Journal of Molecular Endocrinology 32 397–414.
Vladusic EA, Hombsby AE, Guerra-Vladusic FK & Lupa R 1998 Expression of estrogen receptor beta messenger RNA variant in breast cancer. Cancer Research 58 210–214.
Zhao C, Toresson G, Xu L, Koehler KF, Gustafsson JA & Dahlman-Wright K 2005 Mouse estrogen receptor beta isoforms exhibit differences in ligand selectivity and coactivator recruitment. Biochemistry 44 7936–7944.