A molecular glue degrader of the WIZ transcription factor for fetal hemoglobin induction.

Sickle cell disease (SCD) is a prevalent, life-threatening condition attributable to a heritable mutation in ß-hemoglobin. Therapeutic induction of fetal hemoglobin (HbF) can ameliorate disease complications and has been intently pursued. However, safe and effective small-molecule inducers of HbF remain elusive. We report the discovery of dWIZ-1 and dWIZ-2, molecular glue degraders of the WIZ transcription factor that robustly induce HbF in erythroblasts. Phenotypic screening of a cereblon (CRBN)-biased chemical library revealed WIZ as a previously unknown repressor of HbF. WIZ degradation is mediated by recruitment of WIZ(ZF7) to CRBN by dWIZ-1, as resolved by crystallography of the ternary complex. Pharmacological degradation of WIZ was well tolerated and induced HbF in humanized mice and cynomolgus monkeys. These findings establish WIZ degradation as a globally accessible therapeutic strategy for SCD.

Unaltered hepatic wound healing response in male rats with ancestral liver injury.

The possibility that ancestral environmental exposure could result in adaptive inherited effects in mammals has been long debated. Numerous rodent models of transgenerational responses to various environmental factors have been published but due to technical, operational and resource burden, most still await independent confirmation. A previous study reported multigenerational epigenetic adaptation of the hepatic wound healing response upon exposure to the hepatotoxicant carbon tetrachloride (CCl4) in male rats. Here, we comprehensively investigate the transgenerational effects by repeating the original CCl4 multigenerational study with increased power, pedigree tracing, F2 dose-response and suitable randomization schemes. Detailed pathology evaluations do not support adaptive phenotypic suppression of the hepatic wound healing response or a greater fitness of F2 animals with ancestral liver injury exposure. However, transcriptomic analyses identified genes whose expression correlates with ancestral liver injury, although the biological relevance of this apparent transgenerational transmission at the molecular level remains to be determined. This work overall highlights the need for independent evaluation of transgenerational epigenetic inheritance paradigms in mammals.

Wdr5 is essential for fetal erythropoiesis and hematopoiesis.

WDR5 is a highly conserved protein that performs multiple scaffolding functions in the context of chromatin. However, efforts to understand the function of WDR5 in normal tissues physiologically are quite limited so far. In our study, we explored the function of Wdr5 in erythropoiesis and hematopoiesis by using a hematopoietic-specific Wdr5 knockout mouse model. We found that loss of Wdr5 mediated by Vav-iCre leads to embryonic lethality with defective erythropoiesis. In addition, Wdr5-deficiency completely impairs the hematopoietic stem and progenitor cells function and might alter the immunophenotype of these stem cells and progenitors by decreasing c-Kit expression. Collectively, we identified the pivotal role of Wdr5 in fetal hematopoiesis and erythropoiesis as the de novo findings.

Discovery of the Clinical Candidate MAK683: An EED-Directed, Allosteric, and Selective PRC2 Inhibitor for the Treatment of Advanced Malignancies.

Polycomb Repressive Complex 2 (PRC2) plays an important role in transcriptional regulation during animal development and in cell differentiation, and alteration of PRC2 activity has been associated with cancer. On a molecular level, PRC2 catalyzes methylation of histone H3 lysine 27 (H3K27), resulting in mono-, di-, or trimethylated forms of H3K27, of which the trimethylated form H3K27me3 leads to transcriptional repression of polycomb target genes. Previously, we have shown that binding of the low-molecular-weight compound EED226 to the H3K27me3 binding pocket of the regulatory subunit EED can effectively inhibit PRC2 activity in cells and reduce tumor growth in mouse xenograft models. Here, we report the stepwise optimization of the tool compound EED226 toward the potent and selective EED inhibitor MAK683 (compound 22) and its subsequent preclinical characterization. Based on a balanced PK/PD profile, efficacy, and mitigated risk of forming reactive metabolites, MAK683 has been selected for clinical development.

BET bromodomain inhibitors regulate keratinocyte plasticity.

Although most acute skin wounds heal rapidly, non-healing skin ulcers represent an increasing and substantial unmet medical need that urgently requires effective therapeutics. Keratinocytes resurface wounds to re-establish the epidermal barrier by transitioning to an activated, migratory state, but this ability is lost in dysfunctional chronic wounds. Small-molecule regulators of keratinocyte plasticity with the potential to reverse keratinocyte malfunction in situ could offer a novel therapeutic approach in skin wound healing. Utilizing high-throughput phenotypic screening of primary keratinocytes, we identify such small molecules, including bromodomain and extra-terminal domain (BET) protein family inhibitors (BETi). BETi induce a sustained activated, migratory state in keratinocytes in vitro, increase activation markers in human epidermis ex vivo and enhance skin wound healing in vivo. Our findings suggest potential clinical utility of BETi in promoting keratinocyte re-epithelialization of skin wounds. Importantly, this novel property of BETi is exclusively observed after transient low-dose exposure, revealing new potential for this compound class.

Drug-induced chromatin accessibility changes associate with sensitivity to liver tumor promotion.

Liver cancer susceptibility varies amongst humans and between experimental animal models because of multiple genetic and epigenetic factors. The molecular characterization of such susceptibilities has the potential to enhance cancer risk assessment of xenobiotic exposures and disease prevention strategies. Here, using DNase I hypersensitivity mapping coupled with transcriptomic profiling, we investigate perturbations in cis-acting gene regulatory elements associated with the early stages of phenobarbital (PB)-mediated liver tumor promotion in susceptible versus resistant mouse strains (B6C3F1 versus C57BL/6J). Integrated computational analyses of strain-selective changes in liver chromatin accessibility underlying PB response reveal differential epigenetic regulation of molecular pathways associated with PB-mediated tumor promotion, including Wnt/ß-catenin signaling. Complementary transcription factor motif analyses reveal mouse strain-selective gene regulatory networks and a novel role for Stat, Smad, and Fox transcription factors in the early stages of PB-mediated tumor promotion. Mapping perturbations in cis-acting gene regulatory elements provides novel insights into the molecular basis for susceptibility to xenobiotic-induced rodent liver tumor promotion and has the potential to enhance mechanism-based cancer risk assessments of xenobiotic exposures.

IRF2 is a master regulator of human keratinocyte stem cell fate.

Resident adult epithelial stem cells maintain tissue homeostasis by balancing self-renewal and differentiation. The stem cell potential of human epidermal keratinocytes is retained in vitro but lost over time suggesting extrinsic and intrinsic regulation. Transcription factor-controlled regulatory circuitries govern cell identity, are sufficient to induce pluripotency and transdifferentiate cells. We investigate whether transcriptional circuitry also governs phenotypic changes within a given cell type by comparing human primary keratinocytes with intrinsically high versus low stem cell potential. Using integrated chromatin and transcriptional profiling, we implicate IRF2 as antagonistic to stemness and show that it binds and regulates active cis-regulatory elements at interferon response and antigen presentation genes. CRISPR-KD of IRF2 in keratinocytes with low stem cell potential increases self-renewal, migration and epidermis formation. These data demonstrate that transcription factor regulatory circuitries, in addition to maintaining cell identity, control plasticity within cell types and offer potential for therapeutic modulation of cell function.

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 localized 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.

Chromatin dynamics underlying latent responses to xenobiotics.

Pleiotropic xenobiotics can trigger dynamic alterations in mammalian chromatin structure and function but many of these are likely non-adverse and simply reflect short-term changes in DNA transactions underlying normal homeostatic, adaptive and protective cellular responses. However, it is plausible that a subset of xenobiotic-induced perturbations of somatic tissue or germline epigenomes result in delayed-onset and long-lasting adverse effects, in particular if they occur during critical stages of growth and development. These could include reprogramming, dedifferentiation, uncontrolled growth, and cumulative toxicity effects through molecular memory of prior xenobiotic exposures or altered susceptibility to subsequent xenobiotic exposures. Here we discuss the current evidence for epigenetic mechanisms underlying latent responses to xenobiotics, and the potential for identifying molecular epigenetic changes that are prodromal to overt morphologic or functional toxicity phenotypes.

Xenobiotic CAR Activators Induce Dlk1-Dio3 Locus Noncoding RNA Expression in Mouse Liver.

Derisking xenobiotic-induced nongenotoxic carcinogenesis (NGC) represents a significant challenge during the safety assessment of chemicals and therapeutic drugs. The identification of robust mechanism-based NGC biomarkers has the potential to enhance cancer hazard identification. We previously demonstrated Constitutive Androstane Receptor (CAR) and WNT signaling-dependent up-regulation of the pluripotency associated Dlk1-Dio3 imprinted gene cluster noncoding RNAs (ncRNAs) in the liver of mice treated with tumor-promoting doses of phenobarbital (PB). Here, we have compared phenotypic, transcriptional ,and proteomic data from wild-type, CAR/PXR double knock-out and CAR/PXR double humanized mice treated with either PB or chlordane, and show that hepatic Dlk1-Dio3 locus long ncRNAs are upregulated in a CAR/PXR-dependent manner by two structurally distinct CAR activators. We further explored the specificity of Dlk1-Dio3 locus ncRNAs as hepatic NGC biomarkers in mice treated with additional compounds working through distinct NGC modes of action. We propose that up-regulation of Dlk1-Dio3 cluster ncRNAs can serve as an early biomarker for CAR activator-induced nongenotoxic hepatocarcinogenesis and thus may contribute to mechanism-based assessments of carcinogenicity risk for chemicals and novel therapeutics.

Optimization of 3-Pyrimidin-4-yl-oxazolidin-2-ones as Allosteric and Mutant Specific Inhibitors of IDH1.

High throughput screening and subsequent hit validation identified 4-isopropyl-3-(2-((1-phenylethyl)amino)pyrimidin-4-yl)oxazolidin-2-one as a potent inhibitor of IDH1R132H. Synthesis of the four separate stereoisomers identified the (S,S)-diastereomer (IDH125, 1f) as the most potent isomer. This also showed reasonable cellular activity and excellent selectivity vs IDH1wt. Initial structure-activity relationship exploration identified the key tolerances and potential for optimization. X-ray crystallography identified a functionally relevant allosteric binding site amenable to inhibitors, which can penetrate the blood-brain barrier, and aided rational optimization. Potency improvement and modulation of the physicochemical properties identified (S,S)-oxazolidinone IDH889 (5x) with good exposure and 2-HG inhibitory activity in a mutant IDH1 xenograft mouse model.

Loss of Tet1-Associated 5-Hydroxymethylcytosine Is Concomitant with Aberrant Promoter Hypermethylation in Liver Cancer.

Aberrant hypermethylation of CpG islands (CGI) in human tumors occurs predominantly at repressed genes in the host tissue, but the preceding events driving this phenomenon are poorly understood. In this study, we temporally tracked epigenetic and transcriptomic perturbations that occur in a mouse model of liver carcinogenesis. Hypermethylated CGI events in the model were predicted by enrichment of the DNA modification 5-hydroxymethylcytosine (5hmC) and the histone H3 modification H3K27me3 at silenced promoters in the host tissue. During cancer progression, selected CGIs underwent hypo-hydroxymethylation prior to hypermethylation, while retaining H3K27me3. In livers from mice deficient in Tet1, a tumor suppressor involved in cytosine demethylation, we observed a similar loss of promoter core 5hmC, suggesting that reduced Tet1 activity at CGI may contribute to epigenetic dysregulation during hepatocarcinogenesis. Consistent with this possibility, mouse liver tumors exhibited reduced Tet1 protein levels. Similar to humans, DNA methylation changes at CGI in mice did not appear to be direct drivers of hepatocellular carcinoma progression, rather, dynamic changes in H3K27me3 promoter deposition correlated strongly with tumor-specific activation and repression of transcription. Overall, our results suggest that loss of promoter-associated 5hmC in liver tumors licenses reprograming of DNA methylation at silent CGI during progression. Cancer Res; 76(10); 3097-108. ©2016 AACR.

Reciprocal changes in DNA methylation and hydroxymethylation and a broad repressive epigenetic switch characterize FMR1 transcriptional silencing in fragile X syndrome.

BACKGROUND: Fragile X syndrome (FXS) is the most common form of inherited intellectual disability, resulting from the loss of function of the fragile X mental retardation 1 (FMR1) gene. The molecular pathways associated with FMR1 epigenetic silencing are still elusive, and their characterization may enhance the discovery of novel therapeutic targets as well as the development of novel clinical biomarkers for disease status. RESULTS: We have deployed customized epigenomic profiling assays to comprehensively map the FMR1 locus chromatin landscape in peripheral mononuclear blood cells (PBMCs) from eight FXS patients and in fibroblast cell lines derived from three FXS patient. Deoxyribonucleic acid (DNA) methylation (5-methylcytosine (5mC)) and hydroxymethylation (5-hydroxymethylcytosine (5hmC)) profiling using methylated DNA immunoprecipitation (MeDIP) combined with a custom FMR1 microarray identifies novel regions of DNA (hydroxy)methylation changes within the FMR1 gene body as well as in proximal flanking regions. At the region surrounding the FMR1 transcriptional start sites, increased levels of 5mC were associated to reciprocal changes in 5hmC, representing a novel molecular feature of FXS disease. Locus-specific validation of FMR1 5mC and 5hmC changes highlighted inter-individual differences that may account for the expected DNA methylation mosaicism observed at the FMR1 locus in FXS patients. Chromatin immunoprecipitation (ChIP) profiling of FMR1 histone modifications, together with 5mC/5hmC and gene expression analyses, support a functional relationship between 5hmC levels and FMR1 transcriptional activation and reveal cell-type specific differences in FMR1 epigenetic regulation. Furthermore, whilst 5mC FMR1 levels positively correlated with FXS disease severity (clinical scores of aberrant behavior), our data reveal for the first time an inverse correlation between 5hmC FMR1 levels and FXS disease severity. CONCLUSIONS: We identify novel, cell-type specific, regions of FMR1 epigenetic changes in FXS patient cells, providing new insights into the molecular mechanisms of FXS. We propose that the combined measurement of 5mC and 5hmC at selected regions of the FMR1 locus may significantly enhance FXS clinical diagnostics and patient stratification.

Evaluation of 5-methylcytosine and 5-hydroxymethylcytosine as potential biomarkers for characterisation of chemical allergens.

Epigenetic regulation of gene expression plays a pivotal role in the orchestration of immune responses. Chemical allergens form two categories: skin sensitizing chemicals associated with allergic contact dermatitis, and chemicals that cause sensitization of the respiratory tract and occupational asthma. In mice these are characterized by different T helper (Th) cell responses. Changes in DNA methylation in particular have been implicated in the in vivo responses to chemical allergy. As such it was hypothesised that differentially methylated regions (DMR) may provide candidates biomarkers of chemical allergy To examine this, mice were exposed to 2,4-dinitrochlorobenzene (DNCB; a contact allergen) or trimellitic anhydride (TMA; a respiratory allergen). DNA from draining lymph nodes was processed for methylated (5mC) and hydroxymethylated (5hmC) DNA immunoprecipitation (MeDIP/hMeDIP) then selected DMR analysed by qPCR. We describe a number of DMRs which, by combined analysis of 5mC and 5hmC, differentiate between responses induced by DNCB and those by TMA. Furthermore, these changes in methylation are specific to the draining lymph node. The Gmpr DMR is suggested as a possible biomarker for contact allergen-induced immune responses; it is characterised by divergent levels of 5mC and 5hmC DNCB-treated mice only. In contrast, the Nwc DMR was characterised by divergent 5mC and 5hmC specifically in response to TMA, highlighting its possible utility as a biomarker for responses induced by chemical respiratory allergens. These data not only represent novel analysis of 5hmC in response to chemical allergy in vivo, but with further investigation, may also provide a possible basis for differentiation between classes of chemical allergens.

Phenobarbital induces cell cycle transcriptional responses in mouse liver humanized for constitutive androstane and pregnane x receptors.

The constitutive androstane receptor (CAR) and the pregnane X receptor (PXR) are closely related nuclear receptors involved in drug metabolism and play important roles in the mechanism of phenobarbital (PB)-induced rodent nongenotoxic hepatocarcinogenesis. Here, we have used a humanized CAR/PXR mouse model to examine potential species differences in receptor-dependent mechanisms underlying liver tissue molecular responses to PB. Early and late transcriptomic responses to sustained PB exposure were investigated in liver tissue from double knock-out CAR and PXR (CAR(KO)-PXR(KO)), double humanized CAR and PXR (CAR(h)-PXR(h)), and wild-type C57BL/6 mice. Wild-type and CAR(h)-PXR(h) mouse livers exhibited temporally and quantitatively similar transcriptional responses during 91 days of PB exposure including the sustained induction of the xenobiotic response gene Cyp2b10, the Wnt signaling inhibitor Wisp1, and noncoding RNA biomarkers from the Dlk1-Dio3 locus. Transient induction of DNA replication (Hells, Mcm6, and Esco2) and mitotic genes (Ccnb2, Cdc20, and Cdk1) and the proliferation-related nuclear antigen Mki67 were observed with peak expression occurring between 1 and 7 days PB exposure. All these transcriptional responses were absent in CAR(KO)-PXR(KO) mouse livers and largely reversible in wild-type and CAR(h)-PXR(h) mouse livers following 91 days of PB exposure and a subsequent 4-week recovery period. Furthermore, PB-mediated upregulation of the noncoding RNA Meg3, which has recently been associated with cellular pluripotency, exhibited a similar dose response and perivenous hepatocyte-specific localization in both wild-type and CAR(h)-PXR(h) mice. Thus, mouse livers coexpressing human CAR and PXR support both the xenobiotic metabolizing and the proliferative transcriptional responses following exposure to PB.

Chemical allergen induced perturbations of the mouse lymph node DNA methylome.

Epigenetic regulation of gene expression plays a pivotal role in the orchestration of immune responses and may determine the vigor, quality, or longevity of such responses. Chemical allergens can be divided into two categories: skin sensitizing chemicals associated with allergic contact dermatitis, and chemicals that cause sensitization of the respiratory tract and occupational asthma. In mice, these are characterized by different T helper cell responses. To explore the regulation and maintenance of these divergent responses, mice were exposed to 2,4-dinitrochlorobenzene (DNCB, a contact allergen) or trimellitic anhydride (TMA, a respiratory allergen). DNA from draining lymph nodes was processed for methylated DNA immunoprecipitation followed by hybridization to a whole-genome DNA promoter array. 6319 differently methylated regions (DMRs) were identified following DNCB treatment, whereas 2178 DMRs were measured following TMA treatment, with approximately half of the TMA DMRs common to DNCB. When limited to promoter region-associated DMRs, 637 genes were uniquely associated with DNCB-induced DMRs but only 164 genes were unique to TMA DMRs. Promoter-associated DMRs unique to either DNCB or TMA were generally hypomethylated whereas DMRs common to both allergens tended to be hypermethylated. Pathway analyses highlighted a number of immune-related pathways, including chemokine and cytokine signaling. These data demonstrate that chemical allergen exposure results in characteristic patterns of DNA methylation indicative of epigenetic regulation of the allergic response.

Ha-ras and ß-catenin oncoproteins orchestrate metabolic programs in mouse liver tumors.

The process of hepatocarcinogenesis in the diethylnitrosamine (DEN) initiation/phenobarbital (PB) promotion mouse model involves the selective clonal outgrowth of cells harboring oncogene mutations in Ctnnb1, while spontaneous or DEN-only-induced tumors are often Ha-ras- or B-raf-mutated. The molecular mechanisms and pathways underlying these different tumor sub-types are not well characterized. Their identification may help identify markers for xenobiotic promoted versus spontaneously occurring liver tumors. Here, we have characterized mouse liver tumors harboring either Ctnnb1 or Ha-ras mutations via integrated molecular profiling at the transcriptional, translational and post-translational levels. In addition, metabolites of the intermediary metabolism were quantified by high resolution (1)H magic angle nuclear magnetic resonance. We have identified tumor genotype-specific differences in mRNA and miRNA expression, protein levels, post-translational modifications, and metabolite levels that facilitate the molecular and biochemical stratification of tumor phenotypes. Bioinformatic integration of these data at the pathway level led to novel insights into tumor genotype-specific aberrant cell signaling and in particular to a better understanding of alterations in pathways of the cell intermediary metabolism, which are driven by the constitutive activation of the ß-Catenin and Ha-ras oncoproteins in tumors of the two genotypes.

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 hepatocarcinogenesis.

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 5hmC.

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 exposure.