BackgroundInduction 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.ResultsExposure 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.ConclusionsCollectively, 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.

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.

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

No abstract

Evidence suggests that epigenetic perturbations are involved in the adverse effects associated with some drugs and toxicants, including certain classes of non-genotoxic carcinogens. Such epigenetic changes (altered DNA methylation and covalent histone modifications) may take place at the earliest stages of carcinogenesis and their identification holds great promise for biomedical research. Here, we evaluate the sensitivity and specificity of genome-wide epigenomic and transcriptomic profiling in phenobarbital (PB)-treated B6C3F1 mice, a well-characterized rodent model of non-genotoxic liver carcinogenesis. Methylated DNA Immunoprecipitation (MeDIP)-coupled microarray profiling of 17,967 promoter regions and 4,566 intergenic CpG islands was combined with genome-wide mRNA expression profiling to identify liver tissue-specific PB-mediated DNA methylation and transcriptional alterations. Only a limited number of significant anti-correlations were observed between PB-induced transcriptional and promoter-based DNA methylation perturbations. However, the constitutive androstane receptor (CAR) target gene Cyp2b10 was found to be concomitantly hypomethylated and transcriptionally activated in a liver tissue-specific manner following PB treatment. Furthermore, analysis of active and repressive histone modifications using chromatin immunoprecipitation revealed a strong PB-mediated epigenetic switch at the Cyp2b10 promoter. Our data reveal that PB-induced transcriptional perturbations are not generally associated with broad changes in the DNA methylation status at proximal promoters and suggest that the drug-inducible CAR pathway regulates an epigenetic switch from repressive to active chromatin at the target gene Cyp2b10. This study demonstrates the utility of integrated epigenomic and transcriptomic profiling for elucidating early mechanisms and biomarkers of non-genotoxic carcinogenesis.

The molecular events during nongenotoxic carcinogenesis and their temporal order are poorly understood but thought to include long-lasting perturbations of gene expression. Here, we have investigated the temporal sequence of molecular and pathological perturbations at early stages of phenobarbital (PB) mediated liver tumor promotion in vivo. Molecular profiling (mRNA, microRNA [miRNA], DNA methylation, and proteins) of mouse liver during 13 weeks of PB treatment revealed progressive increases in hepatic expression of long noncoding RNAs and miRNAs originating from the Dlk1-Dio3 imprinted gene cluster, a locus that has recently been associated with stem cell pluripotency in mice and various neoplasms in humans. PB induction of the Dlk1-Dio3 cluster noncoding RNA (ncRNA) Meg3 was localized to glutamine synthetase-positive hypertrophic perivenous hepatocytes, suggesting a role for β-catenin signaling in the dysregulation of Dlk1-Dio3 ncRNAs. The carcinogenic relevance of Dlk1-Dio3 locus ncRNA induction was further supported by in vivo genetic dependence on constitutive androstane receptor and β-catenin pathways. Our data identify Dlk1-Dio3 ncRNAs as novel candidate early biomarkers for mouse liver tumor promotion and provide new opportunities for assessing the carcinogenic potential of novel compounds.

This paper reports an analysis of the encoded proteins (the proteome) of the genomes of human, fly, worm, yeast, and representatives of bacteria and archaea in terms of the three-dimensional structures of their globular domains together with a general sequence-based study. We show that 39% of the human proteome can be assigned to known structures. We estimate that for 77% of the proteome, there is some functional annotation, but only 26% of the proteome can be assigned to standard sequence motifs that characterize function. Of the human protein sequences, 13% are transmembrane proteins, but only 3% of the residues in the proteome form membrane-spanning regions. There are substantial differences in the composition of globular domains of transmembrane proteins between the proteomes we have analyzed. Commonly occurring structural superfamilies are identified within the proteome. The frequencies of these superfamilies enable us to estimate that 98% of the human proteome evolved by domain duplication, with four of the 10 most duplicated superfamilies specific for multicellular organisms. The zinc-finger superfamily is massively duplicated in human compared to fly and worm, and occurrence of domains in repeats is more common in metazoa than in single cellular organisms. Structural superfamilies over-and underrepresented in human disease genes have been identified. Data and results can be downloaded and analyzed via web-based applications at http://www.sbg. bio.ic.ac.uk.[Supplemental material is available online at http://www.genome.org.]The interpretation and exploitation of the wealth of biological knowledge that can be derived from the human genome (Lander et al. 2001;Venter et al. 2001) requires an analysis of the three-dimensional structures and the functions of the encoded proteins (the proteome). Comparison of this analysis with those of other eukaryotic and prokaryotic proteomes will identify which structural and functional features are common and which confer species specificity. In this paper, we present an integrated analysis of the proteomes of human and 13 other species considering the folds of globular domains, the presence of transmembrane proteins, and the extent to which the proteomes can be functionally annotated. This integrated approach enables us to consider the relationship between these different aspects of annotation and thereby enhance previous analyses of the human and other proteomes (e.g., Koonin et al. 2000;Frishman et al. 2001;Iliopoulos et al. 2001), including the seminal papers reporting the human genome sequence (Lander et al. 2001;Venter et al. 2001).A widely used first step in a bioinformatics-based functional annotation is to identify known sequence motifs and domains from manually curated databases such as PFAM/ INTERPRO (Bateman et al. 2000) and PANTHER (Venter et al. 2001). This strategy was used in the original analyses of the human proteome (Lander et al. 2001;Venter et al. 2001). These annotations tend to be reliable, as these libraries have been carefully constructed to avoid ...

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.