Transcriptome profiling of human hepatocytes treated with Aroclor 1254 reveals transcription factor regulatory networks and clusters of regulated genes

Differential expression of genes in C. elegans reveals transcriptional responses to indirect-acting xenobiotic compounds and insensitivity to 2,3,7,8-tetrachlorodibenzodioxin

Caenorhabditis elegans is a well-established model organism for toxicity testing of chemical substances. We recently demonstrated its potential for bioanalysis of the toxic potency of chemical contaminants in water. While many detoxification genes are homologues to those in mammalians, C. elegans is reported to be deficient in cytochrome CYP1-like P450 metabolism and that its aryl hydrocarbon receptor (AhR) homolog encoded by ahr-1 purportedly does not interact with dioxins or any other known xenobiotic ligand. This suggests that C. elegans is insensitive for compounds that require bioactivation (indirectly acting compounds) and for dioxins or dioxin-like compounds. This study analysed genome-wide gene expression of the nematode in response to 30 μM of aflatoxin B1 (AFB1), benzo(a)pyrene (B(a)P), Aroclor 1254 (PCB1254), and 10 μM of 2,3,7,8-tetrachlorodibenzodioxin (TCDD). After 24 h of exposure in the early L4 larval stage, microarray analysis revealed 182, 86, and 321 differentially expressed genes in the nematodes treated with 30 μM of AFB1, B(a)P, and PCB1254, respectively. Among these genes, many encode xenobiotic-metabolizing enzymes, and their transcription levels were among the highest-ranked fold-changed genes. Interestingly, only one gene (F59B1.8) was upregulated in the nematodes exposed to 10 μM TCDD. Genes related to metabolic processes and catalytic activity were the most induced by exposure to 30 μM of AFB1, B(a)P, and PCB1254. Despite the genotoxic nature of AFB1 and B(a)P, no differential expression was found in the genes encoding DNA repair and cell cycle checkpoint proteins. Analysis of concentration-response curves was performed to determine the Lowest Observed Transcriptomic Effect Levels (LOTEL) of AFB1, B(a)P, and PCB1254. The obtained LOTEL values showed that gene expression changes in C. elegans are more sensitive to toxicants than reproductive effects. Overall, transcriptional responses of metabolic enzymes suggest that the nematode does metabolize AFB1, B(a)P, and PCB1254. Our findings also support the assumption that the transcription factor AhR homolog in C. elegans does not bind typical xenobiotic ligands, rendering the nematode transcriptionally insensitive to TCDD effects.

Hepatocyte size fractionation allows dissection of human liver zonation

Human hepatocytes show marked differences in cell size, gene expression, and function throughout the liver lobules, an arrangement termed liver zonation. However, it is not clear if these zonal size differences, and the associated phenotypic differences, are retained in isolated human hepatocytes, the "gold standard" for in vitro studies of human liver function. Here, we therefore explored size differences among isolated human hepatocytes and investigated whether separation by size can be used to study liver zonation in vitro. We used counterflow centrifugal elutriation to separate cells into different size fractions and analyzed them with label-free quantitative proteomics, which revealed an enrichment of 151 and 758 proteins (out of 5163) in small and large hepatocytes, respectively. Further analysis showed that protein abundances in different hepatocyte size fractions recapitulated the in vivo expression patterns of previously described zonal markers and biological processes. We also found that the expression of zone-specific cytochrome P450 enzymes correlated with their metabolic activity in the different fractions. In summary, our results show that differences in hepatocyte size matches zonal expression patterns, and that our size fractionation approach can be used to study zone-specific liver functions in vitro.

Aroclor 1254 and BDE-47 inhibit dopaminergic function manifesting as changes in locomotion behaviors in zebrafish embryos

Contamination with polychlorinated biphenyls (PCBs) and polybrominated diphenyl ethers (PBDEs) in the environment is a major concern due to their persistent bioaccumulative toxicity that can disturb neurobehavioral functions including movements. Recently, it was reported that some PBDE including BDE-47 stimulates locomotor activities of zebrafish embryos by unknown mechanism. In this study, motor movements of the zebrafish embryo were used as a model system to evaluate the neuronal toxicity of a non-coplanar PCB-dominant mixture (Aroclor 1254) and BDE-47. Both organohalogens increased tail shaking and rotation of embryos in a concentration-dependent manner. Chemical inhibition and gene knock-down of tyrosine hydroxylase and vesicular monoamine transporter 2 (VMAT2) also induced hyperactivities. Hyperactivities induced by these treatments were all inhibited by supplementation of l-tyrosine and l-dopa, precursors of dopamine synthesis. Both organohalogens reduced dopamine contents and increased the 3,4-dihydroxyphenylacetic acid (DOPAC)/dopamine ratio in whole embryos. The results suggest that functional inhibition of dopaminergic neurons is involved in hyperactivities of zebrafish embryos caused by Aroclor 1254 and BDE-47.

Differential expression profile of MicroRNAs during differentiation of cardiomyocytes exposed to polychlorinated biphenyls

Exposure to persistent environmental pollutants, such as polychlorinated biphenyls (PCBs), is a risk factor for the development of congenital heart defects. MicroRNAs (miRNAs) have been shown to be involved in cardiac development. The objective of this study was to investigate changes in miRNA expression profiles during the differentiation of cardiomyocytes exposed to PCBs. For that purpose, PCBs (Aroclor 1254) at a concentration of 2.5 μmol/L were added on day 0 of differentiation of P19 mouse embryonal carcinoma cells into cardiac myocytes. The relative expression of miRNA genes was determined by miRNA microarray and real-time reverse transcriptase polymerase chain reaction (real-time RT-PCR) analyses. The microarray results revealed that 45 miRNAs, of which 14 were upregulated and 31 were downregulated, were differentially expressed in P19 cells treated with PCBs compared with control cells. The miRNA expression data was validated with real-time RT-PCR. The expression of certain potential target genes (Wnt1) was found to be reduced in P19 cells treated with PCBs, whereas the expression of other potential predicted target genes (GSK3β) was increased. Our results demonstrate a critical role of miRNAs in mediating the effect of PCBs during the differentiation of P19 cells into cardiac myocytes.

Disease-driven detection of differential inherited SNP modules from SNP network

Detection of the synergetic effects between variants, such as single-nucleotide polymorphisms (SNPs), is crucial for understanding the genetic characters of complex diseases. Here, we proposed a two-step approach to detect differentially inherited SNP modules (synergetic SNP units) from a SNP network. First, SNP-SNP interactions are identified based on prior biological knowledge, such as their adjacency on the chromosome or degree of relatedness between the functional relationships of their genes. These interactions form SNP networks. Second, disease-risk SNP modules (or sub-networks) are prioritised by their differentially inherited properties in IBD (Identity by Descent) profiles of affected and unaffected sibpairs. The search process is driven by the disease information and follows the structure of a SNP network. Simulation studies have indicated that this approach achieves high accuracy and a low false-positive rate in the identification of known disease-susceptible SNPs. Applying this method to an alcoholism dataset, we found that flexible patterns of susceptible SNP combinations do play a role in complex diseases, and some known genes were detected through these risk SNP modules. One example is GRM7, a known alcoholism gene successfully detected by a SNP module comprised of two SNPs, but neither of the two SNPs was significantly associated with the disease in single-locus analysis. These identified genes are also enriched in some pathways associated with alcoholism, including the calcium signalling pathway, axon guidance and neuroactive ligand-receptor interaction. The integration of network biology and genetic analysis provides putative functional bridges between genetic variants and candidate genes or pathways, thereby providing new insight into the aetiology of complex diseases.

Transcriptomics identifies differences between ultrapure non-dioxin-like polychlorinated biphenyls (PCBs) and dioxin-like PCB126 in cultured peripheral blood mononuclear cells

Polychlorinated biphenyls (PCBs) remain ubiquitously present in human lipids despite the ban on their production and use. Their presence can be chemically monitored in peripheral blood samples of the general population. We tested whether in vitro exposure to different PCB congeners induced different gene expression profiles in peripheral blood cells. We have isolated peripheral blood mononuclear cells (PBMC) from whole blood of 8 healthy individuals and exposed these cells in vitro to individual non-dioxin-like (NDL)-PCB congeners (PCB52, 138 or 180; 10μM) or dioxin-like (DL)-PCB congener PCB126 (1μM) during 18h. Differential gene expression response was measured using Agilent whole-human genome microarrays. Two-way ANOVA analysis of the data showed that both gender and PCB exposure are important factors influencing gene expression responses in blood cells. Hierarchical cluster analysis of genes influenced by PCB exposure, revealed that DL-PCB126 induced a different gene expression response compared to the NDL-PCBs. Biological interpretation of the results revealed that exposure to PCB126 induced the AhR signaling pathway, whereas the induction of nuclear receptor pathways by the NDL-PCBs was limited in blood cells. Nevertheless, molecular responses of blood cells to individual PCB congeners revealed significantly expressed genes that play a role in biological functions and processes known to be affected by PCB exposure in vivo. Observed gene expression changes in this in vitro model were found to be related to hepatotoxicity, immune and inflammatory response and disturbance of lipid and cholesterol homeostasis.

Genetic architecture of susceptibility to PCB126-induced developmental cardiotoxicity in zebrafish

Variability in risk of developmental defects caused by dioxin-like compounds (DLCs) has been demonstrated within and among several vertebrate species. Beyond our knowledge of the aryl hydrocarbon receptor (AHR) and its role in mediating toxicity for this class of compounds, little else is known concerning precise downstream targets influencing this vulnerability. In the present study, zebrafish with divergent genetic backgrounds were screened for susceptibility to developmental cardiotoxicity caused by the prototypical DLC, 3,3',4,4',5-pentachlorobiphenyl (PCB126); a range up to ∼40-fold differences was observed. Differentially sensitive zebrafish were chosen for a genetic cross, and the recombinant generation was used for genome-wide quantitative trait loci (QTL) mapping. Multiple QTLs were identified--several acting alone, one additively, and two others via epistatic interaction. Together, these QTLs account for 24% of the phenotypic variance observed in cardioteratogenicity resulting from PCB126 exposure (logarithm of the odds = 13.55, p = 1.89 × 10⁻¹⁰). Candidate genes in these QTL regions include the following: ahr2, bcor, and capn1 (Chr 22); e2f1 and pdyn (Chr 23); ctnnt2, plcg1, eno3, tgm1, and tgm2 (interacting on Chr 23); and vezf1 (Chr 15). These data demonstrate that DLC-induced cardiac teratogenicity is a multifactorial complex trait influenced by gene × gene and gene × environment interactions. The identified QTLs harbor many DLC-responsive genes critical to cardiovascular development and provide insight into the genetic basis of susceptibility to AHR-mediated developmental toxicity.

Reactivation of L1 retrotransposon by benzo(a)pyrene involves complex genetic and epigenetic regulation

Benzo(a)pyrene (BaP), is an environmental pollutant present in tobacco smoke and a byproduct of fossil fuel combustion which likely contributes to the tumorigenic processes in human cancers including lung and esophageal. Long Interspersed Nuclear Element-1 (LINE-1) or L1 is a mobile element within the mammalian genome that propagates via a "copy-and-paste" mechanism using reverse transcriptase and RNA intermediates. L1 is strongly expressed during early embryogenesis and then silenced as cells initiate differentiation programming. Although the complex transcriptional control mechanisms of L1 are not well understood, L1 reactivation has been described in several human cancers and following exposure of mouse or human cells to BaP. In this study we investigated the molecular mechanisms and epigenetic events that regulate L1 reactivation following BaP exposure. We show that challenge of HeLa cells with BaP induces early enrichment of the transcriptionally-active chromatin markers histone H3 trimethylated at lysine 4 (H3K4Me3) and histone H3 acetylated at lysine 9 (H3K9Ac), and reduces association of DNA methyltransferase-1 (DNMT1) with the L1 promoter. These changes are followed by proteasome-dependent decreases in cellular DNMT1 expression and sustained reduction of cytosine methylation within the L1 promoter CpG island. Pharmacological inhibition of the proteasome signaling pathway with the inhibitor MG132 blocks degradation of DNMT1 and alters BaP-mediated histone epigenetic modifications. We conclude that genetic reactivation of L1 by BaP involves an ordered cascade of epigenetic events that begin with nucleosomal histone modifications and is completed with alterations in DNMT1 recruitment to the L1 promoter and reduced DNA methylation of CpG islands.

An integrated genomic analysis of aryl hydrocarbon receptor-mediated inhibition of B-cell differentiation

The aryl hydrocarbon receptor (AHR) agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) alters differentiation of B cells and suppresses antibody production. A combination of whole-genome, microarray-based chromatin immunoprecipitation (ChIP-on-chip), and time course gene expression microarray analysis was performed on the mouse B-cell line CH12.LX following exposure to lipopolysaccharide (LPS) or LPS and TCDD to identify the primary and downstream transcriptional elements of B-cell differentiation that are altered by the AHR. ChIP-on-chip analysis identified 1893 regions with a significant increase in AHR binding with TCDD treatment. Transcription factor binding site analysis on the ChIP-on-chip data showed enrichment in AHR response elements. Other transcription factors showed significant coenrichment with AHR response elements. When ChIP-on-chip regions were compared with gene expression changes at the early time points, 78 genes were identified as potential direct targets of the AHR. AHR binding and expression changes were confirmed for a subset of genes in primary mouse B cells. Network analysis examining connections between the 78 potential AHR target genes and three transcription factors known to regulate B-cell differentiation indicated multiple paths for potential regulation by the AHR. Enrichment analysis on the differentially expressed genes at each time point evaluated the downstream impact of AHR-regulated gene expression changes on B-cell-related processes. AHR-mediated impairment of B-cell differentiation occurred at multiple nodes of the B-cell differentiation network and potentially through multiple mechanisms including direct cis-acting effects on key regulators of B-cell differentiation, indirect regulation of B-cell differentiation-related pathways, and transcriptional coregulation of target genes by AHR and other transcription factors.

Comparative analysis of AhR-mediated TCDD-elicited gene expression in human liver adult stem cells

Time course and dose-response studies were conducted in HL1-1 cells, a human liver cell line with stem cell-like characteristics, to assess the differential gene expression elicited by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) compared with other established models. Cells were treated with 0.001, 0.01, 0.1, 1, 10, or 100nM TCDD or dimethyl sulfoxide vehicle control for 12 h for the dose-response study, or with 10nM TCDD or vehicle for 1, 2, 4, 8, 12, 24, or 48 h for the time course study. Elicited changes were monitored using a human cDNA microarray with 6995 represented genes. Empirical Bayes analysis identified 144 genes differentially expressed at one or more time points following treatment. Most genes exhibited dose-dependent responses including CYP1A1, CYP1B1, ALDH1A3, and SLC7A5 genes. Comparative analysis of HL1-1 differential gene expression to human HepG2 data identified 74 genes with comparable temporal expression profiles including 12 putative primary responses. HL1-1-specific changes were related to lipid metabolism and immune responses, consistent with effects elicited in vivo. Furthermore, comparative analysis of HL1-1 cells with mouse Hepa1c1c7 hepatoma cell lines and C57BL/6 hepatic tissue identified 18 and 32 commonly regulated orthologous genes, respectively, with functions associated with signal transduction, transcriptional regulation, metabolism and transport. Although some common pathways are affected, the results suggest that TCDD elicits species- and model-specific gene expression profiles.

Prediction and integration of regulatory and protein-protein interactions

Knowledge of transcriptional regulatory interactions (TRIs) is essential for exploring functional genomics and systems biology in any organism. While several results from genome-wide analysis of transcriptional regulatory networks are available, they are limited to model organisms such as yeast ( 1 ) and worm ( 2 ). Beyond these networks, experiments on TRIs study only individual genes and proteins of specific interest. In this chapter, we present a method for the integration of various data sets to predict TRIs for 54 organisms in the Bioverse ( 3 ). We describe how to compile and handle various formats and identifiers of data sets from different sources and how to predict TRIs using a homology-based approach, utilizing the compiled data sets. Integrated data sets include experimentally verified TRIs, binding sites of transcription factors, promoter sequences, protein subcellular localization, and protein families. Predicted TRIs expand the networks of gene regulation for a large number of organisms. The integration of experimentally verified and predicted TRIs with other known protein-protein interactions (PPIs) gives insight into specific pathways, network motifs, and the topological dynamics of an integrated network with gene expression under different conditions, essential for exploring functional genomics and systems biology.

Reproducible chemical-induced changes in gene expression profiles in human hepatoma HepaRG cells under various experimental conditions

The use of in vitro human liver cell models is an attractive approach in toxicogenomic studies designed to analyze gene expression changes induced by a toxic chemical. However, in such studies, reliability, reproducibility and interlaboratory concordance of microarrays, as well as the choice of the most suitable cell model, remain a matter of debate. This work was aimed at evaluating the robustness of microarray technologies and the suitability of the highly differentiated human HepaRG cell line in the investigation of gene expression changes induced by a toxic compound in human liver. The influence of various experimental conditions including cell cultures grown at different test sites, different generations of microarrays, RNA analysis platforms and softwares, was tested on gene expression profiles induced by a 20h treatment with an 8mM concentration of phenobarbital as the toxic compound. As many as 1099 genes (p-value<0.01 and 1.5-fold-change), representing 74% and 30% of the signature genes detected with Agilent 22 and 44K pangenomic microarrays, respectively, were shown to be modulated in common in six independently performed experiments. The most modulated genes included both those known to be regulated by phenobarbital, such as cytochromes P450 and membrane transporters, and those involved in oxidative stress, inflammation and apoptosis, typifying a toxic insult. These data provide strong support for the use of a toxicogenomic approach for the in vitro prediction of chemical toxicity, and for the choice of human HepaRG cells as a promising model system for human hepatotoxicity testing.

Dose- and time-dependent effects of phenobarbital on gene expression profiling in human hepatoma HepaRG cells

Phenobarbital (PB) induces or represses a wide spectrum of genes in rodent liver. Much less is known about its effects in human liver. We used pangenomic cDNA microarrays to analyze concentration- and time-dependent gene expression profile changes induced by PB in the well-differentiated human HepaRG cell line. Changes in gene expression profiles clustered at specific concentration ranges and treatment times. The number of correctly annotated genes significantly modulated by at least three different PB concentration ranges (spanning 0.5 to 3.2 mM) at 20 h exposure amounted to 77 and 128 genes (p< or =0.01) at 2- and 1.8-fold filter changes, respectively. At low concentrations (0.5 and 1 mM), PB-responsive genes included the well-recognized CAR- and PXR-dependent responsive cytochromes P450 (CYP2B6, CYP3A4), sulfotransferase 2A1 and plasma transporters (ABCB1, ABCC2), as well as a number of genes critically involved in various metabolic pathways, including lipid (CYP4A11, CYP4F3), vitamin D (CYP24A1) and bile (CYP7A1 and CYP8B1) metabolism. At concentrations of 3.2 mM or higher after 20 h, and especially 48 h, increased cytotoxic effects were associated with disregulation of numerous genes related to oxidative stress, DNA repair and apoptosis. Primary human hepatocyte cultures were also exposed to 1 and 3.2 mM PB for 20 h and the changes were comparable to those found in HepaRG cells treated under the same conditions. Taken altogether, our data provide further evidence that HepaRG cells closely resemble primary human hepatocytes and provide new information on the effects of PB in human liver. These data also emphasize the importance of investigating dose- and time-dependent effects of chemicals when using toxicogenomic approaches.

Evolving concepts in liver tissue modeling and implications for in vitro toxicology

The development of human cell models stably expressing functional properties of the in vivo cells they are derived from for predicting toxicity of chemicals is a major challenge. For mimicking the liver, a major target of toxic chemicals, primary hepatocytes represent the most pertinent model. Their use is limited by interdonor functional variability and early phenotypic changes although their lifespan can be extended not only by culturing in a 2D dimension under sophisticated conditions but also by the use of synthetic and natural scaffolds as 3D supporting templates that allow cells to have a more stable microenvironment. Hepatocytes derived from stem cells could be the most appropriate alternative but up to now only liver progenitors/hepatoblasts are obtained in vitro. A few hepatocyte cell lines have retained a variable set of liver-specific functions. Among them are the human hepatoma HepaRG cells that express drug metabolism capacity at levels close to those found in primary hepatocytes making them a suitable model for both acute and chronic toxicity studies. New screening strategies are now proposed based on miniaturized and automated systems; they include the use of microfluidic chips and cell chips coupled with high content imaging analysis. Toxicogenomics technologies (particularly toxicotranscriptomics) have emerged as promising in vitro approaches for better identification and discrimination of cellular responses to chemicals. They should allow to discriminate compounds on the basis of the identification of a set of markers and/specific signaling pathways.

Topoisomerase II inhibition involves characteristic chromosomal expression patterns

BACKGROUND

The phenomenon of co-localization of transcriptionally upregulated genes showing similar expression levels is known across all eukaryotic genomes. We recently mapped the Aroclor 1254-regulated transcriptome back onto the genome and provided evidence for the statistically significant co-localization of regulated genes. They did, however, not always show similar expression levels, and many of the regulated genes were, in fact, repressed.

RESULTS

In this study, we were able to reproduce this observation with microarray data stemming from 1) human hepatocytes treated with the gyrase and potential topoisomerase II inhibitor trovafloxacin, 2) human hepatocytes treated with the topoisomerase II inhibitor doxorubicin and 3) mouse lymphoma cells treated with the topoisomerase II inhibitor etoposide. We found statistically significant co-localization of regulated gene pairs--induced and repressed--within the window size of 0-100 kbp. Notably, by using microarray data stemming from lung tissue of a mouse transgenic line overexpressing the transcription factor c-myc, which served as a negative control, we found regulated genes to be located with regard to each other nearly in the same way as genes distributed randomly all over the genome (0-100 kbp).

CONCLUSION

We suggest topoisomerase II inhibition by Aroclor 1254, trovafloxacin, doxorubicin, and etoposide to be responsible for significant co-localization of regulated genes through the inability of the stabilized enzyme complexes to religate DNA. Within the permanently opened chromatin domains, neighbored genes might be allowed to be regulated. Overexpression of c-myc, however, does not inhibit topoisomerase II activity. Consequently, the enzyme is able to perform its normal function of transiently breaking and rejoining the DNA double strand. As a result, exclusively target genes are regulated.

Tissue-specific expression of aryl hydrocarbon receptor and putative developmental regulatory modules in Baltic salmon yolk-sac fry

The aryl hydrocarbon receptor (AhR) is an ancient protein that is conserved in vertebrates and invertebrates, indicating its important function throughout evolution. AhR has been studied largely because of its role in toxicology-gene expression via AhR is induced by many aromatic hydrocarbons in mammals. Recently, however, it has become clear that AhR is involved in various aspects of development such as cell proliferation and differentiation, and cell motility and migration. The mechanisms by which AhR regulates these various functions remain poorly understood. Across-species comparative studies of AhR in invertebrates, non-mammalian vertebrates and mammals may help to reveal the multiple functions of AhR. Here, we have studied AhR during larval development of Baltic salmon (Salmon salar). Our results indicate that AhR protein is expressed in nervous system, liver and muscle tissues. We also present putative regulatory modules and module-matching genes, produced by chromatin immunoprecipitation (ChIP) cloning and in silico analysis, which may be associated with evolutionarily conserved functions of AhR during development. For example, the module NFKB-AHRR-CREB found from salmon ChIP sequences is present in human ULK3 (regulating formation of granule cell axons in mouse and axon outgrowth in Caernohabditis elegans) and SRGAP1 (GTPase-activating protein involved in the Slit/Robo pathway) promoters. We suggest that AhR may have an evolutionarily conserved role in neuronal development and nerve cell targeting, and in Wnt signaling pathway.

Two polymorphisms of presenilin-2 gene (PSEN2) 5' regulatory region are not associated with Alzheimer's disease (AD) in the Polish population

Presenilin 2 gene (PSEN2) is one of the causative genes for familial Alzheimer's disease. A delA polymorphism located in PSEN2 promoter was proposed to be a risk factor for early-onset AD. We examined association between AD and PSEN2 polymorphisms located in two 5'UTR regions in group of 217 late-onset AD patients, 109 mild cognitive impairment patients, and 225 non-demented control subjects. No significant differences for genotype and allele distributions of a delA and a novel insAC polymorphisms in the studied groups as compared to controls were observed. Univariate and multivariate risk estimation shows that neither delA, insAC alleles nor the genotypes are risk factors for AD. No significant interaction between the APOE4 and PSEN2 polymorphisms was found. A bioinformatic analysis showed that delA polymorphism influences binding sites of transcription factors involved in the cellular processes related to AD. The rare variants identified in exon 3 of the PSEN2 could have a potential influence on PSEN2 transcript splicing.

Role of AHR2 in the expression of novel cytochrome P450 1 family genes, cell cycle genes, and morphological defects in developing zebra fish exposed to 3,3',4,4',5-pentachlorobiphenyl or 2,3,7,8-tetrachlorodibenzo-p-dioxin

Halogenated agonists for the aryl hydrocarbon receptor (AHR), such as 3,3',4,4',5-pentachlorobiphenyl (PCB126) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), cause developmental toxicity in fish. AHR dependence of these effects is known for TCDD but only presumed for PCB126, and the AHR-regulated genes involved are known only in part. We defined the role of AHR in regulation of four cytochrome P450 1 (CYP1) genes and the effect of PCB126 on cell cycle genes (i.e., PCNA and cyclin E) in zebra fish (Danio rerio) embryos. Basal and PCB126-induced expression of CYP1A, CYP1B1, CYP1C1, and CYP1C2 was examined over time as well as in relation to cell cycle gene expression and morphological effects of PCB126 in developing zebra fish. The four CYP1 genes differed in the time for maximal basal and induced expression, i.e., CYP1B1 peaked within 2 days postfertilization (dpf), the CYP1Cs around hatching (3 dpf), and CYP1A after hatching (14-21 dpf). These results indicate developmental periods when the CYP1s may play physiological roles. PCB126 (0.3-100nM) caused concentration-dependent CYP1 gene induction (EC50: 1.4-2.7nM, Lowest observed effect concentration [LOEC]: 0.3-1nM) and pericardial edema (EC50: 4.4nM, LOEC: 3nM) in 3-dpf embryos. Blockage of AHR2 translation significantly inhibited these effects of PCB126 and TCDD. PCNA gene expression was reduced by PCB126 in a concentration-dependent manner, suggesting that PCB126 could suppress cell proliferation. Our results indicate that the four CYP1 genes examined are regulated by AHR2 and that the effect of PCB126 on morphology in zebra fish embryos is AHR2 dependent. Moreover, the developmental patterns of expression and induction suggest that CYP1 enzymes could function in normal development and in developmental toxicity of PCB126 in fish embryos.

Interplay of genes regulated by estrogen and diindolylmethane in breast cancer cell lines

Diindolylmethane (DIM), a biologically active congener of indole-3-carbinol (I3C) derived from cruciferous vegetables, is a promising agent for the prevention of estrogen-sensitive cancers. Both DIM and estrogen affect transcription of genes by binding receptors, such as aryl hydrocarbon receptor (AhR) or estrogen receptors (ER). Gene regulation by DIM and estradiol (E2) can be very complex. While DIM typically binds the AhR, this complex can directly associate with the ER, recruit co-activators that bind to estrogen-responsive promoters, and activate transcription. Alternately, DIM can bind the ER directly. In this study, we have analyzed gene expression using microarray profiling and quantitative real time-polymerase chain reaction in MCF7 breast cancer cells treated with E2 (1 nM) or DIM (25 microM) alone or in combination for 16 h. The interplay of E2 and DIM was reflected in the expression of a subset of genes (<90) in which the combination of E2 and DIM acted either additively or antagonistically to alter gene expression.