| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
or ERß in Human Osteosarcoma Cells: Distinct and Common Target Genes for These Receptors
Departments of Molecular and Integrative Physiology (F.S., J.F., B.S.K.) and Cell and Structural Biology (D.H.B., B.S.K.), University of Illinois at Urbana-Champaign, Urbana, Illinois 61801; and Womens Health Research Institute (B.K., C.R.L.), Wyeth Research, Collegeville, Pennsylvania 19426
Address all correspondence and requests for reprints to: Dr. Benita S. Katzenellenbogen, University of Illinois, Department of Molecular and Integrative Physiology, 524 Burrill Hall, 407 South Goodwin Avenue, Urbana, Illinois 61801-3704. E-mail: katzenel{at}life.uiuc.edu.
| Abstract |
|---|
|
|
|---|
and ß (ER
and ERß). To compare the actions of these receptors, we generated U2OS human osteosarcoma cells stably expressing ER
or ERß, at levels comparable with those in osteoblasts, and we characterized their response to 17ß-estradiol (E2) over time using Affymetrix GeneChip microarrays to determine the expression of approximately 12,000 genes, followed by quantitative PCR verification of the regulation of selected genes. Of the approximately 100 regulated genes we identified, some were stimulated by E2 equally through ER
and ERß, whereas others were selectively stimulated via ER
or ERß. The E2-regulated genes showed three distinct temporal patterns of expression over the 48-h time course studied. Of the functional categories of the E2-regulated genes, most numerous were those encoding cytokines and factors associated with immune response, signal transduction, and cell migration and cytoskeleton regulation, indicating that E2 can exert effects on multiple pathways in these osteoblast-like cell lines. Of note, E2 up-regulated several genes associated with cell motility selectively via ERß, in keeping with the selective E2 enhancement of the motility of ERß-containing cells. On genes regulated equally by E2 via ER
or ERß, the phytoestrogen genistein preferentially stimulated gene expression via ERß. These studies indicate both common as well as distinct target genes for these two ERs, and identify many novel genes not previously known to be under estrogen regulation. | Introduction |
|---|
|
|
|---|
Estrogen exerts its effects on bone and other target tissues by interacting with two different members of the nuclear receptor superfamily of hormone-regulated transcription factors, named estrogen receptor (ER)
and ERß (5, 6, 7). These two receptors are encoded by different genes, on human chromosome 6 and 14, respectively. ER
and ERß have similar modular domain structures and very high amino acid identity in their DNA-binding domains (97%), whereas they are more divergent in their N-terminal A/B domains (only 18% amino acid identity) and in their ligand-binding domains (59% amino acid identity). After the binding of hormone to these receptors, the hormone-receptor complexes bind to specific sequences on the DNA [estrogen response elements (EREs)] or interact with other transcription factors without direct ER
or ERß binding to DNA (i.e. at activator protein-1, Sp1, and other sites) (8, 9, 10, 11, 12, 13). In both cases, liganded ERs recruit coregulator proteins and components of the transcriptional machinery to regulate the transcription of target genes (14, 15, 16). Recently, the importance of membrane-initiated signaling in the actions of estrogens has been highlighted by studies in vitro and in vivo, where membrane/cytoplasmic ERs as well as nuclear ERs seem to be involved in regulating E2 action and gene transcription in bone and other target cells (1, 17, 18, 19).
In bone, ERs are present in osteoblasts and chondrocytes, and at somewhat lower levels in osteoclasts, bone marrow stromal cells, osteocytes, and bone cell precursors (20, 21). The levels of ERs in bone are generally about 10-fold lower than in reproductive tissues, such as uterus, and their levels can be affected by many parameters including cell differentiation state (22, 23, 24). Estrogen seems to act directly on osteoblasts and probably, as well, on osteoclasts and precursors of both cell types, and indirectly by regulating cytokine production in osteoblasts and bone marrow stromal cells, which in turn affects the actions and formation of osteoclasts (1, 25).
The aims of this study were to determine the effects of E2 on gene expression in a human osteoblastic cell line and, via transcriptional profiling, to better understand the comparative roles of ER
and ERß in mediating the effects of this hormone. In particular, our interest was to elucidate whether the two receptors elicit the same and/or different transcriptional responses on a range of endogenous cellular genes. For this purpose, we have generated U2OS human osteosarcoma cells stably expressing either ER
or ERß and examined several clonal cell lines containing similar, relatively low, physiological levels of functional receptor. Using Affymetrix GeneChip microarrays, which allow the examination of >12,000 genes, we observed that the hormone-regulated genes divided into three distinct temporal patterns and into several functional categories, indicating that E2 exerts effects on multiple regulatory pathways in these cells. In this study, we report the identification of genes that were commonly regulated by E2 through ER
and ERß, as well as some that were preferentially or exclusively regulated by one or the other ER subtype, and many novel genes not previously known to be under regulation by this hormone, and we discuss the relationship of some of these genes to the biological effects of E2 we observed in these cells.
| Materials and Methods |
|---|
|
|
|---|
(encoding amino acids 1595) or ERß (encoding amino acids 1530). Each ER cDNA (26) was subcloned into the pcDNA3.1+ expression vector (Invitrogen, Carlsbad, CA), which contains a neomycin resistance gene. Clones were selected with the antibiotic G418 (800 µg/ml) and 20 clones per ER subtype were screened for ER expression and for transactivation ability with an ERE-containing reporter gene [2ERE-tk-CAT (chloramphenicol acetyltransferase)]. Four clones of each (denoted as ER
clones 14 and ERß clones 14) were chosen for further characterization and for the gene expression studies. The U2OS-ER cells were routinely grown in MEM with phenol red (Sigma, St. Louis, MO) supplemented with 15% fetal bovine serum (Hyclone, Logan, UT), 100 U penicillin/ml, 100 µg streptomycin/ml, and 400 µg G418/ml. Before use in experiments, cells were grown in MEM without phenol red and supplemented with 5% charcoal-dextran-treated fetal bovine serum for at least 4 d before the start of E2 treatment. MCF-7 human breast cancer cells were grown as previously described (27).
Hormone binding and Western immunoblot assays
Whole cell extracts were prepared in cell lysis buffer [20 mM Tris (pH 7.4), 0.5 M NaCl, 1 mM dithiothreitol, 10% glycerol, 50 µg/ml leupeptin, 50 µg/ml aprotinin, 2.5 µg/ml pepstatin A, and 0.2 mM phenylmethylsulfonyl fluoride] using freeze/thaw procedure. Total protein concentration was determined using the BCA kit (Pierce Biotechnology, Rockford, IL). Whole cell extracts were incubated in duplicate with a range of 3H-E2 concentrations alone or with 100-fold excess unlabeled E2 for 1 h on ice. Hydroxylapatite slurry was added and incubated for an additional 15 min on ice. The slurry was washed twice and its radioactivity then determined by scintillation counting. The amount of E2 binding in picomoles was calculated per milligram of protein.
Western blots of whole cell extracts used the ER
-specific antibody H222 and the ERß-specific antibody UCG40 (kindly provided by Geoffrey L. Greene, University of Chicago, Chicago, IL) and were done as previously described (28).
Hormone treatments, DNA microarrays, and analysis of microarray data
U2OS-ER
or ERß containing cells were maintained in culture for 48 h before cell harvest and RNA collection. During this period, cells were treated with 10 nM E2 (Sigma) for 4, 8, 24, or 48 h, and three separate samples were collected for each time point. Control cell samples were also treated with 0.1% ethanol control vehicle for 48 h. In this way, cell density was similar in all samples, consistent with our observations (as discussed in Results), that E2 does not affect proliferation of these cells. Additional independent time course experiments were conducted to generate RNA samples for quantitative RT-PCR verification of gene regulation, as detailed in the next section below. Total RNA was prepared using Trizol Reagent (Invitrogen) according to the manufacturers instructions. RNA was further purified using RNeasy columns (QIAGEN, Valencia, CA) and treatment with ribonuclease-free deoxyribonuclease I (QIAGEN).
Total RNA from each sample was used to generate cRNA, which was labeled with biotin as recommended by Affymetrix (Santa Clara, CA). Each cRNA was then hybridized on an Affymetrix human Hu-U95A GeneChip, which contains oligonucleotide probe sets representing approximately 12,500 human genes and expressed sequence tags. After washing, the chips were scanned and analyzed using MicroArray Suite 5.0 software (Affymetrix, Santa Clara, CA). Average intensities for each GeneChip were globally scaled to a target intensity of 150. Data were then analyzed using GeneSpring version 5.0.1 software (Silicon Genetics, San Carlos, CA). Data were first normalized on a per chip basis by dividing each measurement by the 50th percentile of all measurements on that chip, and then E2-treated samples were normalized to the mean of the vehicle-treated control samples.
We applied a confidence score (CS) to evaluate which genes were estrogen regulated (as described in Ref.29 ; adapted from Ref.30). The CS is based on four parameters: fold change (FC), P value (PV), percentage of present calls (PC) and expression level (EL), (CS = FC + PV + PC + EL). For each parameter, arbitrary scores were assigned. For FC, a score of 5 was assigned if the FC was greater than 1.95, 2 if the FC was between 1.5 and 1.95 and a penalty of 0.5 if it was under 1.5. For P value, a score of 3 was assigned if P value was less than 0.05, 2 if was between 0.05 and 0.1, and a penalty of 0.5 if it was greater than 0.1. If present calls were assigned to more than 50% of the samples, the score was 3, between 2550% was 1, and if less than 25% a penalty of 0.5 was applied. For the expression level, a score of 3 was applied if it was greater than 30, a score of 1 if it was between 15 and 30 and a penalty of 0.5 if expression level was less than 15. The confidence scores ranged from 2 to 14, and genes with a CS value of 11 or higher were considered to be significantly E2-regulated genes. To determine ER
or ERß preference in gene regulation, the average FC for one of the receptors had to be 2.0 or greater, whereas FC was less than 1.3-fold for the other receptor at all time points.
Estrogen-regulated genes were assigned to functional categories according to LocusLink, OMIM, PubMed (http://www.ncbi.nlm.nih.gov), GeneCards (http://bioinfo.weizmann.ac.il/cards) and GenMAPP databases. The entire microarray data set will be available through the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo).
Quantitative real-time PCR
Real-time PCR was carried out on approximately 25 genes to verify E2 regulation as assessed by microarray data analysis, and to evaluate mRNA levels of ER
or ERß in U2OS stably transfected cells. Three independent time-course experiments using three separate cell samples each for ER
or ERß containing cells either E2 or control vehicle-treated at each time point, were conducted. The primers used are listed in Table 1
. One microgram of total RNA from each sample was reverse transcribed in a total volume of 20 µl using 200 U reverse transcriptase, 50 pmol random hexamers, and 1 mM deoxynucleotide triphosphates (New England Biolabs, Beverly, MA). The resulting cDNA was then diluted to a total volume of 100 µl. Each real-time PCR consisted of 5 µl of diluted reverse transcription product, 1x SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA), and 50 nM of forward and reverse primers. Reactions were carried out in an ABI Prism 7700 Sequence Detection System (Applied Biosystems) for 40 cycles (95 C for 15 sec, 60 C for 1 min) after an initial 10-min incubation at 95 C. The FC in expression was calculated using the 
Ct comparative threshold cycle method (31) with the ribosomal protein 36B4 mRNA as an internal control. As described in Ref.31 , gene expression is normalized to an endogenous reference gene (36B4) and the FC in gene expression is then determined relative to the vehicle-treated control.
|
Cell adhesion assays were performed in 96-well plates coated for 2 h at 37 C with 2040 µg/ml of different coating substrates (poly-L-lysine, BSA, collagen type I, fibronectin, or laminin). Cells were treated with 10 nM E2 for 24 or 48 h, and after blocking nonspecific sites with 1% BSA for 2 h at 37 C, cells were harvested, resuspended in MEM minus phenol red containing 0.1% BSA at 105 cells/ml, and 100 µl of cell suspension were seeded per well. Cells were allowed to adhere for 30 min and then the medium was gently removed, the wells washed once with MEM minus phenol red and the cells then fixed with methanol. Cells were stained with 0.1% crystal violet for 30 min at room temperature and washed multiple times with double distilled water. Stain was then extracted from the cells with 5% Triton X-100 overnight at room temperature and absorbance monitored at 570 nm.
| Results |
|---|
|
|
|---|
or ERß
or ERß using a pcDNA3.1+ plasmid containing a neomycin resistance gene as a vector for stable integration. Positive clones were selected with G418 (800 µg/ml), and after two rounds of selection, various clones were tested for the presence of ER mRNA by real-time PCR and for ER protein by Western blot with ER
-specific or ERß-specific antibodies (Fig. 1
and ERß mRNA in U2OS-ER cell clones, as determined by real-time quantitative PCR using standard curves with known amounts of ER
and ERß cDNA (Fig. 1A
-containing human breast cancer cell line MCF-7. In the eight clones (ER
14, ERß 14) used for most of the studies reported in this paper, the number of copies of ER
and ERß mRNA were similar: 8410 ± 777 per 10 ng total RNA for ER
and 9112 ± 1304 for ERß. These ER
mRNA levels are approximately 20% of that measured in our MCF-7 breast cancer cells. In addition, as seen in Fig. 1A
-expressing U2OS cells and likewise no ER
mRNA was present in the U2OS ERß-expressing cells. The parental U2OS cells lacked detectable levels of either receptor mRNA.
|
clone 1 and ERß clone 1 cells to ensure the presence of receptor able to bind hormone. These assays demonstrated similar levels of 3H-E2 binding in these U2OS-ER
and U2OS-ERß cells (
0.1 pmol/mg protein, Fig. 1B
and UCG40 polyclonal antibody for ERß showed ER
or ERß of the correct size (
65 and 56 kDa, respectively), and the presence of only ERß in the ERß clones (ß clones 14 shown) and only ER
in the ER
clones (
clone 1 shown). Further characterization by estrogen-responsive ERE-tk-CAT reporter assays demonstrated that U2OS-ER
and ERß cells possessed E2-dependent transcriptionally functional receptors (data not shown). The level of ER
and ERß in these clones is consistent with previous reports indicating the approximately 10-fold lower ER level in bone vs. uterine tissue (23, 24).
We examined whether E2 affected the motility and adhesion properties of the U2OS-ER
and ERß-containing cells. Motility was assessed using a Boyden chamber chemotactic assay and 5% FBS as a general chemotactic stimulus. With E2 treatment for 24 or 48 h, the ERß-containing cells showed a 3-fold increase in motility; in contrast, although the ER
and ERß-containing cells showed similar motility, E2 elicited no change in motility of the ER
-containing cells (Fig. 1D
). We also evaluated whether E2 treatment for 24 or 48 h affected the adhesive properties of these cells and found that it did not. This employed an adhesion assay that allows a comparison of the rapid adhesion (in 30 min) of the cells to various substrates (poly-L-lysine, BSA, collagen type I, fibronectin, or laminin). Adhesion to collagen type I is shown as an example (Fig. 1E
) because it is a major component of bone matrix, and it is seen that E2 had no effect on the adhesion of the cells. This was the case with all of the substrates tested, indicating that this property is largely E2 independent.
Gene expression profiling and time courses of gene regulation using Affymetrix GeneChip microarrays
To characterize changes in gene expression in response to E2 treatment, U2OS-ER
(clone 1) and U2OS-ERß (clone 1) cells were treated with 10 nM E2 for times from 048 h, and RNA was harvested, purified, and derivatized for gene profiling. Affymetrix Hu-U95A GeneChips, which contain probes for approximately 12,000 human genes, were used. As described in Materials and Methods, a CS was applied to discriminate genes robustly and reproducibly regulated by E2. The CS takes into account several parameters for each gene and weights the FC increase or decrease vs. control, the reproducibility of the output (P value), the expression level related to the intensity of the signal on the chip and a fidelity score based on a decision matrix, which compares hybridization of the perfect match probes vs. the mismatch for each probe set.
We considered E2-regulated genes to be those with a confidence score equal to or greater than 11 (of a maximum score of 14), and we included in our further studies only genes up- or down-regulated more than 2-fold. Using these criteria, we identified 105 genes that were E2-regulated through ER
and/or ERß. This list of genes was subjected to gene tree cluster analysis using a standard correlation algorithm available in the GeneSpring software.
In the gene cluster analysis (Fig. 2
), up-regulated genes are shown in red, down-regulated genes in blue, and unchanging genes in yellow. We observed that many more genes were stimulated than inhibited (85 vs. 20) in ER
- or ERß-containing cells, and in further analyses, we focused on the up-regulated genes because these were more numerous, were the most robustly regulated, and gave good reproducibility of validation by real-time PCR.
|
-containing cells, the early and late pattern was the most frequently represented (49% of E2-stimulated genes through ER
), followed by the late (46%) and early (5%) patterns. In U2OS-ERß cells, the distribution was slightly different with 43%, 41%, and 16%, respectively, in these three time patterns.
|
and ERß (Fig. 4
, whereas nine genes were stimulated selectively by ERß. Thus, the two receptors showed substantial overlap in regulated genes (61% of the genes), but both also have some different target genes. The estrogen-stimulated genes regulated selectively through one of the two ER subtypes are listed in Table 3
|
|
|
|
Functional classification of E2-regulated genes
To classify E2 up-regulated genes into functional groups (Fig. 4
), we chose eight categories: cytokines/immune response, signal transduction, cell motility/cytoskeleton regulation, growth factors/hormones, apoptosis/cell proliferation, housekeeping, nucleic acid processing, and other/unknown. We then assigned each gene to one major category representing its primary function (even though many genes encode proteins that have several physiological functions), based on mining of several databases (LocusLink, PubMed, GeneCards, OMIM, GenMapp).
Many of the genes commonly regulated through ER
and ERß encode proteins associated with cytokine/immune response, signal transduction, and cell motility and cytoskeleton regulation (56% of all genes). Several ER
-selective genes were associated with nucleic acid processing. Very few genes regulated by the ERs were related to cell cycle, proliferation, or apoptosis, supporting our observations (not shown) that E2 had no effect on the proliferation rate of these cells.
Hormone-regulated genes encoding cytokines and factors associated with immune response
Three members of the Natural Killer cell lectin-type receptors (NKG family), NKG2-C, -E, and -F (KLRC2, KLRC3, and KLRC4) were robustly regulated by E2 through both receptors (Table 2A
). NKGs are activating receptors of natural killer cells that mediate HLA-E recognition (32, 33, 34). CD34 (common to both ER
and ERß up-regulation; Table 2A
) and sialomucin/CD164 (ER
selective; Table 3
) are cell surface antigens characteristic of human hematopoietic progenitors. CD164 is an adhesion receptor that inhibits proliferation of CD34+ hematopoietic progenitors (35, 36). IL-8, a cytokine known for its angiogenic and chemotactic properties, is thought to be involved in bone cell function under IL-1ß and TNF
control in bone marrow stromal cells and osteoblasts (37, 38). It was also regulated by E2 in both ER
- and ERß-containing U2OS cells (Table 2A
).
Cyclooxygenase-2 (Cox-2), a key enzyme regulating the production of prostaglandins, and also regulating angiogenesis, tumor cell invasion and inflammatory responses (39), was markedly up-regulated by E2 via both ER
and ERß (Table 2A
). Cox-2 has been shown recently to be under estrogenic regulation mediated through the phosphoinositide 3-kinase/Akt pathway (40).
Hormone-regulated genes associated with cell motility and cytoskeletal function
This category is of particular interest because, using cell motility assays, we observed an increase in motility of U2OS-ERß containing cells after E2 treatment (Fig. 1D
). Among the ERß selectively regulated genes, there are two candidates that might be involved in this response. The first is autotaxin (ENPP2), a known tumor autocrine motility factor that has been shown to stimulate the migration of melanoma cells via a G protein-coupled phosphoinositide 3-kinase
and Cdc42- and Rac1-dependent pathway. This regulation possibly involves Paxillin and focal adhesion kinase, two key players in cell adhesion and motility control (41, 42). Autotaxin has 5-nucleotide pyrophosphatase and phosphodiesterase activity, and it acts in the extracellular space as a lysophospholipase D with production of lysophosphatidic acid, which seems to be the tumor cell motility inducer (43). The experiments described in the section below demonstrate that autotaxin is very selectively regulated by ERß. Another candidate associated with the motility increase by ERß might be Rap1GEF (GFR), a guanine nucleotide exchange factor that constitutively activates Rap1 (44) and has been shown to be involved in cell motility by modulating integrin function (45, 46).
Quantitative real-time PCR validation and verification of gene expression regulation in multiple cell clones and using ER subtype-selective ligands
Table 1
shows the primers used, and Tables 2A
and 3
show the approximately 25 genes that we verified for their E2 regulation by quantitative PCR (names in italics). It is of note that many of these are genes newly identified as being regulated by E2. In addition, we verified their E2 regulation in several clones (four each) of U2OS ER
and/or ERß-containing cells (clones 14) and used ER subtype-selective ligands to confirm common regulation by both ERs or selective gene regulation by ER
or ERß. The tables indicate generally good correspondence in gene regulation between the microarray and the quantitative RT-PCR analyses using several ER
- or ERß-containing cell clones. RT-PCR, however, is often more sensitive, with a greater dynamic range, such that FCs found by RT-PCR were often greater than those obtained from microarray analysis.
Figure 5
shows real-time PCR validation of some genes that were commonly regulated by E2 through both ER subtypes. The time course of stimulation for four of these genes (GREB1, oligophrenin-1, IL-24, and carbonic anhydrase XII) is presented. All of these genes showed up-regulation by E2 by 4 h. There were robust, and very similar FCs in mRNA for GREB1 mediated via ER
or ERß over time (Fig. 5
). Oligophrenin-1 also showed a pattern of mRNA change that was similar for ER
and ERß in which a plateau in stimulation was achieved by 4 h for ER
and 8 h for ERß. Interestingly, IL-24 was preferentially up-regulated by ER
, with ERß being a much less effective mediating receptor. In contrast, the stimulation via ERß was substantially greater for the carbonic anhydrase XII gene than was the stimulation via ER
. The ER mediation of all of these E2-induced changes is clearly suggested by the fact that the stimulation is fully reversed by an excess of the antiestrogen ICI 182,780. In addition, and as expected, no stimulation was observed with this pure antiestrogen alone (Fig. 5
).
|
. The ICI antiestrogen reversal of the E2 stimulation of these genes supports the mediation by ER (Fig. 6
|
or ERß, cells were treated with the ER
-selective ligand propylpyrazoletriol [PPT (47, 48, 49, 50)], the ERß-selective ligand diarylpropionitrile [DPN (51, 52)], or the soy phytoestrogen genistein which has a higher affinity for ERß than ER
(53) (Fig. 7
or ERß, showed stimulation by PPT in the ER
-containing cells and by DPN in the ERß-containing cells, as expected. Interestingly, genistein, which has preferential binding affinity for ERß (53), was considerably more effective in stimulating PDZK1 gene expression via ERß (Fig. 7
. These findings confirm the common, ER
and ERß regulation for PDZK1 observed with E2, and the ERß selectivity for E2 in regulation of cystatin D expression.
|
| Discussion |
|---|
|
|
|---|
and ERß. Yet the individual contributions of ER
and ERß to gene regulation and physiological function in bone remain only partially understood. In this study, we have characterized a model system for examining the actions of the ER subtypes, ER
and ERß, in a human osteoblast-like cell type. These studies have enabled us to evaluate changes in the expression of approximately 12,000 genes by E2 through ER
and/or ERß. We have identified novel estrogen-regulated genes and genes regulated commonly as well as exclusively or preferentially through ER
or ERß. Furthermore, the identification of regulated genes allows some insight into the pathways and gene networks regulated by E2 through these two ERs.
Estrogen regulation of genes associated with cytokines, motility, cytoskeleton, and bone remodeling
Because estrogen plays important roles in bone remodeling and the control of osteoblast and osteoclast numbers (1), it is of note that E2 showed significant regulation of cytokines and of genes encoding proteins associated with regulation of the immune response. Our findings are consistent with the reports of estrogen regulation of genes encoding cytokines, growth factors, and bone matrix proteins in wild-type mouse trabecular bone and the importance of ER
and/or ERß in the mediation of these effects, as evidenced by the use of ER knockout mice (54, 55, 56). We observed E2 up-regulation of genes encoding three members of the Natural Killer cell lectin-type receptors (NKG family, KLRC2, KLRC3, and KLRC4). These NKG family genes were all regulated by E2 at early as well as later times. Other factors associated with cytokine/immune function observed to be under estrogen regulation included CD34 antigen, defensin ß1, IL-13 receptor
2, IL-8, phospholipase A2, selenoprotein P1, and sialomucin/CD164. Many of these genes, as indicated in Tables 2A
and 3
, represent new E2-regulated genes, ones not previously reported to be regulated by this hormone.
Another category of genes, in which many members were regulated by E2, were those associated with regulation of cell motility and the cytoskeleton. We observed a 3-fold increase in the motility of the ERß-containing cells in response to E2 treatment. Hence, it is of note that two genes preferentially regulated by E2 via ERß were autotaxin (ENPP2), a known tumor autocrine motility factor and phosphodiesterase, and Rap1 GEF, the guanine nucleotide exchange factor that constitutively activates Rap1 and is known to modulate integrin function (45, 46).
WISP-2, a connective tissue growth factor isolated from osteoblasts (57), was regulated by E2 via both receptor subtypes, and it might play a role in cytoskeletal regulation in these U2OS cells. Its stimulation by E2 has been reported in MCF-7 breast cancer cells (29, 58). In addition, oligophrenin-1, keratin 19, integrin
6, cyritestin 1, connexin 43, and cadherin 19, all of which play roles in cell motility and the cytoskeleton, were for the most part not previously known to be regulated by estrogen in bone, although several are known to be under estrogen regulation in reproductive target cells (59, 60).
BMP6, a potent osteogenic factor believed to play an important role in the bone-protective actions of estrogens, was observed to be up-regulated through both ERs, as was PTH-like hormone (PTHLH/PTHrP), a homolog of PTH that functions as an autocrine growth inhibitor for osteoblast-like cells and as an anabolic agent in osteoporosis (61, 62). Neurotensin was another growth factor/hormone found to be up-regulated by E2 via ER
and ERß. Cystatin D, an E2-responsive gene newly identified in these studies, was stimulated very preferentially via ERß. Because this gene encodes a protein that is a secreted inhibitor of cysteine peptidases (cathepsins S and H), it might function in several roles in bone physiology.
ER regulation of other cellular functions: similarities and differences in estrogen target genes and temporal patterns of gene regulation in osteosarcoma and breast cancer cells
E2 up-regulated the expression of several genes encoding proteins associated with nucleic acid processing, such as the ER coregulators RIP140 (receptor-interacting protein 140) and TRAP240 (thyroid hormone receptor-associated protein, 240-kDa subunit). However, in these U2OS cells, E2 regulated very few genes associated with apoptosis or cell proliferation, consistent with our findings that E2 had no effect on the proliferation rate of either the ER
- or ERß-containing U2OS cells (data not shown). In contrast, in a recently reported study in ER-positive human breast cancer cells, in which E2 stimulates cell proliferation, we observed that E2 stimulated the expression of positive proliferation regulators including multiple growth factors and genes involved in cell cycle progression, and E2 also down-regulated transcriptional repressors and antiproliferative and proapoptotic genes (29). Hence, the categories of genes stimulated by E2 in the U2OS cells and in breast cancer cells are very different, and they reflect the quite different physiological effects of E2 on these target cells. Furthermore, E2 down-regulated the expression of few genes in these ER-containing osteosarcoma cells, whereas E2 down-regulated the expression of more genes in ER-positive MCF-7 breast cancer cells than it up-regulated (29). Despite these marked differences in genes regulated in these two different types of estrogen target cells, some genes were stimulated by E2 in both U2OS and MCF-7 breast cancer cells. These included pS2, keratin 19, NHERF, GREB1, PDZK1, receptor interacting protein 140, WISP2, and connexin 43, all previously shown to be under estrogen regulation, as referenced in Tables 2A
and 3
. Some have important roles in cell cytoarchitecture and cytoskeletal regulation (e.g. NHERF, connexin 43), supporting the known effects of estrogen on the cytoskeleton of both mammary (63) and bone cells.
It is of interest that our time-course analysis in the U2OS cells (Fig. 3
) revealed that genes regulated by E2 divided into three temporal patterns, denoted E, E + L, or L in Tables 2A
and 3
: those regulated at early times only (4 and/or 8 h), those regulated at both early and late times, and those showing increased expression only at late times (24 and/or 48 h). In both MCF-7 and in ER-containing U2OS cells, approximately equal numbers of stimulated genes fell into the early and late, and late-only categories, whereas relatively few genes showed early-only regulation in the U2OS cells compared with the breast cancer cells, where approximately one third of the total stimulated genes showed the early-only expression pattern. We have used the protein synthesis inhibitor, cycloheximide, to begin to distinguish primary vs. secondary effects of E2-ER complexes on regulation of some of the interesting genes we have identified in the U2OS cells. These studies reveal that some of the genes turned on early that show elevated expression over time and hence fall into the "early and late" category (such as GREB1, autotaxin, oligophrenin-1, and KLRC4) likely represent primary estrogen response genes because their increase in mRNA expression is not prevented by cycloheximide. In contrast, cycloheximide blocked the estrogenic stimulation of cystatin D and PDZK1, a late and an early and late gene, respectively, suggesting that these are probably secondary response genes.
Novel estrogen-regulated genes and genes regulated selectively by ER
or ERß
These studies have enabled us to compare the actions of ER
and ERß. Although the majority of E2-stimulated genes were regulated through both ER subtypes, we identified a number of genes selectively regulated either by ER
or by ERß, as well as many novel estrogen-regulated genes (64). Most of the E2-regulated genes involved in cytoskeleton regulation and motility, signal transduction, cytokine, and immune response and growth factors/hormones were commonly regulated through both receptors. These included a marked up-regulation of Cox-2, cyritestin 1, keratin 19, several ILs, and integrin
6. The latter factor, integrin
6, forms laminin receptors when heterodimerized with integrin ß1 or ß4 subunits, and it is one of the factors involved in tumor invasion through basal membranes (65, 66, 67).
Important growth factors and hormones regulated by E2 through both receptors included several that can act in either an autocrine or paracrine manner, including BMP6 and PTHLH, to effect the bone-protective actions of E2 (68). PTHLH, a homolog of PTH, acts as an anabolic agent in osteoporosis (62), and its overproduction in breast and prostate cancers is associated with bone metastasis of these tumors (69, 70). Hence, its up-regulation by E2 via either ER
or ERß is of interest. Of the genes we observed to be selectively regulated by one ER subtype, nearly three fourths were ER
selective. These ER
-regulated genes included several associated with nucleic acid processing, two zinc finger proteins, the coregulator TRAP240, and the cytochrome P450 family member 2B6 (CYP2B6). Of the genes exclusively regulated by ERß were autotaxin and cystatin D, factors associated with motility/cytoskeleton and cytokine/immune function.
A recent report of related studies in U2OS cells identified some similar and some different gene targets for ER
and ERß (71) as we observed in this report. There are several important methodological differences between the two studies that could account for some of the differences in findings. First, the levels of ER
and ERß in the U2OS cells used by Monroe et al. (71) were very high, namely three times that of MCF-7 cells, whereas we chose to use cells that contained ER at levels more comparable to those in osteoblasts, about 20% that of MCF-7 cells. Also, the studies in Ref.71 examined E2 regulation only at 24 h, whereas we examined gene regulation over a time course of E2 treatment up to 48 h, so that we identified both early, early and late, and late-only responding genes; hence, some different sets of genes would likely be identified in the two studies. In addition, we used a confidence score that includes considerations of FC, P values, expression level and present calls in defining estrogen-regulated genes, whereas Monroe et al. (71) used FC only. Furthermore, the studies reported in Ref.71 used an earlier version Affymetrix microarray chip that contains probes for about half of the genes on the U95A array we used, and also uses fewer and different probe sets per gene. Hence, it is not surprising that there are some differences in the observed pattern of gene regulations between our study and that presented in Ref.71 .
The differences in gene targets of ER
vs. ERß, or of progesterone receptor-A vs. progesterone receptor-B (72), may reflect differences in the activation functions of these receptors, in particular their N-terminal activation regions, as well as the complexity of gene regulation by hormone occupied-nuclear receptors. The latter involves a combination of aspects (7, 16, 73): the diverse nature of target gene promoters (with different response elements that may be simple or composite) where receptors may bind directly to the DNA or act through tethering to other DNA-bound transcription factors; the nature of the coregulators recruited to the hormone-receptor-gene promoter; and the receptor protein conformations induced by different ligands in the context of associated coregulators and transcription factors.
These studies, which have allowed a comparison of the gene regulatory activities of ER
and ERß, should assist in the characterization of genes and cellular signaling pathways regulated by E2 through ER
and/or ERß. They highlight commonality but also significant differences in gene targets for these two ERs and have identified many novel genes not previously known to be under estrogen regulation. Several of these are likely to play significant roles in the bone maintenance and antiosteoporosis effects of estrogens. It will be of interest in subsequent studies to use the preferentially regulated genes we have identified to explore the molecular basis for the ER subtype selectivity in the regulation of these genes.
| Acknowledgments |
|---|
| Footnotes |
|---|
F.S. and D.H.B. contributed equally and should be considered equal first authors.
Abbreviations: BMP6, Bone morphogenetic protein 6; CAT, chloramphenicol acetyltransferase; Cox-2, cyclooxygenase-2; CS, confidence score; DPN, diarylpropionitrile; E2, 17ß-estradiol; ERE, estrogen response element; ER, estrogen receptor; FBS, fetal bovine serum; FC, fold change; GREB, gene regulated by estrogen in breast cancer 1; NHERF, Na+/H+ exchanger regulatory factor; NKG, natural killer cell lectin-type receptors; PDZK1, PDZ domain containing 1; PPT, propylpyrazoletriol; PTHLH, PTH-like hormone, also known as PTHrP; WISP2, WNT1-inducible signaling pathway protein 2.
Received December 11, 2003.
Accepted for publication March 9, 2004.
| References |
|---|
|
|
|---|
/ß)-dependent activation at GC-rich (Sp1) promoter elements. J Biol Chem 275:53795387
is developmentally regulated during osteoblast differentiation and contributes to selective responsiveness of gene expression. Endocrinology 139:20482057
and ß mRNA during differentiation of human osteoblast SV-HFO cells. Endocrinology 138:50675070
receptor chimeras. Endocrinology 139:45134522
regulates gene expression in the kidney. Endocrinology 144:701710
C(T)) method. Methods 25:402408[CrossRef][Medline]
-induced interleukin-8 production in human bone marrow stromal and osteoblast-like cells. Calcif Tissue Int 55:1620[CrossRef][Medline]
in human melanoma cells. FEBS Lett 515:137140[CrossRef][Medline]
-selective agonists. J Med Chem 43:49344947[CrossRef][Medline]
and estrogen receptor-ß: correlations with biological character and distinct differences among SRC coactivator family members. Endocrinology 141:35343545