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

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 which 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, CGI underwent hypo-hydroxymethylation prior to hypermethylation, whilst 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 observed 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 reprogramming of DNA methylation at silent CGI during progression

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