Multiple genes exhibit phenobarbital-induced constitutive active/androstane receptor-mediated DNA methylation changes during liver tumorigenesis and in liver tumors

A prognostic model based on DNA methylation-related gene expression for predicting overall survival in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) continues to increase in morbidity and mortality among all types of cancer. DNA methylation, an important epigenetic modification, is associated with cancer occurrence and progression. The objective of this study was to establish a model based on DNA methylation risk scores for identifying new potential therapeutic targets in HCC and preventing cancer progression.

Methods

Transcriptomic, clinical, and DNA methylation data on 374 tumor tissues and 50 adjacent normal tissues were downloaded from The Cancer Genome Atlas-Liver Hepatocellular Carcinoma database. The gene expression profiles of the GSE54236 liver cancer dataset, which contains data on 161 liver tissue samples, were obtained from the Gene Expression Omnibus database. We analyzed the relationship between DNA methylation and gene expression levels after identifying the differentially methylated and expressed genes. Then, we developed and validated a risk score model based on the DNA methylation-driven genes. A tissue array consisting of 30 human hepatocellular carcinoma samples and adjacent normal tissues was used to assess the protein and mRNA expression levels of the marker genes by immunohistochemistry and qRT-PCR, respectively.

Results

Three methylation-related differential genes were identified in our study: GLS, MEX3B, and GNA14. The results revealed that their DNA methylation levels were negatively correlated with local gene expression regulation. The gene methylation levels correlated strongly with the prognosis of patients with liver cancer. This was confirmed by qRT-PCR and immunohistochemical verification of the expression of these genes or proteins in tumors and adjacent tissues. These results revealed the relationship between the level of relevant gene methylation and the prognosis of patients with liver cancer as well as the underlying cellular and biological mechanisms. This allows our gene signature to provide more accurate and appropriate predictions for clinical applications.

Conclusion

Through bioinformatics analysis and experimental validation, we obtained three DNA methylation marker: GLS, MEX3B, and GNA14. This helps to predict the prognosis and may be a potential therapeutic target for HCC patients.

CAR requires Gadd45β to promote phenobarbital-induced mouse liver tumors in early stage

Phenobarbital (PB) is an archetypal substance used as a mouse hepatocellular carcinoma (HCC) promotor in established experimental protocols. Our previous results showed CAR is the essential factor for PB induced HCC promotion. Subsequent studies suggested Gadd45β, which is induced by PB through CAR activation, is collaborating with CAR to repress TNF-α induced cell death. Here, we used Gadd45β null mice (Gadd45β KO) treated with N-diethylnitrosamine (DEN) at 5 weeks of age and kept the mice with PB supplemented drinking water from 7 to 57 weeks old. Compared with wild type mice, Gadd45β KO mice developed no HCC in the PB treated group. Increases in liver weight were more prominent in wild type mice than KO mice. Microarray analysis of mRNA derived from mouse livers found multiple genes specifically up or down regulated in wild type mice but not null mice in DEN + PB groups. Further qPCR analysis confirmed two genes, Tgfbr2 and irisin/Fndc5, were up-regulated in PB treated wild type mice but no significant increase was observed in Gadd45β KO mice. We focused on these two genes because previous reports showed that hepatic Irisin/Fndc5 expression was significantly higher in HCC patients and that irisin binds to TGF-β receptor complex that includes TGFBR2 subunit. Our results revealed irisin peptide in cell culture media increased the growth rate of mouse hepatocyte-derived AML12 cells. Microarray analysis revealed that irisin-regulated genes in AML12 cells showed a significant association with the genes in the TGF-β pathway. Expression of irisin/Fndc5 and Tgfbr2 induced growth of human HCC cell line HepG2. Thus, Gadd45β plays an indispensable role in mouse HCC development regulating the irisin/Fndc5 and Tgfbr2 genes.

Research Progress in Liver-Regenerating Microenvironment and DNA Methylation in Hepatocellular Carcinoma: The Role of Traditional Chinese Medicine

The development, progression, recurrence, and metastasis of hepatocellular carcinoma (HCC) are closely associated with an abnormal liver-regenerating microenvironment (LRM). Therefore, preventing and reversing an abnormal LRM is a potential therapeutic strategy against HCC. Studies are increasingly focusing on the impact of regeneration, fibrosis, angiogenesis, inflammation, immunomodulation, and hepatic stem cells on HCC development and progression. As a key epigenetic mechanism, DNA methylation is extensively involved in regulating physiological and pathological pathways. In this review, we summarize recent findings on the role of DNA methylation in the fibrotic, angiogenic, inflammatory/immune, and stem cell microenvironments of HCC, and discuss new advances in Traditional Chinese Medicine (TCM) on influencing the abnormal LRM, so as to gain new insights into alleviating the abnormal LRM via regulating DNA methylation by TCM.

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.

Case examples of an evaluation of the human relevance of the pyrethroids/pyrethrins-induced liver tumours in rodents based on the mode of action

Rodent carcinogenicity studies are useful for screening for human carcinogens but they are not perfect. Some modes of action (MOAs) lead to cancers in both experimental rodents and humans, but others that lead to cancers in rodents do not do so in humans. Therefore, analysing the MOAs by which chemicals produce tumours in rodents and determining the relevance of such tumour data for human risk are critical. Recently, experimental data were obtained as case examples of an evaluation of the human relevance of pyrethroid (metofluthrin and momfluorothrin)- and pyrethrins-induced liver tumours in rats based on MOA. The MOA analysis, based on the International Programme on Chemical Safety (IPCS) framework, concluded that experimental data strongly support that the postulated MOA for metofluthrin-, momfluorothrin- and pyrethrins-produced rat hepatocellular tumours is mediated by constitutive androstane receptor (CAR) activation. Since metofluthrin and momfluorothrin are close structural analogues, reproducible outcomes for both chemicals provide confidence in the MOA findings. Furthermore, cultured human hepatocyte studies and humanized chimeric mouse liver studies demonstrated species difference between human hepatocytes (refractory to the mitogenic effects of these compounds) and rat hepatocytes (sensitive to their mitogenic effects). These data strongly support the hypothesis that the CAR-mediated MOA for liver tumorigenesis is of low carcinogenic risk for humans. In this research, in addition to cultured human hepatocyte studies, the usefulness of the humanized chimeric liver mouse models was clearly demonstrated. These data substantially influenced decisions in regulatory toxicology. In this review I comprehensively discuss the human relevance of the CAR-mediated MOA for rodent liver tumorigenesis based on published information, including our recent molecular research on CAR-mediated MOA.

Candidate genes responsible for early key events of phenobarbital-promoted mouse hepatocellular tumorigenesis based on differentiation of regulating genes between wild type mice and humanized chimeric mice

Phenobarbital (PB) is a nongenotoxic hepatocellular carcinogen in rodents. PB induces hepatocellular tumors by activating the constitutive androstane receptor (CAR). Some previous research has suggested the possible involvement of epigenetic regulation in PB-promoted hepatocellular tumorigenesis, but the details of its molecular mechanism are not fully understood. In the present study, comprehensive analyses of DNA methylation, hydroxymethylation and gene expression using microarrays were performed in mouse hepatocellular adenomas induced by a single 90 mg kg-1 intraperitoneal injection dose of diethylnitrosamine (DEN) followed by 500 ppm PB in the diet for 27 weeks. DNA modification and expression of hundreds of genes are coordinately altered in PB-induced mouse hepatocellular adenomas. Of these, gene network analysis showed alterations of CAR signaling and tumor development-related genes. Pathway enrichment analysis revealed that differentially methylated or hydroxymethylated genes belong mainly to pathways involved in development, immune response and cancer cells in contrast to differentially expressed genes belonging primarily to the cell cycle. Furthermore, overlap was evaluated between the genes with altered expression levels with 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) alterations in mouse hepatocellular adenoma induced by DEN/PB and the genes with altered expression levels in the liver of CD-1 mice or humanized chimeric mice treated with PB for 7 days. With the integration of transcriptomic and epigenetic approaches, we detected candidate genes responsible for early key events of PB-promoted mouse hepatocellular tumorigenesis. Interestingly, these genes did not overlap with genes altered by the PB treatment of humanized chimeric mice, thus suggesting a species difference between the effects of PB in mouse and human hepatocytes.

Dose-response analysis of epigenetic, metabolic, and apical endpoints after short-term exposure to experimental hepatotoxicants

Identification of sensitive and novel biomarkers or endpoints associated with toxicity and carcinogenesis is of a high priority. There is increasing interest in the incorporation of epigenetic and metabolic biomarkers to complement apical data; however, a number of questions, including the tissue specificity, dose-response patterns, early detection of those endpoints, and the added value need to be addressed. In this study, we investigated the dose-response relationship between apical, epigenetic, and metabolomics endpoints following short-term exposure to experimental hepatotoxicants, clofibrate (CF) and phenobarbital (PB). Male F344 rats were exposed to PB (0, 5, 25, and 100 mg/kg/day) or CF (0, 10, 50, and 250 mg/kg/day) for seven days. Exposure to PB or CF resulted in dose-dependent increases in relative liver weights, hepatocellular hypertrophy and proliferation, and increases in Cyp2b1 and Cyp4a1 transcripts. These changes were associated with altered histone modifications within the regulatory units of cytochrome genes, LINE-1 DNA hypomethylation, and altered microRNA profiles. Metabolomics data indicated alterations in the metabolism of bile acids. This study provides the first comprehensive analysis of the apical, epigenetic and metabolic alterations, and suggests that the latter two occur within or near the dose response curve of apical endpoint alterations following exposure to experimental hepatotoxicants.

Small-molecule modulators of the constitutive androstane receptor

INTRODUCTION

The constitutive androstane receptor (CAR) induces drug-metabolizing enzymes for xenobiotic metabolism.

AREAS COVERED

This review covers recent advances in elucidating the biological functions of CAR and its modulation by a growing number of agonists and inhibitors.

EXPERT OPINION

Extrapolation of animal CAR function to that of humans should be carefully scrutinized, particularly when rodents are used in evaluating the metabolic profile and carcinogenic properties of clinical drugs and environmental chemicals. Continuous efforts are needed to discover novel CAR inhibitors, with extensive understanding of their inhibitory mechanism, species selectivity, and discriminating power against other xenobiotic sensors.

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.

Phenobarbital induces alterations in the proteome of hepatocytes and mesenchymal cells of rat livers

Preceding studies on the mode of action of non-genotoxic hepatocarcinogens (NGCs) have concentrated on alterations induced in hepatocytes (HCs). A potential role of non-parenchymal liver cells (NPCs) in NGC-driven hepatocarcinogenesis has been largely neglected so far. The aim of this study is to characterize NGC-induced alterations in the proteome profiles of HCs as well as NPCs. We chose the prototypic NGC phenobarbital (PB) which was applied to male rats for a period of 14 days. The livers of PB-treated rats were perfused by collagenase and the cell suspensions obtained were subjected to density gradient centrifugation to separate HCs from NPCs. In addition, HCs and NPC isolated from untreated animals were treated with PB in vitro. Proteome profiling was done by CHIP-HPLC and ion trap mass spectrometry. Proteome analyses of the in vivo experiments showed many of the PB effects previously described in HCs by other methods, e.g. induction of phase I and phase II drug metabolising enzymes. In NPCs proteins related to inflammation and immune regulation such as PAI-1 and S100-A10, ADP-ribosyl cyclase 1 and to cell migration such as kinesin-1 heavy chain, myosin regulatory light chain RLC-A and dihydropyrimidinase-related protein 1 were found to be induced, indicating major PB effects on these cells. Remarkably, in vitro treatment of HCs and NPCs with PB hardly reproduced the proteome alterations observed in vivo, indicating differences of NGC induced responses of cells at culture conditions compared to the intact organism. To conclude, the present study clearly demonstrated that PB induces proteome alterations not only in HCs but also in NPCs. Thus, any profound molecular understanding on the mode of action of NGCs has to consider effects on cells of the hepatic mesenchyme.

Epigenetic profiles as defined signatures of xenobiotic exposure

With the advent of high resolution sequencing technologies there has been increasing interest in the study of genome-wide epigenetic modification patterns that govern the underlying gene expression events of a particular cell or tissue type. There is now mounting evidence that perturbations to the epigenetic landscape occur during a host of cellular processes including normal proliferation/differentiation and aberrant outcomes such as carcinogenesis. Furthermore, epigenetic perturbations have been associated with exposure to a range of drugs and toxicants, including non-genotoxic carcinogens (NGCs). Although a variety of epigenetic modifications induced by NGCs have been studied previously, recent genome-wide integrated epigenomic and transcriptomic studies reveal for the first time the extent and dynamic nature of the epigenetic perturbations resulting from xenobiotic exposure. The interrogation and integration of one such epigenetic mark, the newly discovered 5-hydroxymethylcytosine (5hmC) modification, reveals that drug treatment associated perturbations of the epigenome can result in unique epigenetic signatures. This review focuses on how recent advances in the field of epigenetics can enhance our mechanistic understanding of xenobiotic exposure and provide novel safety biomarkers.

Different Pathways of Constitutive Androstane Receptor–mediated Liver Hypertrophy and Hepatocarcinogenesis in Mice Treated with Piperonyl Butoxide or Decabromodiphenyl Ether

The constitutive androstane receptor (CAR) is essential for Cyp2b induction, liver hypertrophy, and hepatocarcinogenesis in response to phenobarbital (PB). Liver hypertrophy with Cyp2b induction is a major mode of action of hepatocarcinogenesis in rodents. However, it remains unclear whether CAR is involved in the response to many other nongenotoxic hepatocarcinogens besides PB. In this study, we investigated CAR involvement in liver hypertrophy and hepatocarcinogenesis of Cyp2b-inducing nongenotoxic hepatocarcinogens, piperonyl butoxide (PBO), and decabromodiphenyl ether (DBDE), using wild-type and CAR knockout (CARKO) male mice. PB was used as the positive control. In the wild-type mice, 4-week treatment with PBO, DBDE, or PB induced hepatocellular hypertrophy with increased Cyp2b10 messenger RNA and Cyp2b protein expression. In CARKO mice, only PBO showed liver hypertrophy with Cyp2b10 and Cyp3a11 induction. After 27-week treatment following diethylnitrosamine initiation, PBO and PB generated many eosinophilic altered foci/adenomas in wild-type mice; however, the lesions were far less frequent in CARKO mice. DBDE increased the multiplicity of basophilic altered foci/adenomas in wild-type and CARKO mice. Our findings indicate that murine CAR plays major roles in hepatocarcinogenesis but not in liver hypertrophy of PBO. DBDE may act via CAR-independent pathways during hepatocarcinogenesis.

Environmental toxicants, epigenetics, and cancer

Tumorigenesis, a complex and multifactorial progressive process of transformation of normal cells into malignant cells, is characterized by the accumulation of multiple cancer-specific heritable phenotypes triggered by the mutational and/or non-mutational (i.e., epigenetic) events. Accumulating evidence suggests that environmental and occupational exposures to natural substances, as well as man-made chemical and physical agents, play a causative role in human cancer. In a broad sense, carcinogenesis may be induced through either genotoxic or non-genotoxic mechanisms; however, both genotoxic and non-genotoxic carcinogens also cause prominent epigenetic changes. This review presents current evidence of the epigenetic alterations induced by various chemical carcinogens, including arsenic, 1,3-butadine, and pharmaceutical and biological agents, and highlights the potential for epigenetic changes to serve as markers for carcinogen exposure and cancer risk assessment.

Phenobarbital mediates an epigenetic switch at the constitutive androstane receptor (CAR) target gene Cyp2b10 in the liver of B6C3F1 mice

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.

Low doses of the carcinogen furan alter cell cycle and apoptosis gene expression in rat liver independent of DNA methylation

BACKGROUND

Evidence of potent rodent carcinogenicity via an unclear mechanism suggests that furan in various foods [leading to an intake of up to 3.5 microg/kg body weight (bw)/day] may present a potential risk to human health.

OBJECTIVES

We tested the hypothesis that altered expression of genes related to cell cycle control, apoptosis, and DNA damage may contribute to the carcinogenicity of furan in rodents. In addition, we investigated the reversibility of such changes and the potential role of epigenetic mechanisms in response to furan doses that approach the maximum estimated dietary intake in humans.

METHODS

The mRNA expression profiles of genes related to cell cycle, apoptosis, and DNA damage in rat liver treated with furan concentrations of 0.1 and 2 mg/kg bw were measured by quantitative polymerase chain reaction (PCR) arrays. We assessed epigenetic changes by analysis of global and gene-specific DNA methylation [methylation-specific PCR, combined bisulfite restriction analysis (COBRA), and methylated DNA immunoprecipitation chip] and microRNA (miRNA) analyses.

RESULTS

The expression profiles of apoptosis-related and cell-cycle-related genes were unchanged after 5 days of treatment, although we observed a statistically significant change in the expression of genes related to cell cycle control and apoptosis, but not DNA damage, after 4 weeks of treatment. These changes were reversed after an off-dose period of 2 weeks. None of the gene expression changes was associated with a change in DNA methylation, although we detected minor changes in the miRNA expression profile (5 miRNA alterations out of 349 measured) that may have contributed to modification of gene expression in some cases.

CONCLUSION

Nongenotoxic changes in gene expression may contribute to the carcinogenicity of furan in rodents. These findings highlight the need for a more comprehensive risk assessment of furan exposure in humans.

Applications of genetically modified tools to safety assessment in drug development

The process of new drug development consists of several stages; after identifying potential candidate compounds, preclinical studies using animal models link the laboratory and human clinical trials. Among many steps in preclinical studies, toxicology and safety assessments contribute to identify potential adverse events and provide rationale for setting the initial doses in clinical trials. Gene modulation is one of the important tools of modern biology, and is commonly employed to examine the function of genes of interest. Advances in new drug development have been achieved by exploding information on target selection and validation using genetically modified animal models as well as those of cells. In this review, a recent trend of genetically modified methods is discussed with reference to safety assessments, and the exemplary applications of gene-modulating tools to the tests in new drug development were summarized.

The constitutive active/androstane receptor facilitates unique phenobarbital-induced expression changes of genes involved in key pathways in precancerous liver and liver tumors

Our overall goal is to elucidate progressive changes, in expression and methylation status, of genes which play key roles in phenobarbital (PB)-induced liver tumorigenesis, with an emphasis on their potential to affect signaling through critical pathways involved in the regulation of cell growth and differentiation. PB-elicited unique expression changes of genes, including some of those identified previously as exhibiting regions of altered DNA methylation, were discerned in precancerous liver tissue and/or individual liver tumors from susceptible constitutive active/androstane receptor (CAR) wild-type (WT) compared with resistant CAR knockout (KO) mice. Many of these function in crucial cancer-related processes, for example, angiogenesis, apoptosis, cell cycle, DNA methylation, Hedgehog signaling, invasion/metastasis, Notch signaling, and Wnt signaling. Furthermore, a subset of the uniquely altered genes contained CAR response elements (CAREs). This included Gadd45b, a coactivator of CAR and inhibitor of apoptosis, and two DNA methyltransferases (Dnmt1, Dnmt3a). The presence of CAREs in Dnmts suggests a potential direct link between PB and altered DNA methylation. The current data are juxtaposed with the effects of PB on DNA methylation and gene expression which occurred uniquely in liver tumor-prone B6C3F1 mice, as compared with the resistant C57BL/6, following 2 or 4 weeks of treatment. Collectively, these data reveal a comprehensive view of PB-elicited molecular alterations (i.e., changes in gene expression and DNA methylation) that can facilitate hepatocarcinogenesis. Notably, candidate genes for initial "fingerprints" of early and late stages of PB-induced tumorigenesis are proposed.

Phenobarbital elicits unique, early changes in the expression of hepatic genes that affect critical pathways in tumor-prone B6C3F1 mice

At 2 and 4 weeks following treatment with phenobarbital (PB), the classical nongenotoxic rodent liver carcinogen, we elucidated unique gene expression changes (both induction and repression) in liver tumor-susceptible B6C3F1 mice, as compared with the relatively resistant C57BL/6. Based on their cancer-related roles, we believe that altered expression of at least some of these genes might underlie PB-induced liver tumorigenesis. Putative constitutive active/androstane (CAR) response elements (CAREs), a subset of PB response elements, were present within multiple genes whose expression was uniquely altered in the B6C3F1 mice, suggesting a role for CAR in their regulation. Additionally, three DNA methyltransferase genes (Dnmt1, Dnmt3a, and Dnmt3b) were repressed uniquely in the tumor-prone B6C3F1 mice, and all possess putative CAREs, providing a potential direct link between PB and expression of key genes that regulate DNA methylation status. Previously, we demonstrated that PB-elicited unique regions of altered methylation (RAMs) in B6C3F1 mice, as compared with the relatively resistant C57BL/6, at 2 and 4 weeks, and annotation of the regions harboring these changes revealed 51 genes. This is extended by the current study, which employed RNA isolated from the same liver tissue used in the earlier investigations. Genes elucidated from both the methylation and expression analyses are involved in identical processes/pathways (e.g., cell cycle, apoptosis, angiogenesis, epithelial-mesenchymal cell transition, invasion/metastasis, and mitogen-activated protein kinase, transforming growth factor-beta, and Wnt signaling). Therefore, these changes might represent very early events that directly contribute to PB-induced tumorigenesis. It is instructive to consider the possibility that, in a hypothesis-driven fashion, these genes are initial candidates that could be utilized to develop a biomarker "fingerprint" of early exposure to PB and PB-like compounds.