help button home button Endocrine Society JCEM
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Journal of Clinical Endocrinology & Metabolism , doi:10.1210/jc.2004-1423
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hanifi-Moghaddam, P.
Right arrow Articles by Blok, L. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hanifi-Moghaddam, P.
Right arrow Articles by Blok, L. J.
Related Collections
Right arrow Female Endocrinology
The Journal of Clinical Endocrinology & Metabolism Vol. 90, No. 2 973-983
Copyright © 2005 by The Endocrine Society

Molecular Portrait of the Progestagenic and Estrogenic Actions of Tibolone: Behavior of Cellular Networks in Response to Tibolone

P. Hanifi-Moghaddam, S. C. J. P. Gielen, H. J. Kloosterboer, M. E. De Gooyer, A. M. Sijbers, A. J. van Gool, M. Smid, M. Moorhouse, F. H. van Wijk, C. W. Burger and L. J. Blok

Departments of Reproduction and Development (P.H.-M., L.J.B.), Obstetrics and Gynecology (S.C.J.P.G., F.H.v.W., C.W.B.), Urology (M.S.), and Bioinformatics (M.M.), Erasmus MC, 3000 CA Rotterdam, The Netherlands; and Research and Development Laboratories, N.V. Organon (H.J.K., M.E.D.G., A.M.S., A.J.v.G.), 5340 BH Oss, The Netherlands

Address all correspondence and requests for reprints to: Dr. P. Hanifi-Moghaddam, Department of Reproduction and Development, Erasmus Medical Center, P.O. Box 1738, 3000 DR Rotterdam, The Netherlands. E-mail: p.hanifi_moghaddam{at}erasmusmc.nl.


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Tibolone is a synthetic steroid with estrogenic effects on brain, vagina, and bone without stimulating the endometrium. During tibolone treatment, it is thought that the progestagenic properties of tibolone stimulate cell differentiation, which effectively counterbalances the growth-stimulating effects of the estrogenic properties of tibolone. The objective of this study was to characterize the expression profile that reflects the endometrial responses to the separated estrogenic (growth-inducing) and progestagenic (growth-inhibiting) actions of tibolone, thus gaining insight into the counteracting effect of these properties of tibolone on the endometrium. The estrogenic action of tibolone was studied in the estrogen-responsive ECC1 cell line (expressing estrogen receptor {alpha}), and the progestagenic action was studied in the progesterone-responsive cell line Ishikawa PRAB-36 (expressing PRA and PRB). The data showed that the progestagenic and estrogenic effects of tibolone produce different expression profiles with a narrow overlap in genes; however, both properties modulate the same biological processes. The final genetic network analysis indicated that the estrogenic effect of tibolone is potentially counterbalanced by the progestagenic metabolite of tibolone via differential regulation of similar cellular processes. For example, both progestagenic and estrogenic properties stimulate proliferation, but they exert the opposite effect on apoptosis. The apoptosis network was stimulated by the progestagenic properties of tibolone; in contrast, the estrogenic effect of tibolone suppressed the apoptosis network. The current results indicate that this differential regulation is realized through modulation of a different group of genes and rarely via contraregulation of the same set of genes.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
TIBOLONE IS A steroidogenic compound, which is used for hormone replacement therapy and prevention of osteoporosis in several countries. During passage through the intestine and liver, tibolone is converted into its 3{alpha}-reduced derivative (3{alpha}-OH-tibolone), its 3ß-reduced derivative (3ß-OH-tibolone), and its {Delta}4-isomer ({Delta}4-tibolone). The 3{alpha}- and 3ß-hydroxy metabolites bind exclusively to the estrogen receptor (ER), whereas the {Delta}4-isomer and tibolone activate the progesterone (PR) and androgen (AR) receptors (1, 2).

The 3{alpha}- and 3ß-OH metabolites of tibolone have a half-life of about 7 h. These metabolites are also present in large amounts as sulfated metabolites and serve as a source of estrogenic activity in tissues where sulfatases are active. Surprisingly, we and others have found that the 3ß-OH metabolite possesses progestagenic activity (3, 4), which may be due to the conversion of the 3ß-OH metabolite back into tibolone (which has progestagenic properties), whereas subsequent metabolism of tibolone into the {Delta}4-isomer is also possible. The {Delta}4-isomer is found in low concentration and is only detectable for 6 h. In the endometrium, tibolone is locally converted into its progestagenic {Delta}4-isomer (5). Progestagens can induce sulfotransferase (3, 6, 7) in the endometrium, and because the 3-OH metabolites are good substrates for at least three isoforms of sulfotransferase, this may also contribute to abrogation of the estrogenic response. It appears that the progestagenic response in the endometrium is dominant, which is confirmed by the induction of many progesterone-sensitive parameters in endometrial cells, as indicated below.

In animal and human studies, it has been shown that tibolone has tissue-specific effects; it acts as an estrogenic compound on bone and on the central nervous system, whereas no stimulation is found on the breast or endometrium (1, 8, 9, 10, 11, 12, 13). In the endometrium, under physiological circumstances, estrogenic cell growth stimulation is counterbalanced by progesterone-induced cell differentiation (3). During tibolone treatment, it is thought that local formation of the {Delta}4-isomer, which has progestagenic and androgenic activities, and sulfation of the estrogenic metabolites (1, 14) effectively counterbalance any putative adverse estrogenic effects of the 3{alpha}- and 3ß-reduced derivatives of tibolone (5, 14, 15, 16). Knowledge about the biological processes through which tibolone exerts its counterbalancing actions on the endometrium, however, is lacking (10, 14).

The main objective of the present study was to characterize the transcription profile that reflects the endometrial response to the estrogenic (growth-inducing) and progestagenic (growth-inhibiting) metabolites of tibolone. Using this information, it was our goal to gain insight into the counteracting effect of the different activities of tibolone. The study of estrogenic and progestagenic actions of tibolone separately in one cell line that expresses both ER and PR would be very difficult due to cross-talk and interactions between the receptors. Therefore, we used two endometrial cell lines, which respond only to progesterone or estrogen: PRAB-36 and ECC1. The Ishikawa PRAB-36 cell line is progesterone responsive, and the ECC1 cell line is estrogen responsive. Like progesterone, tibolone inhibits the growth of PRAB-36 cells, and both estrogen and tibolone are able to stimulate the growth of ECC1 cells. After testing of the cell lines for their ability to convert tibolone into its metabolites, the cell lines were treated with tibolone for 6, 24, and 48 h. The expression profiles of the cell lines at different time points were generated using a 9600 cDNA microarray. The differentially expressed genes were analyzed to determine which biological processes were (significantly) regulated in each cell line and at each time point. Finally, to gain insight into the counterbalancing effect of tibolone, we generated a genetic network, which indicated the behavior of cellular processes in response to the progestagenic and estrogenic activities of tibolone.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Cell lines

Ishikawa (17) and ECC1 cells (18) were derived from two different well differentiated human endometrial adenocarcinomas. The parental Ishikawa cell line (17) does not express detectable levels of PR, ER, or AR (19, 20). Using stable transfection techniques we have generated modified Ishikawa cell lines, which express PRA and PRB (19). The human endometrial carcinoma ECC1 cell line contains ERs (mostly ER{alpha}), and in our laboratory these cells are estrogen responsive in growth. Furthermore, this cell line contains low amounts of PR and AR, but no growth response to progesterone or androgen administration was observed (20). Both cell lines were maintained in DMEM/Ham’s F-12 medium (Invitrogen Life Technologies, Gaithersburg, MD), containing 5% fetal bovine serum (Perbio Science, Helsingborg, Sweden) and antibiotics, in a humid atmosphere in the presence of 5% CO2. Growth studies were performed and conversion of tibolone into its metabolites was measured as described by Blok et al. (4).

RNA isolation and gene expression profiling

Total RNA was isolated according to the method described by Blok et al. (4). Microarray hybridizations were performed with the Human Unigene 1 microarray (Incyte Genomics, Inc., St. Louis, MO), containing 9.600 cDNAs. Both labeling and hybridizations were performed according to previously described protocols (21).

Statistical analysis

Because most subsequent analyses depend on the accuracy of assessment of the expression ratios between the treatment and control groups, one of the most critical steps in microarray profiling experiments is the accurate assessment of these ratios. Spots were selected for additional analysis when the signal to background ratio was greater than 2.5 and the area occupied by the signal or control on the cDNA spot on the array was greater than 40%. After signal quantification, a signal correction algorithm was used to correct for systematic differences between the Cy3 and Cy5 labels. This algorithm applied a second order polynomial regression model to the data by fitting a parabola through log-transformed Cy3 vs. Cy5 intensities. The resulting new data from the regression model were taken as the new gene expression ratios. Genes that were significantly differentially expressed in the treated vs. the control sample were selected by a statistical method derived from box plots, which are widely used to visualize the overall shape of a dataset (22). In our experiments the expression of the majority of genes is not expected to differ between the treatment and control samples, and their ratios will be in the center of the distribution of expression ratios. Therefore, we computed the first (Q1) and third (Q3) quartiles of the distribution of the residuals of the regression model and the interquartile range (IQR) of the distribution as a measure for the variation in the expression ratios of nondifferentially expressed genes. Then, in analogy with box plots, an inner fence was set at Q1 – (1.58 x IQR) and Q3 + (1.58 x IQR). Using these criteria, genes within the inner fence have a probability of P = 0.995 to be nondifferentially expressed, whereas the outlier group will harbor the differentially expressed genes. Genes were selected for additional study if they fell outside the inner fence (fold change, ≥1.3; P ≤ 0.005). The complete list of differentially regulated genes is available on our website (http://www.eur.nl/fgg/rede/hanifi_moghaddam/differentaly_regulated_genelist.xls).

Overlap of differentially expressed genes between two cell lines and different times points

To identify co-occurrence (or overlap) of differentially expressed genes at different time points in both cell lines, Venn Mapper was used (23). Venn Mapper is a program that compares heterologous microarray datasets, based on the number of common, differentially expressed genes. The application loads microarray data (gene expression ratios) and determines which genes are up- or down-regulated by a user-defined ratio cut-off level. For each experiment, lists of differentially expressed genes are computed. Every list will be compared with every other list, and the number of co-occurring genes will be calculated.

Analysis of functional gene categories

For functional classification of genes and identification of enrichment and depletion of gene categories, four programs were used: GO-Mapper (24), Ease (25), GoMiner (26), and FatiGo (27). This analysis was performed for all genes regulated in three time points and at each time point separately.

A custom visualization system was developed to display the GO classification of genes (using the GO-Mapper output) found to be of interest. By means of hierarchical tree, the context of each GO term is maintained relative to both its parents and children (if any). The ratios of the numbers of genes between the two cell lines at each of the three different time points is presented graphically. The absolute numbers of up- and down-regulated genes are included in the graphic: grouped in six pairs, split as two per cell line, three per time point.

Clustering and genetic network building

The set of genes significantly differentially expressed (P < 0.005; fold change, 1.3) was clustered by the similarity of their expression profiles using a Self-Organizing Three Algorithm (28). For the construction of a genetic network, the Pathway Assist version 2.5 (Ariadne Genomics, Inc., Rockville, MD) was used. The database (updated April 2004) of the Pathway Assist contains biological knowledge represented in a formalized form with the focus on how proteins, cellular processes, and small molecules interact, modify, and regulate each other. The Pathway Assist provides a method for searching objects individually by keyword, string, or attributes. These include, for example, type, effect (positive, negative, or unknown), mechanism (transcription or phosphorylation), tissue type, biological process, belonging to cell structure, and others. In addition, the program provides a search method for objects such as proteins, small molecules, and enzymes. The complete databases of Kyoto Encyclopedia of Genes and Genomes (www.genome.ad.jp/kegg), Database of Interacting Proteins (dip.doe-mbi.ucla.edu), Bimolecular Interaction Network Database (bind.mshri.on.ca), and Gene Ontology (www.geneontology.org) were imported into the pathway Assist database.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Evaluation of the model cell lines

Responsiveness to the hormones: cell growth study. Tibolone is a compound that, through metabolic conversion, can act on ER-, AR-, and PR-positive cells. In previous studies we have shown that PRAB-36 cells express PRA and PRB at high levels, and no expression of ER and AR could be detected. In ECC1 cells, high expression of ER{alpha} was observed, whereas expression of PR and AR was very low. Furthermore, in the currently used cell lines, effects of androgens on growth were never observed (data not shown). Therefore, addition of tibolone to the progesterone-responsive Ishikawa cell line (PRAB-36) was expected to result in {Delta}4-isomer-induced growth inhibition, whereas in the estrogen-responsive ECC1 cell line, tibolone treatment was expected to result in 3{alpha}-OH- and/or 3ß-OH-tibolone-induced growth. It was observed that the progestagenic growth-inhibiting effect of tibolone administration was similar to the effect of medroxyprogesterone acetate (MPA) in the PRAB-36 cell line (Fig. 1Go). Furthermore, in ECC1 cells, a comparable growth-inducing effect was measured for estradiol and tibolone (Fig. 1Go).



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 1. Effect of MPA, estradiol (E2), and tibolone administration on the growth of progesterone-responsive Ishikawa cells (PRAB-36) and on estrogen-responsive ECC1 cells. Cells were cultured in the presence of tibolone (100 nM), MPA (1 nM), or E2 (1 nM) for 10 d. The physiological effective doses of tibolone, E2, and MPA on cell growth were determined in a dose-dependent manner, as described previously (4 20 ). Cell growth is expressed as a percentage of control growth. Each bar represents the average of four different incubations ± SD within one representative experiment.

 
Conversion of tibolone into its metabolites. To have an indication which metabolites were present during incubations with tibolone, the conversion of tibolone into its metabolites was measured in both cell lines by HPLC analyses. The results in Table 1Go show that after 3 d of culture, tibolone is converted into the progestagenic/androgenic {Delta}4-isomer (80% in PRAB-36 and 70% in ECC1 cells, respectively) and the estrogenic 3ß-OH metabolite (19% in PRAB-36 and 27% in ECC1 cells, respectively). No 3{alpha}-OH-tibolone was detected.


View this table:
[in this window]
[in a new window]
 
TABLE 1. Conversion of tibolone into its metabolites ({Delta}4-tibolone, 3{alpha}OH-tibolone, and 3ßOH-tibolone) after 3 d of culture

 
Similarity in basal gene expression profile between cell lines. PRAB-36 and ECC1 cell lines were derived from two patients at different points in time. The question was raised of whether the genetic backgrounds of the cell lines, rather than the effects of hormonal treatment, would have an important impact on the generated expression profiles. This question was addressed by calculating the expression ratios between untreated PRAB-36 and ECC1 cells for all genes on the array (calculations were performed on three independent sets of control data from the PRAB-36 cell line and three other independent sets of control data from the ECC1 cell line). If the genetic background does not have a major impact on gene expression, than the expression ratios between genes expressed in PRAB-36 and in ECC1 cells should be approximately 1. It was observed that 99.48% of ratios were between 0.8 and 1.2. In total, only 32 genes that showed (no significant) differential regulation [for 11 genes, fold change between 2 and 3.5 (P > 0.05); for 21 genes, fold change between –2 and –3.5 (P > 0.05)] in the control samples were excluded from subsequent analysis. These data indicate that the two cell lines have a very similar basal expression pattern.

Analyses of gene expression profiles

Identification of up- and down-regulated genes. In Fig. 2Go, all raw data obtained are depicted in so-called box plots to indicate the quality of the array data. Because microarray experiments are prone to systematic variation, the high correlation between logG and logR shows that measured fluorescent intensities arise from biological differences between the RNA samples rather than from systematic biases.



View larger version (25K):
[in this window]
[in a new window]
 
FIG. 2. Box plot data of all six array analysis. The left panels represent PRAB-36 cells cultured for 6, 24, or 48 h with 100 nM tibolone compared with control cells. The right panels represent ECC1 cells cultured for 6, 24, or 48 h with 100 nM tibolone compared with control cells. For each plot the y-axis represents the intensity of the control signal (intensity P2), and the x-axis represents the intensity of the tibolone-stimulated signal (intensity P1).

 
When all data for each cell line were compiled (6, 24, and 48 h of tibolone treatment), 570 genes in total were significantly regulated by tibolone (fold change, ≥1.3; P < 0.005). In both cell lines, the majority of genes were down-regulated by tibolone (Table 2Go), and the number of genes regulated in ECC1 cells was higher than that in PRAB-36 cells [ECC1, 385 (309 + 76); PRAB-36, 261 (185 + 76)].


View this table:
[in this window]
[in a new window]
 
TABLE 2. Number of genes regulated by tibolone treatment in PRAB-36 and/or ECC1 cells at all three time points

 
In Table 3AGo, it is shown that in the PRAB-36 cell line, about 54 (14 + 40) of 112 (19 + 93) genes up- and down-regulated at 6 h were also regulated after 24 h, and 58 (13 + 45) of these 112 genes were also regulated after 48 h. The overlap in regulated genes between 24 and 48 h is high (109 genes). In the ECC1 cell line, about 45% of the regulated genes showed regulation at more than one time point. The large overlap in the number of genes between different time points in a cell line suggests that a number of biological processes are activated over a longer period of time. In Table 3BGo, it is shown that only a small number of genes are regulated by both progestagenic and estrogenic properties of tibolone. Interestingly, only 13% (76 genes; Table 2Go) of the significantly regulated genes were coregulated by both progestagenic and estrogenic properties of tibolone. Opposite regulation was found in only 3.9% of the total of regulated genes (22 genes; Table 2Go).


View this table:
[in this window]
[in a new window]
 
TABLE 3. Matrix of co-occurring genes at three time points and in both cell lines

 
Regulation of ER-, PR-, and AR-sensitive genes by tibolone. To identify which tibolone-regulated genes are known to be regulated by estrogens, progestagens, or androgens, we have composed a list of genes associated with these hormones and their receptors. A search of the literature and public databases (including our own database) for ER-, PR-, and AR-regulated genes resulted in a list containing approximately 1500 genes. From these 1500, 922 genes were present on the Incyte array. When our 570 tibolone-regulated genes were compared with the 922 genes, 229 genes were found to overlap. As indicated in Table 4Go, of these 229 overlapping genes, 58 were regulated in PRAB-36 cells, 131 were regulated in ECC1 cells, and 40 were regulated in both cell lines. The results in Table 4Go show that the {Delta}4-isomer-mediated effect (progestagenic effect on PRAB-36 cells) of tibolone predominantly affects known PR-regulated genes, whereas the 3ß-OH-tibolone-mediated effect (estrogenic effect on ECC1 cells) of tibolone predominantly affects known ER-regulated genes. Furthermore, in the current experiments, androgenic effects of the {Delta}4-isomer of tibolone are scarcely apparent; only four AR-responsive genes were identified. This indicates that, as observed in growth experiments, the androgenic effect of tibolone is negligible in our cell lines.


View this table:
[in this window]
[in a new window]
 
TABLE 4. ER-, PR-, and AR-regulated genes derived from the literature or public databases that are present on the 9600 cDNA array (n = 922) were compared to the 570 tibolone-regulated genes

 
In the current experiments we can distinguish between two groups of genes; 229 of the tibolone-regulated genes are known, from the literature or public databases, to be regulated via ER, PR, or AR, and 341 genes are not known to be regulated by these hormone receptors.

Progestagenic and estrogenic effects of tibolone on cellular processes: analysis of gene categories, clustering, and building genetic networks. As shown in Tables 2Go and 3BGo, only 13% of the genes are coregulated, and only 3.9% of the genes are contraregulated between the two cell lines. This indicates that the progestagenic and estrogenic properties of tibolone regulate very distinct sets of genes. The question was then raised of whether the distinct expression profiles generated by the progestagenic and estrogenic properties of tibolone act at the same biological processes or at different biological processes.

To answer this question, for each cell line the regulated genes at three time points were compiled. Subsequently, the genes were categorized into their associated biological processes using four different programs. Table 5Go data show an example of (GO-Mapper) GO categories and the number and percentage of genes in each category. The percentage of genes that fell into each category of biological processes was compared between the two cell lines to find over- or underrepresentation of gene numbers in different categories. If the progestagenic and estrogenic properties of tibolone affect different biological processes, then one of the cell lines should show depletion or enrichment in a certain category. Significant enrichment or depletion of a category suggests the activation or repression of a biological process. However, no significant enrichment or depletion in categories was found between the two cell lines using the four programs (the complete data of GO-Mapper analysis and the visualization of functional classification are available on our website: www.eur.nl/fgg/rede/hanifi_moghaddam/functional_classification.xls, and www.eur.nl/fgg/rede/hanifi_moghaddam/main.html. In addition, repeating this analysis to compare time points within and between cell lines showed no significant differences between both cell lines or between different time points regarding over- or underrepresentation of gene numbers in different categories. Taking into account the narrow overlap of regulated genes in both cell lines, these results indicate that the progestagenic and estrogenic properties of tibolone modulate similar biological processes and do so through modulation of a different set of genes.


View this table:
[in this window]
[in a new window]
 
TABLE 5. Tibolone-regulated genes in the PRAB-36 and ECC1 cell lines were grouped into biologically functional categories (biological processes)

 
This finding was also demonstrated in the Fig. 3Go, which shows cluster analyses of the genes involved in apoptosis (Fig. 3AGo), proliferation and differentiation (Fig. 3BGo), and transcription (Fig. 3CGo). The drawback of clustering, however, is that the cluster program groups together genes regardless of their function and without reference to the valuable biological information available in public databases. To elucidate how the counterbalancing effect of the estrogenic and progestagenic activities of tibolone modulates cellular processes, a genetic network analysis was performed. Genetic network models are typically constructed to describe which genes regulate each other and how environmental inputs could affect gene expression.



View larger version (34K):
[in this window]
[in a new window]
 
FIG. 3. Clustering of genes classified into the functional categories of apoptosis (A), proliferation and differentiation (B), and transcription (C). From left to right: columns 1–3, PRAB-36 (6, 24, and 48 h of tibolone); columns 4–6, ECC1 (6, 24, and 48 h of tibolone). Black, Unchanged expression; green, down-regulated genes; red, up-regulated genes.

 
Of 570 genes, 212 could be classified in 40 cellular processes. The rest of genes either could not be placed in any cellular processes or they formed small discontented networks. Due to the high complexity of the generated genetic network, we selected four cellular processes that play major roles in cell growth and differentiation. The genetic network represented in Fig. 4Go represents transcription, apoptosis, differentiation, and proliferation (in total, 134 genes could be annotated to these cellular processes). The genetic network approach depicted in Fig. 4Go also confirmed that the PR- and ER-mediated effects of tibolone on cellular processes indeed occur predominantly via modulation of a different set of genes and rarely via contraregulation of the same set of genes.



View larger version (46K):
[in this window]
[in a new window]
 
FIG. 4. Behavior of cellular processes in the endometrial cancer cell lines in response to tibolone administration. A, Behavior of cellular processes in the PRAB-36 cell line in response to tibolone. B, Behavior of cellular processes in the ECC1 cell line in response to tibolone. Four cellular processes were studied: apoptosis, proliferation, differentiation, and transcription. Yellow circles, Unchanged expression; green circle, down-regulated expression; red circle, up-regulated expression. Type of connections: effect on regulation (dashed gray line, gray box), effect on expression (blue, solid blue line, blue box). A high-resolution copy of this figure and a gene list of the network are available on our web site: www.eur.nl/fgg/rede/hanifi_moghaddam/network.htm.

 
Finally, to comprehend the response of the network to the progestagenic and estrogenic effects of tibolone, we calculated the net suppressing or inducing effects of the progestagenic and estrogenic properties of tibolone on apoptosis and proliferation (Table 6Go). In this analysis we have taken into account that genes can either induce or suppress a cellular process. The results in Table 6Go suggest that both progestagenic and estrogenic properties of tibolone stimulate proliferation, but they have an opposite effect on apoptosis. The progestagenic effect of tibolone stimulates apoptosis, whereas the estrogenic effect of tibolone inhibits apoptosis. It should be noted, however, that a certain caution when interpreting these data is needed, because the gene regulatory information used is not based on endometrium-specific data. Furthermore, all genes have been annotated a similar weight, but we did not take into account that some genes are key regulators of a process while others are not. Results from the network analysis will be evaluated in Discussion.


View this table:
[in this window]
[in a new window]
 
TABLE 6. The numbers of up- and down-regulated genes involved in apoptosis and proliferation, and their net inducing or suppressing effect on each process are shown

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Hormone signaling often converges or is intertwined with signaling by other hormones, growth factors, or cytokines. For example, estrogen induces the expression of PR; therefore, many of the in vivo effects attributed to progesterone are, to some degree, dependent on estrogens. In contrast, progesterone attenuates the proliferative effect of estrogens in the uterus partly through down-regulation of the ER (29). In the present cell models, PRA and PRB are predominantly expressed in the PRAB-36 cell line, and ER in the ECC1 cell line; therefore, we did not taken into account any possible direct PR- and ER-mediated cross-talk.

For the present study, we used two well differentiated endometrial adenocarcinoma cell lines, PRAB-36 and ECC1. The PRAB-36 cell line expresses PR and is inhibited in growth by progestagenic compounds, whereas the ECC1 cell line expresses ER and is stimulated in growth by estrogenic compounds. Furthermore, both cell lines were shown to convert tibolone partly into its progestagenic {Delta}4-metabolite and its estrogenic 3ß-OH-metabolite. Because the cell lines are derived from different patients, their genetic backgrounds could possibly interfere with their hormone-stimulated gene expression profiles. Therefore, the basal gene expression profiles of the two cell lines were compared. A very high similarity between the basal expression profiles of the cell lines suggested that the generated expression portrait under effect of hormones is probably hardly affected by differences in the genetic constellation of the cell lines. Overall, these results indicated that the PRAB-36 and ECC1 cell lines are suitable models to study the separated progestagenic and estrogenic properties of tibolone administration to endometrial cells.

Analysis of the gene expression profiles generated in the tibolone-stimulated cell lines revealed that one group of regulated genes had already been recognized in the literature to be regulated via ER or PR (as indicated in Table 4Go), whereas the other group of genes was not known to be regulated by these hormone receptors. One explanation could be that tibolone regulates a tibolone-specific set of genes; however, it is also possible that the information in the public domain is not complete in this respect.

The most striking finding in the current investigations was that despite distinct expression profiles, both estrogenic and progestagenic properties of tibolone regulate similar biological processes. This suggests that the balance between the estrogenic, growth-inducing activities and the progestagenic, differentiative, and consequently growth-inhibiting activities of tibolone in the current model cannot be explained by opposed regulation of the same genes or by modulation of completely different biological processes. Most interestingly, the analysis of genetic networks and overlying the expression data on these networks (as shown in Fig. 4Go) revealed that the counterbalancing effect of tibolone on cellular processes is realized through regulation of a different set of genes, not via contraregulation of the same set of genes. Going one step further, it was shown that the actual modulation of a process is realized by differential regulation of genes with inducing or suppressing effects. Therefore, the number of up- or down-regulated genes may not be enough to analyze whether a pathway is activated or inhibited.

As an example, we discuss below the differential regulation of apoptosis and transcription by PR- and ER-mediated signaling of tibolone. The selection of the discussed genes is based on the number of connections with other genes in the network, and the consistency of their expression. Genes with a high connectivity play a vital role in the stability of a network or any other complex system (30, 31).

Differential modulation of apoptosis

Many physiological growth control mechanisms that govern cell proliferation and tissue homeostasis are linked to apoptosis. The genetic network, as depicted in Fig. 4Go, suggests a central role for apoptosis in the regulation of proliferation and differentiation by tibolone. By calculating the net effect of each gene in the network on apoptosis and proliferation, we observed that both progestagenic and estrogenic properties of tibolone stimulate proliferation, but they have opposite effects on apoptosis. The progestagenic effect of tibolone stimulates apoptosis, whereas the estrogenic effect of tibolone inhibits the apoptosis. In fact, Smid-Koopman (Smid-Koopman, E., L. C. M. Kuhne, E. E. Hanekamp, S. C. J. P. Gielen, P. E. De Ruiter, J. A. Grootegoed, T. J. M. Helmerhorst, C. W. Burger, A. O. Brinkmann, F. J. Huikeshoven, L. J. Blok, submitted for publication) showed that long-term culture (28 d) in the presence of MPA did indeed induce apoptosis and brought cell growth to a halt. A similar experiment in PRAB-36 cells with tibolone also induced apoptosis (not shown). In contrast to this, recently it was shown that tibolone did not affect the rate of apoptosis in the postmenopausal endometrium (32). These results seem to contradict each other, but may be explained when we consider that progestagen-induced ({Delta}4-tibolone) apoptosis may in vivo be counterbalanced by the estrogenic metabolite of tibolone.

This counterbalancing effect may be seen in the regulation of IGF-binding proteins (IGFBPs), which are differently regulated in each cell line. Down-regulation of IGFBP3 was observed only in PRAB-36 cells, whereas in ECC1 cells, IGFBP-4 and IGFBP-5 were up-regulated. The IGFs, their receptors, and their binding proteins play key roles in regulating cell proliferation and apoptosis. Noteworthy, in the cell lines investigated in this study, tibolone did not show any effect on the expression of either IGFs or their receptors. This suggests that any effect resulting from regulation of IGFBPs by tibolone occurs independently of IGFs or their receptors. A series of studies performed over the past years has established that certain actions of the IGFBPs are indeed IGF independent. IGFBP-3 and IGFBP-5, in particular, have been shown to have effects on proliferation, migration, and sensitivity to apoptosis that are independent of IGF signaling per se (33, 34, 35). These data combined with the present findings suggest that tibolone may differentially regulate apoptotic signaling pathways, partly through differential regulation of several IGFBPs (17).

Differential regulation of the transcription machinery

Tibolone metabolites can bind to PR ({Delta}4-isomer-mediated effect) and ER (3ß-OH-tibolone-mediated effect) and modulate PR- and ER-mediated transcription. Both steroid receptors are nuclear transcription factors and can regulate target gene transcription directly by binding to the promoters of these genes or indirectly through regulation of the expression of one or more transcription factors that will interact with the promoter regions of other target genes.

Binding of a receptor to a hormone response element is associated with the recruitment of coregulators and subsequent histone acetylation (36, 37). From the list of known steroid receptor coregulators, nuclear receptor interacting protein 1 (NRIP1 or RIP140) and thymine-guanine-interacting factor (TGIF) were regulated by tibolone (Fig. 4Go, Transcription network). NRIP1 is a transcription coactivator that specifically interacts with the hormone-dependent activation domain AF2 of nuclear receptors, whereas TGIF is a transcription corepressor that interacts with HDAC1 and CBP/P300, thereby inhibiting nuclear hormone receptor-mediated transcription.

Tibolone-induced up-regulation of NRIP1 (in ECC1 cells, at all three time points) and down-regulation of TGIF (in PRAB-36 cells, at all three time points) may point to the fact that these genes are necessary determinants for the estrogenic and progestagenic actions of tibolone. A similar finding was reported in a recent paper showing that the estrogen-like activity of tamoxifen in the uterus and in endometrial cells requires a high level of steroid receptor coactivator 1 (38).

This observation shows that in the endometrium, the specific coregulator requirements for estrogenic and progestagenic properties of tibolone are distinct. Because coregulators are themselves involved in multiple signal transduction pathways, the ability of each tibolone metabolite to modulate the expression of a specific coregulator could be a key factor underlying the counterbalancing character of tibolone. It has to be mentioned that besides regulation of the expression of coregulators, other factors, such as differential regulation of coregulator activity, also play an important role in the determination of the biocharacter of a (synthetic) hormone. This, however, could not be studied using the current experimental set-up.

As mentioned above, PR and ER can indirectly regulate target gene transcription through a mechanism involving other transcription factors (39). The transcription factors Specificity Protein 3 and B-cell lymphoma 6 protein, for example, are both down-regulated in PRAB-36 cells, whereas myeloblastosis protooncogene protein (MYB) was up-regulated in the ECC1 cell line. SP3 is a bifunctional transcription factor (transcription activator or repressor) (39, 40), which is involved in recruitment of other transcription factors to promoter regions (41); BCL6 is a zinc finger transcription factor with cell cycle regulatory functions that can inhibit apoptosis (42). MYB is a transcriptional activator, and recently it was shown that the expression of IGFBP-5 is regulated by MYB (43).

Concluding remarks

The biological effect the tibolone treatment exerts on an endometrial cell is orchestrated by the selective expression of genes in response to tibolone’s estrogenic and progestagenic properties. It seems logical that the counterbalancing effect can be realized through regulation of the same biological processes, but it was intriguing to find that different genes are involved in this process and not counterregulation of the same genes.

The study of the expression dynamics of many genes simultaneously enabled us to dissect the complex genetic network that controls the patterns and rhythms of gene expression in the cells in response to tibolone treatment. However, we are aware that constructing a genetic network, based on experimental data and literature-driven knowledge, is a simplification of the cellular system, and that the processes selected in this study (transcription, apoptosis, differentiation, and proliferation) are only part of the total biochemical network in a cell. Furthermore, not taking into account the progestagenic and estrogenic effects on the expression of each others receptors introduces an additional simplification of tibolone working mechanism. In contrast, we do think that our work helps in understanding at least partly the behavior of an endometrial cellular system in response to the progestogenic and estrogenic activities of tibolone.


    Footnotes
 
This work was supported by N.V. Organon.

First Published Online November 30, 2004

Abbreviations: AR, Androgen receptor; ER, estrogen receptor; IGFBP, IGF-binding protein; IQR, interquartile range; MPA, medroxyprogesterone acetate; NRIP, nuclear receptor interacting protein; 3{alpha}-OH-tibolone, 3{alpha}-reduced derivative of tibolone; 3ß-OH-tibolone, 3ß-reduced derivative of tibolone; PR, progesterone receptor; PRA, progesterone receptor A; PRB, progesterone receptor B; TGIF, thymine-guanine-interacting factor.

Received July 20, 2004.

Accepted November 16, 2004.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Kloosterboer HJ 2001 Tibolone: a steroid with a tissue-specific mode of action. J Steroid Biochem Mol Biol 76:231–238[CrossRef][Medline]
  2. de Gooyer ME, Deckers GH, Schoonen WG, Verheul HA, Kloosterboer HJ 2003 Tibolone is not converted by human aromatase to 7{alpha}-methyl-17{alpha}-ethynylestradiol (7{alpha}-MEE): analyses with sensitive bioassays for estrogens and androgens and with LC-MSMS. Steroids 68:21–30[CrossRef][Medline]
  3. Markiewicz L, Gurpide E 1990 In vitro evaluation of estrogenic, estrogen antagonistic and progestagenic effects of a steroidal drug (Org OD-14) and its metabolites on human endometrium. J Steroid Biochem 35:535–541[CrossRef][Medline]
  4. Blok LJ, De Ruiter PE, Kuhne EC, Hanekamp EE, Grootegoed JA, Smid-Koopman E, Gielen SC, De Gooyer ME, Kloosterboer HJ, Burger CW 2003 Progestogenic effects of tibolone on human endometrial cancer cells. J Clin Endocrinol Metab 88:2327–2334[Abstract/Free Full Text]
  5. Tang B, Markiewicz L, Kloosterboer HJ, Gurpide E 1993 Human endometrial 3ß-hydroxysteroid dehydrogenase/isomerase can locally reduce intrinsic estrogenic/progestagenic activity ratios of a steroidal drug (Org OD 14). J Steroid Biochem Mol Biol 45:345–351[CrossRef][Medline]
  6. Tseng L, Gurpide E 1975 Induction of human endometrial estradiol dehydrogenase by progestins. Endocrinology 97:825–833[Abstract/Free Full Text]
  7. Kohorn EI, Tchao R 1969 Conversion of proliferative endometrium to secretory endometrium by progesterone in organ culture. J Endocrinol 45:401–405[Abstract/Free Full Text]
  8. Tax L, Goorissen EM, Kicovic PM 1987 Clinical profile of Org OD 14. Maturitas Suppl 1:3–13
  9. Ederveen AG, Kloosterboer HJ 2001 Tibolone exerts its protective effect on trabecular bone loss through the estrogen receptor. J Bone Miner Res 16:1651–1657[CrossRef][Medline]
  10. Moore RA 1999 Livial: a review of clinical studies. Br J Obstet Gynaecol 106(Suppl 19):1–21
  11. Sendag F, Terek MC, Ozsener S, Oztekin K 2001 Mammographic density changes in postmenopausal women using tibolone therapy. Int J Gynaecol Obstet 74:63–64[CrossRef][Medline]
  12. Christodoulakos GE, Lambrinoudaki IV, Vourtsi AD, Panoulis KP, Kelekis DA, Creatsas GC 2002 Mammographic changes associated with raloxifene and tibolone therapy in postmenopausal women: a prospective study. Menopause 9:110–116[CrossRef][Medline]
  13. Lundstrom E, Christow A, Kersemaekers W, Svane G, Azavedo E, Soderqvist G, Mol-Arts M, Barkfeldt J, von Schoultz B 2002 Effects of tibolone and continuous combined hormone replacement therapy on mammographic breast density. Am J Obstet Gynecol 186:717–722[CrossRef][Medline]
  14. de Gooyer ME, Kleyn GT, Smits KC, Ederveen AG, Verheul HA, Kloosterboer HJ 2001 Tibolone: a compound with tissue specific inhibitory effects on sulfatase. Mol Cell Endocrinol 183:55–62[CrossRef][Medline]
  15. Rymer J, Fogelman I, Chapman MG 1994 The incidence of vaginal bleeding with tibolone treatment. Br J Obstet Gynaecol 101:53–56[Medline]
  16. von Dadelszen P, Gillmer MD, Gray MD, McEwan HP, Pyper RJ, Rollason TP, Wright A 1994 Endometrial hyperplasia and adenocarcinoma during tibolone (Livial) therapy. Br J Obstet Gynaecol 101:158–161[Medline]
  17. Nishida M, Kasahara K, Oki A, Satoh T, Arai Y, Kubo T 1996 Establishment of eighteen clones of Ishikawa cells. Hum Cell 9:109–16[Medline]
  18. Satyaswaroop PG, Zaino RJ, Mortel R 1983 Human endometrial adenocarcinoma transplanted into nude mice: growth regulation by estradiol. Science 219:58–60[Abstract/Free Full Text]
  19. Smid-Koopman E, Blok LJ, Kuhne LC, Burger CW, Helmerhorst TJ, Brinkmann AO, Huikeshoven FJ 2003 Distinct functional differences of human progesterone receptors A and B on gene expression and growth regulation in two endometrial carcinoma cell lines. J Soc Gynecol Invest 10:49–57[Medline]
  20. Gielen SCJP, Hanekamp EE, Hanifi-Moghaddam P, Sijbers AM, Van Gool AJ, Burger CW, Blok LJ, Huikeshoven FJ, Growth regulation and transcriptional activities of estrogen and progesterone in human endometrial cancer cells. Gynecol Oncol, in press
  21. Schena M, Shalon D, Davis RW, Brown PO 1995 Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 270:467–470[Abstract/Free Full Text]
  22. Venables WN, Ripley DB 1997 Modern applied statistics with S-PLUS, 2nd Ed. New York: Springer
  23. Smid M, Dorssers LC, Jenster G 2003 Venn mapping: clustering of heterologous microarray data based on the number of co-occurring differentially expressed genes. Bioinformatics 19:2065–2071[Abstract/Free Full Text]
  24. Smid M, Lambert CJ, Dorssers LC 2004 GO-Mapper: functional analysis of gene expression data using the expression level as a score to evaluate Gene Ontology terms. Bioinformatics Advance Access published online on May 6
  25. Hosack DA, Dennis Jr G, Sherman BT, Lane HC, Lempicki RA 2003 Identifying biological themes within lists of genes with EASE. Genome Biol 4:R70
  26. Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S, Bussey KJ, Riss J, Barrett JC, Weinstein JN 2003 GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol 4:R28
  27. Herrero J, Al-Shahrour F, Diaz-Uriarte R, Mateos A, Vaquerizas JM, Santoyo J, Dopazo J 2003 GEPAS: a web-based resource for microarray gene expression data analysis. Nucleic Acids Res 31:3461–3467[Abstract/Free Full Text]
  28. Herrero J, Valencia A, Dopazo J 2001 A hierarchical unsupervised growing neural network for clustering gene expression patterns. Bioinformatics 17:126–362[Abstract/Free Full Text]
  29. Hewitt SC, Korach KS 2000 Progesterone action and responses in the {alpha}ERKO mouse. Steroids 65:551–557[CrossRef][Medline]
  30. Albert R, Jeong H, Barabasi AL 2000 Error and attack tolerance of complex networks. Nature 406:378–382[CrossRef][Medline]
  31. Bray D 2003 Molecular networks: the top-down view. Science 301:1864–1865[Abstract/Free Full Text]
  32. Gokaslan H, Oktem O, Durmusoglu F, Eren F 2003 Apoptosis in the endometrium of postmenopausal women receiving tibolone. Gynecol Obstet Invest 56:124–127[CrossRef][Medline]
  33. LeRoith D, Roberts Jr CT 2003 The insulin-like growth factor system and cancer. Cancer Lett 195:127–137[Medline]
  34. Firth SM, Baxter RC 2002 Cellular actions of the insulin-like growth factor binding proteins. Endocr Rev 23:824–854[Abstract/Free Full Text]
  35. Gui Y, Murphy LJ 2003 Interaction of insulin-like growth factor binding protein-3 with latent transforming growth factor-ß binding protein-1. Mol Cell Biochem 250:189–195[CrossRef][Medline]
  36. Halachmi S, Marden E, Martin G, MacKay H, Abbondanza C, Brown M 1994 Estrogen receptor-associated proteins: possible mediators of hormone-induced transcription. Science 264:1455–1458[Abstract/Free Full Text]
  37. Chen H, Tini M, Evans RM 2001 HATs on and beyond chromatin. Curr Opin Cell Biol 13:218–224[CrossRef][Medline]
  38. Shang Y, Brown M 2002 Molecular determinants for the tissue specificity of SERMs. Science 295:2465–2468[Abstract/Free Full Text]
  39. Majello B, De Luca P, Lania L 1997 Sp3 is a bifunctional transcription regulator with modular independent activation and repression domains. J Biol Chem 272:4021–4026[Abstract/Free Full Text]
  40. Ammanamanchi S, Freeman JW, Brattain MG 2003 Acetylated sp3 is a transcriptional activator. J Biol Chem 278:35775–35780[Abstract/Free Full Text]
  41. Sriraman V, Sharma SC, Richards JS 2003 Transactivation of the progesterone receptor gene in granulosa cells: evidence that Sp1/Sp3 binding sites in the proximal promoter play a key role in luteinizing hormone inducibility. Mol Endocrinol 17:436–449[Abstract/Free Full Text]
  42. Kurosu T, Fukuda T, Miki T, Miura O 2003 Transactivation of the progesterone receptor gene in granulosa cells: evidence that Sp1/Sp3 binding sites in the proximal promoter play a key role in luteinizing hormone inducibility. Oncogene 22:4459–4468[CrossRef][Medline]
  43. Tanno B, Negroni A, Vitali R, Pirozzoli MC, Cesi V, Mancini C, Calabretta B, Raschella G 2002 Expression of insulin-like growth factor-binding protein 5 in neuroblastoma cells is regulated at the transcriptional level by c-Myb and B-Myb via direct and indirect mechanisms. J Biol Chem 277:23172–23180[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Cancer Epidemiol. Biomarkers Prev.Home page
G. C. Rodriguez, B.J. Rimel, W. Watkin, J. M. Turbov, C. Barry, H. Du, G. L. Maxwell, and J.M. Cline
Progestin Treatment Induces Apoptosis and Modulates Transforming Growth Factor-{beta} in the Uterine Endometrium
Cancer Epidemiol. Biomarkers Prev., March 1, 2008; 17(3): 578 - 584.
[Abstract] [Full Text] [PDF]


Home page
Mol Hum ReprodHome page
S. Bai, G. Grossman, L. Yuan, B. A. Lessey, F. S. French, S. L. Young, and E. M. Wilson
Hormone control and expression of androgen receptor coregulator MAGE-11 in human endometrium during the window of receptivity to embryo implantation
Mol. Hum. Reprod., February 1, 2008; 14(2): 107 - 116.
[Abstract] [Full Text] [PDF]


Home page
Reproductive SciencesHome page
C. Barbier, H. J. Kloosterboer, and D. G. Kaufman
Effects of Tibolone Metabolites on Human Endometrial Cell Lines in Co-culture
Reproductive Sciences, January 1, 2008; 15(1): 75 - 82.
[Abstract] [PDF]


Home page
Reproductive SciencesHome page
H. A. M. Verheul, L. J. Blok, C. W. Burger, P. Hanifi-Moghaddam, and H. J. Kloosterboer
Levels of Tibolone and Estradiol and their Nonsulfated and Sulfated Metabolites in Serum, Myometrium, and Vagina of Postmenopausal Women Following Treatment for 21 Days With Tibolone, Estradiol, or Estradiol Plus Medroxyprogestrone Acetate
Reproductive Sciences, February 1, 2007; 14(2): 160 - 168.
[Abstract] [PDF]


Home page
J Mol EndocrinolHome page
P Hanifi-Moghaddam, B Sijmons, M C Ott, W F J van IJcken, D Nowzari, E C M Kuhne, P van der Spek, H J Kloosterboer, C W Burger, and L J Blok
The hormone replacement therapy drug tibolone acts very similar to medroxyprogesterone acetate in an estrogen-and progesterone-responsive endometrial cancer cell line
J. Mol. Endocrinol., December 1, 2006; 37(3): 405 - 413.
[Abstract] [Full Text] [PDF]


Home page
Reproductive SciencesHome page
W. Zhang, J. Mazella, H. J. Kloosterboer, and L. Tseng
Progestagenic Effects of Tibolone are Target Gene--Specific In Human Endometrial Cells
Reproductive Sciences, September 1, 2006; 13(6): 459 - 465.
[Abstract] [PDF]


Home page
Eur Heart JHome page
C. Grohe
'Tib the balance': the search for the optimal hormone replacement therapy
Eur. Heart J., July 2, 2005; 26(14): 1345 - 1346.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hanifi-Moghaddam, P.
Right arrow Articles by Blok, L. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hanifi-Moghaddam, P.
Right arrow Articles by Blok, L. J.
Related Collections
Right arrow Female Endocrinology


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals