DNA methylation landscape of hepatoblastomas reveals arrest at early stages of liver differentiation and cancer-related alterations

INTRODUCTION

DNA methylation in mammals plays a key role in embryonic and differentiation processes through the modulation of chromatin compaction states leading to domains of variation in transcriptional activities [1]. The Developmental Origins of Health and Disease (DOHaD) hypothesis postulates that environmental factors acting in utero, together with genetic factors, may greatly influence disease susceptibility and outcomes in both infancy and adulthood [2, 3], particularly modulating alterations in DNA methylation in pediatric pathological conditions [4].

A specific class of pediatric cancers – the embryonal tumors – is defined based on embryonic features of the tumors and their very early age of onset. A relatively small number of genetic mutations have been reported in pediatric tumors compared to adult solid tumors [5, 6], and this paucity of mutations may be partially correlated with the early age of onset [5]. Another possible explanation is that mechanisms other than DNA mutation may disrupt cell differentiation processes during development, a model that has been proposed for embryonal tumors [7]. Hepatoblastomas are embryonal tumors of the liver and are the most common hepatic tumors of early childhood, with an incidence in the United States estimated to be between 5.2 and 10.5 cases per million children <1 and 1-4 years old, respectively [8]. In Brazil, an age-adjusted incidence rate of 2.8 per million children less than 15 years of age has been reported [9]. Due to their rarity, knowledge regarding genetic predisposing factors for hepatoblastoma is limited. Indeed, few genetic alterations have been found in both cytogenetic and point mutation studies of hepatoblastomas. The identified genetic alterations are primarily linked to the Wnt signaling pathway, with a relatively high frequency of CTNNB1 activating mutations [10, 11]; its role was demonstrated after insertion of CTNNB1 mutations in precursor liver cells, resulting in high incidence of hepatocellular carcinomas and hepatoblastomas in mice [12].

A less studied pathway that could partially explain both developmental disruption and cancer development is the occurrence of epigenetic changes [3, 4]. DNA methylation is a stable modification involved in many biological processes, including tissue differentiation, embryogenesis, and disease development [13]. In general, whereas normal cells usually have unmethylated gene promoters and CpG islands, and heavily methylated repetitive sequences, the opposite pattern is found in tumor cells. Anomalous DNA methylation represents an important mechanism driving tumor development and progression [14]. Although the roles of DNA methylation in oncogenesis have not yet been fully elucidated, promoter hypermethylation has been shown to be associated with inappropriate transcriptional repression of tumor suppressor genes, and hypomethylation of repetitive sequences may result in activation of retro-transposon elements and genomic instability [15, 16].

To date, most of studies investigating DNA methylation in hepatoblastomas have focused on specific genes or genomic regions. Site-specific hypermethylation of selected imprinted regions were reported in hepatoblastomas in comparison to normal liver counterparts, suggesting that aberrant DNA methylation is associated to hepatoblastoma formation [17]. Hypermethylation of the promoter regions of RASSF1A [18–20], CASP8 [18], SOCS1 [21–23], APC, CDH1, MT1G [23], HHIP [10], CDKN2A [24], and IGF2 [18, 23] were also reported. Additionally, hypomethylation of IGFBP3 [25] has been described in a subgroup of tumors. To date, only two genome-wide methylation studies in hepatoblastoma samples have been published. In one of them, the comparison of three matched tumor-normal samples revealed genes exhibiting aberrant methylation that were involved with cell adhesion, blood coagulation and nervous system development, in addition to low methylation level near the transcriptional start site of the alpha fetoprotein gene [26]. The second one reported the investigation of hepatoblastomas from two patients, and revealed four genes, GPR180, MST1R, OCIAD2 and PARP6, with methylation status that could be associated with clinical parameters, such as age at diagnosis and poor outcome [27].

In the present study, we carried out an extensive methylation analysis across the entire genome using two independent cohorts of hepatoblastomas (8 tumors paired with their adjacent normal liver tissues from the same patients, and 11 hepatoblastomas and 3 normal liver tissues from different patients). We examined the DNA methylation status of >450,000 CpG sites, as well as LINE-1 sequences throughout the entire genome. The methylation levels of hepatoblastomas were compared on a locus by locus basis with those found in non-tumoral liver samples. Furthermore, we interrogated the extent of the overlap between the differentially methylated sites (DMSs) here detected in hepatoblastomas and the CpG sites reported to exhibit methylation changes through liver development.

RESULTS

A trend toward global DNA hypomethylation was detected in hepatoblastomas

We analyzed 19 hepatoblastomas using a discovery group (Hepatoblastoma set #1) and a validation group (Hepatoblastoma set #2). This strategy was used to lessen the number of false negatives because we studied a small sample size. The clinical features of the patients with hepatoblastomas are summarized in Table 1.

Prior to the differential methylation analyses, several parameters were considered to infer surrogate variables that could impact the comparison of interest. These analyses did not detect any association with the investigated co-variables (age at diagnosis, tumor histology, alpha-fetoprotein levels, treatment protocol, progression and overall survival). The mean and standard deviation of the methylation levels (beta-values ranging from 0 to 1) were very similar within each of the four groups of comparison, Hepatoblastomas set #1 (0.49 +/- 0.37) and set #2 (0.48 +/- 0.37), and control groups of non-tumor differentiated liver set #1 (0.52 +/- 0.36) and set #2 (0.52 +/- 0.36), indicating similar levels of inter- and intra-group heterogeneity, i.e., among the different samples of hepatoblastomas and the different samples of differentiated livers.

The comparison between hepatoblastomas and their matched non-tumoral liver tissues from set #1 revealed 1,399 differentially methylated CpG sites (adjP<0.05) – DMSs -, with 820 hypomethylated CpGs (551 of them related to 457 genes), and 579 hypermethylated CpGs (456 of them related to 363 genes) (Figure 1A, left). A validation set of 12 hepatoblastomas (Hepatoblastoma set #2) was used to replicate the analysis, using as a control group the same non-tumoral liver samples plus additional three samples. The set #2 comparison exhibited 61,384 DMSs (adjP<0.05) with a large proportion of them (86.7%; 53,222 CpGs) showing hypomethylation. There was an overlap of 1,359 DMSs (97% of the CpGs detected in the paired differential methylation analysis) between both analyses (Figure 1A; Supplementary Table 1) affecting 765 genes; common DMSs were either hypermethylated or hypomethylated CpGs in both hepatoblastoma sets, with a maximum beta-value variation of 0.21, and an average deviation of 0.009. Unsupervised hierarchical clustering analysis based on the methylation levels of these common 1,359 DMSs discriminated all control livers from 18 out of 19 hepatoblastoma samples (Supplementary Figure 1), except for tumor sample HB17. We did not obtain a clear discrimination of tumor subgroups based on the clinical characteristics of the patients, but given the small size of the cohort, this observation is not conclusive.

The common set of DMSs displayed a non-random distribution regarding different categories of genomic sequences being mostly hypermethylated at gene bodies and CpG islands, and hypomethylated at intergenic regions (p<0.0001; chi-square distribution test, Table 2). Additionally, we inspected the distribution of the DMSs across the autosomes considering the paired and non-paired analyses separately, seeking to uncover a preferential gain or loss of methylation in a specific chromosome (Figure 1B). The autosomes in which the differences of the DMSs were positive (i.e. meaning that the CpG sites were more methylated in hepatoblastomas compared with non-tumoral livers) were chromosomes 6, 9, 11, 20 and 22. Therefore, most of the DMSs were hypomethylated in the majority of the chromosomes, suggesting a widespread and non-randomic pattern of hypomethylation.

DNA methylation is considered a repressive epigenetic mark that characterizes most CpG sites within mammalian genomes. However, CpG dense regions, mainly CpG islands, are generally unmethylated in normal cells [28]. Additionally, DNA methylation values represent a mix of different cell types, thus, small methylation changes may correspond to an epigenetic drift of a subset of cells. Therefore, trying to identify CpG sites with more impact for the disease as a whole, we used an additional threshold of ≥20% methylation difference in hepatoblastomas compared to control livers. We defined differentially methylated regions (DMRs) based on sequences containing ≥3 significant CpG sites, all in the same direction (hypomethylation or hypermethylation) in the same genomic location, with at least one site displaying methylation differences of ≥20%. Thirty-eight genes were found to harbor sequences that met these parameters (Table 3; Supplementary Table 2 contains the details of the CpG sites which are related to DMRs). Among these 38 genes, 14 exhibited hypermethylated CpGs, and 23 exhibited hypomethylation, and one gene, TSPAN9, showed both gene body hypermethylation and promoter hypomethylation. One of these genes, EIF4E, has already been reported in liver cancer, and other 27 are directly related to liver function.

Differentially methylated genes in hepatoblastomas are involved with liver cell differentiation and cancer

Among the 765 genes associated with the set of common DMSs, there was an enrichment for Metabolic, Cancer, Insulin, MAPK, Hedgehog and Wnt signaling pathways, as well as Regulation of the actin cytoskeleton, Adherens junction, Glycolysis/Gluconeogenesis and Drug Metabolism (Supplementary Table 3). Genes with hypomethylated sites were related to negative regulation of cell differentiation (19 genes), including somatic stem cell maintenance; negative regulation of RNA metabolic processes (32 genes); and negative regulation of DNA-dependent transcription (31 genes). Hypermethylated CpG sites were mapped in genes related to organic acid transport (14 genes), including fatty acid transport; acyl-CoA metabolic processes (8 genes); response to hormone stimulus (28 genes); response to nutrient levels (17 genes); and liver development (8 genes).

Several of the enriched pathways as well as the detected relevant biological processes (liver cell differentiation and stem cell maintenance) might be related to liver development; therefore, we compared the differentially methylated CpG sites detected in hepatoblastomas with the 28,447 CpG sites found to undergo methylation changes during liver development [29]. There were 434 differentially methylated CpG sites associated to 295 genes, this set was found to be significantly enriched during liver development (p<0.0001). Among these 434 CpG sites, 142 (17%) were hypomethylated and 292 (50%) were hypermethylated, with similar patterns of hypomethylation or hypermethylation in both hepatoblastomas and fetal livers compared with adult livers. Functional annotation of the 295 genes revealed enrichment in Metabolic, Cancer, Insulin, MAPK and Wnt signaling pathways, as well as Regulation of the actin cytoskeleton and Adherents junction. Three additional pathways, compared to the whole set of genes, stood out in relevance, namely citrate cycle, tight junction and endocytosis.

Therefore, we sought to verify the performance of the global methylation profile in hepatoblastomas in relation to the methylomes of fetal and adult liver samples using this same study [29] as independent groups of liver samples. Beta-values were averaged across samples (hepatoblastomas, fetal or adult livers) and plotted using the notch boxplot method. This method, although not a formal test, infers that if the notches of two boxplots do not overlap, then there is strong evidence (95% confidence) that their medians differ. The results revealed that hepatoblastomas display a significant global hypomethylation pattern when compared with both fetal and adult differentiated livers (Figure 2).

Upon unsupervised hierarchical clustering of our set of 1,359 DMSs, two clusters were formed: cluster 1, containing all differentiated liver samples, including those from our study except for one sample, which is a paired non-tumoral sample of a congenital hepatoblastoma (HB33) that grouped with tumors; and cluster 2, which further separated in two branches, one containing all hepatoblastoma samples with exception for HB17, which grouped with all fetal liver samples in the other branch (Figure 3).

Furthermore, we performed an in silico gene expression analysis of the 765 genes associated with the 1,359 differentially methylated CpGs using two studies which contain microarray expression data from hepatoblastomas and non-tumoral liver samples [30, 31]. In the first dataset, there were 717 out of 765 genes in common to both studies, and expression levels of these genes discriminated all non-tumoral liver samples from 94% (47 out of 50) of the hepatoblastoma samples (Supplementary Figure 2A). Using the second dataset, we found that 516 out of 765 genes were common to both studies, i.e., 516 genes harboring differentially methylated sites that were present in the expression microarray; the expression levels exhibited by these common set of genes discriminated all non-tumoral liver samples from 88% (22 out of 25) of the hepatoblastoma samples (Supplementary Figure 2B).

Profile of methylation in CpGs sites related to retrotransposon sequences (LINE-1) and X-chromosome inactivation

LINE-1 sequences are not hypomethylated in hepatoblastomas

To clarify whether the hypomethylation in hepatoblastomas extends to repetitive sequences, we used pyrosequencing to analyze the methylation status of four CpG sites mapped within LINE-1 sequences. Generally, the LINE-1 methylation levels were very similar across hepatoblastomas and normal differentiated and fetal liver samples. However, focusing on the analysis of individual CpGs, hepatoblastomas from both sets showed a slightly hypomethylated pattern only at the first CpG of the LINE-1 sequence when compared with the differentiated livers (p = 0.03802 for set #1 and p = 0.05174 for set #2; Figure 4). The copy number alteration profile for all samples of hepatoblastoma set #1 has been previously evaluated by array-CGH [32]; tumor samples with and without chromosomal alterations >100 kb were compared with control differentiated livers regarding the methylation levels of LINE-1. Significant differences in LINE-1 methylation levels were found again only for the first CpG site (Supplementary Figure 3), which is slightly hypomethylated in hepatoblastomas with copy number alterations (p = 0.0057; the Kruskall-Wallis test also indicated a difference of p = 0.0166 between the groups at this CpG).

X-chromosome CpG sites exhibit hypomethylation in male hepatoblastomas

Using all the samples from our study, the methylation levels of CpG sites mapped to the X chromosome were recovered to compare tumor and control liver samples using a sex-matched approach. The methylation profile of the X chromosome was similar in female hepatoblastomas compared with female liver controls, but hypomethylation was detected in male hepatoblastomas (Supplementary Figure 4).

DISCUSSION

Embryonal tumors are believed to evolve from the disruption of normal embryogenesis [33]. We hypothesized that methylation changes probably contribute to this process, and we analyzed DNA methylation in normal liver and hepatoblastoma samples. Although these are certainly not homogeneous, the methylation data showed minimal variance between samples and groups, indicating that the correction we made before the comparison was efficient. The paired comparison between hepatoblastomas and liver samples revealed few differentially methylated sites exhibiting low-level differences, although the non-paired analysis showed higher differences on top of higher number of sites. This second finding may be related to inherent differences between tumor samples and non-matched controls, reinforcing the contribution of the paired comparison. Notably, the smaller DMS set detected in the paired group is almost completely included in the larger DMS set detected in the non-paired analysis.

Our investigation revealed that hepatoblastomas exhibit a global low-level hypomethylation pattern compared with both normal fetal and differentiated livers. This hypomethylation was generally restricted to non-repetitive DNA, a pattern unusual for most adult solid tumors which exhibit a major hypermethylation at promoters and hypomethylation in repetitive sequences, mainly LINE-1 [34, 35]. However, similar global hypomethylation patterns have recently been reported for hepatocellular carcinoma [32], embryonal tumor medulloblastomas [36] and, lately, in three hepatoblastoma samples [26].

Methylome studies in two other embryonal tumors (Wilms tumor and neuroblastoma) reported a methylation fingerprint in tissue development genes, similar to what we observed in hepatoblastomas. Differentially methylated regions detected in Wilms tumor showed hypomethylation of renal development genes and were enriched for bivalent domains in embryonic stem cells [37]. In neuroblastomas, changes in methylation occurred in sequences related to functional chromatin domains of development and cancer-related genes [38].

Recently, a genome-wide methylation analysis on hepatoblastomas disclosed four novel tumor suppressor genes potentially related to progression in hepatoblastoma, with the increase of methylation levels in more than four genes being associated with poorer prognosis [27].

Functional enrichment analysis of the genes affected by abnormal methylation have highlighted the Wnt signaling pathway, typically altered in hepatoblastomas mainly through a high frequency of CTNNB1 activating mutations [10, 39]. Here, we found evidence that methylation changes in the genome also interfere with Wnt signaling. Moreover, genes impacted by DMSs were related to pathways that have been associated with cancer, energy and drug metabolism. Glycolysis/gluconeogenesis (GNG) are metabolic pathways reciprocally regulated in the liver that result in the generation of energy. Interestingly, restoration of gluconeogenesis through the inhibition of mTOR appears to improve therapeutic efficacy against hepatocarcinomas [40]. In addition, pathways in cancer, endocytosis and adherens junction were cellular signaling pathways in common previously reported to be deregulated in hepatoblastomas at methylation levels [26]. Regarding the biological processes likely to be impaired by methylation changes, we identified involvement of genes that maintain the interface between the environment and the phenotype, such as those related to the maintenance of nutrient levels and response to hormone stimuli.

To date, it is not clear to what extent the cancer-related epigenetic alterations detected in this study correspond to natural epigenetic variability in human livers; we are aware that fully differentiated liver samples are not an ideal control for methylation studies of embryonal tumors. Therefore, we also compared our data to a recently published study that analyzed methylation changes occurring during liver differentiation [29]. Even with different approaches to determine the DMSs in the two studies, we found an overlap in the frequency and distribution of hypomethylated and hypermethylated CpGs in hepatoblastomas/fetal livers compared with fully differentiated adult livers. The function of such shared CpG sites delineates a stemness profile for hepatoblastoma cells. The respective methylation changes may be related to the postulated impairment of liver cell differentiation leading to the development of hepatic embryonal tumor. It was reported a trend towards hypomethylation in fetal compared to adult liver samples in all CpG sites [29]. Upon clustering, the set of 1,359 DMSs detected in hepatoblastomas clearly showed that tumors and normal fetal livers clustered together into one branch, and all adult liver tissues into another, a finding that supports the hypothesis that hepatoblastomas exhibit a methylation pattern resembling the early stages of liver development. It is likely that hepatoblastoma cells have become fixed in time and space also due to methylation changes, i.e. they have lost their ability to differentiate further into adult liver cells. Adding to this scenario, hepatoblastomas also were reported to display an expression pattern comparable to fetal livers [30].

In our analysis, the signaling pathways related to liver development that may be disturbed by methylation changes in hepatoblastomas were Metabolism, Cancer, Insulin, MAPK and Wnt signaling pathways, as well as Regulation of the actin cytoskeleton and Adherents junction. Recovering publicly available expression data from two independent studies, we showed that the expression levels of the genes associated with the detected differentially methylated CpG sites were able to discriminate all non-tumoral livers from the great majority of the hepatoblastoma samples [30, 31]. This result supports the impact of our methylation findings on tumor gene expression as well as reinforces the role of the detected pathways in the tumorigenic process of hepatoblastoma. Altogether, these data suggest that DNA methylation has a high impact on the expression profile for hepatoblastoma development.

Five of the genes with altered methylation have been previously reported as altered in liver carcinogenesis. HSD17B13 does not have a known function, but was suggested as a biomarker for hepatocellular carcinomas [41]. DKK3 is part of the Wnt signaling pathway and is overexpressed in hepatoblastomas, regardless of tumor histology [42]. MYC is a transcriptional factor that contributes to liver tumor maintenance [43]. HAL is regulated by glucocorticoids and glucagon, and catalyzes the first reaction in histidine metabolism via preferentially the Protein Kinase A (PKA) pathway [44]. IGF2 is a growth factor frequently overexpressed in hepatoblastomas due to loss of imprinting at 11p15 or paternal uniparental disomy at 11p15 [21].

We found two striking cases in which we believe that exome sequencing would help to further clarify their tumorigenesis processes. One hepatoblastoma (HB17) grouped with the normal liver samples, displaying a similar distribution of methylation levels; one can speculate that in this particular tumor, the occurrence of driver mutations could have triggered the malignant transformation early during liver cells differentiation. Additionally, a non-tumoral liver sample paired with one hepatoblastoma exhibited a methylation pattern more similar to undifferentiated liver cells, grouping with the fetal rather than the differentiated livers, and near tumors in the clustering. This sample was derived from one congenital tumor, suggesting that the methylation changes have been previously established in the non-tumoral liver as a consequence of a highly penetrant germline mutation, that might be related to the epigenetic machinery.

Three genes identified in our hepatoblastoma analyses were also listed in the top 20 differentially methylated genes during liver development: CRP, NNMT (both hypermethylated) and C3P1 (hypomethylated). Among them, NNMT stood out as an enzyme (nicotinamide N-methyltransferase) that regulates hepatic nutrient metabolism [45]. NNMT has also been reported to be overexpressed in a variety of human cancers and to promote epigenetic remodeling in cancer by consuming methyl units from S-adenosyl methionine (SAM), resulting in a metabolic sink [46]. Recently, it has been shown that NNMT and the metabolic state regulate pluripotency in human embryonic stem cells (hESCs) through the consumption of SAM in naive cells, making it unavailable for histone methylation, which represses Wnt and activates the HIF pathway [47]. According to authors, their findings support the hypothesis that the metabolome regulates the epigenetic landscape of the earliest steps in human development. Therefore, we speculate that NNMT dysregulation caused by epigenetic changes may be a suitable mechanism linked to the origin of hepatoblastomas, thus deserving further investigation.

Subsequently, we hypothesized that the remaining DMSs detected in hepatoblastomas might be related to cell transformation and tumor progression rather than the maintenance of an undifferentiated cell state. The disease enrichment analyses revealed that ten of these remaining genes have previously been found in liver neoplasms: THRSP, HEPACAM [48], MACC1 [49], TERT [50], NUAK1 [51], NDRG1 [52], RASSF1 [53], PRDM2 [54], ALDOB [55] and MT1G [23], with MACC1, NDRG1 and RASSF1 already associated with prognosis in liver cancers.

Altogether, we suggest that DNA hypomethylation of intergenic sequences that were primarily detected in hepatoblastomas, as well as the other methylation changes, contributes to the impairment of cell differentiation, which seems to be required for tumorigenesis in embryonal tumors, as reported for meduloblastoma [36], Wilms tumor and neuroblastomas [37, 38]. Alternatively, it may represent an inherent aspect of hepatic tumorigenesis [56]. Both hypotheses are supported by our findings of enrichment for methylation changes occurring in genes relevant for liver differentiation.

Cytosine methylation not only contributes to tissue-specific gene expression but also is associated with the silencing of retrotransposons, genomic imprinting and X-chromosome inactivation in mammals. Therefore, to complement the methylation profile of hepatoblastomas, these aspects were explored. LINE-1 is a retrotransposon repetitive element that constitutes approximately 17% of the human genome. It has been reported to be hypomethylated in many adult tumors, and its methylation level correlates with clinical/pathological features [57, 58]. No differences were found in LINE-1 methylation of hepatoblastomas compared with differentiated livers except for the first LINE-1 CpG, which showed a minor though significant reduction in methylation, similarly to results previously described [17]. Additionally, in our data, low-level hypomethylation of LINE-1 were associated with the occurrence of copy number alterations; when comparing tumors with and without copy number alterations, hypomethylation of the first CpG of LINE-1 was mainly detected in hepatoblastomas harboring chromosomal alterations. However, the absence of significant hypomethylation of repetitive sequences in hepatoblastomas clearly contrasts with most of the solid adult tumors. Therefore, the low-level global hypomethylation reported here is restricted to non-repetitive sequences of the tumors, including intergenic sequences. Such a hypomethylation pattern raises the possibility that intergenic regulatory elements may be affected, as has recently been reported for acute lymphoblastic leukemia [59].

No differences were detected between female hepatoblastomas and control livers regarding the CpG sites located on the X chromosome, suggesting that the X-inactivation, which takes place in the multipotent cell stage during very early embryonic development, was already established when the tumor precursor cells were laid down, and preserved. Conversely, male hepatoblastomas showed a significant loss of methylation at X chromosome CpGs, an observation that needs to be clarified, but that are in accordance with the hypomethylation trend detected in hepatoblastomas in general.

Our results pinpoint that DNA methylation, mostly a low-level hypomethylation, is a key epigenetic mechanism related to hepatoblastoma tumorigenesis, in addition to genetic mutations. It remains to be addressed in future studies whether there is an active process of demethylation working on hepatoblastoma origin, possibly involving TET enzymes and 5hmC, or alternatively, whether a passive pathway of demethylation occurs via inhibition of DNMTs, a possibility previously outlined. Low basal levels of DNMT1 in hepatoblastoma cell lines confirmed an essential role of DNMT1 depletion in the enhancement of cancer stem cell properties [60]. The authors have demonstrated that epigenetic reprogramming induced by transient DNMT1 inhibition influences both malignant properties and the pool of hepatic cancer stem cells.

CONCLUSIONS

Our results showed that hepatoblastomas exhibit a global low-level hypomethylation pattern in non-repetitive sequences, and an apparent arrest at early stages of liver differentiation. This suggests either that embryonal tumors are driven by different oncogenic mechanisms, as has been observed in other embryonal tumors, or that hypomethylation is a major feature of liver carcinogenesis because it has also been reported in hepatocellular carcinoma. The pathway analysis of the differentially methylated genes revealed an enrichment for metabolism, negative regulation of cell differentiation, and cancer. Moreover, confirming previous genetic findings, the Wnt signaling pathway was highlighted, suggesting a non-stochastic mechanism of CpG methylation in hepatoblastomas. In addition to the methylation landscape, we also provided a list of potential candidate genes for hepatoblastomas.

MATERIALS AND METHODS

Human tissue samples

Hepatoblastomas and differentiated non-tumoral liver samples were provided by the Biobank of the A. C. Camargo Cancer Center (ACCCC) and the Pediatric Oncology Institute (GRAACC), which are both cancer hospitals in São Paulo, Brazil. The institutional review board of each institution independently approved the study, and informed consent was obtained from the patients’ legal guardians. All procedures followed the guidelines of the Declaration of Helsinki. All patients received pre-surgical chemotherapy according to SIOPEL (http://www.siopel.org/) or COG (http://www.childrensoncologygroup.org/) protocols. Patients were followed by clinical examination, imaging tests and measurement of alpha-fetoprotein for a minimum period of 18 months.

Sections from tumor tissue blocks stained with hematoxylin and eosin were reviewed by a pathologist. Samples were grouped into a training set (consisting of 8 cases with paired hepatoblastoma and differentiated non-tumoral liver samples; described here as Hepatoblastoma set#1) and a validation set (11 hepatoblastomas and 3 unrelated non-tumoral liver samples - one patient with undisclosed age, the other two patients were 4 and 13 years old; described here as Hepatoblastoma set#2).

Genomic DNA isolation

Samples were macrodissected for enrichment of neoplastic (at least 70% of homogeneity) or normal tissues, followed by DNA extraction performed at the Macromolecule Bank (A. C. Camargo Cancer Center). Genomic DNA was extracted using phenol:chloroform. NanoDrop (Thermo Fisher Scientific) and Qubit (Thermo Fisher Scientific) equipment were used to assess DNA purity and quantity, and DNA quality was checked by electrophoresis in 0.8% agarose gels.

DNA methylation analysis

Bisulfite conversion of 500 ng of DNA was performed using the EZ DNA Methylation kit (Zymo Research) according to the manufacturer’s recommendations. Bisulfite-converted DNA samples were hybridized in the Human Methylation 450 BeadChip microarrays (HM450K, Illumina), following the Illumina Infinium HD methylation protocol. The HM450K platform measures the DNA methylation level of 485,577 loci distributed across the genome at single-nucleotide resolution. Data was submitted to NCBI/GEO and is available as GSE78732.

The Illumina iScan SQ scanner (Illumina) was used to obtain images of the microarrays. The fluorescence signals were interpreted with the GenomeStudio software (v.2011.1) with the methylation module v.1.9.0 (Illumina). Probes were annotated according to the Illumina annotation file using UCSC version hg19 of the human reference genome. The methylation levels for each CpG probe were provided as beta-values ranging from 0 to 1 (0 indicating unmethylated CpGs, and 1 indicating fully methylated CpGs).

The RnBeads package [61] was applied to the dataset. A total of 2,001 unreliable probes were identified using Greedycut (p>0.05), and removed. Probes located in SNPs (4,713 probes) as well as those that do not address CpG methylation levels (non-CpG control probes; n=3,140) were also excluded. The background was corrected using the noob method, which is based on a normal-exponential convolution using out-of-band probes [62]. Signal intensities from type I and II probes were normalized using the SWAN method, which adjusts intensities based on a quantile approach [63]. Probes located on the X and Y chromosomes were excluded from downstream analyses, unless otherwise disclosed, resulting in a total of 422,324 normalized β-values for each of the 29 samples (19 hepatoblastomas and 10 control livers). The SVA package [64], that estimates and removes unwanted sources of variation in -omics data, was applied to infer surrogate variables from patient's characteristics (age at diagnosis, tumor histology, alpha-fetoprotein level, treatment protocol, progression and overall survival). Principal Component Analyses (PCA) showed no association between these variables and the groups of comparison regarding methylation levels. However, adjustments for technical effects, such as batch, were performed using default parameters of RnBeads to avoid bias.

Differential methylation analysis

Beta-values were transformed into M-values before performing comparison between groups, employing an empirical Bayesian framework linear model from limma [65]. CpG sites with adjusted p-values (adjP) <0.05 as defined by Benjamini and Hochberg’s method were considered to be differentially methylated. Comparisons were performed in two steps: a paired analysis in 8 cases (Hepatoblastoma set #1), and unpaired analyses comparing an independent set of 11 hepatoblastomas with 3 unpaired non-tumor differentiated liver samples (Hepatoblastoma set #2) plus the same set of 8 non-tumor differentiated liver samples used in the first comparison. The differentially methylated CpG sites (DMSs) that were common to both comparisons (paired and unpaired) were used for downstream analyses. Notch box plots were used to access the distribution of methylation levels of both the DMSs across the autosomal chromosomes, and violin plots were used to evaluate the methylation levels in the X chromosomes of males and females in both cases and controls.

Comparison between wild-type (n=15) and CTNNB1-mutated (n=4) tumors were performed using the same parameters. We identified a weak correlation between age at diagnosis (0.48) and alpha-fetoprotein levels (-0.53) with group comparison. Accordingly, adjustment for the M-values of each sample were performed before differential methylation analyses at CpG site levels.

Functional enrichment analyses (Biological Processes from Gene Ontology and cellular signaling pathways from KEGG) were performed using WebGestalt [62] with the whole genome as the background. Features with adjP<0.05 (hypergeometric test with Benjamini–Hochberg adjustment) comprising at least five genes were considered significant.

Hepatoblastomas and liver development

Next, we compared our data to CpGs reported to exhibit methylation changes during liver development [29]. In the mentioned study, HM450K arrays were used to identify methylation changes between 14 fetal and 96 adult livers, revealing 28,447 differentially methylated CpG sites. These data were compared with the detected DMSs from our study, taking into consideration whether the sites were hypomethylated and hypermethylated. In addition, using Modular Single-set Enrichment Test (MSET), we assessed the enrichment of our DMSs to the differentially methylated CpGs between fetal and adult liver samples [29], generating an in silico p-value based on 10,000 random datasets [66].

Furthermore, to obtain a methylation profile of fetal and adult livers, the IDAT files from GEO (GSE61279) were downloaded, preprocessed and normalized using the same analytical procedure applied to our study. There were 420,432 common B-values for both datasets, ours and the afore mentioned liver development study. Average B-values from fetal and adult livers were compared with the average B-values of hepatoblastomas sets #1 and #2.

Pearson correlation with complete linkage was used for non-supervised hierarchical clustering of hepatoblastomas and differentiated and fetal liver samples; reliability was assessed by bootstrapping using multiExperiment Viewer (MeV) software [67].

Differential gene expression analysis

The 765 genes associated with the 1,359 differentially methylated CpGs detected in our study were searched in two different datasets of expression microarrays, one of them downloaded from ArrayExpress (http://www.ebi.ac.uk/arrayexpress/experiments/E-MEXP-1851/), which contains data of 25 hepatoblastomas and four non-tumoral liver samples [30], and the second set recovered from GEO (http://www.ncbi.nlm.nih.gov/geo/; GSE75271), corresponding to expression data from 50 hepatoblastomas and five non-tumoral liver samples [31]. Expression levels representing the set of genes common to our study were used for a non-supervised hierarchical clustering of hepatoblastomas and non-tumoral liver samples using Pearson correlation with average linkage.

Authors’ contributions

Authors with substantial contributions to the conception and design; or the acquisition, analysis, or interpretation of the data: MM, TCR, AYK, ESSA, TFMA, LRV, PLP, ACVK, CR.

Authors that drafted the article: MM, TCR, CR, ACVK.

Authors that made a critical revision for important intellectual content: CMLC, IWC, MC, HB, SRCT, DMC.

To be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved: CMLC, IWC, MC, HB, SRCT

All authors read and approved the final version of this manuscript.

ACKNOWLEDGMENTS

We would like to thank the Biobank of the A. C. Camargo Cancer Center and the GRAACC (both from São Paulo, Brazil) for providing tissues and DNA samples. We are very grateful to the patients and their families for their precious collaboration.

CONFLICTS OF INTEREST

The authors disclose that they do not have any conflicts of interest.

FUNDING

The present study was supported by grants from FAPESP (2009/00898-1, 2011/24007-9, 2013/08028-1; 2016/04785-0), and CNPq (470446_2013-7). MM was supported by FAPESP 2015/06281-7.

REFERENCES

1. Reik W. Stability and flexibility of epigenetic gene regulation in mammalian development. Nature. 2007; 447:425–32. https://doi.org/10.1038/nature05918.

2. Barker DJ. The origins of the developmental origins theory. J Intern Med. 2007; 261:412–17. https://doi.org/10.1111/j.1365-2796.2007.01809.x.

3. Hanson MA, Gluckman PD. Developmental origins of health and disease: new insights. Basic Clin Pharmacol Toxicol. 2008; 102:90–93. https://doi.org/10.1111/j.1742-7843.2007.00186.x.

4. Armstrong DA, Lesseur C, Conradt E, Lester BM, Marsit CJ. Global and gene-specific DNA methylation across multiple tissues in early infancy: implications for children’s health research. FASEB J. 2014; 28:2088–97. https://doi.org/10.1096/fj.13-238402.

5. Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr, Kinzler KW. Cancer genome landscapes. Science. 2013; 339:1546–58. https://doi.org/10.1126/science.1235122.

6. Zhang J, Walsh MF, Wu G, Edmonson MN, Gruber TA, Easton J, Hedges D, Ma X, Zhou X, Yergeau DA, Wilkinson MR, Vadodaria B, Chen X, et al. Germline Mutations in Predisposition Genes in Pediatric Cancer. N Engl J Med. 2015; 373:2336–46. https://doi.org/10.1056/NEJMoa1508054.

7. Maris JM, Denny CT. Focus on embryonal malignancies. Cancer Cell. 2002; 2:447–50. https://doi.org/10.1016/S1535-6108(02)00206-4.

8. Howlader N, Noone AM, Krapcho M, Neyman N, Aminou R, Waldron W, Altekruse SF, Kosary CL, Ruhl J, Tatalovich Z, Cho H, Mariotto A, Eisner MP, et al. (eds). SEER Cancer Statistics Review, 1975-2008, National Cancer Institute. Bethesda, MD, http://seer.cancer.gov/csr/1975_2008/, based on November 2010 SEER data submission, posted to the SEER web site, 2011.

9. de Camargo B, de Oliveira Ferreira JM, de Souza Reis R, Ferman S, de Oliveira Santos M, Pombo-de-Oliveira MS. Socioeconomic status and the incidence of non-central nervous system childhood embryonic tumours in Brazil. BMC Cancer. 2011; 11:160. https://doi.org/10.1186/1471-2407-11-160.

10. Eichenmüller M, Trippel F, Kreuder M, Beck A, Schwarzmayr T, Häberle B, Cairo S, Leuschner I, von Schweinitz D, Strom TM, Kappler R. The genomic landscape of hepatoblastoma and their progenies with HCC-like features. J Hepatol. 2014; 61:1312–20. https://doi.org/10.1016/j.jhep.2014.08.009.

11. Udatsu Y, Kusafuka T, Kuroda S, Miao J, Okada A. High frequency of beta-catenin mutations in hepatoblastoma. Pediatr Surg Int. 2001; 17:508–12. https://doi.org/10.1007/s003830000576.

12. Mokkapati S, Niopek K, Huang L, Cunniff KJ, Ruteshouser EC, deCaestecker M, Finegold MJ, Huff V. β-catenin activation in a novel liver progenitor cell type is sufficient to cause hepatocellular carcinoma and hepatoblastoma. Cancer Res. 2014; 74:4515–25. https://doi.org/10.1158/0008-5472.CAN-13-3275.

13. Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, Edsall L, Antosiewicz-Bourget J, Stewart R, et al. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009; 462:315–22. https://doi.org/10.1038/nature08514.

14. Robertson KD. DNA methylation and human disease. Nat Rev Genet. 2005; 6:597–610. https://doi.org/10.1038/nrg1655.

15. Baylin SB, Jones PA. A decade of exploring the cancer epigenome - biological and translational implications. Nat Rev Cancer. 2011; 11:726–34. https://doi.org/10.1038/nrc3130.

16. Chatterjee R, Vinson C. CpG methylation recruits sequence specific transcription factors essential for tissue specific gene expression. Biochim Biophys Acta. 2012; 1819:763–70. https://doi.org/10.1016/j.bbagrm.2012.02.014.

17. Rumbajan JM, Maeda T, Souzaki R, Mitsui K, Higashimoto K, Nakabayashi K, Yatsuki H, Nishioka K, Harada R, Aoki S, Kohashi K, Oda Y, Hata K, et al. Comprehensive analyses of imprinted differentially methylated regions reveal epigenetic and genetic characteristics in hepatoblastoma. BMC Cancer. 2013; 13:608. https://doi.org/10.1186/1471-2407-13-608.

18. Honda S, Haruta M, Sugawara W, Sasaki F, Ohira M, Matsunaga T, Yamaoka H, Horie H, Ohnuma N, Nakagawara A, Hiyama E, Todo S, Kaneko Y. The methylation status of RASSF1A promoter predicts responsiveness to chemotherapy and eventual cure in hepatoblastoma patients. Int J Cancer. 2008; 123:1117–25. https://doi.org/10.1002/ijc.23613.

19. Honda S, Miyagi H, Suzuki H, Minato M, Haruta M, Kaneko Y, Hatanaka KC, Hiyama E, Kamijo T, Okada T, Taketomi A. RASSF1A methylation indicates a poor prognosis in hepatoblastoma patients. Pediatr Surg Int. 2013; 29:1147–52. https://doi.org/10.1007/s00383-013-3371-z.

20. Sugawara W, Haruta M, Sasaki F, Watanabe N, Tsunematsu Y, Kikuta A, Kaneko Y. Promoter hypermethylation of the RASSF1A gene predicts the poor outcome of patients with hepatoblastoma. Pediatr Blood Cancer. 2007; 49:240–49. https://doi.org/10.1002/pbc.21031.

21. Honda S, Arai Y, Haruta M, Sasaki F, Ohira M, Yamaoka H, Horie H, Nakagawara A, Hiyama E, Todo S, Kaneko Y. Loss of imprinting of IGF2 correlates with hypermethylation of the H19 differentially methylated region in hepatoblastoma. Br J Cancer. 2008; 99:1891–99. https://doi.org/10.1038/sj.bjc.6604754.

22. Nagai H, Naka T, Terada Y, Komazaki T, Yabe A, Jin E, Kawanami O, Kishimoto T, Konishi N, Nakamura M, Kobayashi Y, Emi M. Hypermethylation associated with inactivation of the SOCS-1 gene, a JAK/STAT inhibitor, in human hepatoblastomas. J Hum Genet. 2003; 48:65–69. https://doi.org/10.1007/s100380300008.

23. Sakamoto LH, DE Camargo B, Cajaiba M, Soares FA, Vettore AL. MT1G hypermethylation: a potential prognostic marker for hepatoblastoma. Pediatr Res. 2010; 67:387–93. https://doi.org/10.1203/PDR.0b013e3181d01863.

24. Shim YH, Park HJ, Choi MS, Kim JS, Kim H, Kim JJ, Jang JJ, Yu E. Hypermethylation of the p16 gene and lack of p16 expression in hepatoblastoma. Mod Pathol. 2003; 16:430–36. https://doi.org/10.1097/01.MP.0000066799.99032.A7.

25. Regel I, Eichenmüller M, Joppien S, Liebl J, Häberle B, Müller-Höcker J, Vollmar A, von Schweinitz D, Kappler R. IGFBP3 impedes aggressive growth of pediatric liver cancer and is epigenetically silenced in vascular invasive and metastatic tumors. Mol Cancer. 2012; 11:9. https://doi.org/10.1186/1476-4598-11-9.

26. Cui X, Liu B, Zheng S, Dong K, Dong R. Genome-wide analysis of DNA methylation in hepatoblastoma tissues. Oncol Lett. 2016; 12:1529–34. https://doi.org/10.3892/ol.2016.4789.

27. Honda S, Minato M, Suzuki H, Fujiyoshi M, Miyagi H, Haruta M, Kaneko Y, Hatanaka KC, Hiyama E, Kamijo T, Okada T, Taketomi A. Clinical prognostic value of DNA methylation in hepatoblastoma: four novel tumor suppressor candidates. Cancer Sci. 2016; 107:812–19. https://doi.org/10.1111/cas.12928.

28. Deaton AM, Bird A. CpG islands and the regulation of transcription. Genes Dev. 2011; 25:1010–22. https://doi.org/10.1101/gad.2037511.

29. Bonder MJ, Kasela S, Kals M, Tamm R, Lokk K, Barragan I, Buurman WA, Deelen P, Greve JW, Ivanov M, Rensen SS, van Vliet-Ostaptchouk JV, Wolfs MG, et al. Genetic and epigenetic regulation of gene expression in fetal and adult human livers. BMC Genomics. 2014; 15:860. https://doi.org/10.1186/1471-2164-15-860.

30. Cairo S, Armengol C, De Reyniès A, Wei Y, Thomas E, Renard CA, Goga A, Balakrishnan A, Semeraro M, Gresh L, Pontoglio M, Strick-Marchand H, Levillayer F, et al. Hepatic stem-like phenotype and interplay of Wnt/beta-catenin and Myc signaling in aggressive childhood liver cancer. Cancer Cell. 2008; 14:471–84. https://doi.org/10.1016/j.ccr.2008.11.002.

31. Sumazin P, Chen Y, Treviño LR, Sarabia SF, Hampton OA, Patel K, Mistretta TA, Zorman B, Thompson P, Heczey A, Comerford S, Wheeler DA, Chintagumpala M, et al. Genomic analysis of hepatoblastoma identifies distinct molecular and prognostic subgroups. Hepatology. 2017; 65:104–21. https://doi.org/10.1002/hep.28888.

32. Rodrigues TC, Fidalgo F, da Costa CM, Ferreira EN, da Cunha IW, Carraro DM, Krepischi AC, Rosenberg C. Upregulated genes at 2q24 gains as candidate oncogenes in hepatoblastomas. Future Oncol. 2014; 10:2449–57. https://doi.org/10.2217/fon.14.149.

33. Rushton J, López-Terrada D. Molecular and genetic basis of childhood cancer. Cancer Biomark. 2010; 9:211–34. https://doi.org/10.3233/CBM-2011-0199.

34. Heyn H, Esteller M. DNA methylation profiling in the clinic: applications and challenges. Nat Rev Genet. 2012; 13:679–92. https://doi.org/10.1038/nrg3270.

35. Miousse IR, Koturbash I. The Fine LINE: Methylation Drawing the Cancer Landscape. BioMed Res Int. 2015; 2015:131547. https://doi.org/10.1155/2015/131547.

36. Hovestadt V, Remke M, Kool M, Pietsch T, Northcott PA, Fischer R, Cavalli FM, Ramaswamy V, Zapatka M, Reifenberger G, Rutkowski S, Schick M, Bewerunge-Hudler M, et al. Robust molecular subgrouping and copy-number profiling of medulloblastoma from small amounts of archival tumour material using high-density DNA methylation arrays. Acta Neuropathol. 2013; 125:913–16. https://doi.org/10.1007/s00401-013-1126-5.

37. Charlton J, Williams RD, Sebire NJ, Popov S, Vujanic G, Chagtai T, Alcaide-German M, Morris T, Butcher LM, Guilhamon P, Beck S, Pritchard-Jones K. Comparative methylome analysis identifies new tumour subtypes and biomarkers for transformation of nephrogenic rests into Wilms tumour. Genome Med. 2015; 7:11. https://doi.org/10.1186/s13073-015-0136-4.

38. Gómez S, Castellano G, Mayol G, Suñol M, Queiros A, Bibikova M, Nazor KL, Loring JF, Lemos I, Rodríguez E, de Torres C, Mora J, Martín-Subero JI, Lavarino C. DNA methylation fingerprint of neuroblastoma reveals new biological and clinical insights. Epigenomics. 2015; 7:1137–53. https://doi.org/10.2217/epi.15.49.

39. Kosaki R, Takenouchi T, Takeda N, Kagami M, Nakabayashi K, Hata K, Kosaki K. Somatic CTNNB1 mutation in hepatoblastoma from a patient with Simpson-Golabi-Behmel syndrome and germline GPC3 mutation. Am J Med Genet A. 2014; 164A:993–97. https://doi.org/10.1002/ajmg.a.36364.

40. Khan MW, Biswas D, Ghosh M, Mandloi S, Chakrabarti S, Chakrabarti P. mTORC2 controls cancer cell survival by modulating gluconeogenesis. Cell Death Dis. 2015; 1:15016. https://doi.org/10.1038/cddiscovery.2015.16.

41. Xing X, Huang Y, Wang S, Chi M, Zeng Y, Chen L, Li L, Zeng J, Lin M, Han X, Liu X, Liu J. Comparative analysis of primary hepatocellular carcinoma with single and multiple lesions by iTRAQ-based quantitative proteomics. J Proteomics. 2015; 128:262–71. https://doi.org/10.1016/j.jprot.2015.08.007.

42. Pei Y, Kano J, Iijima T, Morishita Y, Inadome Y, Noguchi M. Overexpression of Dickkopf 3 in hepatoblastomas and hepatocellular carcinomas. Virchows Arch. 2009; 454:639–46. https://doi.org/10.1007/s00428-009-0772-4.

43. Kress TR, Pellanda P, Pellegrinet L, Bianchi V, Nicoli P, Doni M, Recordati C, Bianchi S, Rotta L, Capra T, Ravà M, Verrecchia A, Radaelli E, et al. Identification of MYC-Dependent Transcriptional Programs in Oncogene-Addicted Liver Tumors. Cancer Res. 2016; 76:3463–72. https://doi.org/10.1158/0008-5472.CAN-16-0316.

44. Alemán G, Ortíz V, Langley E, Tovar AR, Torres N. Regulation by glucagon of the rat histidase gene promoter in cultured rat hepatocytes and human hepatoblastoma cells. Am J Physiol Endocrinol Metab. 2005; 289:E172–79. https://doi.org/10.1152/ajpendo.00584.2004.

45. Hong S, Moreno-Navarrete JM, Wei X, Kikukawa Y, Tzameli I, Prasad D, Lee Y, Asara JM, Fernandez-Real JM, Maratos-Flier E, Pissios P. Nicotinamide N-methyltransferase regulates hepatic nutrient metabolism through Sirt1 protein stabilization. Nat Med. 2015; 21:887–94. https://doi.org/10.1038/nm.3882.

46. Ferreri AJ, Illerhaus G, Zucca E, Cavalli F, and International Extranodal Lymphoma Study Group. Flows and flaws in primary central nervous system lymphoma. Nat Rev Clin Oncol. 2010; 7:1038. https://doi.org/10.1038/nrclinonc.2010.9-c1.

47. Sperber H, Mathieu J, Wang Y, Ferreccio A, Hesson J, Xu Z, Fischer KA, Devi A, Detraux D, Gu H, Battle SL, Showalter M, Valensisi C, et al. The metabolome regulates the epigenetic landscape during naive-to-primed human embryonic stem cell transition. Nat Cell Biol. 2015; 17:1523–35. https://doi.org/10.1038/ncb3264.

48. Chung Moh M, Hoon Lee L, Shen S. Cloning and characterization of hepaCAM, a novel Ig-like cell adhesion molecule suppressed in human hepatocellular carcinoma. J Hepatol. 2005; 42:833–41. https://doi.org/10.1016/j.jhep.2005.01.025.

49. Xie C, Wu J, Yun J, Lai J, Yuan Y, Gao Z, Li M, Li J, Song L. MACC1 as a prognostic biomarker for early-stage and AFP-normal hepatocellular carcinoma. PLoS One. 2013; 8:e64235. https://doi.org/10.1371/journal.pone.0064235.

50. Kim H, Yoo JE, Cho JY, Oh BK, Yoon YS, Han HS, Lee HS, Jang JJ, Jeong SH, Kim JW, Park YN. Telomere length, TERT and shelterin complex proteins in hepatocellular carcinomas expressing “stemness”-related markers. J Hepatol. 2013; 59:746–52. https://doi.org/10.1016/j.jhep.2013.05.011.

51. Cui J, Yu Y, Lu GF, Liu C, Liu X, Xu YX, Zheng PY. Overexpression of ARK5 is associated with poor prognosis in hepatocellular carcinoma. Tumour Biol. 2013; 34:1913–18. https://doi.org/10.1007/s13277-013-0735-x.

52. Cheng J, Xie HY, Xu X, Wu J, Wei X, Su R, Zhang W, Lv Z, Zheng S, Zhou L. NDRG1 as a biomarker for metastasis, recurrence and of poor prognosis in hepatocellular carcinoma. Cancer Lett. 2011; 310:35–45. https://doi.org/10.1016/j.canlet.2011.06.001.

53. Zhang C, Li J, Huang T, Duan S, Dai D, Jiang D, Sui X, Li D, Chen Y, Ding F, Huang C, Chen G, Wang K. Meta-analysis of DNA methylation biomarkers in hepatocellular carcinoma. Oncotarget. 2016; 7:81255–67. https://doi.org/10.18632/oncotarget.13221.

54. Piao GH, Piao WH, He Y, Zhang HH, Wang GQ, Piao Z. Hyper-methylation of RIZ1 tumor suppressor gene is involved in the early tumorigenesis of hepatocellular carcinoma. Histol Histopathol. 2008; 23:1171–75.

55. Wu J, Dan C, Zhao HB, Xiao CX, Liu YP, Si LJ, Ren JL, Guleng B. ALDOB acts as a novel HBsAg-binding protein and its coexistence inhibits cisplatin-induced HepG2 cell apoptosis. Crit Rev Eukaryot Gene Expr. 2014; 24:181–91. https://doi.org/10.1615/CritRevEukaryotGene Expr.2014010087.

56. Shen J, Wang S, Zhang YJ, Wu HC, Kibriya MG, Jasmine F, Ahsan H, Wu DP, Siegel AB, Remotti H, Santella RM. Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips. Epigenetics. 2013; 8:34–43. https://doi.org/10.4161/epi.23062.

57. Baba Y, Huttenhower C, Nosho K, Tanaka N, Shima K, Hazra A, Schernhammer ES, Hunter DJ, Giovannucci EL, Fuchs CS, Ogino S. Epigenomic diversity of colorectal cancer indicated by LINE-1 methylation in a database of 869 tumors. Mol Cancer. 2010; 9:125. https://doi.org/10.1186/1476-4598-9-125.

58. Beck CR, Collier P, Macfarlane C, Malig M, Kidd JM, Eichler EE, Badge RM, Moran JV. LINE-1 retrotransposition activity in human genomes. Cell. 2010; 141:1159–70. https://doi.org/10.1016/j.cell.2010.05.021.

59. Almamun M, Levinson BT, van Swaay AC, Johnson NT, McKay SD, Arthur GL, Davis JW, Taylor KH. Integrated methylome and transcriptome analysis reveals novel regulatory elements in pediatric acute lymphoblastic leukemia. Epigenetics. 2015; 10:882–90. https://doi.org/10.1080/15592294.2015.1078050.

60. Clark AT. DNA methylation remodeling in vitro and in vivo. Curr Opin Genet Dev. 2015; 34:82–87. https://doi.org/10.1016/j.gde.2015.09.002.

61. Assenov Y, Müller F, Lutsik P, Walter J, Lengauer T, Bock C. Comprehensive analysis of DNA methylation data with RnBeads. Nat Methods. 2014; 11:1138–40. https://doi.org/10.1038/nmeth.3115.

62. Triche TJ Jr, Weisenberger DJ, Van Den Berg D, Laird PW, Siegmund KD. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucleic Acids Res. 2013; 41:e90. https://doi.org/10.1093/nar/gkt090.

63. Maksimovic J, Gordon L, Oshlack A. SWAN: subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol. 2012; 13:R44. https://doi.org/10.1186/gb-2012-13-6-r44.

64. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012; 28:882–83. https://doi.org/10.1093/bioinformatics/bts034.

65. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43:e47. https://doi.org/10.1093/nar/gkv007.

66. Eisinger BE, Saul MC, Driessen TM, Gammie SC. Development of a versatile enrichment analysis tool reveals associations between the maternal brain and mental health disorders, including autism. BMC Neurosci. 2013; 14:147. https://doi.org/10.1186/1471-2202-14-147.

67. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003; 34:374–78.