Agnes L. Forgacs

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.

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.

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.

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 temporal and dose-dependent morphological and transcriptional uterine effects elicited by tamoxifen and ethynylestradiol in immature, ovariectomized mice

BackgroundUterine 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 17α-ethynylestradiol (EE) to identify pathways similarly and uniquely modified by each compound.ResultsThree daily doses of 100 μg/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 μg/kg TAM. cDNA microarray analysis identified 2235 differentially expressed genes following a single dose of 100 μg/kg TAM at 2, 4, 8, 12, 18 and 24 hrs, and at 72 hrs after three daily doses (3 × 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.ConclusionOverall, 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.

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 EC₅₀s 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.

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.

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.

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

BackgroundTamoxifen (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. TAMs 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.ResultsA 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.ConclusionComparative 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.

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.