Hepatic gene downregulation following acute and subchronic exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin

Effect of exposure to 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and polychlorinated biphenyls (PCBs) on mitochondrial DNA (mtDNA) copy number in rats

Mitochondria are intracellular organelles responsible for biological oxidation and energy production. These organelles are susceptible to damage from oxidative stress and compensate for damage by increasing the number of copies of their own genome, mitochondrial DNA (mtDNA). Cancer and environmental exposure to some pollutants have also been associated with altered mtDNA copy number. Since exposures to polychlorinated biphenyls (PCBs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) have been shown to increase oxidative stress, we hypothesize that mtDNA copy number will be altered with exposure to these compounds. mtDNA copy number was measured in DNA from archived frozen liver and lung specimens from the National Toxicology Program (NTP) study of female Harlan Sprague Dawley rats exposed to TCDD (3, 10, or 100 ng/kg/day), dioxin-like (DL) PCB 126 (10, 100, or 1000 ng/kg/day), non-DL PCB 153 (10, 100, or 1000 μg/kg/day), and PCB 126 + PCB 153 (10 ng/kg/day + 10 μg/kg/day, 100 ng/kg/day + 100 μg/kg/day, or 1000 ng/kg/day + 1000 μg/kg/day, respectively) for 13 and 52 weeks. An increase in mtDNA copy number was observed in the liver and lung of rats exposed to TCDD and the lung of rats exposed to the mixture of PCB 126 and PCB 153. A statistically significant positive dose-dependent trend was also observed in the lung of rats exposed to PCB 126 and a mixture of PCB 153 and PCB 126, although in neither case was the control copy number significantly exceeded at any dose level. These exposures produced a range of pathological responses in these organs in the two-year NTP studies. Conversely, there was a significant decrease or no change in mtDNA copy number in the liver and lung of rats exposed to non-DL PCB 153. This is consistent with a general lack of PCB 153 mediated liver or lung injury in the NTP study, with the exception of liver hypertrophy. Together, the results suggest that an increase in mtDNA copy number may serve as a sensitive, early biomarker of mitochondrial injury and oxidative stress that contributes to the development of the toxicity of dioxin-like compounds.

Telomeres as targets for the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and polychlorinated biphenyls (PCBs) in rats

Telomere length (TL) can be affected by various factors, including age and oxidative stress. Changes in TL have been associated with chronic disease, including a higher risk for several types of cancer. Environmental exposure of humans to PCBs and dioxins has been associated with longer or shorter leukocyte TL. Relative telomere length (RTL) may serve as a biomarker associated with neoplastic and/or non-neoplastic responses observed with chronic exposures to TCDD and PCBs. RTL was measured in DNA isolated from archived frozen liver and lung tissues from the National Toxicology Program (NTP) studies conducted in female Harlan Sprague Dawley rats exposed for 13, 30, and 52 weeks to TCDD, dioxin-like (DL) PCB 126, non-DL PCB 153, and a mixture of PCB 126 and PCB 153. RTL was assessed by quantitative polymerase chain reaction (qPCR). Consistent with literature, decreased liver and lung RTL was seen with aging. Relative to time-matched vehicle controls, RTL was increased in both the liver and lung tissues of rats exposed to TCDD, PCB 126, PCB 153, and the mixture of PCB 126 and PCB 153, which is consistent with most epidemiological studies that found PCB exposures were associated with increased leukocyte RTL. Increased RTL was observed at doses and/or time points where little to no pathology was observed. In addition to serving as a biomarker of exposure to these compounds in rats and humans, increases in RTL may be an early indicator of neoplastic and non-neoplastic responses that occur following chronic exposure to TCDD and PCBs.

Identification of Time-Invariant Biomarkers for Non-Genotoxic Hepatocarcinogen Assessment

Non-genotoxic hepatocarcinogens (NGHCs) can only be confirmed by 2-year rodent studies. Toxicogenomics (TGx) approaches using gene expression profiles from short-term animal studies could enable early assessment of NGHCs. However, high variance in the modulation of the genes had been noted among exposure styles and datasets. Expanding from our previous strategy in identifying consensus biomarkers in multiple experiments, we aimed to identify time-invariant biomarkers for NGHCs in short-term exposure styles and validate their applicability to long-term exposure styles. In this study, nine time-invariant biomarkers, namely A2m, Akr7a3, Aqp7, Ca3, Cdc2a, Cdkn3, Cyp2c11, Ntf3, and Sds, were identified from four large-scale microarray datasets. Machine learning techniques were subsequently employed to assess the prediction performance of the biomarkers. The biomarker set along with the Random Forest models gave the highest median area under the receiver operating characteristic curve (AUC) of 0.824 and a low interquartile range (IQR) variance of 0.036 based on a leave-one-out cross-validation. The application of the models to the external validation datasets achieved high AUC values of greater than or equal to 0.857. Enrichment analysis of the biomarkers inferred the involvement of chronic inflammatory diseases such as liver cirrhosis, fibrosis, and hepatocellular carcinoma in NGHCs. The time-invariant biomarkers provided a robust alternative for NGHC prediction.

Assessment of the short-term toxicity of TiO2 nanofiber in Sprague Dawley rats

Synthetic nanomaterials have many unique chemical and physical properties, mainly due to their high specific surface area and quantum confinement effect. Specifically, titanium dioxide (TiO₂ ) nanomaterial has high stability, anticorrosive, and photocatalytic properties. However, there are concerns over adverse biological effects resulting from bioeffects. This study was to investigate adverse effects associated with acute ingestion of TiO₂ nanofiber (TDNF). TDNF was fabricated via electrospinning method, followed by dissolution in water. Six- to seven-week-old male Sprague Dawley rats were exposed to a total of 0, 40, and 60 ppm of TDNF for 2 weeks via oral gavage. Serum total protein and weight gain during the course of this study displayed marginal concentration-dependent alterations. These findings were followed by a global gene expression analysis to identify which transcripts might be responsive to TNDF toxicity. Differentially expressed mRNA levels were dose-dependently higher in animals exposed to TNDF. The majority of the affected genes were biochemically involved in immune response and inflammation. We believe this is due to the fact that TNDF is unable to penetrate the cell and forms phagocytosis sites that trigger inflammatory and immune response. All results taken together, short-term ingestion of TNDF produced marginal effects indicative of inflammation. Finally, the broad gene expression data were validated through quantification of immunoglobulin heavy chain alpha (Igha). Igha gene was upregulated in treated groups, showing similar expression patterns to the global gene expression data.

Transcript variations, phylogenetic tree and chromosomal localization of porcine aryl hydrocarbon receptor (AhR) and AhR nuclear translocator (ARNT) genes

Aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor best known for mediating xenobiotic-induced toxicity. AhR requires aryl hydrocarbon receptor nuclear translocator (ARNT) to form an active transcription complex and promote the activation of genes which have dioxin responsive element in their regulatory regions. The present study was performed to determine the complete cDNA sequences of porcine AhR and ARNT genes and their chromosomal localization. Total RNA from porcine livers were used to obtain the sequence of the entire porcine transcriptome by next-generation sequencing (NGS; lllumina HiSeq2500). In addition, both, in silico analysis and fluorescence in situ hybridization (FISH) were used to determine chromosomal localization of porcine AhR and ARNT genes. In silico analysis of nucleotide sequences showed that there were two transcript variants of AhR and ARNT genes in the pig. In addition, computer analysis revealed that AhR gene in the pig is located on chromosome 9 and ARNT on chromosome 4. The results of FISH experiment confirmed the localization of porcine AhR and ARNT genes. In the present study, for the first time, the full cDNAs of AhR and ARNT were demonstrated in the pig. In future, it would be interesting to determine the tissue distribution of AhR and ARNT transcript variants in the pig and to test whether these variants are associated with different biological functions and/or different activation pathways.

The adverse outcome pathway for rodent liver tumor promotion by sustained activation of the aryl hydrocarbon receptor

An Adverse Outcome Pathway (AOP) represents the existing knowledge of a biological pathway leading from initial molecular interactions of a toxicant and progressing through a series of key events (KEs), culminating with an apical adverse outcome (AO) that has to be of regulatory relevance. An AOP based on the mode of action (MOA) of rodent liver tumor promotion by dioxin-like compounds (DLCs) has been developed and the weight of evidence (WoE) of key event relationships (KERs) evaluated using evolved Bradford Hill considerations. Dioxins and DLCs are potent aryl hydrocarbon receptor (AHR) ligands that cause a range of species-specific adverse outcomes. The occurrence of KEs is necessary for inducing downstream biological responses and KEs may occur at the molecular, cellular, tissue and organ levels. The common convention is that an AOP begins with the toxicant interaction with a biological response element; for this AOP, this initial event is binding of a DLC ligand to the AHR. Data from mechanistic studies, lifetime bioassays and approximately thirty initiation-promotion studies have established dioxin and DLCs as rat liver tumor promoters. Such studies clearly show that sustained AHR activation, weeks or months in duration, is necessary to induce rodent liver tumor promotion--hence, sustained AHR activation is deemed the molecular initiating event (MIE). After this MIE, subsequent KEs are 1) changes in cellular growth homeostasis likely associated with expression changes in a number of genes and observed as development of hepatic foci and decreases in apoptosis within foci; 2) extensive liver toxicity observed as the constellation of effects called toxic hepatopathy; 3) cellular proliferation and hyperplasia in several hepatic cell types. This progression of KEs culminates in the AO, the development of hepatocellular adenomas and carcinomas and cholangiolar carcinomas. A rich data set provides both qualitative and quantitative knowledge of the progression of this AOP through KEs and the KERs. Thus, the WoE for this AOP is judged to be strong. Species-specific effects of dioxins and DLCs are well known--humans are less responsive than rodents and rodent species differ in sensitivity between strains. Consequently, application of this AOP to evaluate potential human health risks must take these differences into account.

Cross-species transcriptomic analysis elucidates constitutive aryl hydrocarbon receptor activity

Background

Research on the aryl hydrocarbon receptor (AHR) has largely focused on variations in toxic outcomes resulting from its activation by halogenated aromatic hydrocarbons. But the AHR also plays key roles in regulating pathways critical for development, and after decades of research the mechanisms underlying physiological regulation by the AHR remain poorly characterized. Previous studies identified several core genes that respond to xenobiotic AHR ligands across a broad range of species and tissues. However, only limited inferences have been made regarding its role in regulating constitutive gene activity, i.e. in the absence of exogenous ligands. To address this, we profiled transcriptomic variations between AHR-active and AHR-less-active animals in the absence of an exogenous agonist across five tissues, three of which came from rats (hypothalamus, white adipose and liver) and two of which came from mice (kidney and liver). Because AHR status alone has been shown sufficient to alter transcriptomic responses, we reason that by contrasting profiles amongst AHR-variant animals, we may elucidate effects of the AHR on constitutive mRNA abundances.

Results

We found significantly more overlap in constitutive mRNA abundances amongst tissues within the same species than from tissues between species and identified 13 genes (Agt, Car3, Creg1, Ctsc, E2f6, Enpp1, Gatm, Gstm4, Kcnj8, Me1, Pdk1, Slc35a3, and Sqrdl) that are affected by AHR-status in four of five tissues. One gene, Creg1, was significantly up-regulated in all AHR-less-active animals. We also find greater overlap between tissues at the pathway level than at the gene level, suggesting coherency to the AHR signalling response within these processes. Analysis of regulatory motifs suggests that the AHR mostly mediates transcriptional regulation via direct binding to response elements.

Conclusions

These findings, though preliminary, present a platform for further evaluating the role of the AHR in regulation of constitutive mRNA levels and physiologic function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1053) contains supplementary material, which is available to authorized users.

Gender-specific transcriptomic response to environmental exposure in Flemish adults

Flanders, Belgium, is one of the most densely populated areas in Europe. The Flemish Environment and Health Survey (2002-2006) aimed at determining exposure to pollutants of neonates, adolescents, and older adults and to assess associated biological and health effects. This study investigated genome wide gene expression changes associated with a range of environmental pollutants, including cadmium, lead, PCBs, dioxin, hexachlorobenzene, p,p'-DDE, benzene, and PAHs. Gene expression levels were measured in peripheral blood cells of 20 adults with relatively high and 20 adults with relatively low combined internal exposure levels, all non-smokers aged 50-65. Pearson correlation was used to analyze associations between pollutants and gene expression levels, separately for both genders. Pollutant- and gender-specific correlation analysis results were obtained. For organochlorine pollutants, analysis within genders revealed that genes were predominantly regulated in opposite directions in males and females. Significantly modulated pathways were found to be associated with each of the exposure biomarkers measured. Pathways and/or genes related to estrogen and STAT5 signaling were correlated to organochlorine exposures in both genders. Our work demonstrates that gene expression in peripheral blood is influenced by environmental pollutants. In particular, gender-specific changes are associated with organochlorine pollutants, including gender-specific modulation of endocrine related pathways and genes. These pathways and genes have previously been linked to endocrine disruption related disorders, which in turn have been associated with organochlorine exposure. Based on our results, we recommend that males and females be considered separately when analyzing gene expression changes associated with exposures that may include chemicals with endocrine disrupting properties.

Downregulation of mouse hepatic CYP3A protein by 3-methylcholanthrene does not require cytochrome P450-dependent metabolism

The aryl hydrocarbon receptor (AHR)-dependent induction of cytochromes P450 (P450) such as CYP1A1 by 3-methylcholanthrene (MC) and related polycyclic aromatic hydrocarbons is well characterized. We reported previously that MC treatment triggers a pronounced downregulation, particularly at the protein level, of mouse hepatic Cyp3a11, a counterpart of the key human drug-metabolizing enzyme CYP3A4. To determine whether this effect of MC requires hepatic microsomal P450 activity, we studied liver Cpr-null (LCN) mice with hepatocyte-specific conditional deletion of the NADPH-cytochrome P450 oxidoreductase gene. In vehicle-treated animals, basal levels of CYP3A11 mRNA and CYP3A protein immunoreactivity were elevated by approximately 9-fold in LCN mice compared with wild-type (WT) mice, whereas CYP3A catalytic activity was profoundly compromised in LCN mice. MC treatment caused suppression of CYP3A11 mRNA, CYP3A protein immunoreactivity, and CYP3A catalytic activity in WT mice, and the MC effects at the mRNA and protein levels were maintained in LCN mice. Flavin-containing monooxygenase-3 (Fmo3) induction by MC was suggested previously to occur via an AHR-dependent mechanism requiring conversion of the parent compound to DNA-damaging reactive metabolites; however, hepatic FMO3 mRNA levels were dramatically increased by MC in both WT and LCN mice. MC did not function as a mechanism-based inactivator of CYP3A enzymes in hepatic microsomes prepared from untreated WT mice, under conditions in which 1-aminobenzotriazole caused marked NADPH-dependent loss of total P450 content and CYP3A catalytic activity. These results indicate that MC downregulates mouse hepatic CYP3A protein via a pretranslational mechanism that does not require hepatic microsomal P450-dependent activity.

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.

Hepatic transcriptomic responses to TCDD in dioxin-sensitive and dioxin-resistant rats during the onset of toxicity

The dioxin congener 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) causes a wide range of toxic effects in rodent species, all of which are mediated by a ligand-dependent transcription-factor, the aryl hydrocarbon receptor (AHR). The Han/Wistar (Kuopio) (H/W) strain shows exceptional resistance to many TCDD-induced toxicities; the LD₅₀ of > 9600 μg/kg for H/W rats is higher than for any other wild-type mammal known. We previously showed that this resistance primarily results from H/W rats expressing a variant AHR isoform that has a substantial portion of the AHR transactivation domain deleted. Despite this large deletion, H/W rats are not entirely refractory to the effects of TCDD; the variant AHR in these animals remains fully competent to up-regulate well-known dioxin-inducible genes. TCDD-sensitive (Long-Evans, L-E) and resistant (H/W) rats were treated with either corn-oil (with or without feed-restriction) or 100 μg/kg TCDD for either four or ten days. Hepatic transcriptional profiling was done using microarrays, and was validated by RT-PCR analysis of 41 genes. A core set of genes was altered in both strains at all time points tested, including CYP1A1, CYP1A2, CYP1B1, Nqo1, Aldh3a1, Tiparp, Exoc3, and Inmt. Outside this core, the strains differed significantly in the breadth of response: three-fold more genes were altered in L-E than H/W rats. At ten days almost all expressed genes were dysregulated in L-E rats, likely reflecting emerging toxic responses. Far fewer genes were affected by feed-restriction, suggesting that only a minority of the TCDD-induced changes are secondary to the wasting syndrome.

Toxicogenomic analysis of exposure to TCDD, PCB126 and PCB153: identification of genomic biomarkers of exposure to AhR ligands

Background

Two year cancer bioassays conducted by the National Toxicology Program have shown chronic exposure to dioxin-like compounds (DLCs) to lead to the development of both neoplastic and non-neoplastic lesions in the hepatic tissue of female Sprague Dawley rats. Most, if not all, of the hepatotoxic effects induced by DLC's are believed to involve the binding and activation of the transcription factor, the aryl hydrocarbon receptor (AhR). Toxicogenomics was implemented to identify genomic responses that may be contributing to the development of hepatotoxicity in rats.

Results

Through comparative analysis of time-course microarray data, unique hepatic gene expression signatures were identified for the DLCs, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) (100 ng/kg/day) and 3,3',4,4',5-pentachlorobiphenyl (PCB126) (1000 ng/kg/day) and the non-DLC 2,2',4,4',5,5',-hexachlorobiphenyl (PCB153) (1000 μg/kg/day). A common time independent signature of 41 AhR genomic biomarkers was identified which exhibited at least a 2-fold change in expression following subchronic (13-wk) and chronic (52-wk) p.o. exposure to TCDD and PCB126, but not the non DLC, PCB153. Real time qPCR analysis validated that 30 of these genes also exhibited at least a 2-fold change in hepatic expression at 24 hr following a single exposure to TCDD (5 μg/kg, po). Phenotypic anchoring was conducted which identified forty-six genes that were differently expressed both following chronic p.o. exposure to DLCs and in previously reported studies of cholangiocarcinoma or hepatocellular adenoma.

Conclusions

Together these analyses provide a comprehensive description of the genomic responses which occur in rat hepatic tissue with exposure to AhR ligands and will help to isolate those genomic responses which are contributing to the hepatotoxicity observed with exposure to DLCs. In addition, the time independent gene expression signature of the AhR ligands may assist in identifying other agents with the potential to elicit dioxin-like hepatotoxic responses.

Polymorphisms and functional differences in aryl hydrocarbon receptors (AhR) in Japanese field mice, Apodemus speciosus

Dioxins, which are unintentionally generated toxic pollutants, exert a variety of adverse effects on organisms. The majority of these effects, which include teratogenesis, immunosuppression, tumor promotion, and endocrine disruption, are mediated through aryl hydrocarbon receptor (AhR), a ligand-activated transcription factor. Genetic variations in AhR result in different survivability under exposure to dioxin contamination, which might affect the genetic structure of wildlife populations through differential susceptibility to dioxin exposure. The aim of this study was to clarify the polymorphisms of AhR in Japanese field mice, Apodemus speciosus, and their functional differences in order to develop a molecular indicator for dioxin sensitivity. Wild Japanese field mice had abundant polymorphisms in AhR coding region. Seventy-one single nucleotide polymorphisms, 27 of which occur amino acid substitutions, and consequently 49 alleles were identified in 63 individuals. In the functional analysis of AhR variants using transient reporter assays, a Gln to Arg mutation at amino acid 799 exhibited a significant decrease in the level of transactivational properties (p=0.015) which might modify the dioxin susceptibility of an individual.

Aryl hydrocarbon receptor (AHR)-regulated transcriptomic changes in rats sensitive or resistant to major dioxin toxicities

BACKGROUND

The major toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) appear to result from dysregulation of mRNA levels mediated by the aryl hydrocarbon receptor (AHR). Dioxin-like chemicals alter expression of numerous genes in liver, but it remains unknown which lie in pathways leading to major toxicities such as hepatotoxicity, wasting and lethality. To identify genes involved in these responses we exploited a rat genetic model. Rats expressing an AHR splice-variant lacking a portion of the transactivation domain are highly resistant to dioxin-induced toxicities. We examined changes in hepatic mRNA abundances 19 hours after TCDD treatment in two dioxin-resistant rat strains/lines and two dioxin-sensitive rat strains/lines.

RESULTS

Resistant rat strains/lines exhibited fewer transcriptional changes in response to TCDD than did rats with wildtype AHR. However, well-known AHR-regulated and dioxin-inducible genes such as CYP1A1, CYP1A2, and CYP1B1 remained fully responsive to TCDD in all strains/lines. Pathway analysis indicated that the genes which respond differently to TCDD between sensitive and resistant rats are mainly involved in lipid metabolism, cellular membrane function and energy metabolism. These pathways previously have been shown to respond differently to dioxin treatment in dioxin-sensitive versus dioxin-resistant rats at a biochemical level and in the differential phenotype of toxicologic responses.

CONCLUSION

The transactivation-domain deletion in dioxin-resistant rats does not abolish global AHR transactivational activity but selectively interferes with expression of subsets of genes that are candidates to mediate or protect from major dioxin toxicities such as hepatotoxicity, wasting and death.

Isoform distinct time-, dose-, and castration-dependent alterations in flavin-containing monooxygenase expression in mouse liver after 2,3,7,8-tetrachlorodibenzo-p-dioxin treatment

Flavin-containing monooxygenase (FMO) expression in male mouse liver is altered after 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) exposure or castration. Because TCDD is slowly eliminated from the body, we examined hepatic Fmo mRNA alterations for up to 32 days following 10 or 64 microg/kg TCDD exposure by oral gavage in male C57BL/6J mice. Fmo2 mRNA was significantly induced at 1, 4, and 8 days whereas Fmo3 mRNA was also induced at 32 days relative to controls. Fmo3 mRNA levels exhibited a dose-dependent increase at 4, 8, and 32 days after exposure; Fmo1, Fmo4, and Fmo5 mRNA did not exhibit clear trends. Because castration alone also increased Fmo2, Fmo3, and Fmo4 mRNA we examined the combined effects of castration and TCDD treatment on FMO expression. A greater than additive effect was observed with Fmo2 and Fmo3 mRNA expression. Fmo2 mRNA exhibited a 3-5-fold increase after castration or 10 microg/kg TCDD exposure by oral gavage, whereas an approximately 20-fold increase was observed between the sham-castrated control and castrated TCDD-treated mice. Similarly, treatment with 10 microg/kg TCDD alone increased Fmo3 mRNA 130- and 180-fold in the sham-castrated and castrated mice compared to their controls respectively, whereas, Fmo3 mRNA increased approximately 1900-fold between the sham control and castrated TCDD-treated mice. An increase in hepatic Fmo3 protein in TCDD-treated mice was observed by immunoblotting and assaying methionine S-oxidase activity. Collectively, these results provide evidence for isoform distinct time-, dose-, and castration-dependent effects of TCDD on FMO expression and suggest cross-talk between TCDD and testosterone signal transduction pathways.

Divergent transcriptomic responses to aryl hydrocarbon receptor agonists between rat and human primary hepatocytes

Toxicogenomics has great potential for enhancing our understanding of environmental chemical toxicity, hopefully leading to better informed human health risk assessments. This study employed toxicogenomic technology to reveal species differences in response to two prototypical aryl hydrocarbon receptor (AHR) agonists 2,3,7,8-tetrachlorodibenzo-p-dioxin and the polychlorinated biphenyl (PCB) congener PCB 126. Dose-responses of primary cultures of rat and human hepatocytes were determined using species-specific microarrays sharing over 4000 gene orthologs. Forty-seven human and 79 rat genes satisfied dose-response criteria for both chemicals and were subjected to further analysis including the calculation of the 50% effective concentration and the relative potency (REP) of PCB 126 for each gene. Only five responsive orthologous genes were shared between the two species; yet, the geometric mean of the REPs for all rat and human modeled responsive genes were 0.06 (95% confidence interval [CI]; 0.03-0.1) and 0.002 (95% CI; 0.001-0.005), respectively, suggesting broad species differences in the initial events that follow AHR activation but precede toxicity. This indicates that there are species differences in both the specific genes that responded and the agonist potency and REP for those genes. This observed insensitivity of human cells to PCB 126 is consistent with more traditional measurements of AHR activation (i.e., cytochrome P450 1A1 enzyme activity) and suggests that the species difference in PCB 126 sensitivity is likely due to certain aspects of AHR function. That a species divergence also exists in this expanded AHR-regulated gene repertoire is a novel finding and should help when extrapolating animal data to humans.

The aryl hydrocarbon receptor: a perspective on potential roles in the immune system

The aryl hydrocarbon receptor (AHR) is a protein best known for its role in mediating toxicity. Over 30 years of research has uncovered additional roles for the AHR in xenobiotic metabolism and normal vascular development. Activation of the AHR has long been known to cause immunotoxicity, including thymic involution. Recent data suggesting a role for the AHR in regulatory T-cell (Treg) and T-helper 17 (Th17) cell development have only added to the excitement about this biology. In this review, we will attempt to illustrate what is currently known about AHR biology in the hope that data from fields as diverse as evolutionary biology and pharmacology will help elucidate the mechanism by which AHR modifies immune responses. We also will discuss the complexities of AHR pharmacology and genetics that may influence future studies of AHR in the immune system.

Unsupervised assessment of microarray data quality using a Gaussian mixture model

BACKGROUND

Quality assessment of microarray data is an important and often challenging aspect of gene expression analysis. This task frequently involves the examination of a variety of summary statistics and diagnostic plots. The interpretation of these diagnostics is often subjective, and generally requires careful expert scrutiny.

RESULTS

We show how an unsupervised classification technique based on the Expectation-Maximization (EM) algorithm and the naïve Bayes model can be used to automate microarray quality assessment. The method is flexible and can be easily adapted to accommodate alternate quality statistics and platforms. We evaluate our approach using Affymetrix 3' gene expression and exon arrays and compare the performance of this method to a similar supervised approach.

CONCLUSION

This research illustrates the efficacy of an unsupervised classification approach for the purpose of automated microarray data quality assessment. Since our approach requires only unannotated training data, it is easy to customize and to keep up-to-date as technology evolves. In contrast to other "black box" classification systems, this method also allows for intuitive explanations.

Cytochrome P450 2C11 5'-flanking region and promoter mediate in vivo suppression by 3-methylcholanthrene

Aromatic hydrocarbons such as 3-methylcholanthrene (MC) elicit toxic and adaptive responses through the aryl hydrocarbon receptor (AHR). Aromatic hydrocarbons act via an unknown mechanism to suppress the transcription of CYP2C11, a growth hormone-regulated gene encoding the male-specific rat hepatic cytochrome P450 2C11. We hypothesize that suppression of CYP2C11 by aromatic hydrocarbons is mediated by the gene's promoter and 5'-flank. Using hydrodynamics-based injections to deliver plasmid DNA to the liver of live rats, we studied the MC responsiveness of luciferase constructs containing 10.1, 5.6, and 2.4 kilobases (kb) of the CYP2C11 5'-flank. MC suppressed CYP2C11-luciferase activity of the 10.1- and 5.6-kb constructs to less than 50% of vehicle levels by 24 and 72 h. Luciferase activity of the 2.4-kb CYP2C11 construct was decreased to 63% of vehicle levels 24 h after MC treatment, but no suppression was detected by 72 h. Negative regulatory element(s) responsible for CYP2C11 reporter suppression by MC exist in the proximal 2.4 kb of the 5'-flank; however, additional cis-acting elements located between -5.6 and -2.4 kb mediate persistent reporter suppression. As a positive control for AHR activation, MC dramatically induced the luciferase activity of a Cyp1a1-driven luciferase plasmid under AHR control. Modulation of reporter gene activity by MC was accompanied by induction of endogenous CYP1A1 and suppression of endogenous CYP2C11 mRNA/protein. This is the first demonstration of aromatic hydrocarbon-mediated suppression of a CYP2C11-luciferase construct, and this finding suggests that the 5'-flanking region and promoter mediate down-regulation of this gene in the intact rat.

Patterns of dioxin-altered mRNA expression in livers of dioxin-sensitive versus dioxin-resistant rats

Dioxins exert their major toxicologic effects by binding to the aryl hydrocarbon receptor (AHR) and altering gene transcription. Numerous dioxin-responsive genes previously were identified both by conventional biochemical and molecular techniques and by recent mRNA expression microarray studies. However, of the large set of dioxin-responsive genes the specific genes whose dysregulation leads to death remain unknown. To identify specific genes that may be involved in dioxin lethality we compared changes in liver mRNA levels following exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in three strains/lines of dioxin-sensitive rats with changes in three dioxin-resistant rat strains/lines. The three dioxin-resistant strains/lines all harbor a large deletion in the transactivation domain of the aryl hydrocarbon receptor (AHR). Despite this deletion, many genes exhibited a "Type-I" response-that is, their responses were similar in dioxin-sensitive and dioxin-resistant rats. Several genes that previously were well established as being dioxin-responsive or under AHR regulation emerged as Type-I responses (e.g. CYP1A1, CYP1A2, CYP1B1 and Gsta3). In contrast, a relatively small number of genes exhibited a Type-II response-defined as a difference in responsiveness between dioxin-sensitive and dioxin-resistant rat strains. Type-II genes include: malic enzyme 1, ubiquitin C, cathepsin L, S-adenosylhomocysteine hydrolase and ferritin light chain 1. In silico searches revealed that AH response elements are conserved in the 5'-flanking regions of several genes that respond to TCDD in both the Type-I and Type-II categories. The vast majority of changes in mRNA levels in response to 100 microg/kg TCDD were strain-specific; over 75% of the dioxin-responsive clones were affected in only one of the six strains/lines. Selected genes were assessed by quantitative RT-PCR in dose-response and time-course experiments and responses of some genes were assessed in Ahr-null mice to determine if their response was AHR-dependent. Type-II genes may lie in pathways that are central to the difference in susceptibility to TCDD lethality in this animal model.

Alteration of gene expression profiles in the brain of Japanese medaka (Oryzias latipes) exposed to KC-400 or PCB126

Polychlorinated biphenyls (PCBs) are known as neurotoxic chemicals and possibly alter animal behavior. We previously reported that PCB-exposure induced abnormal schooling behavior in Japanese medaka (Oryzias latipes). This abnormal behavior might be caused by the functional alteration of central or terminal nervous system. To understand the mechanism(s) of behavioral change by PCB-exposure, we analyzed the gene expression profiles in the brain of medaka exposed to 3,3',4,4',5-pentachlorobiphenyl (PCB126) or a PCB mixture (Kanechlor-400: KC-400) using a cDNA microarray that we constructed. Twelve FLF-II strain medaka (six individuals per treatment) were dietary exposed to PCB126 (0.01 microg/g b.w./day) or KC-400 (1 microg/g b.w./day) for three weeks. For the control, six fish were fed a control diet. After the exposure period, fish were dissected, and the brain samples were collected. The samples from control fish were pooled and used as a common reference in the microarray experiment. Microarray data were normalized by the LOWESS method, and we screened the genes whose expression levels were altered more than 1.5-fold. Gene expression profiling showed 97 down-regulated and 379 up-regulated genes in the brain of medaka exposed to PCB126. KC-400 exposure suppressed 15 genes and induced 266 genes in medaka brain. Among these genes, the expression levels of 7 and 188 genes were commonly down- or up-regulated, respectively in both treatment groups. On the other hand, 31 gene expressions were significantly different between PCB126 and KC-400 treatment groups, and three out of 31 genes were received opposite effects. In addition, the microarray data showed that thyroid hormone-responsive genes were up-regulated by PCB-exposure, which may imply that PCBs or their metabolites mimic thyroid hormone effects in the brain of PCB-exposed medaka.

Cytochrome P450 2C11 5'-flanking region and promoter: regulation by aromatic hydrocarbons in vitro

Aromatic hydrocarbons elicit toxic and adaptive responses via the aryl hydrocarbon receptor (AHR). Aromatic hydrocarbons suppress the transcription of the growth hormone-regulated, male-specific rat hepatic cytochrome P450 2C11 gene (CYP2C11) in vivo via an unknown mechanism. We hypothesize that the suppression of CYP2C11 by aromatic hydrocarbons is mediated by the gene's promoter and 5'-flanking region. Following bioinformatic analysis of putative transcription factor (TF) binding sites, we cloned extended lengths of the CYP2C11 5'-flanking region into a promoterless luciferase plasmid. Suppression of CYP2C11 constructs was not observed upon treatment of transfected rat 5L, BP8 or mouse Hepa-1 cells with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) or 3-methylcholanthrene. In human HepG2 cells, the 10.1-kb construct displayed a pronounced 6- to 8-fold induction by TCDD. Deletion analysis localized the paradoxical induction response to a region between -1.8 kb and -1.3 kb, which contains a dioxin-responsive element (DRE) previously shown by us to be capable of binding activated AHR. This was confirmed by site-directed mutagenesis of the DRE. Induction of the 10.1-kb construct by TCDD in HepG2 cells was blocked by alpha-naphthoflavone, an AHR antagonist/partial agonist. The AHR is likely involved in the induction of CYP2C11-luciferase activity by TCDD in HepG2 cells and this response is at least partly DRE-mediated. Although CYP2C11 is suppressed by aromatic hydrocarbons in vivo, CYP2C11-luciferase constructs display a potentially misleading paradoxical induction in vitro that is cell-specific. Regulation of CYP2C11-luciferase plasmids is being studied in vivo in rat liver, where an intact endocrine system and the full complement of TFs needed for CYP2C11 suppression are present.

microRNAs in Adult Rodent Liver Are Refractory to Dioxin Treatment

Dioxin-like chemicals are well known for their ability to upregulate expression of numerous genes via the AH receptor (AHR). However, recent transcriptomic analyses in several laboratories indicate that dioxin-like chemicals or AHR genotype itself also can downregulate levels of mRNAs encoded by numerous genes. The mechanism responsible for such downregulation is unknown. We hypothesized that microRNAs (miRNAs), which have emerged as powerful negative regulators of mRNA levels in several systems, might be responsible for mRNA downregulation in dioxin/AHR pathways. We used two miRNA array platforms as well as quantitative reverse transcriptase-polymerase chain reaction to measure miRNA levels in wild-type (WT) versus Ahr-null mice, in dioxin-sensitive Long-Evans (L-E; Turku/AB) rats versus dioxin-resistant Han/Wistar (H/W; Kuopio) rats and in rat 5L and mouse Hepa-1 hepatoma cells in culture. Treatment with 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in vivo caused few changes in miRNA levels in mouse or rat livers, and those changes that were statistically significant were of modest magnitude. Hepatoma cells in culture also exhibited few changes in miRNA levels in response to TCDD. AHR genotype had little effect on hepatic miRNA levels, either in constitutive expression or in response to TCDD-only a few miRNAs differed in expression between Ahr-null mice compared to mice with WT AHR or between L-E rats (that have WT AHR) compared to H/W rats (whose AHR has a large deletion in the transactivation domain). It is unlikely that mRNA downregulation by dioxins is mediated by miRNAs, nor are miRNAs likely to play a significant role in dioxin toxicity in adult rodent liver.