|
|
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Departments of Biochemistry & Molecular Biology (D.R.B., L.D.B., T.R.Z.) and Pathobiology & Diagnostic Investigation (K.J.W.), Center for Integrative Toxicology (D.R.B., L.D.B., T.R.Z.), National Food Safety & Toxicology Center, Michigan State University, East Lansing, Michigan
Received July 31, 2007; accepted October 17, 2007
| Abstract |
|---|
|
|
|---|
Like the ER, the aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor but is a member of the basic helix-loop-helix-PER/ARNT/SIM (periodicity/aryl hydrocarbon receptor nuclear translocator/simple-minded) family of transcription factors. The AhR is responsible for mediating many, if not all, of the toxic and biochemical responses to TCDD and related compounds. These include a wasting syndrome, tumor promotion, teratogenesis, hepatotoxicity, immunotoxicity, and modulation of endocrine systems, which are mediated in a tissue-, sex-, age-, and species-specific manner (Poland and Knutson, 1982
; Denison and Heath-Pagliuso, 1998
). The proposed mechanism involves ligand binding to the cytoplasmic AhR and translocation to the nucleus, where it heterodimerizes with the aryl hydrocarbon receptor nuclear translocator (ARNT), another member of the basic helix-loop-helix-PER/ARNT/SIM family. This complex then binds specific DNA elements, termed dioxin response elements (DREs), in the regulatory regions of target genes, leading to changes in gene expression (Hankinson, 1995
). Although the mechanisms of AhR/ARNT-mediated changes in gene expression are well established, how changes in gene expression results in the subsequent physiological and toxicological effects remains poorly understood.
|
Only a limited number of estrogen responsive genes that are inhibited by TCDD have been identified to support the ER-AhR gene expression cross-talk mechanism (Krishnan et al., 1995
; Gillesby et al., 1997
; Duan et al., 1999
; Porter et al., 2001
). However, a comprehensive assessment of inhibitory gene expression responses and their relationship to in vivo antiestrogenic endpoints has not been investigated. Therefore, to further examine the inhibitory effects of TCDD on estrogen-mediated uterine gene expression, temporal gene expression responses to EE, TCDD, and EE plus TCDD were investigated. Results indicate that the inhibitory effect of TCDD on EE-induced uterotrophy is associated with the selective inhibition of EE-mediated gene expression responses.
| Materials and Methods |
|---|
|
|
|---|
Animals (n = 5/treatment group/time point) were orally gavaged at time 0 with sesame oil vehicle (Sigma Chemical, St. Louis, MO), TCDD (provided by S. Safe, Texas A&M University, College Station, TX), 17
-ethynylestradiol (EE; Sigma Chemical) or a combination of EE plus TCDD followed by additional doses of vehicle (vehicle and TCDD groups) or EE (EE and EE + TCDD groups) at 24 and 48 h as per the uterotrophic assay (Fig. 1). Doses of 10 and 30 µg/kg EE and TCDD, respectively, were empirically determined to elicit an optimal inhibitory effect on the EE-mediated induction of uterine weight in cotreatment studies (data not shown). Mice were sacrificed by cervical dislocation 4, 12, 24, or 72 h after dosing. Uterine weights were recorded before (wet) and after (blotted) blotting with absorbent tissue. A section of the left uterine horn was removed for histologic examination and fixed in 10% neutral buffered formalin (Sigma). The remaining tissue was subsequently snap-frozen in liquid nitrogen and stored at -80°C. All doses were calculated based on average weights of the animals before dosing. All procedures were performed with the approval of the Michigan State University All-University Committee on Animal Use and Care.
Histological Processing and Assessment. Fixed uteri were embedded in paraffin according to standard histological techniques. Five-micrometer cross-sections were mounted on glass slides and stained with hematoxylin and eosin. Embedding, mounting and staining were performed at the Histology/Immunohistochemistry Laboratory, Michigan State University (http://humanpathology.msu.edu/histology/index.html). Histological slides were evaluated according to standardized National Toxicology Program (NTP) pathology codes. Morphometric analyses were performed for each sample using image analysis software (Scion Image; Scioncorp, Frederick, MD) and standard morphometric techniques. The length of basal lamina underlying the luminal epithelium (LE) and corresponding area of the luminal epithelial cells (LECs) was quantified for multiple representative sectors of each section to calculate LEC height.
RNA Isolation. Total RNA was isolated from uteri using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's protocol. 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 resuspended in RNA Storage Solution (Ambion, Austin, TX). RNA concentrations were calculated by spectrophotometric methods (A260) and purity assessed by the A260:A280 ratio and visual inspection of 1 µgona denaturing gel.
Microarray Experimental Design and Protocols. Spotted mouse cDNA microarrays were prepared in-house and consist of 13,361 features, representing 7,948 unique genes (Unigene Build 144). 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. In breef, PCR-amplified DNA was robotically arrayed onto epoxy-coated glass slides (Schott-Nexterion, Duryea, PA) using an Omnigrid arrayer (GeneMachines, San Carlos, CA) equipped with 48 (12 x 4) Chipmaker 2 pins (Telechem) at the Genomics Technology Support Facility at Michigan State University (http://www.genomics.msu.edu). Changes in uterine gene expression were assessed using a 2 x 2 factorial design (Fig. 2) (Yang and Speed, 2002
). In this design, arrow bases represent samples labeled with Cy3 and arrowheads represent samples labeled with Cy5. Within each replicate, a sample is labeled and hybridized on three independent arrays for a total of six arrays/replicate/time point. Three biological replicates were completed at each time point for a total of 72 microarrays. A 3DNA Array 900 Expression Array Detection Kit (Genisphere, Hatsfield, PA) using 1.0 µg of total RNA was used for probe labeling in all microarray experiments, according to manufacturer's specifications. Samples were hybridized for 18 to 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 an Affymetrix 428 Array Scanner (Santa Clara, CA). Images were analyzed for feature and background intensities using GenePix Pro 5.0 (Molecular Devices, Sunnyvale, CA).
|
Array Data Normalization and Statistical Analysis. Data were normalized using a semiparametric approach (Eckel et al., 2005
). Model-based t values were calculated from normalized data, comparing treated and vehicle 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 (Eckel et al., 2004
). Gene lists were filtered for activity based on the P1(t)-value which indicates a greater likelihood of activity as the value approaches 1.0. A conservative P1(t) cutoff of 0.9999 combined with a differential expression of ± 1.5-fold relative to time-matched vehicle controls (TMVC) was used to filter the expression data and to define active gene lists. All arrays were subjected to quality control assessment to ensure assay performance and data consistency (Burgoon et al., 2005
). Data are stored within dbZach (http://dbzach.fst.msu.edu), a MIAME supportive relational database that ensures proper data management and facilitates data analysis (Burgoon et al., 2006
). Complete data sets with annotation and P1(t) values are available in Supplementary Table 1. Gene expression patterns that passed the established threshold criteria were visualized using hierarchical clustering (GeneSpring 6.0; Silicon Genetics, Redwood City, CA).
Quantitative Real-Time PCR 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 forward and reverse gene-specific primers, designed using Primer3 (Rozen and Skaletsky, 2000
), 3 mM MgCl2, 1.0 mM dNTPs, 0.025 IU of AmpliTaq Gold, and 1x SYBR Green PCR buffer (Molecular Devices). Gene names, accession numbers, forward and reverse primer sequences, and amplicon sizes are listed in Supplementary Table 2. PCR amplification was conducted in MicroAmp Optical 96-well reaction plates (Molecular Devices) on an PRISM 7000 Sequence Detection System (Molecular Devices) 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 6 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 were also included on each plate. Samples with a Ct value within 2 SD of the mean Ct values for the No template controls were considered below the limits of detection. The copy number of each unknown sample for each gene was standardized to Rpl7 to control for differences in RNA loading, quality, and cDNA synthesis (Couse et al., 1995
). Statistical significance of differentially expressed genes was determined using two-way analysis of variance followed by a Tukey's post hoc test (SAS 9.1; SAS Institute, Cary, NC). For graphing purposes, the relative expression levels were scaled such that the expression level of the time-matched vehicle control was equal to 1.
|
| Results |
|---|
|
|
|---|
Histopathology and Morphometry. Treatment of mice with EE resulted in the expected complex uterine histologic examination results, consisting of minimal stromal edema at 4 h that progressed to moderate severity by 12 h. At 24 h, moderate epithelial cell hypertrophy and hyperplasia with moderate stromal edema were observed that progressed to marked epithelial and stromal hypertrophy and hyperplasia with mild stromal edema at 72 h (Fig. 4). Cotreatment of EE and TCDD resulted in comparable histological effects compared with EE treatment alone with the exception of reduced stromal edema at 12, 24, and 72 h, subnuclear vacuolization in epithelial cells at 24 h, and reduced stromal hypertrophy and hyperplasia with marked luminal epithelial cell (LEC) apoptosis at 72 h (Fig. 4). The LEC layer effects are consistent with previous reports of TCDD in the murine uterus (Gallo et al., 1986
). These alterations in uterine histologic results may contribute to the associated decreases in uterine weight observed with TCDD cotreatment.
|
Increased luminal epithelial cell height (LECH) is a well recognized marker of estrogen exposure and has been used to assess the estrogenicity of a number of structurally diverse ligands (O'Connor et al., 1996
; Nakagawa and Tayama, 2001
). EE-induced LECH at 24 and 72 h, as reported previously (Kwekel et al., 2005
). Cotreatment with TCDD did not inhibit this response. The inability to detect antiestrogenic effects on LECH may indicate that TCDD does not influence this response. Alternatively, it may be attributed technical difficulties in the measurement of such a change because of the complex pseudostratified nature of the proliferating LE cells combined with the histopathological alterations induced by TCDD on this cell layer.
Microarray Data Filtering and Clustering. Microarray analyses were preformed using a 2 x 2 factorial design that allowed for the identification of differentially expressed genes after EE and TCDD treatment alone as well as the interactive effects of EE + TCDD compared with the single treatment groups (Yang and Speed, 2002
). A conservative statistical P1(t) cutoff of 0.9999 combined with a differential expression of ±1.5 fold relative to TMVCs was used to identify lists of differentially active genes. Gene expression responses to EE alone displayed the expected complex transcriptional profile as reported previously (Fertuck et al., 2003
; Kwekel et al., 2005
) with a total of 3,746 features, representing 2,753 unique genes, identified as differentially expressed at one or more time points. The number and magnitude of uterine gene expression responses elicited by TCDD were modest compared with EE. Seven hundred ninety-three features representing 628 unique genes were found to be differentially expressed in response to TCDD. EE + TCDD cotreatment resulted in an overall gene expression response similar to that of EE alone, with a total of 3,631 features representing 2,647 unique genes identified as differentially expressed at one or more time points.
To compare the global gene expression responses of EE, TCDD, and EE + TCDD, hierarchical clustering was performed on features that were differentially expressed in any of the treatment groups at any time point relative to the TMVCs. Visualization of the global responses for each treatment group revealed that the temporal expression pattern of EE + TCDD was essentially indistinguishable from EE alone (Fig. 5A). In addition, TCDD displayed similarities to EE, consistent with our previous studies describing the estrogen-like gene expression profile of TCDD (Boverhof et al., 2006
). Clustering by treatment and time point further revealed the temporal similarity of the EE and EE + TCDD treatment groups as each of their gene expression time points clustered with one another (Fig. 5B). The 12-h gene expression responses to TCDD, EE, and EE + TCDD also clustered together, further demonstrating the estrogen-like response to TCDD at this time point.
|
|
Genes Differentially Regulated by EE and EE + TCDD Treatments. Genes identified as differentially expressed by EE + TCDD cotreatment compared with EE treatment alone were identified and investigated further. To be considered in this category, two successive criteria were required to be met. First, these genes needed to exhibit differential expression after EE treatment relative to TM-VCs; second, these genes needed to exhibit differential gene expression after EE + TCDD cotreatment relative to EE. This approach identified genes regulated by EE that were subsequently modulated upon cotreatment with TCDD and identified 163 features representing 133 EE-regulated genes at one or more time points. In many cases, EE-mediated gene expression responses were inhibited by more than 80% upon cotreatment with TCDD. The data also indicate that only a select number of EE-mediated gene expression responses experience inhibition after cotreatment with TCDD, most EE responsive genes being unaffected by the cotreatment. On a per-time-point basis, 9, 23, 32, and 130 features representing 5, 21, 28, and 106 genes were inhibited by TCDD cotreatment at 4, 12, 24, and 72 h, respectively. This indicates a time-dependent increase in the inhibitory effects of TCDD on EE-mediated gene expression responses, suggesting that direct early primary responses may subsequently mediate more extensive secondary and tertiary indirect inhibitory responses.
A small number of gene expression responses were differentially expressed between the EE + TCDD and EE groups but were not EE-regulated responses; rather, they were differentially expressed due to TCDD alone. This included the well characterized induction of Cyp1a1 as well as the induction of inhibitor of growth 1 (Ing1), karyopherin alpha 6 (Kpna6), and replication protein A2 (Rpa2). The induction of these genes cannot be dismissed as a contributing factor to the antiestrogenic effects of TCDD because the induction of inhibitory factors is a previously proposed mechanism (Rogers and Denison, 2002
).
Functional Categorization of Microarray Data. Functional annotation of gene expression responses was performed using data extracted from public databases and published literature. The functions of EE differentially expressed genes have been associated with transcription factors, mRNA and protein synthesis, cell cycle regulation, cellular proliferation, energetics and structural constituents (Fertuck et al., 2003
; Moggs et al., 2004
; Kwekel et al., 2005
). Functional annotation of EE-mediated gene expression responses that were inhibited upon cotreatment with TCDD were associated with the regulation of cell proliferation and growth, water/ion transport, and the maintenance of cellular structural architecture (Table 1). Cellular growth and proliferation genes included branched chain aminotransferase 1 (Bcat), serine proteinase inhibitor B5 (Serpinb5), sestrin 1 (Sesn1), stratifin (Sfn), and trefoil factor 1 (Tff1). Inhibition of this functional category is consistent with previous reports of decreased cellular growth responses in breast and endometrial cancer cell lines (Wang et al., 1998
; Castro-Rivera et al., 1999
; Puga et al., 2000
; Wormke et al., 2000
) and uterine tissue (Buchanan et al., 2002
). TCDD-inhibited water and ion transport genes included aquaporins 1 and 3, (Aqp1 and 3) solute carriers 4a2, 38a3, and 40a1 (Slc4a2, 38a3 and 40a1), and FXYD ion transport regulator 4 (Fxyd4). Inhibition of these responses may contribute to TCDD-mediated decreases in stromal edema and uterine wet weight. Desmocollin 2 (Dsc2), keratins 4, 7, 14 and 19 (Krt2-4, Krt2-7, Krt1-14 and Krt1-19), macrophage receptor with collagenous structure (Marco), TP53 apoptosis effector (Perp), and small prolinerich protein 2A (Sprr2a) represent structural genes inhibited by TCDD. Collectively, the inhibition of these EE-mediated responses may contribute to the antiestrogenic effects of TCDD on uterine histology, growth, and LEC integrity.
|
Verification of Microarray Results. QRTPCR was used to verify changes in transcript levels for a selected subset of EE-inducible genes inhibited by TCDD (Fig. 7). There was good agreement between the microarray and QRTPCR results, although compression of the gene expression response was observed in the microarray data, which has been previously reported comparing microarray analysis to other methods (Yuen et al., 2002
). QRTPCR revealed that Tff1 transcripts were induced more than 400-fold by EE at 12 and 24 h, whereas cotreatment with TCDD inhibited this response by more than 90%. Similar confirmatory responses were noted for Dsc2, Krt1-14, Sprr2a, and Sfn, which were maximally induced 20-, 171-, 206-, and 6.-fold by EE treatment and inhibited 95, 90, 83, and 93%, respectively, by TCDD. QRTPCR was also used to verify the EE induction of proliferating cell nuclear antigen (PCNA) and solute carrier family 25, member 5 (Slc25a5), because microarray analysis indicated that these genes were not affected by TCDD cotreatment. Induction levels of PCNA and Slc25a5 by EE and EE + TCDD were comparable, verifying that TCDD inhibited select EE induced responses.
|
| Discussion |
|---|
|
|
|---|
Inhibition of Cellular Growth and Proliferation Responses. Estrogen induction of uterine weight involves a coordinated proliferative response that is mediated through a well orchestrated series of changes in gene expression (Fertuck et al., 2003
; Moggs et al., 2004
; Kwekel et al., 2005
). Cotreatment with TCDD disrupted several EE-induced genes with important functions in cell cycle regulation, growth, and proliferation. For example, EE-mediated induction of Bcat and Sfn, important regulators of cell cycle progression, was inhibited by TCDD between 12 and 72 h. Bcat regulates G1-to-S phase transition, and cells with reduced expression exhibit faster growth rates, a shorter G1 stage, and an increased frequency of mutations (Schuldiner et al., 1996
). Sfn serves as a G2 checkpoint component as a positive mediator of growth-factor-induced cell cycle progression (Hermeking et al., 1997
; Zhang et al., 2004
). TCDD also inhibited EE induction of Serpinb5 at 12, 24, and 72 h, which plays an essential role in development, as exhibited by embryonic lethality in knockout mice (Gao et al., 2004
), whereas decreased expression results in reduced cellular proliferation and adhesion (Gao et al., 2004
). Sesn1 is a positive regulator of cell growth and protects against apoptosis (Velasco-Miguel et al., 1999
; Budanov et al., 2004
) and was also inhibited by TCDD cotreatment. Additional EE-induced genes implicated in cellular growth, proliferation, and development included retinol binding protein 2, Tnfsf8, and Vezf1, which were also inhibited upon cotreatment between 12 and 72 h. The inhibition of genes involved in regulating cell cycle progression is consistent with that reported in a previous study; however, the inhibition of estrogen-induced cyclin transcripts was not detected, which may be attributed to different experimental treatments and time points (Buchanan et al., 2002
). Overall, the alteration of these responses may be a contributing factor to the observed reduction in stromal cell hypertrophy and hyperplasia as well as the marked LEC apoptosis.
One of the most dramatic TCDD-inhibited responses was that of Tff1, which plays a fundamental role in epithelial maintenance, protection, and regeneration (Lefebvre et al., 1996
; Playford et al., 1996
). Tffs block p53-dependent and independent pathways of apoptosis and promote growth and regeneration by allowing cells to break attachments with the basement membrane to replace epithelial defects without cell death (Hoffmann et al., 2001
). Tff peptides also have anti-inflammatory actions and protect the epithelial mucous layers (Hoffmann et al., 2001
; Vieten et al., 2005
). Therefore, TCDD's inhibition of EE-induced Tff1 by more than 90% may play an important role in the increased LEC degeneration and apoptosis.
Tff1 is also a prognostic marker in human breast cancer and is an estrogen-responsive gene in breast cancer cells and the human endometrium (Gillesby and Zacharewski, 1999
; Punyadeera et al., 2005
). The inhibitory effect of TCDD on estrogen induction of Tff1 has been characterized in human MCF-7 breast cancer cells and depends on an inhibitory DRE that interferes with AP-1- and ERE-mediated transcriptional activation (Gillesby et al., 1997
). Tff1 induction and inhibition by TCDD in human MCF-7 breast cancer cells and the mouse uterus indicates that this mechanism may be conserved across sensitive species and tissues. Examination of the mouse promoter region for Tff1 identified a variant ERE at -475 as well as an AP-1 site at -998 relative to the transcriptional start site. Although a DRE does not overlap with the AP-1 site, two putative DREs are located further upstream at -1920 and -2637 and may play a role in mediating the inhibition.
Water/Ion Transport Responses. EE induced stromal edema, which was significantly inhibited by TCDD. Several EE-regulated genes involved in water and ion transport were inhibited by TCDD, including the EE-mediated down-regulation of Aqp1 and up-regulation of Aqp3. Isoform-specific regulation of aquaporins plays an integral role in mediating the water imbibition in the uterus (Richard et al., 2003
). EE-regulated transcripts involved in sodium and chloride transport were also inhibited by TCDD, including Slc4a2, Slc38a3, and Fxyd4 (Garty et al., 2003
; Quentin et al., 2004
). Collectively, the inhibition of these EE-mediated gene expression responses may have contributed to reductions in uterine wet weight after cotreatment with TCDD.
Inhibition of Structural Constituents. The uterus undergoes extensive cytoarchitectural changes to accommodate the dramatic proliferation and growth response to estrogen, which involves numerous structural, adhesion, and extracellular matrix genes, including a number of keratins, actins, procollagens, tubulins, desmocollins, and small proline-rich proteins in mice, rats, and humans (Watanabe et al., 2002
; Hewitt et al., 2003
; Watanabe et al., 2003
; Moggs et al., 2004
). TCDD cotreatment inhibited a number of these EE-mediated responses, which probably contributed to its anti-uterotrophic effects as well as increases in apoptosis.
In this study, the EE induction of keratins 4, 7, 14, and 19 was significantly inhibited by TCDD cotreatment. Keratins are involved in the formation of the cytoskeleton, which consists of an extensive array of filamentous networks. Their disruption results in epithelial cell fragility and lysis (Sørensen et al., 2003
; Wong et al., 2005
). Keratins 18 and 19 are estrogen-inducible transcripts, the induction of which is blocked by TCDD in MCF-7 cells (Chen et al., 2001
). Inhibition of uterine keratin 18 and 19 suggests that this may represent a conserved response between rodents and humans. Moreover, the inhibition of multiple keratin genes suggests that TCDD may disrupt signaling at a common regulatory region as the basic (Krt2-2 through -8) and acidic (Krt1-9 through -19) keratin genes are encoded in a tandem array on chromosomes 15 and 11, respectively (Chu and Weiss, 2002
).
Desmocollin 2 (Dsc2), a component of desmosomes, is expressed primarily in epithelial cells, serves an integral role in cell adhesion by forming links with the intermediate filament network (Marsden et al., 1997
), and was inhibited by TCDD. TCDD also inhibited Perp induction, which promotes desmosomal complex assembly (Ihrie et al., 2005
). EE induction of small proline rich protein 2a (Sprr2a) was also inhibited by cotreatment with TCDD. The Sprr2 family consists of 11 genes (Sprr2a-2k), which are important structural components of epithelial cells because of their ability to form extensive cross-links (Hong et al., 2004
). Several Sprr2 genes are up-regulated in the luminal epithelial cells of the uterus in response to estrogen, where they are important for cytoarchitectural changes (Hong et al., 2004
). TCDD also inhibited the induction of other structural molecules including Marco, procollagen 6a2, troponin T1, and tubulin β6. Together, the inhibition of these structural constituents could compromise the rapid proliferation and growth induced by EE, resulting in altered results on histologic examination, increased apoptosis, and overall decreased uterine growth.
TCDD as Estrogen and Antiestrogen. Reports have indicated that TCDD elicits an estrogen-like, ER-dependent gene expression profile in the uterus (Ohtake et al., 2003
; Watanabe et al., 2004
; Boverhof et al., 2006
). The regulation of similar genes by EE and TCDD suggests that these responses may represent targets for inhibition. However, these genes were largely unaffected by TCDD, including well characterized estrogen-responsive genes, such as PCNA, Slc25a5, cell division cycle 2 homolog A (Cdc2a), and ornithine decarboxylase (Odc). Instead, many of the inhibited responses were unaffected by TCDD treatment alone, consistent with previous reports of estrogen/TCDD gene expression cross-talk (Zacharewski et al., 1994
; Porter et al., 2001
; Safe and Wormke, 2003
). These data suggest that the anti-uterotrophic effects of TCDD are independent of its weak estrogenic activity.
Decreased ER levels (Romkes et al., 1987
) and increased estrogen metabolism (Spink et al., 1994
) have also been proposed as mechanisms for the antiestrogenic effects of TCDD. However, several reports indicate that these mechanisms do not account for the antiestrogenic effects of TCDD. TCDD does not increase estrogen metabolism in vivo (DeVito et al., 1992
; Petroff and Mizinga, 2003
), and uterine ER levels were unaffected by TCDD (DeVito et al., 1994
; White et al., 1995
). In the present study, TCDD inhibited only a select subset of EE responses, further indicating that increased metabolism and decreased ER levels are not the primary mechanisms, because these would be expected to elicit a more global effect on gene expression. The gene-specific inhibitory effects may involve an inhibitory DRE mechanism, alteration in the activity or availability of gene-specific cofactors, or induction/inhibition of an upstream regulatory protein that are likely to be tissue and species-specific.
Summary. The present study has identified a small subset of EE-induced uterine gene expression responses that are inhibited by TCDD. Moreover, the repressed functional categories can be related to the observed histological and physiological responses and therefore represent potential mediators of TCDD's anti-uterotrophic effect. Select responses, including Tff1, are also in agreement with in vitro studies, indicating the potential conservation of these responses among different models and species. Furthermore, the results indicate that the estrogenic and antiestrogenic gene expression effects of TCDD are independent. The importance of these genes in the uterotrophic response as targets for the antiestrogenicity of TCDD can be addressed through comparative studies with rats and mice and through receptor- and gene-specific null mice. Additional research is also required to fully delineate the mechanism of the gene-specific inhibitory response and should employ molecular approaches including receptor interaction (fluorescence resonance energy transfer), chromatin immunoprecipitation, promoter dissection/reporter gene assays and siRNA approaches in appropriate in vitro models (e.g., uterine based cell lines, primary cells or cocultured primary systems) that closely simulate in vivo conditions. These approaches will more comprehensively elucidate the dual nature of TCDD as an estrogenic and antiestrogenic compound.
| Acknowledgements |
|---|
| Footnotes |
|---|
ABBREVIATIONS: ER, estrogen receptor; ERE, estrogen response element; AP-1, activator protein-1; AhR, aryl hydrocarbon receptor; TCDD, 2,3,7,8-tetrachlorodibenzo-p-dioxin; ARNT, aryl hydrocarbon receptor nuclear translocator; DRE, dioxin response element; EE, 17
-ethynylestradiol; LE, luminal epithelium; LEC, luminal epithelial cell; PCR, polymerase chain reaction; TMVC, time-matched vehicle controls; QRTPCR, quantitative real-time PCR; LECH, luminal epithelial cell height; PCNA, proliferating cell nuclear antigen.
The online version of this article (available at http://molpharm.aspetjournals.org) contains supplemental material. ![]()
Address correspondence to: Dr. Timothy R. Zacharewski, Michigan State University, Department of Biochemistry and Molecular Biology, 224 Biochemistry Building, Wilson Road, East Lansing, MI, USA 48824-1319. E-mail: tzachare{at}msu.edu
| References |
|---|
|
|
|---|
Astroff B, and Safe S (1988) Comparative antiestrogenic activities of 2,3,7,8-tetrachlorodibenzo-p-dioxin and 6-methyl-1,3,8-trichlorodibenzofuran in the female rat. Toxicol Appl Pharmacol 95: 435-443.[CrossRef][Medline]
Boverhof DR, Kwekel JC, Humes DG, Burgoon LD, and Zacharewski TR (2006) Dioxin induces an estrogen-like, estrogen receptor-dependent gene expression response in the murine uterus. Mol Pharmacol 69: 1599-1606.
Buchanan DL, Ohsako S, Tohyama C, Cooke PS, and Iguchi T (2002) Dioxin inhibition of estrogen-induced mouse uterine epithelial mitogenesis involves changes in cyclin and transforming growth factor-beta expression. Toxicol Sci 66: 62-68.
Budanov AV, Sablina AA, Feinstein E, Koonin EV, and Chumakov PM (2004) Regeneration of peroxiredoxins by p53-regulated sestrins, homologs of bacterial AhpD. Science 304: 596-600.
Burgoon LD, Boutros PC, Dere E, and Zacharewski TR (2006) dbZach: A MIAME-compliant toxicogenomic supportive relational database. Toxicol Sci 90: 558-568.
Burgoon LD, Eckel-Passow JE, Gennings C, Boverhof DR, Burt JW, Fong CJ and Zacharewski TR (2005) Protocols for the assurance of microarray data quality and process control. Nucleic Acids Res 33: e172-.
Castro-Rivera E, Wormke M and Safe S (1999) Estrogen and aryl hydrocarbon responsiveness of ECC-1 endometrial cancer cells. Mol Cell Endocrinol 150: 11-21.[CrossRef][Medline]
Chen I, Hsieh T, Thomas T and Safe S (2001) Identification of estrogen-induced genes downregulated by AhR agonists in MCF-7 breast cancer cells using suppression subtractive hybridization. Gene 262: 207-214.[CrossRef][Medline]
Chu PG and Weiss LM (2002) Keratin expression in human tissues and neoplasms. Histopathology 40: 403-439.[CrossRef][Medline]
Couse J, Curtis S, Washburn T, Lindzey J, Golding T, Lubahn D, Smithies O, and Korach K (1995) Analysis of transcription and estrogen insensitivity in the female mouse after targeted disruption of the estrogen receptor gene. Mol Endocrinol 9: 1441-1454.
Denison MS, and Heath-Pagliuso S (1998) The Ah receptor: a regulator of the biochemical and toxicological actions of structurally diverse chemicals. Bull Environ Contam Toxicol 61: 557-568.[CrossRef][Medline]
DeVito MJ, Ma X, Babish JG, Menache M, and Birnbaum LS (1994) Dose-response relationships in mice following subchronic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin: CYP1A1, CYP1A2, estrogen receptor, and protein tyrosine phosphorylation. Toxicol Appl Pharmacol 124: 82-90.[CrossRef][Medline]
DeVito MJ, Thomas T, Martin E, Umbreit TH, and Gallo MA (1992) Antiestrogenic action of 2,3,7,8-tetrachlorodibenzo-p-dioxin: tissue-specific regulation of estrogen receptor in CD1 mice. Toxicol Appl Pharmacol 113: 284-292.[CrossRef][Medline]
Diel P, Schmidt S and Vollmer G (2002) In vivo test systems for the quantitative and qualitative analysis of the biological activity of phytoestrogens. J Chromatogr B Analyt Technol Biomed Life Sci 777: 191-202.[Medline]
Duan R, Porter W, Samudio I, Vyhlidal C, Kladde M, and Safe S (1999) Transcriptional activation of c-fos protooncogene by 17beta-estradiol: mechanism of aryl hydrocarbon receptor-mediated inhibition. Mol Endocrinol 13: 1511-1521.
Eckel JE, Gennings C, Chinchilli VM, Burgoon LD, and Zacharewski TR (2004) Empirical bayes gene screening tool for time-course or dose-response microarray data. J Biopharm Stat 14: 647-670.[CrossRef][Medline]
Eckel JE, Gennings C, Therneau TM, Burgoon LD, Boverhof DR, and Zacharewski TR (2005) Normalization of two-channel microarray experiments: a semiparametric approach. Bioinformatics 21: 1078-1083.
Fertuck KC, Eckel JE, Gennings C, and Zacharewski TR (2003) Identification of temporal patterns of gene expression in the uteri of immature, ovariectomized mice following exposure to ethynylestradiol. Physiol Genomics 15: 127-141.
Gallo MA, Hesse EJ, Macdonald GJ and Umbreit TH (1986) Interactive effects of estradiol and 2,3,7,8-tetrachlorodibenzo-p-dioxin on hepatic cytochrome P-450 and mouse uterus. Toxicol Lett 32: 123-132.[CrossRef][Medline]
Gao F, Shi HY, Daughty C, Cella N, and Zhang M (2004) Maspin plays an essential role in early embryonic development. Development 131: 1479-1489.
Garty H, Lindzen M, Fuzesi M, Aizman R, Goldshleger R, Asher C, and Karlish SJ (2003) A specific functional interaction between CHIF and Na, K-ATPase: role of FXYD proteins in the cellular regulation of the pump. Ann N Y Acad Sci 986: 395-400.[Medline]
Gillesby BE, Stanostefano M, Porter W, Safe S, Wu ZF, and Zacharewski TR (1997) Identification of a motif within the 5' regulatory region of pS2 which is responsible for AP-1 binding and TCDD-mediated suppression. Biochemistry 36: 6080-6089.[CrossRef][Medline]
Gillesby BE, and Zacharewski TR (1999) pS2 (TFF1) levels in human breast cancer tumor samples: correlation with clinical and histological prognostic markers. Breast Cancer Res Treat 56: 253-265.[CrossRef][Medline]
Hall JM, Couse JF, and Korach KS (2001) The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem 276: 36869-36872.
Hankinson O (1995) The aryl hydrocarbon receptor complex. Annu Rev Pharmacol Toxicol 35: 307-340.[CrossRef][Medline]
Hermeking H, Lengauer C, Polyak K, He TC, Zhang L, Thiagalingam S, Kinzler KW, and Vogelstein B (1997) 14-3-3 sigma is a p53-regulated inhibitor of G2/M progression. Mol Cell 1: 3-11.[CrossRef][Medline]
Hewitt SC, Deroo BJ, Hansen K, Collins J, Grissom S, Afshari CA, and Korach KS (2003) Estrogen receptor dependent genomic responses in the uterus mirror the biphasic physiological response to estrogen. Mol Endocrinol 17: 2070-2083.
Hewitt SC, Harrell JC, and Korach KS (2005) Lessons in estrogen biology from knockout and transgenic animals. Annu Rev Physiol 67: 285-308.[CrossRef][Medline]
Hoffmann W, Jagla W, and Wiede A (2001) Molecular medicine of TFF-peptides: from gut to brain. Histol Histopathol 16: 319-334.[Medline]
Hong SH, Nah HY, Lee JY, Lee YJ, Lee JW, Gye MC, Kim CH, Kang BM, and Kim MK (2004) Estrogen regulates the expression of the small proline-rich 2 gene family in the mouse uterus. Mol Cells 17: 477-484.[Medline]
Ihrie RA, Marques MR, Nguyen BT, Horner JS, Papazoglu C, Bronson RT, Mills AA, and Attardi LD (2005) Perp is a p63-regulated gene essential for epithelial integrity. Cell 120: 843-856.[CrossRef][Medline]
Klinge CM (2001) Estrogen receptor interaction with estrogen response elements. Nucleic Acids Res 29: 2905-2919.
Klinge CM, Bowers JL, Kulakosky PC, Kamboj KK and Swanson HI (1999) The aryl hydrocarbon receptor (AHR)/AHR nuclear translocator (ARNT) heterodimer interacts with naturally occurring estrogen response elements. Mol Cell Endocrinol 157: 105-119.[CrossRef][Medline]
Kociba RJ, Keyes DG, Beyer JE, Carreon RM, Wade CE, Dittenber DA, Kalnins RP, Frauson LE, Park CN, Barnard SD, et al. (1978) Results of a two-year chronic toxicity and oncogenicity study of 2,3,7,8-tetrachlorodibenzo-p-dioxin in rats. Toxicol Appl Pharmacol 46: 279-303.[CrossRef][Medline]
Krishnan V, Porter W, Santostefano M, Wang X, and Safe S (1995) Molecular mechanism of inhibition of estrogen-induced cathepsin D gene expression by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in MCF-7 cells. Mol Cell Biol 15: 6710-6719.[Abstract]
Kwekel JC, Burgoon LD, Burt JW, Harkema JR, and Zacharewski TR (2005) A cross-species analysis of the rodent uterotrophic program: elucidation of conserved responses and targets of estrogen signaling. Physiol Genomics 23: 327-342.
Lefebvre O, Chenard MP, Masson R, Linares J, Dierich A, LeMeur M, Wendling C, Tomasetto C, Chambon P, and Rio MC (1996) Gastric mucosa abnormalities and tumorigenesis in mice lacking the pS2 trefoil protein. Science 274: 259-262.
Marsden MD, Collins JE, Greenwood MD, Adams MJ, Fleming TP, Magee AI, and Buxton RS (1997) Cloning and transcriptional analysis of the promoter of the human type 2 desmocollin gene (DSC2). Gene 186: 237-247.[CrossRef][Medline]
Moggs JG, Tinwell H, Spurway T, Chang HS, Pate I, Lim FL, Moore DJ, Soames A, Stuckey R, Currie R, et al. (2004) Phenotypic anchoring of gene expression changes during estrogen-induced uterine growth. Environ Health Perspect 112: 1589-1606.[Medline]
Nakagawa Y, and Tayama K (2001) Estrogenic potency of benzophenone and its metabolites in juvenile female rats. Arch Toxicol 75: 74-79.[CrossRef][Medline]
Nilsson S, Makela S, Treuter E, Tujague M, Thomsen J, Andersson G, Enmark E, Pettersson K, Warner M, and Gustafsson J-A (2001) Mechanisms of estrogen action. Physiol Rev 81: 1535-1565.
O'Connor JC, Cook JC, Craven SC, Van Pelt CS, and Obourn JD (1996) An in vivo battery for identifying endocrine modulators that are estrogenic or dopamine regulators. Fundam Appl Toxicol 33: 182-195.[CrossRef][Medline]
Ohtake F, Takeyama K, Matsumoto T, Kitagawa H, Yamamoto Y, Nohara K, Tohyama C, Krust A, Mimura J, Chambon P, et al. (2003) Modulation of oestrogen receptor signalling by association with the activated dioxin receptor. Nature 423: 545-550.[CrossRef][Medline]
Petroff BK, and Mizinga KM (2003) Pharmacokinetics of ovarian steroids in Sprague-Dawley rats after acute exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Reprod Biol 3: 131-141.[Medline]
Playford RJ, Marchbank T, Goodlad RA, Chinery RA, Poulsom R, and Hanby AM (1996) Transgenic mice that overexpress the human trefoil peptide pS2 have an increased resistance to intestinal damage. Proc Natl Acad Sci U S A 93: 2137-2142.
Poland A, and Knutson JC (1982) 2,3,7,8-Tetrachlorodibenzo-p-dioxin and related halogenated aromatic hydrocarbons: examination of the mechanism of toxicity. Annu Rev Pharmacol Toxicol 22: 517-554.[CrossRef][Medline]
Porter W, Wang F, Duan R, Qin C, Castro-Rivera E, Kim K, and Safe S (2001) Transcriptional activation of heat shock protein 27 gene expression by 17betaestradiol and modulation by antiestrogens and aryl hydrocarbon receptor agonists. J Mol Endocrinol 26: 31-42.[Abstract]
Puga A, Barnes SJ, Dalton TP, Chang C, Knudsen ES, and Maier MA (2000) Aromatic hydrocarbon receptor interaction with the retinoblastoma protein potentiates repression of E2F-dependent transcription and cell cycle arrest. J Biol Chem 275: 2943-2950.
Punyadeera C, Dassen H, Klomp J, Dunselman G, Kamps R, Dijcks F, Ederveen A, de Goeij A, and Groothuis P (2005) Oestrogen-modulated gene expression in the human endometrium. Cell Mol Life Sci 62: 239-250.[CrossRef][Medline]
Quentin F, Eladari D, Frische S, Cambillau M, Nielsen S, Alper SL, Paillard M, and Chambrey R (2004) Regulation of the Cl-/HCO3-exchanger AE2 in rat thick ascending limb of Henle's loop in response to changes in acid-base and sodium balance. J Am Soc Nephrol 15: 2988-2997.
Richard C, Gao J, Brown N, and Reese J (2003) Aquaporin water channel genes are differentially expressed and regulated by ovarian steroids during the periimplantation period in the mouse. Endocrinology 144: 1533-1541.
Rogers JM, and Denison MS (2002) Analysis of the antiestrogenic activity of 2,3,7,8-tetrachlorodibenzo-p-dioxin in human ovarian carcinoma BG-1 cells. Mol Pharmacol 61: 1393-1403.
Romkes M, Piskorska-Pliszczynska J, and Safe S (1987) Effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin on hepatic and uterine estrogen receptor levels in rats. Toxicol Appl Pharmacol 87: 306-314.[CrossRef][Medline]
Rozen S and Skaletsky HJ (2000) Primer3 on the WWW for general users and for biologist programmers, in Bioinformatics Methods and Protocols: Methods in Molecular Biology (Krawetz S MS ed) pp 365-386, Humana Press, Totowa, NJ.
Safe S, and Wormke M (2003) Inhibitory aryl hydrocarbon receptor-estrogen receptor alpha cross-talk and mechanisms of action. Chem Res Toxicol 16: 807-816.[Medline]
Schuldiner O, Eden A, Ben-Yosef T, Yanuka O, Simchen G, and Benvenisty N (1996) ECA39, a conserved gene regulated by c-Myc in mice, is involved in G1/S cell cycle regulation in yeast. Proc Natl Acad Sci U S A 93: 7143-7148.
Sørensen CB, Andresen BS, Jensen UB, Jensen TG, Jensen PK, Gregersen N, and Bolund L (2003) Functional testing of keratin 14 mutant proteins associated with the three major subtypes of epidermolysis bullosa simplex. Exp Dermatol 12: 472-479.[CrossRef][Medline]
Spink DC, Hayes CL, Young NR, Christou M, Sutter TR, Jefcoate CR and Gierthy JF (1994) The effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin on estrogen metabolism in MCF-7 breast cancer cells: evidence for induction of a novel 17 beta-estradiol 4-hydroxylase. J Steroid Biochem Mol Biol 51: 251-258.[CrossRef][Medline]
Umbreit TH, Hesse EJ, Macdonald GJ, and Gallo MA (1988) Effects of TCDD-estradiol interactions in three strains of mice. Toxicol Lett 40: 1-9.[CrossRef][Medline]
Velasco-Miguel S, Buckbinder L, Jean P, Gelbert L, Talbott R, Laidlaw J, Seizinger B, and Kley N (1999) PA26, a novel target of the p53 tumor suppressor and member of the GADD family of DNA damage and growth arrest inducible genes. Oncogene 18: 127-137.[CrossRef][Medline]
Vieten D, Corfield A, Carroll D, Ramani P, and Spicer R (2005) Impaired mucosal regeneration in neonatal necrotising enterocolitis. Pediatr Surg Int 21: 153-160.[CrossRef][Medline]
Wang W, Smith I, Roger and Safe S (1998) Aryl hydrocarbon receptor-mediated antiestrogenicity in MCF-7 cells: modulation of hormone-induced cell cycle enzymes. Arch Biochem Biophys 356: 239-248.[CrossRef][Medline]
Watanabe H, Suzuki A, Goto M, Ohsako S, Tohyama C, Handa H, and Iguchi T (2004) Comparative uterine gene expression analysis after dioxin and estradiol administration. J Mol Endocrinol 33: 763-771.
Watanabe H, Suzuki A, Kobayashi M, Takahashi E, Itamoto M, Lubahn DB, Handa H, and Iguchi T (2003) Analysis of temporal changes in the expression of estrogen-regulated genes in the uterus. J Mol Endocrinol 30: 347-358.[Abstract]
Watanabe H, Suzuki A, Mizutani T, Khono S, Lubahn DB, Handa H, and Iguchi T (2002) Genome-wide analysis of changes in early gene expression induced by oestrogen. Genes Cells 7: 497-507.[Abstract]
White TE, Rucci G, Liu Z, and Gasiewicz TA (1995) Weanling female Sprague-Dawley rats are not sensitive to the antiestrogenic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Toxicol Appl Pharmacol 133: 313-320.[CrossRef][Medline]
Wong P, Domergue R, and Coulombe PA (2005) Overcoming functional redundancy to elicit pachyonychia congenita-like nail lesions in transgenic mice. Mol Cell Biol 25: 197-205.
Wormke M, Castro-Rivera E, Chen I, and Safe S (2000) Estrogen and aryl hydrocarbon receptor expression and crosstalk in human Ishikawa endometrial cancer cells. J Steroid Biochem Mol Biol 72: 197-207.[CrossRef][Medline]
Yang YH, and Speed T (2002) Design issues for cDNA microarray experiments. Nat Rev Genet 3: 579-588.[CrossRef][Medline]
Yuen T, Wurmbach E, Pfeffer RL, Ebersole BJ, and Sealfon SC (2002) Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res 30: e48.
Zacharewski TR, Bondy KL, McDonell P, and Wu ZF (1994) Antiestrogenic effect of 2,3,7,8-tetrachlorodibenzo-p-dioxin on 17 beta-estradiol-induced pS2 expression. Cancer Res 54: 2707-2713.
Zhang Y, Karas M, Zhao H, Yakar S, and LeRoith D (2004) 14-3-3sigma mediation of cell cycle progression is p53-independent in response to insulin-like growth factor-I receptor activation. J Biol Chem 279: 34353-34360.
This article has been cited by other articles:
![]() |
G. P. Daston Gene Expression, Dose-Response, and Phenotypic Anchoring: Applications for Toxicogenomics in Risk Assessment Toxicol. Sci., October 1, 2008; 105(2): 233 - 234. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||