Dynamic changes in 5-hydroxymethylation signatures underpin early and late events in drug exposed liver

Aberrant DNA methylation is a common feature of neoplastic lesions, and early detection of such changes may provide powerful mechanistic insights and biomarkers for carcinogenesis. Here, we investigate dynamic changes in the mouse liver DNA methylome associated with short (1 day) and prolonged (7, 28 and 91 days) exposure to the rodent liver non-genotoxic carcinogen, phenobarbital (PB). We find that the distribution of 5mC/5hmC is highly consistent between untreated individuals of a similar age; yet, changes during liver maturation in a transcriptionally dependent manner. Following drug treatment, we identify and validate a series of differentially methylated or hydroxymethylated regions: exposure results in staged transcriptional responses with distinct kinetic profiles that strongly correlate with promoter proximal region 5hmC levels. Furthermore, reciprocal changes for both 5mC and 5hmC in response to PB suggest that active demethylation may be taking place at each set of these loci via a 5hmC intermediate. Finally, we identify potential early biomarkers for non-genotoxic carcinogenesis, including several genes aberrantly expressed in liver cancer. Our work suggests that 5hmC profiling can be used as an indicator of cell states during organ maturation and drug-induced responses and provides novel epigenetic signatures for non-genotoxic carcinogen exposur

Computational modeling identifies key gene regulatory interactions underlying phenobarbital-mediated tumor promotion

Gene regulatory interactions underlying the early stages of non-genotoxic carcinogenesis are poorly understood. Here, we have identified key candidate regulators of phenobarbital (PB)-mediated mouse liver tumorigenesis, a well-characterized model of non-genotoxic carcinogenesis, by applying a new computational modeling approach to a comprehensive collection of in vivo gene expression studies. We have combined our previously developed motif activity response analysis (MARA), which models gene expression patterns in terms of computationally predicted transcription factor binding sites with singular value decomposition (SVD) of the inferred motif activities, to disentangle the roles that different transcriptional regulators play in specific biological pathways of tumor promotion. Furthermore, transgenic mouse models enabled us to identify which of these regulatory activities was downstream of constitutive androstane receptor and β-catenin signaling, both crucial components of PB-mediated liver tumorigenesis. We propose novel roles for E2F and ZFP161 in PB-mediated hepatocyte proliferation and suggest that PB-mediated suppression of ESR1 activity contributes to the development of a tumor-prone environment. Our study shows that combining MARA with SVD allows for automated identification of independent transcription regulatory programs within a complex in vivo tissue environment and provides novel mechanistic insights into PB-mediated hepatocarcinogenesi

Comparative analysis of affinity-based 5-hydroxymethylation enrichment techniques

The epigenetic modification of 5-hydroxymethylcytosine (5hmC) is receiving great attention due to its potential role in DNA methylation reprogramming and as a cell state identifier. Given this interest, it is important to identify reliable and cost-effective methods for the enrichment of 5hmC marked DNA for downstream analysis. We tested three commonly used affinity-based enrichment techniques; (i) antibody, (ii) chemical capture and (iii) protein affinity enrichment and assessed their ability to accurately and reproducibly report 5hmC profiles in mouse tissues containing high (brain) and lower (liver) levels of 5hmC. The protein-affinity technique is a poor reporter of 5hmC profiles, delivering 5hmC patterns that are incompatible with other methods. Both antibody and chemical capture-based techniques generate highly similar genome-wide patterns for 5hmC, which are independently validated by standard quantitative PCR (qPCR) and glucosyl-sensitive restriction enzyme digestion (gRES-qPCR). Both antibody and chemical capture generated profiles reproducibly link to unique chromatin modification profiles associated with 5hmC. However, there appears to be a slight bias of the antibody to bind to regions of DNA rich in simple repeats. Ultimately, the increased specificity observed with chemical capture-based approaches makes this an attractive method for the analysis of locus-specific or genome-wide patterns of 5hm

The Need to Decide If All Estrogens Are Intrinsically Similar

We used gene expression profiling to investigate whether the molecular effects induced by estrogens of different provenance are intrinsically similar. In this article we show that the physiologic estrogen 17β-estradiol, the phytoestrogen genistein, and the synthetic estrogen diethylstilbestrol alter the expression of the same 179 genes in the intact immature mouse uterus under conditions where each chemical has produced an equivalent gravimetric and histologic uterotrophic effect, using the standard 3-day assay protocol. Data are also presented indicating the limitations associated with comparison of gene expression profiles for different chemicals at times before the uterotrophic effects are fully realized. We conclude that the case has yet to be made for regarding synthetic estrogens as presenting a unique human hazard compared with phytoestrogens and physiologic estrogens

Constitutive androstane receptor 1 is constitutively bound to chromatin and ‘primed’ for transactivation in hepatocytes

The constitutive androstane receptor (CAR) is a xenobiotic sensor expressed in hepatocytes that activates genes involved in drug metabolism, lipid homeostasis, and cell proliferation. Much progress has been made in understanding the mechanism of activation of human CAR by drugs and xenobiotics. However, many aspects of the activation pathway remain to be elucidated. In this report, we have used viral constructs to express human CAR, its splice variants, and mutant CAR forms in hepatocytes from Car-/- mice in vitro and in vivo. We demonstrate CAR expression rescued the ability of Car-/- hepatocytes to respond to a wide range of CAR activators including phenobarbital. Additionally, two major splice isoforms of human CAR, CAR2 and CAR3, were inactive with almost all the agents tested. In contrast to the current model of CAR activation, ectopic CAR1 is constitutively localised in the nucleus and is loaded onto Cyp2b10 gene in the absence of an inducing agent. In studies to elucidate the role of threonine T38 in CAR regulation, we found that the T38D mutant was inactive even in the presence of CAR activators. However, the T38A mutant was activated by CAR inducers, showing that T38 is not essential for CAR activation. Also, using the inhibitor erlotinib, we could not confirm a role for the epidermal growth factor receptor in CAR regulation. Our data suggest that CAR is constitutively bound to gene regulatory regions and is regulated by exogenous agents through a mechanism which involves protein phosphorylation in the nucleus

Genomic and protein expression analysis reveals flap endonuclease 1 (FEN1) as a key biomarker in breast and ovarian cancer

FEN1 has key roles in Okazaki fragment maturation during replication, long patch base excision repair, rescue of stalled replication forks, maintenance of telomere stability and apoptosis. FEN1 may be dysregulated in breast and ovarian cancers and have clinicopathological significance in patients. We comprehensively investigated FEN1 mRNA expression in multiple cohorts of breast cancer [training set (128), test set (249), external validation (1952)]. FEN1 protein expression was evaluated in 568 oestrogen receptor (ER) negative breast cancers, 894 ER positive breast cancers and 156 ovarian epithelial cancers. FEN1 mRNA overexpression was highly significantly associated with high grade (p= 4.89 x 10 - 57) , high mitotic index (p= 5.25 x 10 - 28), pleomorphism (p= 6.31 x 10-19), ER negative (p= 9.02 x 10-35 ), PR negative (p= 9.24 x 10-24 ), triple negative phenotype (p= 6.67 x 10-21) , PAM50.Her2 (p=5.19 x 10-13 ), PAM50.Basal (p=2.7 x 10-41), PAM50.LumB (p=1.56 x 10-26), integrative molecular cluster 1 (intClust.1) ( p=7.47 x 10-12), intClust.5 (p=4.05 x 10-12) and intClust. 10 (p=7.59 x 10-38 ) breast cancers. FEN1 mRNA overexpression is associated with poor breast cancer specific survival in univariate (p=4.4 x 10-16) and multivariate analysis (p=9.19 x 10-7). At the protein level, in ER positive tumours , FEN1 overexpression remains significantly linked to high grade, high mitotic index and pleomorphism (ps< 0.01). In ER negative tumours, high FEN1 is significantly associated with pleomorphism, tumour type, lymphovascular invasion, triple negative phenotype, EGFR and HER2 expression (ps<0.05). In ER positive as well as in ER negative tumours, FEN1 protein over expression is associated with poor survival in univariate and multivariate analysis (ps<0.01). In ovarian epithelial cancers , similarly, FEN1 overexpression is associated with high grade, high stage and poor survival (ps<0.05). We conclude that FEN1 is a promising biomarker in breast and ovarian epithelial cancer

Non-genotoxic carcinogen exposure induces defined changes in the 5-hydroxymethylome.

BACKGROUND: Induction and promotion of liver cancer by exposure to non-genotoxic carcinogens coincides with epigenetic perturbations, including specific changes in DNA methylation. Here we investigate the genome-wide dynamics of 5-hydroxymethylcytosine (5hmC) as a likely intermediate of 5-methylcytosine (5mC) demethylation in a DNA methylation reprogramming pathway. We use a rodent model of non-genotoxic carcinogen exposure using the drug phenobarbital. RESULTS: Exposure to phenobarbital results in dynamic and reciprocal changes to the 5mC/5hmC patterns over the promoter regions of a cohort of genes that are transcriptionally upregulated. This reprogramming of 5mC/5hmC coincides with characteristic changes in the histone marks H3K4me2, H3K27me3 and H3K36me3. Quantitative analysis of phenobarbital-induced genes that are involved in xenobiotic metabolism reveals that both DNA modifications are lost at the transcription start site, while there is a reciprocal relationship between increasing levels of 5hmC and loss of 5mC at regions immediately adjacent to core promoters. CONCLUSIONS: Collectively, these experiments support the hypothesis that 5hmC is a potential intermediate in a demethylation pathway and reveal precise perturbations of the mouse liver DNA methylome and hydroxymethylome upon exposure to a rodent hepatocarcinogen.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Phenotypic Anchoring of Gene Expression Changes during Estrogen-Induced Uterine Growth

A major challenge in the emerging field of toxicogenomics is to define the relationships between chemically induced changes in gene expression and alterations in conventional toxicologic parameters such as clinical chemistry and histopathology. We have explored these relationships in detail using the rodent uterotrophic assay as a model system. Gene expression levels, uterine weights, and histologic parameters were analyzed 1, 2, 4, 8, 24, 48, and 72 hr after exposure to the reference physiologic estrogen 17β-estradiol (E(2)). A multistep analysis method, involving unsupervised hierarchical clustering followed by supervised gene ontology–driven clustering, was used to define the transcriptional program associated with E(2)-induced uterine growth and to identify groups of genes that may drive specific histologic changes in the uterus. This revealed that uterine growth and maturation are preceded and accompanied by a complex, multistage molecular program. The program begins with the induction of genes involved in transcriptional regulation and signal transduction and is followed, sequentially, by the regulation of genes involved in protein biosynthesis, cell proliferation, and epithelial cell differentiation. Furthermore, we have identified genes with common molecular functions that may drive fluid uptake, coordinated cell division, and remodeling of luminal epithelial cells. These data define the mechanism by which an estrogen induces organ growth and tissue maturation, and demonstrate that comparison of temporal changes in gene expression and conventional toxicology end points can facilitate the phenotypic anchoring of toxicogenomic data