- Research article
- Open Access
Species-specific regulation of PXR/CAR/ER-target genes in the mouse and rat liver elicited by o, p'-DDT
BMC Genomics volume 9, Article number: 487 (2008)
Dichlorodiphenyltrichloroethane (DDT) is a persistent estrogenic organochlorine pesticide that is a rodent hepatic tumor promoter, with inconclusive carcinogenicity in humans. We have previously reported that o, p'-DDT elicits primarily PXR/CAR-mediated activity, rather than ER-mediated hepatic responses, and suggested that CAR-mediated effects, as opposed to ER-mediated effects, may be more important in tumor promotion in the rat liver. To further characterize species-specific hepatic responses, gene expression analysis, with complementary histopathology and tissue level analyses were investigated in immature, ovariectomized C57BL/6 mice treated with 300 mg/kg o, p'-DDT, and compared to Sprague-Dawley rat data.
Rats and mice exhibited negligible histopathology with rapid o, p'-DDT metabolism. Gene expression profiles were also similar, exhibiting PXR/CAR regulation with the characteristic induction of Cyp2b10 and Cyp3a11. However, PXR-specific target genes such as Apoa4 or Insig2 exhibited more pronounced induction compared to CAR-specific genes in the mouse. In addition, mouse Car mRNA levels decreased, possibly contributing to the preferential activation of mouse PXR. ER-regulated genes Cyp17a1 and Cyp7b1 were also induced, suggesting o, p'-DDT also elicits ER-mediated gene expression in the mouse, while ER-mediated effects were negligible in the rat, possibly due to the inhibitory effects of CAR on ER activities. In addition, o, p'-DDT induced Gadd45a, Gadd45b and Cdkn1, suggesting DNA damage may be an additional risk factor. Furthermore, elevated blood DHEA-S levels at 12 h after treatment in the mouse may also contribute to the endocrine-related effects of o, p'-DDT.
Although DDT is known to cause rodent hepatic tumors, the marked species differences in PXR/CAR structure, expression patterns and ligand preference as well as significant species-specific differences in steroidogenesis, especially CYP17A1 expression and activity, confound the extrapolation of these results to humans. Nevertheless, the identification of potential modes of action as well as species-specific responses may assist in the selection and further development of more appropriate models for assessing the toxicity of DDT to humans and wildlife.
Dichlorodiphenyltrichloroethane (DDT) is a persistent organochlorine pesticide that is a hepatic tumor promoter in the mouse and rat , with inconclusive carcinogenicity in the human liver [2, 3]. Technical grade DDT is a mixture of p,p'-DDT, o, p'-DDT and their metabolites such as 1,1-dichloro-2,2-bis(p-chlorophenyl) ethylene (DDE) and l,l-dichloro-2,2-bis(p,p-chlorophenylethane) (DDD) . o, p'-DDT exerts estrogenic effects in responsive tissues such as the uterus  by binding to estrogen receptor (ER)α and β subtypes . In addition, o, p'-DDT is an agonist for pregnane X receptor (PXR) and constitutive androstane receptor (CAR) [7, 8]. Thus, it may elicit complex responses through multiple nuclear receptors.
Several structurally diverse estrogenic compounds elicit ER-mediated responses in the liver [9, 10]. We have previously reported that o, p'-DDT elicits primarily PXR/CAR-mediated, rather than ER-mediated hepatic responses, and suggested that CAR-mediated effects, as opposed to ER-mediated effects, may be more important in DDT-induced tumor promotion in the rat liver .
Although the rat is a preferred toxicological model because of the accumulated background knowledge regarding chemical-induced toxicity compared to mice, several studies report significant species-specific differences in response to chemical exposures. For example, 1,4-bis-[2-(3,5,-dichloropyridyloxy)] benzene (TCPOBOP) acts as potent phenobarbital-type enzyme inducer in the mouse liver but not in the rat or human. This is due to a substitution of Thr350 in mice with Met in the rat and human CAR [12–14]. Alternatively, another phenobarbital-type enzyme inducer 2,4,6-triphenyldioxane-1,3 induces hepatic CYP2B in rats but not in mice . Such differences may affect CAR-mediated hepatic drug metabolism and disposition following xenobiotic exposure. In addition, rat genome annotation lags behind the mouse and human in maturity, confounding a more comprehensive assessment of rat gene expression data. Furthermore, targeted gene-disruption can be used in mice to further elucidate the roles of PXR, CAR and ER in o, p'-DDT elicited effects. Thus, we were motivated to investigate hepatic gene expression in the mouse following o, p'-DDT treatment to identify species-specific and conserved responses. Overall, the hepatic gene expression profiles were comparable, however there were marked species-specific differences in genes regulated by ER, which may involve species-difference in CAR activation. Furthermore, blood levels of DHEA-S, a precursor of androgen and estrogen, exhibited mouse-specific elevation, which may also contribute to endocrine system disruption. Thus, the present study further elucidates the dynamics of o, p'-DDT elicited responses in the rodent, and identifies species-specific responses that may be important to o, p'-DDT exposure for humans and wildlife.
Female C57BL/6 mice, ovariectomized on postnatal day 20 were obtained from Charles River Laboratories (Raleigh, NC) on day 25. Mice were housed in polycarbonate cages containing cellulose fiber 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 h 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.
Treatments and Necropsy
Mice (n = 5) were orally gavaged once or once daily for three consecutive days with 300 mg/kg b.w. o, p'-DDT (99.2% purity, Sigma-Aldrich, St Louis, MO) in 0.1 ml of sesame oil (Sigma-Aldrich) vehicle. An equal number of time-matched vehicle control animals (n = 5) were also treated in the same manner. Mice receiving one dose were sacrificed 2, 4, 8, 12, 18, and 24 h after treatment. Mice receiving three daily doses were sacrificed 24 h after the third treatment (72 h). 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. Whole liver weights were recorded and sections of the left lateral lobe (approximately 0.1 g) were snap-frozen in liquid nitrogen and stored at -80°C. The right lateral lobe was placed in 10% neutral buffered formalin (NBF, VWR International, West Chester, PA) for histopathology and stored at room temperature.
Following fixation of the right lateral lobe for at least 24 h in 10% NBF, the samples were embedded in paraffin according to standard techniques. Five μm sections were mounted on glass slides and stained with hematoxylin and eosin. All embedding, mounting and staining of tissues were performed at the Histology Laboratory, (Department of Physiology, Michigan State University). The histopathology of each liver section was scored according to the NTP Pathology guidelines.
Measurement of blood DHEA-S and androstenedione
Rat plasma samples were obtained from a previous study . Mouse serum and rat plasma samples were used for determining dehydroepiandrosterone sulfate (DHEA-S) and androstenedione levels using a DHEA-S ELISA kit (Calbiotech Inc., Spring Valley, CA) and Androstenedione ELISA kit (Genway Biotech Inc., San Diego, CA), respectively. The detection limits, defined by the vendors, are 0.02 μg/ml of DHEA-S and 0.043 ng/ml of androstenedione, respectively.
Protein preparation and Western analysis
Total protein was extracted from liver samples from all the time points using the Total Protein Extraction kit (Chemicon International Inc., Temecula, CA). Protein concentrations were determined using Bicinchoninic Acid Protein Assay Kit (Sigma-Aldrich). Total protein (100 μg) was resolved on a 11% denaturing SDS-polyacrylamide gel, and transferred to Hybond-ECL membrane (GE Healthcare, Waukesha, WI). Blots were incubated with blocking buffer (0.1% Tris-buffered saline pH 7.4 containing 1% low-fat dry milk) for 10 min at room temperature. Goat anti-CYP17A1, goat anti-actin antibodies and donkey horseradish peroxidase-conjugated anti-goat IgG were purchased from Santa Cruz Biotechnology Inc. (Santa Cruz, CA). Immunochemical staining was performed as described previously , with dilution of 1:500 (anti-CYP17A1), 1:500 (anti-actin) and 1:10000 (anti-goat IgG) using blocking buffer, respectively. SuperSignal West Dura substrate (Thermo Fisher Scientific, Inc., Waltham, MA) was used for signal detection. The Western analysis was performed on three independent biological replicates.
Total RNA was isolated from left lateral liver sections using TRIZOL Reagent (Invitrogen, Carlsbad, CA) and resuspended in The RNA Storage Solution (Ambion, Austin, TX). RNA concentrations were determined by spectrophotometry (A260) and purity was assessed by the A260:A280 ratio and by visual inspection of 3 μg on a denaturing gel.
Temporal changes in gene expression were assessed using an independent reference design in which o, p'-DDT-treated samples were co-hybridized with time-matched vehicle controls using 3 biological replicates and 2 independent labelings of each sample (i.e. dye swap) for each time point. Whole Mouse Genome 4 × 44 K Oligo Microarray Kit (Agilent Technologies, Inc, Santa Clara, CA) was used for global gene expression analysis. All the reagents and enzymes were provided by Agilent Technologies, and the microarray analysis was performed according to the vendor's protocol. The microarrays were scanned at 635 nm (Cy5) and 532 nm (Cy3) using 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). All data were managed in the toxicogenomic information management system dbZach relational database .
Microarray data normalization and statistical analysis
Data were normalized using a semi-parametric approach . 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 (p1 [t]-value) of activity on a per gene and time-point basis using the model-based t-value . Genes were filtered for activity based on the p1(t)-value. p1(t) values approaching one indicate changes in gene expression which are more robust. In this study, unique genes with a p1(t) > 0.999 for a minimum of two time points and absolute value fold change ≥ 1.5-fold compared to time-matched vehicle control for at least one time point were considered differentially expressed.
A total of 996 orthologous rat and mouse genes were in common on the rat cDNA, mouse cDNA, and Agilent mouse oligonucleotide microarrays based on HomoloGene http://www.ncbi.nlm.nih.gov/sites/entrez?db=homologene. Differentially expressed orthologs (|fold change| ≥ 1.5 for at least one time point in either species) for ethynylestradiol (EE)-treated mouse liver, o, p'-DDT treated rat liver and o, p'-DDT-treated mouse liver, were hierarchical clustered using Cluster 3.0 and TreeView software http://rana.lbl.gov/eisen/. In addition, differentially expressed genes known to be regulated by PXR, CAR and ER agonists [20–23], were identified from the literature and their expression profiles were subjected to hierarchical clustering.
Correlation analysis was performed using relaxed criteria to allow for more inclusive comparison between mouse and rat data sets. Genes with a p1(t) > 0.99 and absolute fold change ≥ 1.5-fold at one or more time points in the o, p'-DDT-treated mouse liver samples were selected and used for correlation analysis. This analysis involved a multivariate correlation-based visualization application  that has been used to compare two independent gene expression profile datasets [11, 25]. It calculates the temporal correlations between gene expression (fold change) and significance values (p1 [t]-value) for orthologous o, p'-DDT-treated mouse (this study) and o, p'-DDT-treated rat genes , and summarizes the results in a scatter plot.
Quantitative Real-Time PCR (QRT-PCR)
Expression levels of mouse Car (Nr1i3), Cyp17a1, Cyp2b10, Cyp3a11, Cyp7b1, Gadd45a, Gadd45b, Gapdh, Gclm, Pxr (Nr1i2), and Srebf1, and rat Gapdh and Cyp7b1 were measured by QRT-PCR. For each sample, 2 μ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 , 3 mM MgCl2, 1.0 mM dNTPs, 0.025 IU AmpliTaq Gold and 1 × SYBR Green PCR buffer (Applied Biosystems, Foster City, CA). PCR amplification was conducted in MicroAmp Optical 96-well reaction plates (Applied Biosystems) on an Applied Biosystems PRISM 7000 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 eight orders of magnitude to interpolate relative template concentrations of the samples from the standard curves of log copy number versus threshold cycle. The copy number of each unknown sample for each gene was standardized to that of Gapdh gene to control for differences in RNA loading, quality and cDNA synthesis. Primer sequences and amplicon sizes are provided as Additional file 1.
Body weight, relative liver weight, hepatic concentrations of o, p'-DDT, o, p'-DDD, and o, p'-DDE and QRT-PCR data are presented as the mean ± SE. Statistical analysis was performed with two-way ANOVA followed by pairwise comparisons using Tukey's Honestly Significant Difference post hoc test to control Type I error (α = 0.05). For QRT-PCR data, the relative expression levels of target genes were scaled such that the standardized expression level of the time-matched vehicle control group was equal to 1 for graphing purposes. All statistics were performed using SAS 9.1.3 software (SAS Institute Inc., Cary, NC).
Body Weight, Relative Liver Weight, and Histopathology
Relative liver weight was significantly (p < 0.05) increased at 72 h compared to vehicle-treated mice (Fig. 1A). No effect on body weight was observed. These results are similar to those previously observed in the rat o, p'-DDT study . There were no significant histological changes by o, p'-DDT treatment, except for two animals exhibiting cell death of small numbers of hepatocytes at 4 h or 72 h, where the nuclei were shrunken and hyperchromatic.
Hepatic o, p'-DDT, o, p'-DDD and o, p'-DDE concentration
Animals treated with o, p'-DDT showed significantly greater levels of o, p'-DDT, from 2 to 18 h compared to vehicle controls (Fig. 1B). The levels of the metabolites o, p'-DDD and o, p'-DDE were significantly greater than controls from 2 to 12 h, and at 72 h (Figs. 1B and 1C), indicating that o, p'-DDT is readily metabolized in the liver, and that much of the o, p'-DDT and its metabolites remain primarily distributed to the liver from 2 to 12 h, similar to the rat . Note that there was a statistical outlier in the 72 h treated group, which lead to significant o, p'-DDE and o, p'-DDD levels, but was still included in the analysis as there was no technical or biological justification for its removal.
Gene expression was assessed using the 4 × 44 K Mouse Agilent array containing approximately 44,000 oligonucleotide probes, representing 34,204 annotated genes including approximately 21,000 unique genes. Model-based t-values that compared treated and vehicle responses on a per time-point basis followed by Empirical Bayes analysis identified 1,206 differentially expressed genes across the time course based on a p1(t) > 0.999 for a minimum of two time points and absolute fold change ≥ 1.5-fold compared to time-matched vehicle control for at least one time point (Fig. 2A). Overall, 859 genes were induced while 382 genes were repressed. Changes in gene expression ranged from 36.4-fold induction for Cyp2c55 gene to -13.0-fold repression for gene Car3. The 12 h time point exhibited the greatest number of differentially expressed genes (Fig. 2B). All of the gene expression data as well as the identified differentially regulated genes are provided in Additional files 2 and 3, respectively.
A number of phase I drug metabolizing enzyme genes were induced such as Cyp1a1 (~2.2-fold), Cyp2b9 (~5.0-fold), Cyp2b10 (~18.4-fold), Cyp2b13 (~4.1-fold), Cyp2c39 (~2.1-fold), and Cyp2c55 (~36.4-fold) (Table 1). The induction of Cyp2b10 and Cyp3a11 were verified by QRT-PCR (Fig. 2C). The steroid metabolism genes, Cyp17a1 (~6.2-fold) and Cyp7b1 (~5.8-fold), were also induced, while Cyp51 and Cyp7a1 were down-regulated. Steroidogenic Cyp11a1 mRNA level was not detected by microarray or QRT-PCR (data not shown). Note that the hepatic induction of Cyp17a1 and Cyp7b1 by o, p'-DDT was observed in the mouse but not in the rat . Several phase II genes (Gsta1 (~3.4-fold), Gstm2 (~3.2-fold), Ugt2b35 (~3.4-fold)) were also induced.
A number of DNA damage-related genes such as Cdkn1a (~15.4-fold), Gadd45a (~13.5-fold), Gadd45b (~8.4-fold) and Tnfrsf19 (~5.9-fold) were induced following o, p'-DDT treatment, while Fas and Casp9 were down-regulated (~2.0-fold). The oxidative stress-responsive genes Gclm (~1.6-fold), Gsr (~2.5-fold), Hmox2 (~1.6-fold) and Nqo1 (~2.6-fold) exhibited induction. There was no clear induction for the cell proliferation-related genes such as Ccnd1, Ccnb2, Mdm2 or Stmn1, while their orthologs were induced by o, p'-DDT in the rat liver . mRNA levels for Ahr and Pxr (Nr1i2) were slightly elevated, whereas Car (Nr1i3) mRNA levels decreased following o, p'-DDT treatment. Srebf1 which was previously shown to be down-regulated by o, p'-DDT in the rat liver , was also down-regulated in the mouse liver. Overall, the QRT-PCR results closely paralleled the expression pattern seen in the microarray data (Fig. 2C).
We also examined the expression of genes regulated by CAR, PXR or ER that were selected based on the null animal studies [20–23]. While CAR-specific targets such as Cyp1a1 (~2.2-fold), Fmo5 (~2.3-fold), Sult1d1 (~2.8-fold) and Abcc2 (~2.7-fold) showed relatively weak induction, several PXR-regulated genes were induced more strongly including ApoA4 (~5.8-fold), Ces2 (~3.4-fold), Gstm2 (~3.2-fold) or Insig2 (~5.0-fold) genes (Fig. 2D). However, other PXR genes, such as Hmgcs1 and Hmgcs2, were unexpectedly down-regulated.
Comparative Analysis of Gene Expression
Hepatic gene expression profiles of EE-treated mice, o, p'-DDT-treated rats and o, p'-DDT-treated mice were compared by hierarchical clustering focusing on the 538 orthologs represented on both microarray platforms that satisfied the |fold change| ≥ 1.5 for at least one time point in any data set. The filtering criteria used here were relaxed compared to those used in Figs. 2A and 2B to include more orthologs in the correlation analysis, and thus be more informative of the overall similarity between the two data sets. Despite the differences in microarray platform (i.e., cDNA vs. oligonucleotide), o, p'-DDT elicited gene expression profiles in the rat and mouse that were more similar to each other than to the EE-treated mouse profile (Fig. 3A).
One hundred and forty orthologs were identified as differentially expressed on both rat cDNA and mouse oligo microarrays (Fig. 3B), which were further examined by correlation analysis to determine if their gene expression profiles were comparable. Again, the filtering criteria used here were relaxed compared to those used in Figs. 2A and 2B to include more orthologs in the correlation analysis. The scatter plot shows that a majority of the spots localized to either the upper-right (32.9%) or lower-right (47.9%) quadrants (Fig. 3C), indicating that the gene expression profiles between o, p'-DDT-treated rat and mouse were highly similar in terms of their gene expression pattern. Temporal correlation of the p1(t)-value (statistical significance) showed relatively low correlation between the two data sets, with 60% of spots located in the lower-right or lower-left quadrants. All correlation analyses results are provided in Additional file 4.
For example, Cyp17a1, a key enzyme gene for steroidogenesis (Fig. 4A) fell into the lower left quadrant (Fig. 3C) suggesting divergent gene expression profiles in the mouse and rat. Fig. 4B clearly demonstrates that Cyp17a1 is significantly induced in the mouse and is non-responsive in the rat in both microarray and QRT-PCR data. Cyp7b1, another steroid metabolism gene, exhibited a similar species-specific gene expression profile (Fig 4B). Moreover, hepatic CYP17A1 protein induction was evident at 18 and 24 h in the mouse liver (Fig. 4C).
Effects on DHEA-S Levels
Species-specific induction of CYP17A1 mRNA and protein levels suggests that steroid metabolism may also be differentially affected. More specifically, since CYP17A1 metabolizes pregnenolone and 17α-hydroxypregnenolone to produce DHEA (Fig. 4A), blood levels of DHEA-S were measured using an enzyme-linked immunosorbent assay. Serum DHEA-S levels exhibited a significant (p < 0.05) increase relative to vehicle at 12 h in the mouse, (Fig. 4D). In contrast, plasma DHEA-S was not significantly different between treated and control rats. Mouse serum androstenedione levels were below the detection level (data not shown).
o, p'-DDT elicits phenobarbital-type activities in the mouse liver, characterized by the increase in relative liver weight or induction of Cyp2b10, Cyp3a11 and GST mRNA levels [27, 28]. In the rat, we previously reported that o, p'-DDT elicits PXR/CAR-mediated responses, with negligible ER-mediated activity , despite the responsive of the liver to estrogens [9, 29], and the association of estrogenic activity with hepatic tumorigenesis . Consequently, PXR/CAR-mediated responses to DDT may have a more significant role in hepatic tumor promotion . In order to further investigate o, p'-DDT elicited gene expression mediated by PXR/CAR, we investigated its hepatic expression profile using a comparable mouse model, study design and data analysis . Comparative analysis indicated that the hepatic gene expression profiles were similar between rat and mouse. However, correlation analysis of orthologous genes revealed that hepatic Cyp17a1 mRNA and protein levels as well as serum DHEA-S were only elevated in the mouse following o, p'-DDT treatment.
Cyp17a1 exhibits dose-dependent induction by EE in the mouse liver , suggesting ER-mediated induction. In addition, induction of Cyp7b1 and down-regulation of Cyp7a1 by EE are abolished in ER alpha-null mice but not in ER beta-, Fxr-, Pxr- or Car-null mice . Collectively, these results indicate that Cyp7a1 mRNA levels are regulated by ERα in the mouse liver. Considering the strong inductions of Cyp17a1 and Cyp7b1 observed only in the mouse liver as well as down-regulation of Cyp7a1 in the mouse liver, o, p'-DDT may elicit ERα-mediated activity in the mouse liver but not in the rat liver.
In addition, the mouse profile suggests that o, p'-DDT elicited gene expression was predominately PXR regulated based on genes known to be regulated by PXR-, CAR- or PXR/CAR [20–23]. In part, this could be due to the repression of mouse Car mRNA levels by o, p'-DDT compared to its induction in the rat liver . Furthermore, several genes associated with cell proliferation (Ccnb1, Ccnb2, Mdm2 and Stmn1) were induced in the rat liver , but not in the mouse. Ccnd1, a known CAR-regulated cell proliferation-related gene , was not significantly induced in the mouse, further supporting the hypothesis that o, p'-DDT preferentially activated PXR.
CAR also has an inhibitory effect on ER-mediated gene expression . Therefore, the availability of CAR may inhibit ER-mediated effects elicited by o, p'-DDT in the rat liver , while the down regulation of Car mRNA in the mouse and the preferential activation of PXR may facilitate more ER-mediated gene expression. Consequently, ER-mediated responses may be important when assessing DDT-elicited hepatic responses in the mouse in addition to PXR/CAR-mediated response when compared to the rat.
CYP17A1 is an important enzyme in steroid biosynthesis that metabolizes pregnenolone and progesterone to produce DHEA and androstenedione, respectively [33, 34]. Therefore, induction of CYP17A1 mRNA and protein levels may affect DHEA metabolism. o, p'-DDT treatment increased DHEA-S levels in the mouse whereas no change was detected in the rat. Although the liver is usually not a major steroidogenesis organ, microsomes are capable of participating in steroidogenesis . In addition, hepatic microsomal CYP17A1 exhibits age-dependent expression, with higher expression and activity in immature rats [35, 36]. o, p'-DDT elicited changes in other steroidogenic enzymes (e.g., Cyp11a1) and metabolites (e.g., androstenedione) were below the level of detection. Collectively, alterations in steroid levels following o, p'-DDT treatment could affect steroid levels, since peripheral tissues can use circulating DHEA and DHEA-S to produce androgens and estrogens . Furthermore, the 7α-hydroxylation of DHEA and pregnenolone by o, p'-DDT induced CYP7B1 could also facilitate elimination of DHEA or the synthesis of neuro-active hormones [38, 39].
In addition to PXR-, CAR- or ER-mediated gene expression changes, o, p'-DDT induced Gadd45a, Gadd45b and Cdkn1, all of which are DNA damage-responsive genes [40, 41]. Consistent with this are the reports of the DNA damaging potential of DDT [42–44]. Consequently, DNA damage may be an additional risk factor for tumor initiation/promotion following o, p'-DDT exposure in addition to PXR/CAR- and ER-mediated activities. Considering that the induction of DNA damage-responsive genes precedes Cyp2b10 or Cyp3a11 induction, the DNA damage may not be caused by oxidative stress derived from enzyme induction. Moreover, Gclm and Hmox, both known oxidative stress-responsive genes [45, 46], exhibited relatively weak induction compared to rats , suggesting that oxidative stress was not strongly induced.
In conclusion, o, p'-DDT elicits a broad spectrum of species-conserved and specific effects. This includes PXR/CAR- and ER-mediated responses, altered steroidogenesis, oxidative stress, and DNA damage (Fig. 5). Although DDT is known to cause liver tumors in both mice and rats, the marked species differences in PXR/CAR structure, expression patterns and ligand preference as well as significant species-specific differences in steroidogenesis, especially CYP17A1 expression and activity, confound the extrapolation of these results to humans. Nevertheless, the identification of potential modes of action as well as species-specific responses may assist in the development or selection of more appropriate models for assessing the toxicity of DDT to humans and wildlife.
constitutive androstane receptor (Nr1i3)
cyclin-dependent kinase inhibitor 1A
farnesoid X receptor
glutamate-cysteine ligase modifier subunit
growth arrest and DNA-damage-inducible 45
heme oxygenase (decycling) 2
3-hydroxy-3-methylglutaryl-Coenzyme A synthase
hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase
insulin induced gene 2
transformed mouse 3T3 cell double minute 2
NAD(P)H dehydrogenase, quinone 1
pregnane X receptor (Nr1i2)
Quantitative Real-Time PCR
sterol regulatory element binding factor 1
signal transducer and activator of transcription
tumor necrosis factor receptor superfamily member 19
Ito N, Tsuda H, Hasegawa R, Imaida K: Comparison of the promoting effects of various agents in induction of preneoplastic lesions in rat liver. Environ Health Perspect. 1983, 50: 131-138. 10.2307/3429543.
McGlynn KA, Abnet CC, Zhang M, Sun XD, Fan JH, O'Brien TR, Wei WQ, Ortiz-Conde BA, Dawsey SM, Weber JP: Serum concentrations of 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane (DDT) and 1,1-dichloro-2,2-bis(p-chlorophenyl)ethylene (DDE) and risk of primary liver cancer. J Natl Cancer Inst. 2006, 98: 1005-1010.
Cocco P, Fadda D, Billai B, D'Atri M, Melis M, Blair A: Cancer mortality among men occupationally exposed to dichlorodiphenyltrichloroethane. Cancer Res. 2005, 65: 9588-9594. 10.1158/0008-5472.CAN-05-1487.
DDT and its derivatives. [http://www.inchem.org/documents/ehc/ehc/ehc009.htm]
Shelby MD, Newbold RR, Tully DB, Chae K, Davis VL: Assessing environmental chemicals for estrogenicity using a combination of in vitro and in vivo assays. Environ Health Perspect. 1996, 104 (12): 1296-1300. 10.2307/3432965.
Kojima H, Katsura E, Takeuchi S, Niiyama K, Kobayashi K: Screening for estrogen and androgen receptor activities in 200 pesticides by in vitro reporter gene assays using Chinese hamster ovary cells. Environ Health Perspect. 2004, 112 (5): 524-531.
Kretschmer XC, Baldwin WS: CAR and PXR: xenosensors of endocrine disrupters?. Chem Biol Interact. 2005, 155 (3): 111-128. 10.1016/j.cbi.2005.06.003.
Medina-Diaz IM, Arteaga-Illan G, de Leon MB, Cisneros B, Sierra-Santoyo A, Vega L, Gonzalez FJ, Elizondo G: Pregnane X receptor-dependent induction of the CYP3A4 gene by o, p'-1,1,1,-trichloro-2,2-bis (p-chlorophenyl)ethane. Drug Metab Dispos. 2007, 35 (1): 95-102. 10.1124/dmd.106.011759.
Ciana P, Raviscioni M, Mussi P, Vegeto E, Que I, Parker MG, Lowik C, Maggi A: In vivo imaging of transcriptionally active estrogen receptors. Nat Med. 2003, 9 (1): 82-86. 10.1038/nm809.
Boverhof DR, Fertuck KC, Burgoon LD, Eckel JE, Gennings C, Zacharewski TR: Temporal- and dose-dependent hepatic gene expression changes in immature ovariectomized mice following exposure to ethynyl estradiol. Carcinogenesis. 2004, 25 (7): 1277-1291. 10.1093/carcin/bgh114.
Kiyosawa N, Kwekel JC, Burgoon LD, Williams KJ, Tashiro C, Chittim B, Zacharewski TR: o, p'-DDT Elicits PXR/CAR-, Not ER-, Mediated Responses in the Immature Ovariectomized Rat Liver. Toxicol Sci. 2008, 101 (2): 350-363. 10.1093/toxsci/kfm275.
Poland A, Mak I, Glover E: Species differences in responsiveness to 1,4-bis[2-(3,5-dichloropyridyloxy)]-benzene, a potent phenobarbital-like inducer of microsomal monooxygenase activity. Mol Pharmacol. 1981, 20 (2): 442-450.
Pustylnyak VO, Gulyaeva LF, Lyakhovich VV: Induction of cytochrome P4502B: role of regulatory elements and nuclear receptors. Biochemistry (Mosc). 2007, 72 (6): 608-617. 10.1134/S000629790706003X.
Suino K, Peng L, Reynolds R, Li Y, Cha JY, Repa JJ, Kliewer SA, Xu HE: The nuclear xenobiotic receptor CAR: structural determinants of constitutive activation and heterodimerization. Mol Cell. 2004, 16 (6): 893-905.
Pustylnyak VO, Lebedev AN, Gulyaeva LF, Lyakhovich VV, Slynko NM: Comparative study of CYP2B induction in the liver of rats and mice by different compounds. Life Sci. 2007, 80 (4): 324-328. 10.1016/j.lfs.2006.09.015.
Williams CE, Crawford RB, Holsapple MP, Kaminski NE: Identification of functional aryl hydrocarbon receptor and aryl hydrocarbon receptor nuclear translocator in murine splenocytes. Biochem Pharmacol. 1996, 52 (5): 771-780. 10.1016/0006-2952(96)00360-7.
Burgoon LD, Zacharewski TR: dbZach toxicogenomic information management system. Pharmacogenomics. 2007, 8 (3): 287-291. 10.2217/14622418.104.22.1687.
Eckel JE, Gennings C, Therneau TM, Burgoon LD, Boverhof DR, Zacharewski TR: Normalization of two-channel microarray experiments: a semiparametric approach. Bioinformatics. 2005, 21 (7): 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 (3): 647-670. 10.1081/BIP-200025656.
Ueda A, Hamadeh HK, Webb HK, Yamamoto Y, Sueyoshi T, Afshari CA, Lehmann JM, Negishi M: Diverse roles of the nuclear orphan receptor CAR in regulating hepatic genes in response to phenobarbital. Mol Pharmacol. 2002, 61 (1): 1-6. 10.1124/mol.61.1.1.
Rosenfeld JM, Vargas R, Xie W, Evans RM: Genetic profiling defines the xenobiotic gene network controlled by the nuclear receptor pregnane X receptor. Mol Endocrinol. 2003, 17 (7): 1268-1282. 10.1210/me.2002-0421.
Maglich JM, Stoltz CM, Goodwin B, Hawkins-Brown D, Moore JT, Kliewer SA: Nuclear pregnane x receptor and constitutive androstane receptor regulate overlapping but distinct sets of genes involved in xenobiotic detoxification. Mol Pharmacol. 2002, 62 (3): 638-646. 10.1124/mol.62.3.638.
Yamamoto Y, Moore R, Hess HA, Guo GL, Gonzalez FJ, Korach KS, Maronpot RR, Negishi M: Estrogen receptor alpha mediates 17alpha-ethynylestradiol causing hepatotoxicity. J Biol Chem. 2006, 281 (24): 16625-16631. 10.1074/jbc.M602723200.
Burgoon LD, Boutros PC, Dere E, Zacharewski TR: dbZach: A MIAME-compliant toxicogenomic supportive relational database. Toxicol Sci. 2006, 90 (2): 558-568. 10.1093/toxsci/kfj097.
Boverhof DR, Kwekel JC, Humes DG, Burgoon LD, Zacharewski TR: Dioxin induces an estrogen-like, estrogen receptor-dependent gene expression response in the murine uterus. Mol Pharmacol. 2006, 69 (5): 1599-1606. 10.1124/mol.105.019638.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
Sueyoshi T, Kawamoto T, Zelko I, Honkakoski P, Negishi M: The repressed nuclear receptor CAR responds to phenobarbital in activating the human CYP2B6 gene. J Biol Chem. 1999, 274 (10): 6043-6046. 10.1074/jbc.274.10.6043.
Wagner M, Halilbasic E, Marschall HU, Zollner G, Fickert P, Langner C, Zatloukal K, Denk H, Trauner M: CAR and PXR agonists stimulate hepatic bile acid and bilirubin detoxification and elimination pathways in mice. Hepatology. 2005, 42 (2): 420-430. 10.1002/hep.20784.
Nagel SC, Hagelbarger JL, McDonnell DP: Development of an ER action indicator mouse for the study of estrogens, selective ER modulators (SERMs), and Xenobiotics. Endocrinology. 2001, 142 (11): 4721-4728. 10.1210/en.142.11.4721.
Giannitrapani L, Soresi M, La Spada E, Cervello M, D'Alessandro N, Montalto G: Sex hormones and risk of liver tumor. Ann N Y Acad Sci. 2006, 1089: 228-236. 10.1196/annals.1386.044.
Columbano A, Ledda-Columbano GM, Pibiri M, Cossu C, Menegazzi M, Moore DD, Huang W, Tian J, Locker J: Gadd45beta is induced through a CAR-dependent, TNF-independent pathway in murine liver hyperplasia. Hepatology. 2005, 42 (5): 1118-1126. 10.1002/hep.20883.
Min G, Kim H, Bae Y, Petz L, Kemper JK: Inhibitory cross-talk between estrogen receptor (ER) and constitutively activated androstane receptor (CAR). CAR inhibits ER-mediated signaling pathway by squelching p160 coactivators. J Biol Chem. 2002, 277 (37): 34626-34633. 10.1074/jbc.M205239200.
Liu Y, Yao ZX, Papadopoulos V: Cytochrome P450 17alpha hydroxylase/17,20 lyase (CYP17) function in cholesterol biosynthesis: identification of squalene monooxygenase (epoxidase) activity associated with CYP17 in Leydig cells. Mol Endocrinol. 2005, 19 (7): 1918-1931. 10.1210/me.2004-0271.
Brock BJ, Waterman MR: Biochemical differences between rat and human cytochrome P450c17 support the different steroidogenic needs of these two species. Biochemistry. 1999, 38 (5): 1598-1606. 10.1021/bi9821059.
Katagiri M, Tatsuta K, Imaoka S, Funae Y, Honma K, Matsuo N, Yokoi H, Ishimura K, Ishibashi F, Kagawa N: Evidence that immature rat liver is capable of participating in steroidogenesis by expressing 17alpha-hydroxylase/17,20-lyase P450c17. J Steroid Biochem Mol Biol. 1998, 64: 1-2. 10.1016/S0960-0760(97)00164-7.
Imaoka S, Kamataki T, Funae Y: Purification and characterization of six cytochromes P-450 from hepatic microsomes of immature female rats. J Biochem. 1987, 102 (4): 843-851.
Labrie F, Luu-The V, Labrie C, Simard J: DHEA and its transformation into androgens and estrogens in peripheral target tissues: intracrinology. Front Neuroendocrinol. 2001, 22 (3): 185-212. 10.1006/frne.2001.0216.
Muller C, Hennebert O, Morfin R: The native anti-glucocorticoid paradigm. J Steroid Biochem Mol Biol. 2006, 100 (1-3): 95-105. 10.1016/j.jsbmb.2006.03.001.
Tang W, Eggertsen G, Chiang JY, Norlin M: Estrogen-mediated regulation of CYP7B1: a possible role for controlling DHEA levels in human tissues. J Steroid Biochem Mol Biol. 2006, 100 (1–3): 42-51. 10.1016/j.jsbmb.2006.02.005.
Taylor WR, Stark GR: Regulation of the G2/M transition by p53. Oncogene. 2001, 20 (15): 1803-1815. 10.1038/sj.onc.1204252.
Sheikh MS, Hollander MC, Fornance AJ: Role of Gadd45 in apoptosis. Biochem Pharmacol. 2000, 59 (1): 43-45. 10.1016/S0006-2952(99)00291-9.
Hassoun E, Bagchi M, Bagchi D, Stohs SJ: Comparative studies on lipid peroxidation and DNA-single strand breaks induced by lindane, DDT, chlordane and endrin in rats. Comp Biochem Physiol C. 1993, 104 (3): 427-431. 10.1016/0742-8413(93)90013-B.
Yanez L, Borja-Aburto VH, Rojas E, de la Fuente H, Gonzalez-Amaro R, Gomez H, Jongitud AA, Diaz-Barriga F: DDT induces DNA damage in blood cells. Studies in vitro and in women chronically exposed to this insecticide. Environ Res. 2004, 94 (1): 18-24. 10.1016/S0013-9351(03)00047-1.
Harada T, Yamaguchi S, Ohtsuka R, Takeda M, Fujisawa H, Yoshida T, Enomoto A, Chiba Y, Fukumori J, Kojima S: Mechanisms of promotion and progression of preneoplastic lesions in hepatocarcinogenesis by DDT in F344 rats. Toxicol Pathol. 2003, 31 (1): 87-98.
Bauer M, Bauer I: Heme oxygenase-1: redox regulation and role in the hepatic response to oxidative stress. Antioxid Redox Signal. 2002, 4 (5): 749-758. 10.1089/152308602760598891.
Kiyosawa N, Ito K, Sakuma K, Niino N, Kanbori M, Yamoto T, Manabe S, Matsunuma N: Evaluation of glutathione deficiency in rat livers by microarray analysis. Biochem Pharmacol. 2004, 68 (7): 1465-1475. 10.1016/j.bcp.2004.05.053.
The authors thank Anna Kopec for helpful discussion and critical review of this manuscript. This study was supported by National Institute of General Medical Sciences (GM075838), U.S. Environmental Protection Agency (RD83184701). T.R.Z. is partially supported by the Michigan Agricultural Experiment Station
NK performed all microarray, ELISA and Western assays, conducted and statistically analyzed QRT-PCR assays, as well as compiled and interpreted the data collected and generated the primary draft manuscript. JCK designed/conducted animal studies and participated in interpretation of microarray data. LDB provided database support for the microarray data including quality control assessments and statistical analysis. ED participated in microarray data analysis. KJW carried out histology examination. CT and BC carried out tissue level analysis. TRZ conceived the study and its design and supervised its completion. All authors read and approved the final draft of the manuscript.
Electronic supplementary material
Additional File 2: o, p' -DDT-treated mouse liver time-course study, full data set. Blank gene expression ratios or p1(t) values are due to lack of or inapproprate statistical data for computation. (XLS 11 MB)
Additional File 3: Differentially regulated genes in the o , p '-DDT-treated mouse liver.a Highligted values indicate |fold change| > 1.5. b Highlighted values indicate p1(t) > 0.999. (XLS 412 KB)
About this article
Cite this article
Kiyosawa, N., Kwekel, J.C., Burgoon, L.D. et al. Species-specific regulation of PXR/CAR/ER-target genes in the mouse and rat liver elicited by o, p'-DDT. BMC Genomics 9, 487 (2008). https://doi.org/10.1186/1471-2164-9-487
- Constitutive Androstane Receptor
- Relative Liver Weight
- Constitutive Androstane Receptor Activation
- Human Constitutive Androstane Receptor