- Research article
- Open Access
Tamoxifen-elicited uterotrophy: cross-species and cross-ligand analysis of the gene expression program
BMC Medical Genomics volume 2, Article number: 19 (2009)
Tamoxifen (TAM) is a well characterized breast cancer drug and selective estrogen receptor modulator (SERM) which also has been associated with a small increase in risk for uterine cancers. TAM's partial agonist activation of estrogen receptor has been characterized for specific gene promoters but not at the genomic level in vivo.Furthermore, reducing uncertainties associated with cross-species extrapolations of pharmaco- and toxicogenomic data remains a formidable challenge.
A comparative ligand and species analysis approach was conducted to systematically assess the physiological, morphological and uterine gene expression alterations elicited across time by TAM and ethynylestradiol (EE) in immature ovariectomized Sprague-Dawley rats and C57BL/6 mice. Differential gene expression was evaluated using custom cDNA microarrays, and the data was compared to identify conserved and divergent responses. 902 genes were differentially regulated in all four studies, 398 of which exhibit identical temporal expression patterns.
Comparative analysis of EE and TAM differentially expressed gene lists suggest TAM regulates no unique uterine genes that are conserved in the rat and mouse. This demonstrates that the partial agonist activities of TAM extend to molecular targets in regulating only a subset of EE-responsive genes. Ligand-conserved, species-divergent expression of carbonic anhydrase 2 was observed in the microarray data and confirmed by real time PCR. The identification of comparable temporal phenotypic responses linked to related gene expression profiles demonstrates that systematic comparative genomic assessments can elucidate important conserved and divergent mechanisms in rodent estrogen signalling during uterine proliferation.
The estrogen receptor (ER) is a master transcriptional regulator involved in the proliferation and differentiation of many tissues, most notably the female reproductive tract. It functions as a ligand-dependent transcription factor with two activation functions (AF-1 and AF-2) which interact with cofactors (SRC-1, GRIP1, CBP/p300) to modulate transcription using different mechanisms [1, 2]. The genomic activities of ER are mediated via direct DNA binding at estrogen response elements (EREs) or through indirect tethering mechanisms involving AP-1, Sp1, or Nf-κB . It also has been shown to use non-genomic mechanisms via membrane associated ERs which activate various protein kinase cascades . Furthermore, two distinct ER isoforms exist which have divergent functionality as well as tissue- and cell type-specific expression. For example, whereas mammary tissue expresses both ERα and ERβ  at comparable levels in different cells, the uterus predominantly expresses ERα.
Tamoxifen (TAM) is a well characterized breast cancer drug and prophylactic that is a selective estrogen receptor modulator (SERM). SERMs are structurally diverse compounds that bind the ER and elicit ligand- and tissue-specific effects [6, 7], such as inhibiting the proliferative effects of estrogen in ER-positive breast cancers, while maintaining partial agonistic activity in other tissues . SERM binding causes a unique conformation in the ER ligand binding domain which alters the apposition of ER helix 12 and function of AF-2 relative to 17β-estradiol (E2), while still allowing for the ligand-independent functionality of AF-1 .
TAM-induced ER-mediated gene expression has been characterized on a promoter/gene specific basis [10, 11]. However, the effect of TAM on global uterine gene expression has not been comprehensively examined. Whether or not the modulation of helix 12/AF-2 by TAM results in merely decreased efficacy or specificity for typical ER gene targets or potentiates unique cofactor interactions and thus novel genomic targets has not been fully examined. Therefore, elucidating the genomic targets of TAM is important to the understanding of SERM-ER proliferative activities, especially due to its association with uterine cancer following continuous treatment , which is currently understood to be a function of its partial agonist activity.
The classic rodent uterotrophic assay to assess the estrogenic potential of a chemical examines both the physiological and histomorphological endpoints in the uterus . Uterine wet weight, water imbibition and luminal epithelial cell height induction are typically evaluated after three consecutive daily treatments . This assay can be extended by including temporal gene expression analysis, which can be anchored to these more apical endpoints in order to provide a more comprehensive assessment of the molecular and physiological effects of treatment [15, 16].
Furthermore, this enriched assessment can also be used to evaluate the ability of surrogate models to accurately predict human responses to drugs, industrial chemicals, natural products and environmental contaminants. In the case of tamoxifen metabolism rats have been shown to exhibit a metabolic profile more similar to humans than mice [17, 18]. Thus elucidating the dose-responsive and temporal effects of tamoxifen on gene expression will inform the relative importance of this divergence in rodent metabolic profile. Furthermore, elucidating species-specific differences in both rats and mice for either gene function or regulation is also an important factor for reducing uncertainties associated with cross-species extrapolation of data in risk assessment. Therefore, a cross-species comparative method was employed to examine temporal gene expression in rodents during TAM and ethynylestradiol (EE) induced uterotrophy in order to better understand the early molecular events associated with tamoxifen-related uterine cancers in humans.
Uterine Wet Weight (UWW) and Water Content
Increases in UWW were used to evaluate the responsiveness to EE and TAM. Dose response studies were performed using 0.01, 0.1, 1, 10, 100, 300 μg/kg b.w. EE, and 3, 10, 30, 100, 300, 1000, 3000 μg/kg b.w. TAM. In each case, 100 μg/kg approached the maximum uterotrophic response and this dose was subsequently used for the time course studies (Figure 1). Whereas TAM was equipotent to EE in eliciting uterotrophy in C57BL/6 mice and Sprague-Dawley rats, it was 43% less efficacious in both species eliciting only 4- and 5-fold induction (in the rat and mouse, respectively) compared to the ~9-fold increase induced by EE. In the time course studies, 100 μg/kg of EE or TAM was orally administered once daily for three consecutive days. UWW and water content changes were measured at each time point [see Additional File 1 and 2]. EE induced UWW increases only at 24 and 72 hrs in the mouse, while in the rat the classic water imbibition response occurred between 4 and 12 hrs followed by maximum induction at 72 hrs. TAM induction of water imbibition was delayed approximately 8 hrs in the rat and subdued in both (45% and 65% of EE in rat and mouse, respectively). The increases in uterine water content suggest that early increases in UWW are due at least in part to this water imbibition. As in the wet weight, the changes in water content after TAM treatment temporally lagged behind EE and were notably less efficacious. Therefore, the large difference in wet weight between EE and TAM at 72 hrs is possibly due to early differences in gene expression responsible for water imbibition.
Induction of luminal epithelial height (LEH) is a classic estrogen response [19, 20] in the rodent uterus. No increases in LEH were observed in any group between 2 and 12 hrs (Figure 2). However, EE induced LEH as early as 18 hrs in the mouse and 24 hrs in the rat. While EE and TAM produced comparable levels of LEH induction (~2.6-fold) in the mouse at 72 hrs, TAM treatment in the rat elicited a significantly greater LEH increase (p < 0.05, 4.4-fold for TAM and 2.6-fold for EE). EE induced moderate to marked stromal edema beginning as early as 4 and 8 hrs while TAM induced severe stromal and myometrial edema at 12 hrs. Proliferative indices (number of mitotic bodies) were noted in EE treated uteri at 18 and 24 hrs but not detected in TAM treated samples. Moderate to severe hyperplasia and hypertrophy in the stroma and epithelium were present at 72 hrs in all studies. There was mild apoptotic cell death in the epithelium in TAM studies at 72 hrs while mild to marked apoptosis was noted in EE samples. The species difference in uterine architecture (luminal invaginations seen in the mouse but not rat) previously observed  at 72 hrs was not as pronounced in TAM treated samples.
Comparative Species and Ligand Analysis of Gene Expression
Global gene expression changes with respect to VEH controls were measured and compared across time for MmEE, MmTAM, RnEE, and RnTAM treatments. Pair-wise comparisons between compounds per species and between species per compound were made to investigate conserved responses. A two-tiered, bipartite (P1(t) and fold change) approach was used to screen for conserved differentially expressed genes (Figure 3).
In the mouse, 3,663 and 2,821 genes were differentially regulated by EE and TAM respectively, with 2,631 common differentially expressed genes across all time points (Figure 3B). These 2,631 common, differentially expressed genes were then further examined to assess the similarity of their EE- and TAM-elicited temporal expression profiles by comparing their expression profiles in order to assess similarity between the treatments (Figure 3C). This comparison is comparable to a correlation analysis in assessing the similarity of statistically significant differential gene expression relative to a control elicited by EE and TAM. Each differentially expressed gene was designated as either coactive-similar direction (CAS), coactive-divergent direction (CAD), displaced active-similar direction (DAS), or displaced active-divergent direction (DAD) based on the relationship between the time and direction of differential regulation, and the significance (P1(t)) of the expression profile relative to the VEH control. For example, Igf1 is designated as CAS in the rat because it was up regulated in a similar temporal pattern by both EE and TAM between 12–24 hrs, while Junb was designated DAS because although it was up regulated by both compounds, EE induction was at 4 hrs but TAM treatment temporally shifted the induction to 12 hrs. These designations assist in developing gene lists and provide a preliminary assessment of the similarity in gene expression profile on a qualitative level that takes into account the statistical difference from VEH which typically is not considered in correlation analyses. Using this approach, > 85% of EE and TAM induced genes were designated as CAS (2255 genes). Only a few were categorized as CAD or "divergently regulated" (28 genes, 1% of the overlap).
Similarly in the rat, 2,284 and 2,087 genes were differentially expressed by EE and TAM, respectively, resulting in 1,950 commonly regulated genes (Figure 3B). In both species, roughly 93% of TAM-regulated genes overlapped with EE. Of the 1,950 genes differentially expressed by both EE and TAM in the rat, 86% (1,686 genes) were designated as CAS while only 23 were CAD. Comparison of the 28 mouse CAD genes to the 23 rat CAD genes found no overlapping orthologs (Figure 3D), suggesting that there are no conserved differentially regulated genes for which EE and TAM elicit divergent uterine responses in the rat and mouse.
Comparison of Orthologous Gene Expression
A comparable approach was used to examine the cross-species differential gene expression effects of EE and TAM on orthologous rat and mouse genes. 3,417 unique orthologous genes were represented on the rat and mouse array platforms as determined by NCBI's HomoloGene database (Figure 4A). In response to EE treatment, 2,095 and 2,181 of the 3,417 orthologous genes were differentially expressed in the mouse and rat, respectively, with 1,634 orthologs differentially expressed in both species (~75%) (Figure 4B). These 1,634 differentially expressed orthologs were analyzed for coactivity of which 1,116 genes (68%) were designated as CAS and 206 CAD. Similarly, 1,252 mouse and 1,441 rat genes were differentially expressed following TAM treatment with 891 orthologs differentially expressed in both species. 705 (79%) of these were designated as CAS genes while 63 genes were designated CAD. A comparison of the EE and TAM CAD orthologs identified 35 genes that are divergently regulated between the mouse and rat in response to both EE and TAM. Of these 35 genes, 9 were represented by only a single feature on the microarray for both species, 12 were represented by 2 or more features in at least one species that had poor internal correlation, and 14 were represented by 2 or more features in at least one species that were internally consistent. Thus, a strength of evidence approach was taken for pursuing genes exhibiting putative divergent regulation with quantitative real time PCR.
QRT-PCR Verification of Microarray Data
The expression profiles of 25 genes (including known estrogen responsive genes: Fos, C3, Calb3, Igf1, Sult1a1, Aqp5, and Vegf) accounting for the spectrum of temporal patterns and functional categories were verified by QRT-PCR in either species or compound. In general there was a good correlation between the QRT-PCR and microarray results (data not shown).
Fourteen orthologs that exhibited putative divergent regulation in the mouse compared to the rat following EE and TAM treatment were further investigated. However, QRT-PCR indicated that these genes exhibited either a CAS relationship between species or were inconclusive, except for carbonic anhydrase 2 (designated Car2 in the mouse and Ca2 in the rat and human, hereafter referred to as Ca2 to represent all orthologs). QRT-PCR analysis of Ca2 confirmed the microarray data (correlation coefficient r = 0.95, 0.98, 0.82, and 0.78 for MmEE, MmTAM, RnEE and RnTAM, respectively; Figure 5). These results were confirmed by two different sets of primers, with one set querying species specific regions of the mouse and rat mRNA, and the second set designed to amplify the same region in both species using the same primers.
Four-Way Venn Analysis
Following pair-wise comparisons by chemical and species, an integrated comparison of all four data sets was performed where differentially expressed orthologous genes (NCBI, HomoloGene) for each data set were entered into the 4-way Venn Diagram Generator. The number and identity of unique genes populating the overlaps of each data set was determined for common orthologs (Figure 6). The comparison of all four data sets resulted in 902 genes being commonly regulated at any time point (Venn subset I). Of the 902 overlapping genes, 398 exhibited a similar temporal profile (CAS) across all four data sets, suggesting comparable modes of regulation which we refer to as "orthologous expression". These genes represented functional categories associated with cellular proliferation and differentiation; hallmarks of estrogen induced uterotrophy. Selected conserved examples exhibiting orthologous expression in terms of either fold-change alone (Fos, C3, Igf1, Ca3, Sepp1) or large increase in copy number, regardless of fold change (Cd24a, Slc25a5, Krt13, Dcn, Itm2b) are provided in Table 1.
Of interest are the 128 and 12 genes populating subsets B and K, respectively. Subset B represents putative species-conserved EE-specific responses, of which, less than 60% exhibited similar direction of regulation (CAS or DAS designation) by TAM. Most of these genes were regulated by TAM but did not meet the fold change or P1(t) selection criteria, suggesting a potency issue rather than unique regulation. For example S100a8, a calcium binding protein, was induced 1.5-fold (Mm) and 3.5-fold (Rn) by EE while TAM only induced a 1.3-fold change in the mouse that did not meet the statistical cutoff while in the rat the statistical cutoff was achieved but induction was only 1.2-fold. Likewise, further analysis of the 12 species-conserved TAM-specific responses revealed that only Myh6 (myosin, heavy chain 6), a component of muscle fibre that is likely expressed in the myometrium, exhibited an overlap in temporal expression between species. However, comparable expression is only observed at 72 hrs where secondary and tertiary effects are likely to overlap. In addition, there is insufficient data to exclude the possibility that Myh6 is regulated by TAM via an ER-independent mechanism.
Collectively, these results suggest there are no species-conserved, TAM-specific or -divergent genes which are regulated in a separate manner when compared to EE in the rodent uterus. Approximately 400 genes exhibited highly conserved and correlated expression across all four data sets, defining a robust gene set that is predictable and integral to ER-mediated uterotrophy in the rodent. However, species-specific gene expression profiles are difficult to assess and confirm given the complex orthologous relationships between the species specific array probe sets. Many of the putative differences suggested may be attributed to differing probe sensitivities or representations within respective genes. The lone identified exception is the species-divergent regulation of Ca2 that is conserved in response to both ER ligands and confirmed by homologous and species-specific QRT-PCR primer pairs.
The uterotrophic response has been used as an acute assay to assess the estrogenicity of a compound using both physiological and cytoarchitectural endpoints. Our extended design incorporates early monitoring of global gene expression (2–24 hrs) to capture the ER-mediated effects as well as subsequent secondary and tertiary responses. Early gene expression was also examined for conserved responses across ligand and species. The highly parallel nature of the EE and TAM study designs and subsequent analyses facilitated more robust comparisons, thus increasing confidence in the identification of conserved uterine responses between species and ligands.
EE and TAM doses were selected based on UWW dose response studies which identified 100 μg/kg as the most efficacious dose for each compound. Although differing pharmacokinetics may account for discrepancies between rodents and humans, 100 μg/kg TAM is below the pharmacological dose of Nolvadex® (tamoxifen citrate, 20–40 mg/day) prescribed to women (300–400 μg/kg), assuming an average weight of 70 kg. Moreover, the ovariectomized model provides increased sensitivity to estrogen action in the immature rodent and thus allows for a more comprehensive assessment of the partial-agonistic activity of TAM in the uterus.
The physiological effects were comparable between species and ligand with EE eliciting the classic uterotrophic responses while TAM showed decreased efficacy with only partial agonist activity, as previously reported . However, the more pronounced induction of luminal epithelial cell height in the rat compared to the mouse was unexpected. Previous rat studies have not observed increased efficacy of TAM in inducing luminal epithelial cell height relative to estrogen [23, 24], which may due to differences in animal age, ovariectomy status or other study design issues. However, a comparable increase in TAM LEH induction in rats relative to estrogen was previously reported as not statistically different . Furthermore, the uterotrophic effects of EE and TAM in ovariectomized Cynomolgus Macaque monkeys  showed a noticeable increase in TAM induced LEH relative to EE that was not statistically significant. The similarity in LEH response to higher order mammals/primates suggests that the rat may be a better model of the human response, but further study is required.
The partial agonist effects of TAM with regard to its ER-complex conformation and differential interactions at specific canonical response elements (including EREs and AP-1 sites) as well as gene specific promoters (Vitellogenin, Complement component 3, prolactin) have been well studied. However, the global gene targets of TAM-ER complex have not been comprehensively elucidated in the rodent uterus. Although direct primary ligand-ER effects can not be deduced from expression profiling alone, our comparative approach provides evidence of several conserved responses elicited by two structurally diverse ligands. While EE is absorbed and attains peak plasma concentration before TAM does [27, 28], the difference is only a matter of hours and the more likely limiting step is the hydroxylation of TAM to 4-OH-tamoxifen or endoxifen, which could account for the 10–12 h delay observed in gene expression response between EE and TAM. However, as primary gene responses give way to secondary responses, gene expression cascades are sequentially propagated across time. If the two ligand-receptor complexes elicit different behaviours at the primary response genes, it is likely that subsequent differences in secondary and tertiary responses would also be propagated over time. However, this is not the case as EE and TAM elicited comparable expression profiles that are maintained throughout the time course in both species. This suggests that TAM elicits parallel uterine gene expression behavior when compared EE, despite the temporal shift of expression due to pharmacokinetic differences .
An inclusive, multi-step approach to identifying differentially expressed genes was used to ensure conservation of gene expression between ligand or species could not be attributed to strict cut offs applied during screening. Comparison of gene expression profiles revealed highly similar sets of genes between species and compounds (Figure 3 and 4). When profiles were aligned and compared on a quantitative scale, TAM typically elicited a lower fold change when compared to EE, consistent with its partial agonist activity. However, four notable species-conserved orthologs (Hsp90aa1, Slc9a3r1, Acsbg1, Col5a1) were identified where TAM induction exceeded that of EE. For example, Slc9a3r1, a cytoskeletal-membrane protein binding-protein that functions to maintain epithelial cell structure and polarity , was induced 4.1-fold (Rn) and 3.6-fold (Mm) by TAM but only 2.7-fold (Rn) and 3-fold (Mm) by EE, respectively. Few genes exhibited unique differential gene expression suggesting that TAM elicits uterotrophic affects through the same target genes affected by EE (Figure 6). This is significant, as little information is known about how conformational changes in the SERM-ER complex affect the number or types of uterine genes modulated after treatment.
A large number of regulated genes were differentially expressed across all four data sets, of which, 398 exhibited comparable patterns of expression [see Additional Files 3, 4, 5, and 6], and likely represent conserved ER-mediated uterine responses. They were associated with several functional categories including cellular energetics, chaperone and protein folding, cell structure, cell signaling, transcription and RNA processing, and translation and protein turnover (Table 1). Although most genes were up regulated, those associated with the cell cycle and mitogenic activity were down regulated, consistent with an overall proliferative response [31, 32]. Several of these conserved early (e.g., Fos and Vegf), mid (e.g., Ccnd1) and late (e.g., C3) responders responses are known ER targets regulated via estrogen response elements (EREs), AP-1, and Sp1 sites. The temporal differences in expression (between 2 and 72 hrs after treatment) of genes known to be ER regulated affirms the bi- or multi-phasic activities of activated ER during uterotrophy . Although the TAM profiles exhibited a temporal shift for early regulated genes (Fos and Vegf), genes differentially expressed at 18, 24 and 72 hrs (Ccnd1, Tk1) exhibited little temporal shift. This suggests a pharmacokinetic difference which is not applied equally to all responses or is masked by the sustained nature of mid to late gene responses.
TAM perturbs calcium homeostasis causing intracellular influx of Ca2+ ions which is related to cytotoxic, thrombotic and systemic effects [34, 35]. This is believed to be ER independent, and is suggestive of alternative mechanisms activated only by TAM. However, no differential gene expression was observed regarding genes involved in calcium homeostasis, including Ca2+ transporters Atp2a1, -2a2, -2b1, and -2b2, calcium transporters that are targets of the calcineurin/Crz1 signaling pathway which regulates intracellular Ca2+ homeostasis. Moreover, the increased LEH responsiveness in rat is not explained by our current data, although this may be due to differences in genome coverage represented on our mouse and rat cDNA microarrays. It is also likely that there are significant species-specific post-transcriptional differences that are not assessed using gene expression approaches.
In summary, our parallel study design afforded a robust comparative analysis of EE and TAM elicited responses in the rat and mouse uterus. TAM induced UWW and LEH in a manner consistent with its partial-agonist activity at a dose equipotent to EE. This suggests that EE and TAM elicit comparable uterine gene expression profiles despite conformation differences in the liganded ER complexes, and that the transcriptional activity via AF-1, which is activated by both EE and TAM, is sufficient to elicit uterine differential gene expression. This is consistent with mice studies with targeted disruption of the ER DNA binding domain , which did not exhibit acute uterotrophy but still elicited epithelial proliferation and increased LEH. It also supports the sufficiency of the AF-1 domain alone (via TAM binding) to mediate uterotrophy. Although differences in LEH induction were observed between EE and TAM, collectively our data suggests that TAM does not elicit a unique uterine gene expression profile when compared to EE. However, more comprehensive studies are warranted that would examine the differential expression of all orthologous rat and mouse genes.
Experimental designs and methods for the rat tamoxifen data parallel previously published rat and mouse EE and mouse TAM studies [21, 37]. Briefly, female Sprague-Dawley rats and C57BL/6 mice, ovariectomized on PND 20 and all within 10% of the average body weight were obtained from Charles River Laboratories (Raleigh, NC) on day 25. Animals were housed in polycarbonate cages containing cellulose fibre chip bedding (Aspen Chip Laboratory Bedding, Northeastern Products, Warrensberg, NY) and maintained at 40–60% humidity and 23°C in a room with a 12 hrs dark/light cycle (7 am–7 pm). Animals were allowed free access to de-ionized water and Harlan Teklad 22/5 Rodent Diet 8640 (Madison, WI), and acclimatized for 4 days prior to dosing.
In dose response studies, animals received three consecutive, daily oral treatments of EE (0.01 to 300 μg/kg) or TAM (3 to 3000 μg/kg b.w.). Animals were sacrificed 72 hrs after the initial dose. In time course studies animals were treated once or once daily for three consecutive days via oral gavage with 100 μg/kg b.w. EE or TAM in 0.1 mL of sesame oil vehicle (Sigma Chemical, St Louis, MO), (EE time course treatment in mouse, MmEE; EE time course treatment in rat, RnEE; EE in mouse, MmEE; TAM in mouse, MmTAM;). This oral dose was empirically derived and chosen because it elicited a maximal uterotrophic response in both species while showing no acute toxic effects. Animals receiving a single dose were sacrificed 2, 4, 8, 12, 18, or 24 hrs after treatment. Animals receiving three consecutive, daily doses were sacrificed 72 hrs post initial dose. An equal number of time-matched vehicle control (VEH) animals were treated in the same manner, (n = 8 for MmEE, n = 5 for all others). Doses of EE and TAM were calculated based on average weights of each treatment and VEH group prior to dosing. All procedures were performed with the approval of the Michigan State University All-University Committee on Animal Use and Care.
Animals were sacrificed by cervical dislocation and animal body weights were recorded. The uterine body was dissected at the border of the cervix and whole uteri were harvested and stripped of extraneous connective tissue and fat. Whole uterine weights were recorded before (wet) and after (blotted) being blotted under pressure with absorbent tissues and were subsequently snap-frozen in liquid nitrogen and stored at -80°C. Weight due to water was calculated as the difference between the wet and blotted weights. A small (~5 mm) section of the right, distal uterine horn was placed in 10% neutral buffered formalin (NBF; VWR, West Chester, PA) and stored at RT for at least 24 hrs prior to further processing. Statistical analysis of wet weight and water content were conducted using a two-way ANOVA with a Tukey's Honestly Significant Difference post hoc test, p < 0.05 (SAS 9.1, SAS Institute Inc. Cary, NC).
Total RNA was isolated from whole uteri (~20 mg/rat, ~3 mg/mouse) using Trizol® Reagent (Invitrogen, Carlsbad, CA). Uteri were removed from -80°C storage and immediately homogenized in 1 mL Trizol® Reagent using a Mixer Mill 300 tissue homogenizer (Retsch, Germany). Total RNA was isolated according to manufacturer's protocol and resuspended in The RNA Storage Solution (Ambion, Austin, TX). Concentration was calculated by spectrophotometric methods (A260) and purity assessed by the A260:A280 ratio and by visual inspection of 1 μg on a denaturing gel.
NBF-fixed uterine sections were routinely processed and embedded in paraffin according to standard histological techniques. Five μm cross-sections were mounted on glass slides and stained with hematoxylin and eosin. All embedding, mounting and staining of tissues were performed at the Histology/Immunohistochemistry Laboratory, located in the Department of Physiology, of Michigan State University.
Histopathological and Morphometric Assessment
Histological slides were scored according to standardized National Toxicology Program (NTP) pathology codes. Morphometric analyses were performed on cross sections for each animal using image analysis software (Scion Image, Scion Corporation, Frederick, MD) and standard morphometric techniques. Briefly, the contour length of basal lamina underlying the luminal epithelium (LE) and corresponding areas of LE was quantified for multiple, representative sectors of each section. Statistical analyses on all morphometry data were performed using a two-way ANOVA with a Tukey's Honestly Significant Difference post hoc test, n = 5, p < 0.05 (SAS 9.1, SAS Institute Inc.).
Rat cDNA arrays were produced in-house using a LION Bioscience's Rat cDNA library (LION Bioscience, Heidelberg Germany) consisting of 8,567 clones representing 5,684 unique genes (Unigene Build #48). Clones were selected based on their level of annotation as well as sequence similarity to well annotated human and mouse genes. Detailed protocols for microarray construction, labeling of the cDNA probe, sample hybridization and slide washing can be found at http://dbzach.fst.msu.edu/interfaces/microarray.html. Briefly, PCR amplified DNA was robotically arrayed onto epoxy coated glass slides (SCHOTT Louisville, KY), using an Omnigrid arrayer (GeneMachines, San Carlos, CA) equipped with 32 (8 × 4) or 48 (12 × 4) Chipmaker 2 pins (Telechem, Atlanta, GA) for the rat and mouse arrays, respectively, at the Genomics Technology Support Facility at Michigan State University http://www.genomics.msu.edu.
Array Experimental Design
Temporal changes in gene expression were assessed using an independent reference design in which samples from EE and TAM treated animals were co-hybridized with VEH. Comparisons were performed on 3 biological replicates × 2 independent labelings of each sample (incorporating a dye swap) for each time point. Total RNA (15 μg) was reverse transcribed in the presence of Cy3- or Cy5-labeled dUTP (Amersham, Piscataway, NJ) to create fluor-labeled cDNA, which was purified using QIAquick PCR purification kit (Qiagen, Valencia, CA). In contrast, note that mouse array experiments used a 3DNA Array 900 Expression Array Detection Kit (Genisphere, Hatsfield, PA) using 1 μg of total RNA, according to manufacturer's specifications, for probe labeling. Cy3- and Cy5-labeled samples were mixed, vacuum concentrated (~1–2 μl) and resuspended in 32 μl of hybridization buffer (40% formamide, 4× SSC, 1% SDS) with 20 μg polydA and 20 μg of mouse COT-1 DNA (Invitrogen) as a competitor. This probe mixture was heated at 95°C for 2 min and was then hybridized to the array under a 22 × 40 mm lifterslip (Erie Scientific Company, Portsmouth, NH) in a light protected and humidified hybridization chamber (Corning Inc., Corning, NY). Samples were hybridized for 18–24 h at 42°C in a water bath. Slides were then washed, dried by centrifugation and scanned at 635 (Cy5) and 532 nm (Cy3) on a GenePix 4000B microarray scanner (Molecular Devices, Union City, CA). Images were analyzed for feature and background intensities using GenePix Pro 6.0 (Molecular Devices).
Array Data Normalization and Statistical Analysis
Data sets for rat and mouse EE and mouse TAM have been previously published [21, 37] and have been integrated into the current comparative analysis with rat TAM data. As previously described with these data sets, rat TAM microarray data were first examined using a quality assurance protocol prior to further analysis to ensure consistent, high quality data throughout the dose-response and time course studies prior to normalization and further analysis . Data were normalized using a semi-parametric approach . Model-based t-values were calculated from normalized data, comparing treated from VEH responses per time-point. Empirical Bayes analysis was used to calculate posterior probabilities of activity (P1(t)-value) on a per gene and time-point basis using the model-based t-value . Genes were filtered for differential expression based on the P1(t)-value which indicates increasing activity as the value approaches 1.0. A two-tiered set of criteria including a statistical P1(t) > 0.999 and |fold-change| > 1.5 was used as an initial selection filter of the expression data and defined initial differentially expressed gene lists in both the rat and mouse. All data was deposited into dbZach http://dbzach.fst.msu.edu, a MIAME compliant relational database that ensures proper data management and facilitates data analysis. Complete data sets with annotation and P1(t) values are available in Additional Files 3, 4, 5, and 6.
Expression Data Annotation and Coactivity Analysis
Features refer to unique cDNA clones spotted on the array and are assigned a GenBank accession number, gene name, symbol, and Entrez Gene ID where annotation is available. For the sake of brevity and consistency, genes are referenced by their official gene symbol as defined by NCBI. Rat and mouse orthologous gene pairs were derived from the publicly available Ensembl and HomoloGene databases. Comparing gene expression between species or ligand treatment involved examining the time, direction, and statistical significance of the change in expression on a gene by gene basis. Similarities and differences in gene expression patterns were designated as coactive-similar direction (CAS), coactive-divergent direction (CAD), displaced active-similar direction (DAS), and displaced active-divergent direction (DAD). Comparative analysis was conducted using a multivariate correlation-based visualization application developed in-house. Additional analysis was performed using the 4-way Venn Diagram Generator (4VDG, http://www.pangloss.com/seidel/Protocols/venn4.cgi, which involved applying an initial relaxed filtering criteria (|fold change| > 1.3 and P1(t) > 0.99) to each data set, then entering respective Entrez Gene IDs (mouse IDs were used for rat where HomoloGene indicated orthology) into the 4VDG. The output categories were then filtered for only those genes which met more stringent criteria (|fold change| > 1.5 and P1(t) > 0.999) at any one time point. Gene lists for each 4-way Venn category are available [see Additional File 7].
QRT-PCR (Real Time) Analysis
For each sample, 1.0 μg of total RNA was reverse transcribed by SuperScript II using an anchored oligo-dT primer as described by the manufacturer (Invitrogen). The resultant cDNA (1.0 μl) was used as the template in a 30 μl PCR reaction containing 0.1 μM each of forward and reverse gene-specific primers designed using Primer3 (15), 3 mM MgCl2, 1.0 mM dNTPs, 0.025 IU AmpliTaq Gold and 1× SYBR Green PCR buffer (Applied Biosystems, Foster City, CA). Gene names, accession numbers, forward and reverse primer sequences and amplicon sizes are available [see Additional File 8]. PCR amplification was conducted in MicroAmp Optical 96-well reaction plates (Applied Biosystems) on an Applied Biosystems PRISM 7500 Sequence Detection System using the following conditions: initial denaturation and enzyme activation for 10 min at 95°C, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. A dissociation protocol was performed to assess the specificity of the primers and the uniformity of the PCR generated products. Each plate contained duplicate standards of purified PCR products of known template concentration covering six orders of magnitude to interpolate relative template concentrations of the samples from the standard curves of log copy number versus threshold cycle (Ct). No template controls (NTC) were also included on each plate. Samples with a Ct value within 2 SD of the mean Ct values for the NTCs were considered below the limits of detection. The copy number of each unknown sample for each gene was standardized to the geometric mean of house-keeping gene, Rpl7 to control for differences in RNA loading, quality and cDNA synthesis. Statistical significance of induced or repressed genes was determined using two-way ANOVA followed by t-test for VEH treatment comparisons (SAS 9.1, SAS Institute Inc.). For graphing purposes, the relative expression levels were scaled such that the expression level of the time-matched control group was equal to one.
Hall JM, Couse JF, Korach KS: The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem. 2001, 276: 36869-36872. 10.1074/jbc.R100029200.
O'Lone R, Frith MC, Karlsson EK, Hansen U: Genomic targets of nuclear estrogen receptors. Mol Endocrinol. 2004, 18: 1859-1875. 10.1210/me.2003-0044.
Bjornstrom L, Sjoberg M: Mechanisms of estrogen receptor signaling: convergence of genomic and nongenomic actions on target genes. Mol Endocrinol. 2005, 19: 833-842. 10.1210/me.2004-0486.
Losel RM, Falkenstein E, Feuring M, Schultz A, Tillmann HC, Rossol-Haseroth K, Wehling M: Nongenomic steroid action: controversies, questions, and answers. Physiol Rev. 2003, 83: 965-1016.
Saji S, Jensen EV, Nilsson S, Rylander T, Warner M, Gustafsson JA: Estrogen receptors alpha and beta in the rodent mammary gland. Proc Natl Acad Sci USA. 2000, 97: 337-342. 10.1073/pnas.97.1.337.
Gajdos C, Jordan VC: Selective estrogen receptor modulators as a new therapeutic drug group: concept to reality in a decade. Clin Breast Cancer. 2002, 2: 272-281. 10.3816/CBC.2002.n.002.
Jordan VC: The past, present, and future of selective estrogen receptor modulation. Ann N Y Acad Sci. 2001, 949: 72-79.
Valverde P: Pharmacotherapies to manage bone loss-associated diseases: a quest for the perfect benefit-to-risk ratio. Curr Med Chem. 2008, 15: 284-304. 10.2174/092986708783497274.
Pearce ST, Liu H, Jordan VC: Modulation of estrogen receptor alpha function and stability by tamoxifen and a critical amino acid (Asp-538) in helix 12. J Biol Chem. 2003, 278: 7630-7638. 10.1074/jbc.M211129200.
Webb P, Nguyen P, Valentine C, Lopez GN, Kwok GR, McInerney E, Katzenellenbogen BS, Enmark E, Gustafsson JA, Nilsson S, Kushner PJ: The estrogen receptor enhances AP-1 activity by two distinct mechanisms with different requirements for receptor transactivation functions. Mol Endocrinol. 1999, 13: 1672-1685. 10.1210/me.13.10.1672.
Schultz JR, Petz LN, Nardulli AM: Cell- and ligand-specific regulation of promoters containing activator protein-1 and Sp1 sites by estrogen receptors alpha and beta. J Biol Chem. 2005, 280: 347-354.
Lavie O, Barnett-Griness O, Narod SA, Rennert G: The risk of developing uterine sarcoma after tamoxifen use. Int J Gynecol Cancer. 2008, 18: 352-356. 10.1111/j.1525-1438.2007.01025.x.
Odum J, Lefevre PA, Tittensor S, Paton D, Routledge EJ, Beresford NA, Sumpter JP, Ashby J: The rodent uterotrophic assay: critical protocol features, studies with nonyl phenols, and comparison with a yeast estrogenicity assay. Regul Toxicol Pharmacol. 1997, 25: 176-188. 10.1006/rtph.1997.1100.
Kanno J, Onyon L, Peddada S, Ashby J, Jacob E, Owens W: The OECD program to validate the rat uterotrophic bioassay. Phase 2: coded single-dose studies. Environ Health Perspect. 2003, 111: 1550-1558.
Diel P, Geis RB, Caldarelli A, Schmidt S, Leschowsky UL, Voss A, Vollmer G: The differential ability of the phytoestrogen genistein and of estradiol to induce uterine weight and proliferation in the rat is associated with a substance specific modulation of uterine gene expression. Mol Cell Endocrinol. 2004, 221: 21-32. 10.1016/j.mce.2004.04.006.
Diel P, Schmidt S, Vollmer G: In vivo test systems for the quantitative and qualitative analysis of the biological activity of phytoestrogens. J Chromatogr B Analyt Technol Biomed Life Sci. 2002, 777: 191-202. 10.1016/S1570-0232(02)00494-4.
DeGregorio MW, Coronado E, Osborne CK: Tumor and serum tamoxifen concentrations in the athymic nude mouse. Cancer Chemother Pharmacol. 1989, 23: 68-70. 10.1007/BF00647244.
Robinson SP, Langan-Fahey SM, Jordan VC: Implications of tamoxifen metabolism in the athymic mouse for the study of antitumor effects upon human breast cancer xenografts. Eur J Cancer Clin Oncol. 1989, 25: 1769-1776. 10.1016/0277-5379(89)90347-7.
Kang JS, Lee BJ, Ahn B, Kim DJ, Nam SY, Yun YW, Nam KT, Choi M, Kim HS, Jang DD, et al: Expression of estrogen receptor alpha and beta in the uterus and vagina of immature rats treated with 17-ethinyl estradiol. J Vet Med Sci. 2003, 65: 1293-1297. 10.1292/jvms.65.1293.
Padilla-Banks E, Jefferson WN, Newbold RR: The immature mouse is a suitable model for detection of estrogenicity in the uterotropic bioassay. Environ Health Perspect. 2001, 109: 821-826. 10.2307/3454825.
Kwekel JC, Burgoon LD, Burt JW, Harkema JR, Zacharewski TR: A cross-species analysis of the rodent uterotrophic program: elucidation of conserved responses and targets of estrogen signaling. Physiol Genomics. 2005, 23: 327-342. 10.1152/physiolgenomics.00175.2005.
MacNab MW, Tallarida RJ, Joseph R: An evaluation of tamoxifen as a partial agonist by classical receptor theory – an explanation of the dual action of tamoxifen. Eur J Pharmacol. 1984, 103: 321-326. 10.1016/0014-2999(84)90493-X.
Carthew P, Edwards RE, Nolan BM, Tucker MJ, Smith LL: Compartmentalized uterotrophic effects of tamoxifen, toremifene, and estradiol in the ovariectomized Wistar (Han) rat. Toxicol Sci. 1999, 48: 197-205. 10.1093/toxsci/48.2.197.
Stygar D, Muravitskaya N, Eriksson B, Eriksson H, Sahlin L: Effects of SERM (selective estrogen receptor modulator) treatment on growth and proliferation in the rat uterus. Reprod Biol Endocrinol. 2003, 1: 40-10.1186/1477-7827-1-40.
Nephew KP, Osborne E, Lubet RA, Grubbs CJ, Khan SA: Effects of oral administration of tamoxifen, toremifene, dehydroepiandrosterone, and vorozole on uterine histomorphology in the rat. Proc Soc Exp Biol Med. 2000, 223: 288-294. 10.1046/j.1525-1373.2000.22341.x.
Sikoski P, Register TC, Lees CJ, Lundeen S, Hutchison J, Brown KH, Cline JM: Effects of two novel selective estrogen receptor modulators, raloxifene, tamoxifen, and ethinyl estradiol on the uterus, vagina and breast in ovariectomized cynomolgus monkeys (Macaca fascicularis). Am J Obstet Gynecol. 2007, 196 (75): e71-77.
Shin SC, Choi JS, Li X: Enhanced bioavailability of tamoxifen after oral administration of tamoxifen with quercetin in rats. Int J Pharm. 2006, 313: 144-149. 10.1016/j.ijpharm.2006.01.028.
Twaddle NC, Churchwell MI, Newbold RR, Delclos KB, Doerge DR: Determination using liquid-chromatography-electrospray tandem mass spectroscopy of ethinylestradiol serum pharmacokinetics in adult Sprague-Dawley rats. J Chromatogr B Analyt Technol Biomed Life Sci. 2003, 793: 309-315. 10.1016/S1570-0232(03)00331-3.
Robinson SP, Langan-Fahey SM, Johnson DA, Jordan VC: Metabolites, pharmacodynamics, and pharmacokinetics of tamoxifen in rats and mice compared to the breast cancer patient. Drug Metab Dispos. 1991, 19: 36-43.
Morales FC, Takahashi Y, Kreimann EL, Georgescu MM: Ezrin-radixin-moesin (ERM)-binding phosphoprotein 50 organizes ERM proteins at the apical membrane of polarized epithelia. Proc Natl Acad Sci USA. 2004, 101: 17705-17710. 10.1073/pnas.0407974101.
Boland CR, Goel A: The silence of the genes: matching mismatch repair defects with tumors. Cancer. 2003, 98: 2091-2094. 10.1002/cncr.11769.
Kishioka Y, Thomas M, Wakamatsu J, Hattori A, Sharma M, Kambadur R, Nishimura T: Decorin enhances the proliferation and differentiation of myogenic cells through suppressing myostatin activity. J Cell Physiol. 2008, 215: 856-867. 10.1002/jcp.21371.
Hewitt SC, Deroo BJ, Hansen K, Collins J, Grissom S, Afshari CA, Korach KS: Estrogen receptor-dependent genomic responses in the uterus mirror the biphasic physiological response to estrogen. Mol Endocrinol. 2003, 17: 2070-2083. 10.1210/me.2003-0146.
Bollig A, Xu L, Thakur A, Wu J, Kuo TH, Liao JD: Regulation of intracellular calcium release and PP1alpha in a mechanism for 4-hydroxytamoxifen-induced cytotoxicity. Mol Cell Biochem. 2007, 305: 45-54. 10.1007/s11010-007-9526-2.
Dobrydneva Y, Weatherman RV, Trebley JP, Morrell MM, Fitzgerald MC, Fichandler CE, Chatterjie N, Blackmore PF: Tamoxifen stimulates calcium entry into human platelets. J Cardiovasc Pharmacol. 2007, 50: 380-390. 10.1097/FJC.0b013e31811ec748.
O'Brien JE, Peterson TJ, Tong MH, Lee EJ, Pfaff LE, Hewitt SC, Korach KS, Weiss J, Jameson JL: Estrogen-induced proliferation of uterine epithelial cells is independent of estrogen receptor alpha binding to classical estrogen response elements. J Biol Chem. 2006, 281: 26683-26692. 10.1074/jbc.M601522200.
Fong CJ, Burgoon LD, Williams KJ, Forgacs AL, Zacharewski TR: Comparative temporal and dose-dependent morphological and transcriptional uterine effects elicited by tamoxifen and ethynylestradiol in immature, ovariectomized mice. BMC Genomics. 2007, 8: 151-10.1186/1471-2164-8-151.
Burgoon LD, Eckel-Passow JE, Gennings C, Boverhof DR, Burt JW, Fong CJ, Zacharewski TR: Protocols for the assurance of microarray data quality and process control. Nucleic Acids Res. 2005, 33: e172-10.1093/nar/gni167.
Eckel JE, Gennings C, Therneau TM, Burgoon LD, Boverhof DR, Zacharewski TR: Normalization of two-channel microarray experiments: a semiparametric approach. Bioinformatics. 2005, 21: 1078-1083. 10.1093/bioinformatics/bti105.
Eckel JE, Gennings C, Chinchilli VM, Burgoon LD, Zacharewski TR: Empirical bayes gene screening tool for time-course or dose-response microarray data. J Biopharm Stat. 2004, 14: 647-670. 10.1081/BIP-200025656.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/19/prepub
The authors would like to thank Ed Dere and Ania Kopeć for helpful discussions and critical review of this manuscript. JCK is supported by a fellowship from the Center for Integrated Toxicology through National Institute of Environmental Health Sciences (NIEHS) Training Grant T32 ES-07255-16. TRZ is partially supported by the Michigan Agricultural Experimental Station.
JK supervised and participated in all in-life studies, carried out gene expression experiments and analysis and drafted the manuscript. AF coordinated and participated in the in-life studies and histological processing. LB carried out statistical analysis and provided bioinformatics support in data management. KW carried out histopathological assessments and analysis. TZ conceived the study, participated in design of the study and participated in final draft of manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Uterine Relative Wet Weights. data presented represents averaged weights of unblotted whole rodent uteri (XLS 18 KB)
Additional file 2: Uterine Relative Water Content. data presented represents average difference between wet and blotted uterine weight (XLS 20 KB)
Additional file 3: Mouse EE Gene Expression Data. data presented represents averaged fold changes and statistical values per time point for each feature on our array platform for mouse ethynylestradiol. (XLS 5 MB)
Additional file 4: Mouse TAM Gene Expression Data. data presented represents averaged fold changes and statistical values per time point for each feature on our array platform for mouse tamoxifen. (XLS 5 MB)
Additional file 5: Rat EE Gene Expression Data. data presented represents averaged fold changes and statistical values per time point for each feature on our array platform for rat ethynylestradiol. (XLS 4 MB)
Additional file 6: Rat TAM Gene Expression Data. data presented represents averaged fold changes and statistical values per time point for each feature on our array platform for rat tamoxifen (XLS 3 MB)
Additional file 7: 4-Way Venn Analysis Gene Sets. data presented represents gene lists for each sub group in the 4-way Venn analysis. (XLS 140 KB)
Additional file 8: QRT-PCR Primer Sequences. data presented represents the DNA sequences of gene-specific primers used for quantitative real-time PCR analysis. (XLS 15 KB)
Authors’ original submitted files for images
About this article
Cite this article
Kwekel, J.C., Forgacs, A.L., Burgoon, L.D. et al. Tamoxifen-elicited uterotrophy: cross-species and cross-ligand analysis of the gene expression program. BMC Med Genomics 2, 19 (2009). https://doi.org/10.1186/1755-8794-2-19
- Estrogen Receptor
- Express Gene List
- Uterine Response
- Rodent Uterus