o-p'-DDT-mediated uterotrophy and gene expression in immature C57BL/6 mice and Sprague-Dawley rats.

1,1,1-Trichloro-2,2-bis(2-chlorophenyl-4-chlorophenyl)ethane (o,p'-DDT) is an organochlorine pesticide and endocrine disruptor known to activate the estrogen receptor. Comprehensive ligand- and species-comparative dose- and time-dependent studies were conducted to systematically assess the uterine physiological, morphological and gene expression responses elicited by o,p'-DDT and ethynyl estradiol (EE) in immature ovariectomized C57BL/6 mice and Sprague-Dawley rats. Custom cDNA microarrays were used to identify conserved and divergent differential gene expression responses. A total of 1256 genes were differentially expressed by both ligands in both species, 559 of which exhibited similar temporal expression profiles suggesting that o,p'-DDT elicits estrogenic effects at high doses when compared to EE. However, 51 genes exhibited species-specific uterine expression elicited by o,p'-DDT. For example, carbonic anhydrase 2 exhibited species- and ligand-divergent expression as confirmed by quantitative real-time PCR. The identification of comparable temporal phenotypic responses linked to gene expression demonstrates that systematic comparative gene expression assessments are valuable for elucidating conserved and divergent estrogen signaling mechanisms in rodent uterotrophy.

Comparisons of differential gene expression elicited by TCDD, PCB126, ßNF, or ICZ in mouse hepatoma Hepa1c1c7 cells and C57BL/6 mouse liver.

The aryl hydrocarbon receptor (AhR) is a promiscuous receptor activated by structurally diverse synthetic and natural compounds. AhR activation may lead to ligand-specific changes in gene expression despite similarities in mode of action. Therefore, differential gene expression elicited by four structurally diverse, high affinity AhR ligands (2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; 10nM, 30 µg/kg), 3,3',4,4',5-pentachlorobiphenyl (PCB126; 100nM, 300µg/kg), ß-naphthoflavone (ßNF; 10 µM, 90 mg/kg), and indolo[3,2-b]carbazole (ICZ; 1µM)) in mouse Hepa1c1c7 hepatoma cells and C57BL/6 mouse liver samples were compared. A total of 288, 183, 119, and 131 Hepa1c1c7 genes were differentially expressed (|fold-change|≥ 1.5, P1(t)≥ 0.9999) by TCDD, ßNF, PCB126, and ICZ, respectively. Only ∼35% were differentially expressed by all 4 ligands in Hepa1c1c7 cells. In vivo, 661, 479, and 265 hepatic genes were differentially expressed following treatment with TCDD, ßNF, and PCB126, respectively. Similar to Hepa1c1c7 cells, ≤ 34% of gene expression changes were common across all ligands. Principal components analysis identified time-dependent gene expression divergence. Comparisons of ligand-elicited expression between Hepa1c1c7 cells and mouse liver identified only 11 common gene expression changes across all ligands. Although metabolism may explain some ligand-specific gene expression changes, PCB126, ßNF, and ICZ also elicited divergent expression compared to TCDD, suggestive of selective AhR modulation.

Triazine herbicides and their chlorometabolites alter steroidogenesis in BLTK1 murine Leydig cells.

The triazine herbicides, atrazine (ATR), simazine (SIM), propazine (PRO), terbuthylazine (TBA), and their chlorinated metabolites have been implicated in the etiology of testicular dysgenesis by altering steroidogenesis. To further investigate their effects on testosterone biosynthesis, BLTK1 cells were used to evaluate steroid hormone levels and genome-wide gene expression. BLTK1 cells are a novel murine Leydig cell line possessing an intact steroidogenic pathway with constitutive low basal testosterone (T) levels that can be induced by recombinant human chorionic gonadotropin (rhCG). Triazines (ATR, SIM, PRO, and TBA) and their chlorometabolites (DEA, DIA, and DACT) induced concentration-dependent (1, 3, 10, 30, 100, 300, and 600 µM) increases in progesterone (P) and T levels relative to solvent control at 24h. Temporal analysis (300 µM at 1, 2, 4, 8, 12, 24, or 48h) elicited comparable P and T profiles by all compounds with varying efficacies (ATR > TBA > PRO > DEA > DIA > DACT >> SIM) that were similar to rhCG. ATR and TBA elicited time- and concentration-dependent induction of Star, Hsd3b6, and Hsd17b3 mRNA levels, whereas Hsd3b1, Cyp17a1, and Srd5a1 mRNA expression was repressed. PRO elicited similar albeit weaker effects, whereas SIM had negligible effects consistent with their induction of P and T levels. Whole-genome microarrays identified 797 differentially regulated genes elicited by 300 µM ATR, occurring primarily at later time points (> 12h) with overrepresented functions associated with steroidogenesis and cholesterol metabolism. These results indicate that changes in P and T levels can be partially attributed to triazine-elicited alterations in steroidogenic gene expression.

Comparative analysis of temporal and dose-dependent TCDD-elicited gene expression in human, mouse, and rat primary hepatocytes.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD)-elicited time- and dose-dependent differential gene expression was compared in human, mouse, and rat primary hepatocytes. Comprehensive time course (10 nM TCDD or dimethyl sulfoxide vehicle control for 1, 2, 4, 8, 12, 24, and 48h) studies identified 495, 2305, and 711 differentially expressed orthologous genes in human, mouse, and rat hepatocytes, respectively. However, only 16 orthologs were differentially expressed across all three species, with the majority of orthologs exhibiting species-specific expression (399 human, 2097 mouse, and 533 rat), consistent with species-specific expression reported in other in vitro and in vivo comparative studies. TCDD also elicited the dose-dependent induction of 397 human, 100 mouse, and 443 rat genes at 12h and 615 human, 426 mouse, and 314 rat genes at 24h. Comparable EC50 values were obtained for AhR battery genes including Cyp1a1 (0.1 nM human, 0.05 nM mouse, 0.08 nM rat at 24h) and Tiparp (0.97 nM human, 0.63 nM mouse, 0.14 nM rat at 12h). Overrepresented functions and pathways included amino acid metabolism in humans, immune response in mice, and energy homeostasis in rats. Differentially expressed genes functionally associated with lipid transport, processing, and metabolism were overrepresented in all three species but exhibited species-specific expression consistent with the induction of hepatic steatosis in mice but not in rats following a single oral gavage of TCDD. Furthermore, human primary hepatocytes showed lipid accumulation following 48h of treatment with TCDD, suggesting that AhR-mediated steatosis in mice more closely resembles human hepatic fat accumulation compared with that in rats. Collectively, these results suggest that species-specific gene expression profiles mediate the species-specific effects of TCDD despite the conservation of the AhR and its signaling mechanism.

Comparative toxicogenomic analysis of oral Cr(VI) exposure effects in rat and mouse small intestinal epithelia.

Continuous exposure to high concentrations of hexavalent chromium [Cr(VI)] in drinking water results in intestinal tumors in mice but not rats. Concentration-dependent gene expression effects were evaluated in female F344 rat duodenal and jejunal epithelia following 7 and 90 days of exposure to 0.3-520 mg/L (as sodium dichromate dihydrate, SDD) in drinking water. Whole-genome microarrays identified 3269 and 1815 duodenal, and 4557 and 1534 jejunal differentially expressed genes at 8 and 91 days, respectively, with significant overlaps between the intestinal segments. Functional annotation identified gene expression changes associated with oxidative stress, cell cycle, cell death, and immune response that were consistent with reported changes in redox status and histopathology. Comparative analysis with B6C3F1 mouse data from a similarly designed study identified 2790 differentially expressed rat orthologs in the duodenum compared to 5013 mouse orthologs at day 8, and only 1504 rat and 3484 mouse orthologs at day 91. Automated dose-response modeling resulted in similar median EC50s in the rodent duodenal and jejunal mucosae. Comparative examination of differentially expressed genes also identified divergently regulated orthologs. Comparable numbers of differentially expressed genes were observed at equivalent Cr concentrations (µg Cr/g duodenum). However, mice accumulated higher Cr levels than rats at ≥ 170 mg/L SDD, resulting in a ~2-fold increase in the number of differentially expressed genes. These qualitative and quantitative differences in differential gene expression, which correlate with differences in tissue dose, likely contribute to the disparate intestinal tumor outcomes.

BLTK1 murine Leydig cells: a novel steroidogenic model for evaluating the effects of reproductive and developmental toxicants.

Leydig cells are the primary site of androgen biosynthesis in males. Several environmental toxicants target steroidogenesis resulting in both developmental and reproductive effects including testicular dysgenesis syndrome. The aim of this study was to evaluate the effect of several structurally diverse endocrine disrupting compounds (EDCs) on steroidogenesis in a novel BLTK1 murine Leydig cell model. We demonstrate that BLTK1 cells possess a fully functional steroidogenic pathway that produces low basal levels of testosterone (T) and express all the necessary steroidogenic enzymes including Star, Cyp11a1, Cyp17a1, Hsd3b1, Hsd17b3, and Srd5a1. Recombinant human chorionic gonadotropin (rhCG) and forskolin (FSK) elicited concentration- and time-dependent induction of 3',5'-cyclic adenosine monophosphate, progesterone (P), and T, as well as the differential expression of Star, Hsd3b6, Hsd17b3, and Srd5a1 messenger RNA levels. The evaluation of several structurally diverse male reproductive toxicants including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), atrazine, prochloraz, triclosan, monoethylhexyl phthalate (MEHP), glyphosate, and RDX in BLTK1 cells suggests different modes of action perturb steroidogenesis. For example, prochloraz and triclosan antifungals reduced rhCG induction of T, consistent with published in vivo data but did not alter basal T levels. In contrast, atrazine and MEHP elicited modest induction of basal T but antagonized rhCG-mediated induction of T levels, whereas TCDD, glyphosate, and RDX had no effect on basal or rhCG induction of T in BLTK1 cells. These results suggest that BLTK1 cells maintain rhCG-inducible steroidogenesis and are a viable in vitro Leydig cell model to evaluate the effects of EDCs on steroidogenesis. This model can also be used to elucidate the different mechanisms underlying toxicant-mediated disruption of steroidogenesis.

Genome-wide gene expression effects in B6C3F1 mouse intestinal epithelia following 7 and 90days of exposure to hexavalent chromium in drinking water.

Chronic administration of high doses of hexavalent chromium [Cr(VI)] as sodium dichromate dihydrate (SDD) elicits alimentary cancers in mice. To further elucidate key events underlying tumor formation, a 90-day drinking water study was conducted in B6C3F1 mice. Differential gene expression was examined in duodenal and jejunal epithelial samples following 7 or 90days of exposure to 0, 0.3, 4, 14, 60, 170 or 520mg/L SDD in drinking water. Genome-wide microarray analyses identified 6562 duodenal and 4448 jejunal unique differentially expressed genes at day 8, and 4630 and 4845 unique changes, respectively, in the duodenum and jejunum at day 91. Comparative analysis identified significant overlap in duodenal and jejunal differential gene expression. Automated dose-response modeling identified >80% of the differentially expressed genes exhibited sigmoidal dose-response curves with EC(50) values ranging from 10 to 100mg/L SDD. Only 16 genes satisfying the dose-dependent differential expression criteria had EC(50) values <10mg/L SDD, 3 of which were regulated by Nrf2, suggesting oxidative stress in response to SDD at low concentrations. Analyses of differentially expressed genes identified over-represented functions associated with oxidative stress, cell cycle, lipid metabolism, and immune responses consistent with the reported effects on redox status and histopathology at corresponding SDD drinking water concentrations. Collectively, these data are consistent with a mode of action involving oxidative stress and cytotoxicity as early key events. This suggests that the tumorigenic effects of chronic Cr(VI) oral exposure likely require chronic tissue damage and compensatory epithelial cell proliferation.

Comparative metabolomic and genomic analyses of TCDD-elicited metabolic disruption in mouse and rat liver.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) elicits a broad spectrum of species-specific effects that have not yet been fully characterized. This study compares the temporal effects of TCDD on hepatic aqueous and lipid metabolite extracts from immature ovariectomized C57BL/6 mice and Sprague-Dawley rats using gas chromatography-mass spectrometry and nuclear magnetic resonance-based metabolomic approaches and integrates published gene expression data to identify species-specific pathways affected by treatment. TCDD elicited metabolite and gene expression changes associated with lipid metabolism and transport, choline metabolism, bile acid metabolism, glycolysis, and glycerophospholipid metabolism. Lipid metabolism is altered in mice resulting in increased hepatic triacylglycerol as well as mono- and polyunsaturated fatty acid (FA) levels. Mouse-specific changes included the induction of CD36 and other cell surface receptors as well as lipases- and FA-binding proteins consistent with hepatic triglyceride and FA accumulation. In contrast, there was minimal hepatic fat accumulation in rats and decreased CD36 expression. However, choline metabolism was altered in rats, as indicated by decreases in betaine and increases in phosphocholine with the concomitant induction of betaine-homocysteine methyltransferase and choline kinase gene expression. Results from these studies show that aryl hydrocarbon receptor-mediated differential gene expression could be linked to metabolite changes and species-specific alterations of biochemical pathways.

Identification of aryl hydrocarbon receptor binding targets in mouse hepatic tissue treated with 2,3,7,8-tetrachlorodibenzo-p-dioxin.

Genome-wide, promoter-focused ChIP-chip analysis of hepatic aryl hydrocarbon receptor (AHR) binding sites was conducted in 8-week old female C57BL/6 treated with 30 µg/kg/body weight 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) for 2 h and 24 h. These studies identified 1642 and 508 AHR-bound regions at 2h and 24h, respectively. A total of 430 AHR-bound regions were common between the two time points, corresponding to 403 unique genes. Comparison with previous AHR ChIP-chip studies in mouse hepatoma cells revealed that only 62 of the putative target genes overlapped with the 2 h AHR-bound regions in vivo. Transcription factor binding site analysis revealed an over-representation of aryl hydrocarbon response elements (AHREs) in AHR-bound regions with 53% (2 h) and 68% (24 h) of them containing at least one AHRE. In addition to AHREs, E2f-Myc activator motifs previously implicated in AHR function, as well as a number of other motifs, including Sp1, nuclear receptor subfamily 2 factor, and early growth response factor motifs were also identified. Expression microarray studies identified 133 unique genes differentially regulated after 4 h treatment with TCDD. Of which, 39 were identified as AHR-bound genes at 2 h. Ingenuity Pathway Analysis on the 39 AHR-bound TCDD responsive genes identified potential perturbation in biological processes such as lipid metabolism, drug metabolism, and endocrine system development as a result of TCDD-mediated AHR activation. Our findings identify direct AHR target genes in vivo, highlight in vitro and in vivo differences in AHR signaling and show that AHR recruitment does not necessarily result in changes in target gene expression.

Genome-wide computational analysis of dioxin response element location and distribution in the human, mouse, and rat genomes.

The aryl hydrocarbon receptor (AhR) mediates responses elicited by 2,3,7,8-tetrachlorodibenzo-p-dioxin by binding to dioxin response elements (DRE) containing the core consensus sequence 5'-GCGTG-3'. The human, mouse, and rat genomes were computationally searched for all DRE cores. Each core was then extended by 7 bp upstream and downstream, and matrix similarity (MS) scores for the resulting 19 bp DRE sequences were calculated using a revised position weight matrix constructed from bona fide functional DREs. In total, 72318 human, 70720 mouse, and 88651 rat high-scoring (MS ≥ 0.8437) putative DREs were identified. Gene-encoding intragenic DNA regions had ∼1.6 times more putative DREs than the noncoding intergenic DNA regions. Furthermore, the promoter region spanning ±1.5 kb of a TSS had the highest density of putative DREs within the genome. Chromosomal analysis found that the putative DRE densities of chromosomes X and Y were significantly lower than the mean chromosomal density. Interestingly, the 10 kb upstream promoter region on chromosome X of the genomes were significantly less dense than the chromosomal mean, while the same region in chromosome Y was the most dense. In addition to providing a detailed genomic map of all DRE cores in the human, mouse, and rat genomes, these data will further aid the elucidation of AhR-mediated signal transduction.

Effects of tamoxifen and ethynylestradiol cotreatment on uterine gene expression in immature, ovariectomized mice.

Tamoxifen (TAM), the primary treatment for estrogen receptor (ER)-positive breast cancer, has been associated with an increased incidence of endometrial cancer in postmenopausal, but not premenopausal women. TAM elicits a partial ER-mediated uterotrophic response in immature rodents when compared with ethynylestradiol (EE), a potent ER agonist. However, cotreatment with 1000 microg/kg TAM antagonizes the uterotrophic effect induced by 30 microg/kg EE. To further investigate the anti-uterotrophic activity of TAM, immature, ovariectomized C57BL/6 mice were treated with a single oral dose of EE, TAM, EE+TAM, or vehicle, and harvested at 2, 4, 8, 12, 18, and 24 h or after three daily treatments at 72 h. Significant increases in uterine wet weight (UWW) were observed at 18 h for EE, TAM, and the mixture. However, mixture induction of UWW was significantly lower when compared with EE-induced uterotrophy at 72 h. This inhibitory effect is also reflected in decreases in luminal circumference, yet EE-induced luminal epithelial cell height was unaffected by cotreatment with TAM. Gene expression analysis using a 2 x 2 factorial cDNA microarray study design identified 2518 differentially expressed genes following EE treatment alone. However, only 290 EE-elicited gene expression changes were affected by TAM cotreatment, in a manner consistent with the anti-estrogenic response. These data suggest that TAM antagonism of EE-induced UWW increase involves the selective inhibition of EE-induced genes.

Effects of TCDD on the expression of nuclear encoded mitochondrial genes.

Generation of mitochondrial reactive oxygen species (ROS) can be perturbed following exposure to environmental chemicals such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Reports indicate that the aryl hydrocarbon receptor (AhR) mediates TCDD-induced sustained hepatic oxidative stress by decreasing hepatic ATP levels and through hyperpolarization of the inner mitochondrial membrane. To further elucidate the effects of TCDD on the mitochondria, high-throughput quantitative real-time PCR (HTP-QRTPCR) was used to evaluate the expression of 90 nuclear genes encoding mitochondrial proteins involved in electron transport, oxidative phosphorylation, uncoupling, and associated chaperones. HTP-QRTPCR analysis of time course (30 microg/kg TCDD at 2, 4, 8, 12, 18, 24, 72, and 168 h) liver samples obtained from orally gavaged immature, ovariectomized C57BL/6 mice identified 54 differentially expressed genes (/fold change/ > 1.5 and P-value < 0.1). Of these, 8 exhibited a sigmoidal or exponential dose-response profile (0.03 to 300 microg/kg TCDD) at 4, 24 or 72 h. Dose-responsive genes encoded proteins associated with electron transport chain (ETC) complexes I (NADH dehydrogenase), III (cytochrome c reductase), IV (cytochrome c oxidase), and V (ATP synthase) and could be generally categorized as having proton gradient, ATP synthesis, and chaperone activities. In contrast, transcript levels of ETC complex II, succinate dehydrogenase, remained unchanged. Putative dioxin response elements were computationally found in the promoter regions of all 8 dose-responsive genes. This high-throughput approach suggests that TCDD alters the expression of genes associated with mitochondrial function which may contribute to TCDD-elicited mitochondrial toxicity.

Tamoxifen-elicited uterotrophy: cross-species and cross-ligand analysis of the gene expression program.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

Comparative temporal and dose-dependent morphological and transcriptional uterine effects elicited by tamoxifen and ethynylestradiol in immature, ovariectomized mice.

BACKGROUND: Uterine temporal and dose-dependent histopathologic, morphometric and gene expression responses to the selective estrogen receptor modulator tamoxifen (TAM) were comprehensively examined to further elucidate its estrogen receptor-mediated effects. These results were systematically compared to the effects elicited by the potent estrogen receptor ligand 17alpha-ethynylestradiol (EE) to identify pathways similarly and uniquely modified by each compound. RESULTS: Three daily doses of 100 microg/kg TAM elicited a dose-dependent increase in uterine wet weight (UWW) in immature, ovariectomized C57BL/6 mice at 72 hrs with concurrent increases in luminal epithelial cell height (LECH), luminal circumference and glandular epithelial tubule number. Significant UWW and LECH increases were detected at 24 hrs after a single dose of 100 microg/kg TAM. cDNA microarray analysis identified 2235 differentially expressed genes following a single dose of 100 microg/kg TAM at 2, 4, 8, 12, 18 and 24 hrs, and at 72 hrs after three daily doses (3 x 24 hrs). Functional annotation of differentially expressed genes was associated with cell growth and proliferation, cytoskeletal organization, extracellular matrix modification, nucleotide synthesis, DNA replication, protein synthesis and turnover, lipid metabolism, glycolysis and immunological responses as is expected from the uterotrophic response. Comparative analysis of TAM and EE treatments identified 1209 common, differentially expressed genes, the majority of which exhibited similar profiles despite a temporal delay in TAM elicited responses. However, several conserved and treatment specific responses were identified that are consistent with proliferation (Fos, Cdkn1a, Anapc1), and water imbibition (Slc30a3, Slc30a5) responses elicited by EE. CONCLUSION: Overall, TAM and EE share similar gene expression profiles. However, TAM responses exhibit lower efficacy, while responses unique to EE are consistent with the physiological differences elicited between compounds.