Research article | Open | Open Peer Review | Published:
The DNA methylation profile of liver tumors in C3H mice and identification of differentially methylated regions involved in the regulation of tumorigenic genes
BMC Cancervolume 18, Article number: 317 (2018)
C3H mice have been frequently used in cancer studies as animal models of spontaneous liver tumors and chemically induced hepatocellular carcinoma (HCC). Epigenetic modifications, including DNA methylation, are among pivotal control mechanisms of gene expression leading to carcinogenesis. Although information on somatic mutations in liver tumors of C3H mice is available, epigenetic aspects are yet to be clarified.
We performed next generation sequencing-based analysis of DNA methylation and microarray analysis of gene expression to explore genes regulated by DNA methylation in spontaneous liver tumors of C3H mice. Overlaying these data, we selected cancer-related genes whose expressions are inversely correlated with DNA methylation levels in the associated differentially methylated regions (DMRs) located around transcription start sites (TSSs) (promoter DMRs). We further assessed mutuality of the selected genes for expression and DNA methylation in human HCC using the Cancer Genome Atlas (TCGA) database.
We obtained data on genome-wide DNA methylation profiles in the normal and tumor livers of C3H mice. We identified promoter DMRs of genes which are reported to be related to cancer and whose expressions are inversely correlated with the DNA methylation, including Mst1r, Slpi and Extl1. The association between DNA methylation and gene expression was confirmed using a DNA methylation inhibitor 5-aza-2′-deoxycytidine (5-aza-dC) in Hepa1c1c7 cells and Hepa1-6 cells. Overexpression of Mst1r in Hepa1c1c7 cells illuminated a novel downstream pathway via IL-33 upregulation. Database search indicated that gene expressions of Mst1r and Slpi are upregulated and the TSS upstream regions are hypomethylated also in human HCC. These results suggest that DMRs, including those of Mst1r and Slpi, are involved in liver tumorigenesis in C3H mice, and also possibly in human HCC.
Our study clarified genome wide DNA methylation landscape of C3H mice. The data provide useful information for further epigenetic studies of mice models of HCC. The present study particularly proposed novel DNA methylation-regulated pathways for Mst1r and Slpi, which may be applied not only to mouse HCC but also to human HCC.
C3H mice have been employed for numerous studies of carcinogenesis. The mice, particularly their males, are predisposed to spontaneously develop liver tumors in adulthood and have been used as an animal model of spontaneous hepatocellular carcinoma (HCC) [1, 2]. They are also often used as chemically induced HCC models, such as diethylnitrosamine (DEN)- or DEN and phenobarbital (DEN/PB)-induced HCC models [1, 3,4,5]. A recent study established an HCC model with liver cirrhosis, a key feature of human HCC, by administering carbon tetrachloride (CCl4) to C3H mice .
To investigate the process of carcinogenesis and therapeutic methods in proper animal models, information on the genetic and epigenetic background of their tumors should greatly help understanding the mechanisms involved. A well-known genetic feature of spontaneous liver tumors of C3H mice is the occurrence of somatic mutations in the proto-oncogene Ha-ras in approximately 10–60% of the tumors . Somatic point mutations of ras genes (Ha-ras, K-ras, and N-ras) were also detected in a variety of human and animal tumors7. These somatic mutations activate Ras proteins and lead to activation of tumor augmenting pathways, including the RAF/MEK/ERK kinase cascade . Another somatic mutation was reported in B-raf, a member of the RAF kinase family in about 20% of liver tumors in C3H mice administered DEN at 2 weeks of age . B-raf mutation also augments the kinase activity leading to activation of the MEK/ERK cascade [5, 8].
On the other hand, some of the spontaneous liver tumors of C3H mice do not have Ha-ras or B-raf mutations [8, 9]. Epigenetic modifications, such as DNA methylation and histone modifications, are pivotal posteriori control mechanisms of gene expression and are closely involved in tumorigenesis . A recent study performed next-generation sequencing of DNA methylation of human non-cirrhotic HCC and fibrolamellar hepatocellular carcinoma (FLC) and found distinctive epigenetic signatures of the two types of tumors . While we previously studied methylation status in the spontaneous liver tumors of control C3H mice and liver tumors of gestationally arsenite exposed C3H mice by a methylated DNA immunoprecipitation (MeDIP)–CpG island microarray method , precise single base resolution analysis of DNA methylation has not been performed for the liver tumors of C3H mice.
Recent progress in genomic study using next generation sequencing enabled DNA methylation analysis at one-base pair resolution and revealed more precise modes of action of DNA methylation, beyond gene silencing by canonical CpG islands. DNA methylation is involved in suppression/activation of gene expression, transcription factor binding, splicing and nucleosome posioning . In the present study, we clarified differentially methylated CpGs (DMCs) and differentially methylated regions (DMRs), where DNA in the spontaneous liver tumors of C3H mice is hyper- or hypomethylated compared to the normal livers, by reduced representation bisulfite sequencing (RRBS) method, a sequencing-based genome-wide DNA methylation analysis [14, 15]. According to the recent reports showing that DNA methylation levels within ± 1000 bp regions  or ± 2000 bp regions [11, 16, 17] of TSS showed a strong correlation with gene repression, we focused on DMRs within ± 2000 bp of TSS (promoter DMRs). Using the data of promoter DMRs and by in vitro experiments and PubMed search, we identified genes whose expressions are closely associated with the methylation of proximal promoter DMRs and which are involved in tumorigenesis. We also assessed the gene expression and DNA methylation status of these genes in human HCC using TCGA database to explore similarities in tumorigenic regulation in human.
Male C3H/HeN mice around 74 weeks of age were obtained from F2 pups of the control group of a project for investigating the effects of gestational arsenic exposure on the F2. The design of animal breeding was the same as described previously [9, 18]. Briefly, pregnant F0 C3H/HeN mice were purchased from CLEA Japan (Tokyo, Japan) and F1 and F2 pups in the control group were given free access to a standard diet (CA-1; CLEA Japan) and tap water (the control group). The mice were handled in a humane manner in accordance with the National Institute for Environmental Studies (NIES) guidelines for animal experiments.
Cell lines and treatment
Hepa1c1c7 cells were kindly provided by Dr. Y. Fujii-Kuriyama (University of Tsukuba) in 2005  and cultured in the Dulbecco’s Modified Eagle’s Medium (DMEM, SIGMA-ALDRICH, D5796) containing 10% FBS and 1% penicillin-streptomycin. Hepa1-6 cells provided by RIKEN BRC CELL BANK were cultured in DMEM containing 10% FBS, 1% penicillin-streptomycin and 5% sodium pyruvate. 5-Aza-dC was purchased from Santa Cruz (CA, USA). After attaching to the dish, cells were incubated in a culture medium containing 5-aza-dC for 72 h. For overexpression study, mouse ORF clones Mst1r (NM_001287261) and Slpi (NM_011414), and pCMV6-Kan/Neo (pCMV6KN) as an empty vector were purchased from ORIGENE Inc. (MD, USA). Hepa1c1c7 cells were transfected with 2 μg of vectors using HilyMax (Dojindo, Japan) according to the manufacturer’s instructions and gene expression changes were examined after culturing for 24 h.
Genomic DNA extraction
Genomic DNA was prepared as previously described  from the normal livers and macroscopic liver tumor tissues of C3H mice and cell lines. Briefly tissues or cells were lysed in lysis buffer, treated with RNase, and purified with a phenol-chloroform mixture.
RRBS libraries were prepared from genomic DNA of 3 normal liver tissues and 3 tumor tissues having a Ha-ras mutation according to the protocol reported by Boyle et al. . with some modifications. Briefly, 100 ng of genomic DNA was digested with MspI, subjected to gap-filling and A-tailing, and ligated with TruSeq adaptors included in the TruSeq DNA Sample Prep Kit (FC-121-1001, Illumina). After DNA size selection and ligation efficiency were checked, adaptor-ligated DNA was subjected to bisulfite conversion and amplified by PCR. The PCR products were cleaned-up using Agencourt AMPure XP beads and sequenced on an Illumina HiSeq2500. The sequence data have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE111420 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111420).
Alignment of RRBS data and identification of DMCs and DMRs
Bisulfite-treated sequencing reads were aligned by a paired-end alignment method for a unique best hit to the mouse reference genome (NCBI/mm10) using the Bismark program , and adapter trimming and quality control was performed using Trim Galore (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The aligned data were used for generating and analyzing DNA methylation profiles using the methylKit package . Briefly, %methylation scores were determined for CpG sites using the read.bismark function of methylKit as SAM format alignments files created with the Bismark aligner as input files. For each CpG site, the read.bismark function counts the numbers of reads with C (methylated and remained) and with T (unmethylated and converted) in both strands separately, and calculates %methylation. The obtained %methylation profiles were evaluated for their overall statistics and extent of overall similarities among samples using methyKit functions such as getMethylationStats and getCorrelation .
DMCs were extracted using methylKit for the cytosines in CpG sites with < 0.01 q-value by the logistic regression method, ≥10 reads in coverage and ≥ 25% methylation differences between normal and tumor tissues. Under the condition, 2,331,643 CpG sites were identified. DMRs were detected for regions containing ≥3 CpG sites and at least 1 DMC and with ≥30% absolute mean CpG methylation difference by using an eDMR package on R. eDMR implements statistical analysis . Promoter DMRs, which were defined as regions within ±2000 bp from TSSs for Refseq transcripts [11, 16, 17], were extracted using bedtools 2.25.0 (bedtools.readthedocs.io). A summary of the detailed sequencing and alignment characteristics for the samples used in this study will be published elsewhere with the data of arsenite-exposed F2 mice (Okamura et al., in preparation).
Bisulfite sequencing analysis was carried out as previously described with minor modifications . Genomic DNA (1 μg) was digested with EcoRI and subjected to bisulfite modification using EZ DNA Methylation-Gold Kit (Zymo Research, CA, USA). The bisulfite-treated DNA (2 ng) was amplified with the primers shown in Additional file 1: Table S1. The PCR products were purified with the Wizard SV Gel and PCR Clean-Up System (Promega) and cloned into pGEM-T Easy Vector (Promega). Obtained clones were cycle sequenced with M13RV primers (5′- CAGGAAACAGCTATGAC-3′) and a BigDye Terminator version 3.1 Cycle Sequencing Kit (Applied Biosystems, CA, USA), and analyzed with an Applied Biosystems 3730 DNA analyzer.
Microarray analysis of gene expression
Total RNA was prepared from the 3 normal liver tissues, 3 tumor tissues having Ha-ras mutation and Hepa1c1c7 cells with an RNeasy Mini Kit (Qiagen, Valencia, CA, USA). After the quality of total RNA was checked using a 2100 Bioanalyzer (Agilent Technologies), the gene expression profile was analyzed using SurePrint G3 Mouse GE arrays (Agilent Technologies) according to the manufacturer’s protocol for one-color microarray-based gene expression analysis. The scanned images were analyzed with Feature Extraction Software 9.1 to obtain background subtracted and spatially detrended Processed Signal intensities. The microarray data have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE104627 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104627) and GSE104626 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104626).
The expression data from each chip were then normalized to the 75the percentile of all measurement using Agilent Gene-Spring GX. Genes with an expression value of more than 1.0, in either normal tissues or tumor tissues or in both, were considered to be expressed above the threshold and utilized for analysis. The genes with a two- or more fold increase or decrease in mean expression values (n = 3) were judged to have changed. For Hepa1c1c7 cells with or without overexpression of Mst1r or Slpi, genes with raw expression values of more than 50, either after upregulation or before downregulation or in both cases, were utilized for further analysis.
Real-time PCR analysis
Total RNA was prepared from tissues and cells with an RNeasy Mini Kit (Qiagen). cDNA synthesis and real-time PCR was carried out as described previously . The primers used are shown in Additional file 1: Table S1.
Prediction of the transcription factor binding site
Consensus regions of transcription factor binding sites were detected using the JASPAR database (http://jaspar.genereg.net/)  with the relative score threshold set at > 90% and also by using published literature.
Forty one paired datasets of normal and HCC tissues of human livers for RNA-seq and DNA methylation levels were obtained from the TCGA database (http://cancergenome.nih.gov/). They were downloaded and used to find genes with significant expression differences and differences in DNA methylation within ±2000 bp from TSS between normal and HCC tissues.
The differences in gene expression between the two groups were analyzed by a two-tailed paired Student’s t-test. A p-value < 0.05 was considered to be statistically significant. The differences in gene expression among three groups were analyzed by one-way ANOVA followed by the Turkey-Kramer test as a post hoc comparison. Spearman correlation analysis was used to determine the relationship between DNA methylation differences and gene expression ratios of normal and liver tumors in C3H mice. For multiple comparisons to identify DNA methylation differences, paired t-test with Bonferroni correction was performed.
Genome wide analysis of DNA methylation and gene expression in normal and tumor livers of C3H mice
We conducted genome wide DNA methylation analysis for normal liver tissues and tumor tissues of male C3H mice by RRBS. Approximately 22–55 million sequence read pairs were obtained for the RRBS libraries and 63–71% of the reads were aligned to unique genomic locations. Hierarchical clustering analysis showed that DNA methylation status was distinctly separated between the normal and tumor tissues (Fig. 1a). Scatter plots of methylation values show high correlation coefficient value for each pair (Additional file 2: Figure S1).
We identified 103,965 DMCs (4.5% in total CpGs) with higher ratios of hypo-methylated CpGs in tumor tissues (hypo-DMCs), such as 60–80% in total DMCs in autosomal chromosomes and more than 90% in sex chromosomes (Fig. 1b). Annotation analysis showed that the majority of DMCs were in intron and intergenic regions (Fig. 1c) and in the region other than CpG island (CGI) or CGI shore (Fig. 1d). Higher ratio of hyper-DMCs, compared with hypo-DMCs, were found in the region between − 5000 bp upstream and TSS (Fig. 1c) and in CGI and CGI shore (Fig. 1d).
We then identified 3337 DMRs (Additional file 3: Table S2) including 225 promoter DMRs. Among the 225 genes driven by the differentially methylated promoters, 118 genes are expressed (see Methods). The methylation status of the 118 promoter DMRs was inversely associated with gene expression [13, 16, 17]. The scatter plot of DNA methylation differences of DMRs and gene expression ratios between the normal and tumor tissues of the 118 genes showed moderate inverse correlation (Fig. 1e). From the 118 genes, we shortlisted 27 genes which show more than 2-fold expression differences (p < 0.05) between normal and tumor tissues and inverse association with methylation status of the DMR, and which appear in a PubMed search for the term “cancer” (Table 1 and gray square in Fig. 1e).
Expressions of Mst1r, Slpi and Extl1 are associated with DNA methylation alterations in the liver tumors of C3H mice
Among the genes shown in Table 1, we focused on Mst1r and Slpi, which are greatly up- regulated (67- and 408-fold, p-value 0.0005 and 0.009, respectively) and have relatively large numbers of DMCs (13 and 9) in the DMRs, and Extl1, which have a large number of DMCs (44 DMCs) in the DMR. Mst1r is a receptor tyrosine kinase activated by the binding of the ligand, Mst1 (macrophage-stimulating 1, also known as hepatocyte growth factor-like; HGFL) [26, 27]. Overexpression of Mst1r activates the molecule without the need of the ligands . Such overexpression and activation of Mst1r are found in many types of human cancers and often correlate with poor prognosis in various cancers . Slpi has activity as serine protease inhibitor and is overexpressed in various types of cancers , including mouse liver tumors induced by CCl4 or DEN . Extl1 is a putative tumor-suppressor gene found in neuroblastoma patients .
The DNA methylation status of these genes are visualized in Figs. 2a, 3a, and 4a. Mst1r and Slpi, which are upregulated in the tumors, have hypo-DMRs starting from + 1192 and + 697 bps from TSS, respectively. Extl1, whose expression is downregulated in the tumors, has hyper-DMR starting from + 929 bp from TSS. We validated the DNA methylation status of these DMRs by bisulfite sequencing (Figs. 2b, 3b, and 4b) and confirmed that the RRBS data correlated well with the bisulfite sequencing data. We also measured gene expression levels of the normal and tumor tissues by real-time PCR and confirmed up-regulation of Mst1r and Slpi and down-regulation of Extl1 in the tumor tissues (Figs. 2c, 3c, and 4c). The expression levels of these genes did not differ between tumor tissues with and without Ha-ras mutation (Additional file 2: Figure S2).
Expressions of Mst1r, Slpi, and Extl1 are upregulated in cancer cell lines by 5-aza-dC treatment
To determine whether DNA methylation alterations in the DMRs of Mst1r, Slpi, and Extl1 are involved in the expression of these genes, we assessed the effects of the DNA methylation inhibitor 5-aza-dC on the DNA methylation statuses and gene expressions in mouse liver cancer cell lines (Fig. 5). Mst1r and Slpi were up-regulated and hypomethylated at DMRs by culturing Hepa1c1c7 cells with 50 μM and 100 μM of 5-aza-dC (Fig. 5a, c). Expression of Extl1 was not detected in Hepa1c1c7 cells. In Hepa1-6 cells, the three genes were up-regulated and their DMRs were hypomethylated by culturing with 0.1 and 1 μM of 5-aza-dC (Fig. 5b, d, e). These data support the notion that expressions of Mst1r, Slpi and Extl1 are upregulated by hypomethylation of DMRs.
Interestingly, susceptibility against 5-aza-dC largely differed between the two cell lines. Measurement of gene expression levels of DNA methyltransferase Dnmt1/3a/3b and DNA demethylation enzymes Tet1/2/3, the major enzymes determining DNA methylation level, clarified that Hepa1-6 cells have significantly lower expression levels of Dnmts in comparison with Hepa1c1c7 cells (Additional file 2: Figure S3). Expression levels of Tet2 was also significantly lower in Hepa1-6 cells compared to Hepa1c1c7 cells (p < 0.01). Tet1 expression was not detected in either cell line. These results suggest that the lowered expression of DNA methyltransferases is involved in the different sensitivity to 5-aza-dC between Hepa1-6 and Hepa1c1c7 cells.
Transcription factor (TF) binding sites in the DMRs of Mst1r and Slpi
DNA methylation of TF binding sites or their proximity near the TSS can affect the binding of transcription factors and regulate transcription of genes [13, 30]. We searched for TF binding sites in the DMRs of Mst1r, Slpi, and Extl1 using the JASPAR database  and literature. We found a CpG site overlapping NF-kB binding site in the DMR of Mst1r (Fig. 2a) and one overlapping AhR binding site in the DMR of Slpi (Fig. 3a). Methylation of CpG in the binding sequence for NF-kB  and AhR [32, 33] have been reported to inhibit their binding affinities. No CpG site overlapping the TF binding site was detected in the DMR of Extl1 (Fig. 4a).
Overexpression of Mst1r induces robust expression of Il33
We further assessed whether upregulation of Mst1r and Slpi play causal roles in tumorigenesis by overexpressing these genes and analyzing gene expression changes using microarrays in hepa1c1c7 cells.
We selected genes which are statistically significantly up or down regulated more than 10 fold by overexpression of Mst1r or Slpi, also show the same direction of significant expression changes in the liver tumors of C3H mice, and are found to be related with “cancer” by a PubMed search.
Overexpression of Mst1r induced IL33 which meets all the above criteria (Fig. 6a, b). Real-time PCR confirmed up-regulation of IL33 in the liver tumors of C3H mice (Fig. 6c). No gene was identified under the criteria in the cells overexpressing Slpi.
Mst1r and Slpi are upregulated in human HCC with altered DNA methylation status around the TSS
We analyzed the association between gene expression and DNA methylation status of human HCC by use of TCGA database. Data from 41 pairs of normal livers and HCC tissues were used for each analysis (Additional file 4: Table S3). Expression of Mst1r and Slpi are upregulated in the human HCC dataset we selected (Fig. 7a and b), as reported previously [27, 28]. Similar to the methylation status of Mst1r and Slpi in C3H mice, the downstream regions of TSS in both genes were weakly but significantly hypomethylated in human HCC compared to the normal liver tissues (Fig. 7c). These results showed the possibility that both the genes are regulated by DNA methylation, not only in the liver tumors of C3H mice, but also in human HCC. Furthermore, the expression of IL-33 strongly (r = 0.807), although not significantly (p < 0.1), correlated with MST1R expression in human HCC tissues highly expressing MST1R (Fig. 7d). Correlation was not observed between the expressions of the two genes when all the 41 pairs were analyzed (r = 0.095).
In the present study, we investigated the DNA methylation landscape of liver tumors of C3H mice by RRBS analysis. Hypo-DNA methylation is a well-known feature of tumors . The present study clarified a large increase in the number of hypo-DMCs in the liver tumors of C3H mice (Fig. 1b), well illustrating the features of tumors. The present study further identified tumor-specific DMRs, including promoter DMRs for cancer-related genes.
Among the genes with promoter DMRs, we focused on three genes, Mst1r, Slpi and Extl1 (Table 1). Mst1r activation is involved in tumor augmenting signaling pathways, such as the PI3K/AKT and RAS/ERK pathways . Slpi is overexpressed in various types of cancers . While there are studies reporting anti-metastatic activity of Slpi, a large body of studies reported overexpression of Slpi augments or associates with tumor proliferation and invasion in cancers of many origins [28, 35]. Extl1 is reported as a putative tumor-suppressor gene . The present study showed hypo-DMR around + 1000 bp from TSS and robust upregulation in Mst1r and Slpi in the tumors (Figs. 2 and 3) and hyper-DMR around + 1000 bp from TSS and downregulation of Extl1 in the tumors (Fig. 4). These DMRs, which are not overlapped with CpG island, were not detected in our previous study on the spontaneous liver tumors of C3H mice using MeDIP-CpG microarrays . Furthermore, this study firstly reported an association between expression changes of these genes and DMRs by an in vitro experiment using 5-aza-dC (Fig. 5).
We analyzed TCGA data for 41 pairs of normal and HCC tissues from patients and found overexpression of MST1R and SLPI and lowered DNA methylation of TSS down-stream region also in human HCC (Fig. 7). On the other hand, a recent study reported downregulation of MST1R and hypermethylation of the TSS upstream region in human hepatoblastomas . The TSS upstream region of MST1R in the human HCC we analyzed was hypermethylated with transcription upregulation (Fig. 7c). Thus, MST1R expression and the association of DNA methylation status differ between HCC and hepatoblastomas. Aiming to target MST1R for cancer therapy, drugs such as tyrosine kinase inhibitors against the tyrosine kinase activity of MST1R and monoclonal antibodies against MST1R are being developed and some of these are under clinical trials or in research phases . The results of the present study encourage further study on DNA methylation status around TSS of Mst1r to develop epigenetic therapies.
In the present study, we further identified IL33 as the putative down-stream targets of Mst1r by overexpression of Mst1r into hepa1c1c7 cells (Fig. 6c). IL33 was also upregulated in the tumor tissues of C3H mice (Fig. 6e) and its expression is associated with MST1R expression in higher MST1R expressing group (Fig. 7d). IL33 activates NF-kB and ERK signaling pathways and is implicated in HCC . The present study newly demonstrates the possibility that IL33 is a down-stream target of Mst1r and is involved in the liver tumors of C3H mice and in human HCC.
DNA methylation exerts a variety of modes of action. DNA methylation of TF binding sequences affects the affinity of TFs to DNA [13, 30]. Hypermethylation of the Fas promoter, which contains 3 putative NF-kB binding sites inhibits NF-kB binding . NF-kB pathway is deeply involved in HCC in human  as well as C3H mice [39, 40]. The DMR of Mst1r in the C3H mice tumors contained the putative NF-kB consensus sequence (5′- GTGGAGCCCC(G)-3′) and the CpG in the sequence is hypomethylated (Fig. 2b). Thus, hypomethylation of CpGs in the putative NF-kB site in the DMR of Mst1r may take a part in augmenting NF-kB binding and activating Mst1r. The DMR of Slpi includes the core sequence (5’-CGTG-3′) of xenobiotic-responsible element (XRE) and its hypomethylation inhibits the binding of AhR/Arnt complex [32, 33]. Thus, hypermethylation of CpG in the XRE core sequence in Slpi (Fig. 3b) may be involved in Slpi activation through AhR/Arnt binding, while the involvement of AhR/Arnt in spontaneous hepatic tumors is unknown. Contribution of the methylation status in TF binding sites require further studies.
During this study, we found that the sensitivity of Hepa1c1c7 cells and Hepa1-6 cells to 5-aza-dC is largely different (Fig. 5). 5-aza-dC is now utilized as epigenetic medicine to recover the expression of tumor suppressor genes by reducing DNA methylation . 5-Aza-dC is transported into cells by nucleoside transporters, incorporated into DNA after being triphosphorylated by deoxy-cytidine kinase, and inhibits DNMT1 activity by forming a covalent bond with the enzyme or augmenting DNMT1 degradation. We measured the expression of DNA methylation and demethylation enzymes and found significantly lower expression of Dnmt1, Dnmt3a and Dnmt3b in Hepa1-6 cells compared to Hepa1c1c7 cells. Thus the suppression of Dnmts may contribute to the difference in the sensitivity of the cells to 5-aza-dC.
The present study clarified the DNA methylation landscape in association with gene expression in the liver tumors of C3H mice. The data provide useful information for further genome-wide studies of mice models of HCC. The present study particularly proposed novel pathways regulated by DNA methylation for Mst1r and Slpi, which may be applied not only to mouse liver tumors but also to human HCC.
Differentially methylated cytosine
Dulbecco’s modified Eagle’s medium
Differentially methylated region
Fetal bovine serum
Fibrolamellar hepatocellular carcinoma
Methylated DNA immunoprecipitation
Open reading frame
Reduced representation bisulfite sequencing
The Cancer Genome Atlas
Transcription start site
He L, Tian DA, Li PY, He XX. Mouse models of liver cancer: progress and recommendations. Oncotarget. 2015;6(27):23306–22.
Maronpot RR, Fox T, Malarkey DE, Goldsworthy TL. Mutations in the ras proto-oncogene: clues to etiology and molecular pathogenesis of mouse liver tumors. Toxicology. 1995;101(3):125–56.
Delgado E, Okabe H, Preziosi M, Russell JO, Alvarado TF, Oertel M, Nejak-Bowen KN, Zhang Y, Monga SP. Complete response of Ctnnb1-mutated tumours to beta-catenin suppression by locked nucleic acid antisense in a mouse hepatocarcinogenesis model. J Hepatol. 2015;62(2):380–7.
Chen X, Yamamoto M, Fujii K, Nagahama Y, Ooshio T, Xin B, Okada Y, Furukawa H, Nishikawa Y. Differential reactivation of fetal/neonatal genes in mouse liver tumors induced in cirrhotic and non-cirrhotic conditions. Cancer Sci. 2015;106(8):972–81.
Buchmann A, Karcier Z, Schmid B, Strathmann J, Schwarz M. Differential selection for B-raf and ha-ras mutated liver tumors in mice with high and low susceptibility to hepatocarcinogenesis. Mutat Res. 2008;638(1–2):66–74.
Reiberger T, Chen Y, Ramjiawan RR, Hato T, Fan C, Samuel R, Roberge S, Huang P, Lauwers GY, Zhu AX, et al. An orthotopic mouse model of hepatocellular carcinoma with underlying liver cirrhosis. Nat Protoc. 2015;10(8):1264–74.
Pylayeva-Gupta Y, Grabocka E, Bar-Sagi D. RAS oncogenes: weaving a tumorigenic web. Nat Rev Cancer. 2011;11(11):761–74.
Jaworski M, Buchmann A, Bauer P, Riess O, Schwarz M. B-raf and ha-ras mutations in chemically induced mouse liver tumors. Oncogene. 2005;24(7):1290–5.
Nohara K, Okamura K, Suzuki T, Murai H, Ito T, Shinjo K, Takumi S, Michikawa T, Kondo Y, Hata K. Augmenting effects of gestational arsenite exposure of C3H mice on the hepatic tumors of the F(2) male offspring via the F(1) male offspring. J Appl Toxicol. 2016;36(1):105–12.
Hassler MR, Egger G. Epigenomics of cancer - emerging new concepts. Biochimie. 2012;94(11):2219–30.
Malouf GG, Tahara T, Paradis V, Fabre M, Guettier C, Yamazaki J, Long H, Lu Y, Raynal NJ, Jelinek J, et al. Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features. Epigenetics. 2015;10(9):872–81.
Suzuki T, Yamashita S, Ushijima T, Takumi S, Sano T, Michikawa T, Nohara K. Genome-wide analysis of DNA methylation changes induced by gestational arsenic exposure in liver tumors. Cancer Sci. 2013;104(12):1575–85.
Tirado-Magallanes R, Rebbani K, Lim R, Pradhan S, Benoukraf T. Whole genome DNA methylation: beyond genes silencing. Oncotarget. 2017;8(3):5629–37.
Meissner A, Gnirke A, Bell GW, Ramsahoye B, Lander ES, Jaenisch R. Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Res. 2005;33(18):5868–77.
Boyle P, Clement K, Gu H, Smith ZD, Ziller M, Fostel JL, Holmes L, Meldrim J, Kelley F, Gnirke A, et al. Gel-free multiplexed reduced representation bisulfite sequencing for large-scale DNA methylation profiling. Genome Biol. 2012;13(10):R92.
De Souza AP, Planello AC, Marques MR, De Carvalho DD, Line SR. High-throughput DNA analysis shows the importance of methylation in the control of immune inflammatory gene transcription in chronic periodontitis. Clin Epigenetics. 2014;6(1):15.
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, et al. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41(2):178–86.
Nohara K, Tateishi Y, Suzuki T, Okamura K, Murai H, Takumi S, Maekawa F, Nishimura N, Kobori M, Ito T. Late-onset increases in oxidative stress and other tumorigenic activities and tumors with a ha-ras mutation in the liver of adult male C3H mice gestationally exposed to arsenic. Toxicol Sci. 2012;129(2):293–304.
Suzuki T, Nohara K. Regulatory factors involved in species-specific modulation of arylhydrocarbon receptor (AhR)-dependent gene expression in humans and mice. J Biochem. 2007;142(4):443–52.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.
Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A. Mason CE: methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13(10):R87.
Li S, Garrett-Bakelman FE, Akalin A, Zumbo P, Levine R, To BL, Lewis ID, Brown AL, D'Andrea RJ, Melnick A, et al. An optimized algorithm for detecting and annotating regional differential methylation. BMC Bioinformatics. 2013;14(Suppl 5):S10.
Takumi S, Aoki Y, Sano T, Suzuki T, Nohmi T, Nohara K. In vivo mutagenicity of arsenite in the livers of gpt delta transgenic mice. Mutat Res Genet Toxicol Environ Mutagen. 2014;760:42–7.
Nohara K, Ao K, Miyamoto Y, Ito T, Suzuki T, Toyoshiba H, Tohyama C. Comparison of the 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-induced CYP1A1 gene expression profile in lymphocytes from mice, rats, and humans: most potent induction in humans. Toxicology. 2006;225(2–3):204–13.
Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32(Database issue):D91–4.
Camp ER, Liu W, Fan F, Yang A, Somcio R, Ellis LM. RON, a tyrosine kinase receptor involved in tumor progression and metastasis. Ann Surg Oncol. 2005;12(4):273–81.
Yao HP, Zhou YQ, Zhang R, Wang MH. MSP-RON signalling in cancer: pathogenesis and therapeutic potential. Nat Rev Cancer. 2013;13(7):466–81.
Nukiwa T, Suzuki T, Fukuhara T, Kikuchi T. Secretory leukocyte peptidase inhibitor and lung cancer. Cancer Sci. 2008;99(5):849–55.
Mathysen D, Van Roy N, Van Hul W, Laureys G, Ambros P, Speleman F, Wuyts W. Molecular analysis of the putative tumour-suppressor gene EXTL1 in neuroblastoma patients and cell lines. Eur J Cancer. 2004;40(8):1255–61.
Zhu H, Wang G, Qian J. Transcription factors as readers and effectors of DNA methylation. Nat Rev Genet. 2016;17(9):551–65.
Ettou S, Humbrecht C, Benet B, Billot K, d'Allard D, Mariot V, Goodhardt M, Kosmider O, Mayeux P, Solary E, et al. Epigenetic control of NF-kappaB-dependent FAS gene transcription during progression of myelodysplastic syndromes. Mol Cancer Res. 2013;11(7):724–35.
Shen ES, Whitlock JP Jr. The potential role of DNA methylation in the response to 2,3,7,8-tetrachlorodibenzo-p-dioxin. J Biol Chem. 1989;264(30):17754–8.
Takahashi Y, Suzuki C, Kamataki T. Silencing of CYP1A1 expression in rabbits by DNA methylation. Biochem Biophys Res Commun. 1998;247(2):383–6.
Ehrlich M. DNA hypomethylation in cancer cells. Epigenomics. 2009;1(2):239–59.
Zheng D, Gui B, Gray KP, Tinay I, Rafiei S, Huang Q, Sweeney CJ, Kibel AS, Jia L. Secretory leukocyte protease inhibitor is a survival and proliferation factor for castration-resistant prostate cancer. Oncogene. 2016;35(36):4807–15.
Honda S, Minato M, Suzuki H, Fujiyoshi M, Miyagi H, Haruta M, Kaneko Y, Hatanaka KC, Hiyama E, Kamijo T, et al. Clinical prognostic value of DNA methylation in hepatoblastoma: four novel tumor suppressor candidates. Cancer Sci. 2016;107(6):812–9.
Yang Y, Wang JB, Li YM, Zhao YU, Wang R, Wu Q, Zheng RS, Ou YR. Role of IL-33 expression in oncogenesis and development of human hepatocellular carcinoma. Oncol Lett. 2016;12(1):429–36.
Luedde T, Schwabe RF. NF-kappaB in the liver--linking injury, fibrosis and hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2011;8(2):108–18.
Hartwell HJ, Petrosky KY, Fox JG, Horseman ND, Rogers AB. Prolactin prevents hepatocellular carcinoma by restricting innate immune activation of c-Myc in mice. Proc Natl Acad Sci U S A. 2014;111(31):11455–60.
Chen GG, Leung J, Liang NC, Li L, Wu K, Chan UP, Leung BC, Li M, Du J, Deng YF, et al. Ent-11alpha-hydroxy-15-oxo-kaur-16-en-19-oic-acid inhibits hepatocellular carcinoma in vitro and in vivo via stabilizing IkBalpha. Investig New Drugs. 2012;30(6):2210–8.
Pleyer L, Greil R. Digging deep into “dirty” drugs - modulation of the methylation machinery. Drug Metab Rev. 2015;47(2):252–79.
We wish to thank Drs Noriyuki Suzuki, E Koike and A Furuyama (NIES) for their useful discussions, H Murai for her excellent technical assistance, and K Kaneko, M Ikezawa and Y Hayakawa for their helpful secretarial assistances. The authors also wish to thank the staff of Animal Care Company (Tokyo, Japan) for their excellent assistance in the maintenance of mice.
This study was partly supported by the National Institute for Environmental Studies (1115AA082; 1315AT001; 1620AA041), and Grant-in-Aid for Scientific Research (B) (26293154, 15 K15246, KN) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.
Availability of data and materials
The datasets generated during the current study are available in the NCBI’s Gene Expression Omnibus, GSE111420 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111420), GSE104627 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104627) and GSE104626 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104626).
Ethics approval and consent to participate
Animal studies were permitted by the Animal Care and Use Committee of National Institute for Environmental Studies (NIES) and were performed in accordance with guideline for the Care and Use of Laboratory Animals of NIES.
The results on human HCC shown in Fig. 7 are based on existing data generated by The Cancer Genome Atlas project established by the NCI and NHGRI http://cancergenome.nih.gov. New experimental material was not used in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers for bisulfite sequencing and qRT-PCR. (XLSX 10 kb)
Figure S1. Scatter plots for sample pairs. Scatter plots of % methylation values for each pair in normal tissues (n = 3) and tumor tissues (n = 3) were obtained using methylKit. Numbers on upper right corner denote pair-wise Pearson’s correlation scores. The histograms on the diagonal are methylation distribution of CpG sites for each sample. Figure S2/. Expression levels of Mst1r, Slpi and Extl1 in tumor tissues with (n = 6) and without (n = 5) Ha-ras mutation. Figure S3. Hepa1-6 cells show lower expression of Dnmts compared to Hepa1c1c7. cDNA was prepared from Hepa1c1c7 cells and Hepa1-6 cells and the expression was measured by real-time PCR (n = 3). All genes were normalized to the expression of rRNA.. *** p < 0.001. (PPTX 296 kb)
Table S2. List of DMRs position and methylation difference between normal and tumor liver tissues. (XLSX 317 kb)
Table S3. Datasets of DNA methylation and gene expression downloaded from TCGA database. (XLSX 13 kb)