Skip to content

Advertisement

  • Research article
  • Open Access
  • Open Peer Review

Comprehensive analyses of imprinted differentially methylated regions reveal epigenetic and genetic characteristics in hepatoblastoma

  • 1, 2,
  • 1,
  • 3,
  • 4,
  • 1,
  • 5,
  • 1,
  • 1,
  • 4,
  • 6,
  • 7,
  • 7,
  • 5,
  • 4,
  • 3,
  • 8,
  • 1Email author and
  • 1Email author
BMC Cancer201313:608

https://doi.org/10.1186/1471-2407-13-608

  • Received: 13 March 2013
  • Accepted: 20 December 2013
  • Published:
Open Peer Review reports

Abstract

Background

Aberrant methylation at imprinted differentially methylated regions (DMRs) in human 11p15.5 has been reported in many tumors including hepatoblastoma. However, the methylation status of imprinted DMRs in imprinted loci scattered through the human genome has not been analyzed yet in any tumors.

Methods

The methylation statuses of 33 imprinted DMRs were analyzed in 12 hepatoblastomas and adjacent normal liver tissue by MALDI-TOF MS and pyrosequencing. Uniparental disomy (UPD) and copy number abnormalities were investigated with DNA polymorphisms.

Results

Among 33 DMRs analyzed, 18 showed aberrant methylation in at least 1 tumor. There was large deviation in the incidence of aberrant methylation among the DMRs. KvDMR1 and IGF2-DMR0 were the most frequently hypomethylated DMRs. INPP5Fv2-DMR and RB1-DMR were hypermethylated with high frequencies. Hypomethylation was observed at certain DMRs not only in tumors but also in a small number of adjacent histologically normal liver tissue, whereas hypermethylation was observed only in tumor samples. The methylation levels of long interspersed nuclear element-1 (LINE-1) did not show large differences between tumor tissue and normal liver controls. Chromosomal abnormalities were also found in some tumors. 11p15.5 and 20q13.3 loci showed the frequent occurrence of both genetic and epigenetic alterations.

Conclusions

Our analyses revealed tumor-specific aberrant hypermethylation at some imprinted DMRs in 12 hepatoblastomas with additional suggestion for the possibility of hypomethylation prior to tumor development. Some loci showed both genetic and epigenetic alterations with high frequencies. These findings will aid in understanding the development of hepatoblastoma.

Keywords

  • Hepatoblastoma
  • Genomic imprinting
  • Differentially methylated region
  • DNA methylation

Background

Hepatoblastoma is the most common primary liver tumor in children, accounting for just over 1% of pediatric cancers and 79% of liver cancers in children under the age of 15 [1]. Most of these tumors are purely derived from epithelium composed exclusively of immature hepatocytic elements, known as fetal and embryonal types. The fetal type consists of smaller than normal hepatocytes that are arranged in irregular laminae, recapitulating those of the fetal liver. The embryonal type is comprised of smaller cells as compared to the fetal type. It has a more immature appearance and pattern of growth. Some of the tumors, referred to as mixed type, are characterized by epithelial patterns and spindled mesenchymal cells. A much rarer variant of such mixed type tumor harbors teratoid features, which contains foci of mature cartilage, intestinal-type or keratinized epithelium, melanin pigment, or skeletal muscle in addition to the elements mentioned above. To date, several genetic and epigenetic features have been observed in hepatoblastoma (reviewed in [2]). The most recurrent cytogenetic abnormalities include the presence of extra copies of chromosomes 2, 8, 20, and the loss of chromosome 4. Mutations or upregulation of the genes involved in embryonic development have been reported. For example, APC, CTNNB1, AXIN1, and AXIN2 (key factors involved in the Wnt signaling pathway) are frequently mutated, suggesting that aberration of this pathway occurs as an early event during tumorigenesis. Mutation of PIK3CA, amplification of PIK3C2B, and upregulation of hedgehog ligands and their target genes have also been reported. Epigenetic silencing by promoter hypermethylation occurs at several tumor suppressor genes, such as SFRP1, APC, HHIP, SOCS1, CASP8, and RASSF1A. In addition, several imprinted genes, including IGF2, DLK1, PEG3, PEG10, MEG3, and NDN, have been reported to be overexpressed in hepatoblastoma [2].

Imprinted genes are expressed in a parent-of-origin-specific manner. They are usually clustered in subchromosomal regions called imprinting domains. The human genome contains more than 30 imprinting domains (http://www.geneimprint.com). Imprinting domains have at least one DMR that are characterized by DNA methylation on one of the two parental alleles. There are maternally methylated DMRs and paternally methylated DMRs. In addition, two classes of imprinted DMRs, gametic and somatic, have been described. Gametic DMRs acquire methylation during gametogenesis and the methylation is maintained from zygote to somatic cells during all the developmental stages. Most gametic DMRs are known as imprinting control regions (ICRs) that regulate the imprinted expression of the genes in the domain. By contrast, methylations of somatic DMRs are established during early embryogenesis after fertilization under the control of nearby ICRs [3]. Somatic DMRs also regulate the expression of the imprinted genes.

Many imprinted genes regulate cell growth and differentiation, and, thus, disruption of imprinting, mainly due to aberrant DNA methylation at the responsible DMR, is implicated in pre- and/or post-natal growth disorders and in the pathogenesis of cancers [4]. For example, hypermethylation of H19-DMR, which is the ICR of the IGF2/H19 imprinting domain at the 11p15.5 locus, is a cause of Beckwith-Wiedemann syndrome (BWS), the most common overgrowth syndrome characterized by occasional development of embryonal tumors, including hepatoblastoma [5]. The hypermethylation leading to biallelic expression of IGF2 is seen in a range of tumors, also including hepatoblastoma [6, 7]. The LOH of 11p15.5, especially the loss of the maternal allele, is found in approximately 20% of hepatoblastoma cases, and it is reported to be a risk factor for the relapse of this tumor [7, 8]. Furthermore, several imprinted genes are overexpressed in hepatoblastoma as mentioned above. Thus, it is speculated that aberrant DNA methylation at imprinted DMRs is a key mechanism during malignant transformation of progenitor cells in a variety of tissues, including the liver [2, 9]. However, the methylation status of imprinted DMRs scattered through the human genome has yet to be analyzed comprehensively in hepatoblastoma.

In this study, we performed comprehensive methylation analyses and polymorphism analyses of 33 imprinted DMRs in hepatoblastoma. We therefore describe some epigenetic and genetic characteristics of hepatoblastoma. These findings collectively aid in the understanding of the development of hepatoblastoma.

Methods

Samples

Twelve hepatoblastomas and their paired adjacent normal liver tissues were analyzed. Eleven sporadic hepatoblastoma samples (HB01 - HB11) were obtained from the Department of Pediatric Surgery, Faculty of Medicine, Kyushu University, Japan. One hepatoblastoma developed in a BWS patient (BWS109) was obtained from Toho University, Omori Medical Centre, Japan. Histochemical analyses of the tumor tissues indicated that the average of the tumor cell contents was approximately 70%. Ten of the patients were treated based on the Japanese Study Group for Pediatric Liver Tumor-2 (JPLT-2) protocol (HB08 and HB09 were not). Clinical information of the hepatoblastoma cases is shown in Table 1. Three livers (CL7, CL16, CBD1) were used as normal controls. CL7 (a 7-year-old who died from spinal muscular atrophy type I-C with chronic respiratory insufficiency) and CL16 (a 16-year-old who died after head trauma) were provided by the non-profit organization, Human & Animal Bridging Research Organization (Chiba, Japan). CBD1 (a 7-month-old who had congenital biliary dilatation) was obtained from the Department of Pediatric Surgery, Faculty of Medicine, Kyushu University. Written informed consents were obtained from the parents or the guardians of the participants, because the participants were children or dead. This study was approved by the Ethical Committee for Human Genome and Gene Analyses of the Faculty of Medicine, Saga University.
Table 1

Clinical information of hepatoblastoma cases

Case

Sex/agea

Histology

PRETEXT

Preoperative chemotherapyb

POSTTEXT

Outcome

Other information

HB01

F/1y3m

Combined fetal and embryonal type

III

CITA4

III

Alive

 

HB02

F/3y2m

Fetal typec

III

CITA4

III

Aive

 

HB03

F/7y11m

Hepatoblastoma (NOS)d

III

CITA5

III

Alive

Small for gestational age

HB04

M/1y4m

Mixed epithelial and mesenchymal with teratoid featurec

IV

CITA4 + ITEC2

IV

Alive

 

HB05

M/1y2m

Mixed epithelial and mesenchymal with teratoid feature

III

CITA5

II

Alive

 

HB06

M/10m

Mixed epithelial and mesenchymal with teratoid feature

III

CITA4

III

Alive

 

HB07

M/8m

Combined fetal and embryonal type

II

CITA2

II

Alive

 

HB08

F/28d

Combined fetal and embryonal type

II

  

Alive

 

HB09

M/1y6m

Combined fetal and embryonal type

II

  

Treatment related death

Small for gestational age

HB10

F/6y6m

Fetal type

II

CITA2

II

Alive

 

HB11

F/3m

Combined fetal and embryonal type

IV

CITA7

III

Treatment related death

 

BWS109

F/1y0m

Hepatoblastoma (NOS)d

IV,M(+)

CITA7 + ITEC1

IV

Alive

Beckwith-Wiedemann syndrome, liver transplantation at 1 year old

aage at diagnosis, bCITA: cisplatin-pirarubicin, ITEC: Ifosfamide, pirarubicin, etoposide, and carboplatin. The numerals indicate the cycle numbers of the chemotherapy. cdifficult to diagnose due to chemotherapy, dnot otherwise specified.

DNA isolation and bisulphite conversion

Genomic DNA was extracted from each sample using the QIAamp DNA Mini Kit (Qiagen, Germany) according to the manufacturer’s instructions. One microgram of genomic DNA was subjected to bisulfite conversion using the EZ DNA Methylation KitTM (Zymo Research, CA), and then the converted DNA was eluted in 100 μl of water.

MALDI-TOF MS analysis

The methylation status of imprinted DMRs was screened by MALDI-TOF MS analysis with a MassARRAY system (Sequenom, CA) [10], according to the manufacturer’s instructions. MALDI-TOF MS analysis produced signal pattern pairs indicative of non-methylated and methylated DNA. Epityper software analysis of the signals yielded the methylation index which ranged from 0 (no methylation) to 1 (full methylation) of each CpG unit, which contained one or more CpG sites measured as one unit in the MALDI-TOF MS analysis. Aberrant methylation of a CpG unit was defined as when the difference of methylation indexes between two samples exceeded 0.15, which was based on the fact that we have previously found that the differences of H19-DMR hypermethylation or KvDMR1 hypomethylation in BWS patients were at least more than 0.15 (data not shown). Since analyzed DMRs included several CpG units, aberrant methylation of a DMR was defined as when more than 60% of total number of analyzed CpG units showed aberrant methylation (with the difference exceeding 0.15). We used CL7 and CBD1 as normal controls in MALDI-TOF MS analysis.

Pyrosequencing

Pyrosequencing was conducted using QIAGEN PyroMark Q24 according to the manufacturer’s instruction (Qiagen, Germany). Some of the primers for DMR analysis were described by Woodfine et al. [11]. We designed other primers by using PyroMark Assay Design 2.0 (Qiagen, Germany). The primers for LINE-1 (GenBank accession no. X58075) analyses were described by Bollati et al. [12]. The criterion for MALDI-TOF MS analysis was also employed to define the aberrant methylation of each CpG site and an analyzed region. We used three livers, i.e. CL7, CL16, and CBD1, as normal controls in pyrosequencing. The control livers were analyzed in triplicate for LINE-1 and once for DMRs.

DNA Polymorphism analysis

LOH, UPD, and copy number abnormalities were investigated with DNA polymorphisms. For quantitative analyses, tetranucleotide repeat markers near the imprinted DMRs were amplified and separated by electrophoresis on an Applied Biosystems 3130 genetic analyzer. Data were then quantitatively analyzed with GeneMapper software (Applied Biosystems, CA). The peak height ratios of two parental alleles were calculated. A single nucleotide polymorphism (SNP) of KCNQ1DN (rs229897) was also analyzed.

All primers used in this study are shown in Additional file 1: Table S1.

Statistical analysis

The methylation statuses of the samples were compared in three pairs: adjacent normal liver tissue (A) and control livers (C), denoted as AxC; tumors (T) and control livers, denoted as TxC; tumors and adjacent normal liver tissue, denoted as TxA. The binomial distribution test was performed to compare aberrant hypomethylation and aberrant hypermethylation within each comparison pair (AxC, TxC, and TxA). The Chi squared test or Fisher’s exact test was used for comparison between maternally methylated and paternally methylated DMRs and between gametic and somatic DMRs in aberrant hypomethylation and aberrant hypermethylation cases for each comparison pair. For LINE-1 methylation, a paired t-test was used to compare tumor and adjacent normal liver tissue, and an independent t-test (Welch’s t-test) was used for comparing tumor or adjacent normal liver tissue with control liver. A p-value less than 0.05 was considered to be statistically significant. Bonferroni correction was performed when needed.

Results

Clinical information of hepatoblastoma cases

Clinical information of the 12 hepatoblastoma cases analyzed in this study are shown in Table 1. Eleven tumors were sporadic (HB01-HB11), and one was associated with BWS (BWS109). The ratio of males to females was 5:7. The mean age at diagnosis was 25.7 months, ranging from 28 days to 7 years and 11 months. In terms of histological features, 5 cases had combined fetal and embryonal types, 3 cases had mixed epithelial and mesenchymal features with teratoid features, 2 cases were fetal types, and 2 cases were hepatoblastomas (not otherwise specified). Using PRETEXT staging [13], 4 cases were stage II, 5 cases were stage III, and 3 cases were stage IV. Ten of twelve cases were undergoing chemotherapy based on the JPLT-2 protocol, but only two cases (HB05 and HB11) regressed to a lower stage after chemotherapy.

Analyses of aberrant methylation and genetic alterations at imprinted DMRs

We selected 33 regions reported previously as imprinted DMRs in the human genome [11, 14] (refer to http://www.geneimprint.com/). Our strategy in this study involved screening the methylation levels of DMRs in tumors, their paired adjacent normal tissues, and normal control livers by MALDI-TOF MS. The samples that showed aberrant methylation were analyzed again by pyrosequencing to confirm the result. These two methods are the most reliable methods of methylation analysis at present [10, 15, 16]. First, we analyzed the methylation level of these regions in two normal livers, i.e. CL7 and CBD1, by MALDI-TOF MS (Additional file 2: Figure S1). A total of 20 DMRs showed almost 50% methylation, however, 8 DMRs (IGF2R-DMR2, IGF2-DMR0, IGF2-DMR2, IG-DMR-CG4, IG-DMR-CG6, TCEB3C-DMR, USP29-DMR, and NNAT-DMR) showed mostly full methylation, and 5 DMRs (TP73-DMR, SPTBN1-DMR, WT1-AS-DMR, DLK1-DMR, and GNASXL-DMR) showed mostly no methylation. It is highly possible that these regions were not differentially methylated in the normal liver, probably due to tissue-specific and/or age-related features of differential methylation, because most of the regions were also analyzed by pyrosequencing and their methylation statuses were confirmed (Additional file 2: Figure S1).

Next, we screened the methylation status of the 33 DMRs in 12 hepatoblastomas and their paired adjacent normal liver tissue by MALDI-TOF MS. We found aberrant methylation in tumors and also in adjacent liver tissue by comparing the methylation between tumors and normal controls (TxC), tumors and adjacent liver tissue (TxA), and adjacent liver tissue and normal controls (AxC). The definition of aberrant methylation is described in the Methods section. After excluding samples harboring chromosomal abnormalities as described later, we confirmed the aberrant methylation using pyrosequencing, except in the case of H19-DMR (representative data is shown in Figure 1 and all data in Additional file 3: Figure S2). Additional normal control liver, CL16, was used in pyrosequencing analyses. The methylation status of H19-DMR was analyzed by hot-stop combined bisulfite restriction analysis (COBRA) [17] or bisulfite sequencing because of the difficulty in the primer-design for pyrosequencing (Additional file 4: Figure S3).
Figure 1
Figure 1

Representative results of methylation analyses by MALDI-TOF MS and pyrosequencing. (A and B) Representative samples of aberrant hypomethylation. KvDMR1 of HB08 and PEG10 of HB11 are shown. KvDMR1 was hypomethylated in both adjacent normal liver and tumor tissues, whereas PEG10-DMR was hypomethylated only in tumor tissue. (C and D) Representative samples of aberrant hypermethylation. INPP5Fv2-DMR in HB05 and RB1-DMR in HB11 are shown. Only tumors showed hypermethylation at these DMRs. Aberrant methylation of a DMR was defined as when more than 60% of total CpG units or CpG sites were aberrantly methylated. Aberrant methylation of a CpG unit or CpG site was defined as occurring when the difference of its methylation indexes in two samples exceeded 0.15. The vertical axis represents the methylation index; the horizontal axis represents CpG units (MALDI-TOF MS) or CpG sites (pyrosequencing) analyzed. Green line: average of normal control livers; blue line: adjacent normal liver; red line: tumor (hepatoblastoma).

In order to exclude aberrantly methylated DMRs, as associated with chromosome abnormalities, such as UPD or copy number abnormality, DNA polymorphism analyses using microsatellites and SNPs were performed on all regions showing aberrant methylation in the MALDI-TOF MS analyses. We found seven genetic alterations in four tumors resulting in aberrant methylation: abnormal allelic copy number of 11p13-p15.5 in HB01, 20q11-q13 in HB05, and 19q13 and 20q13 in HB11; LOH of 7q32.2 and 11p15.5 in HB11; and paternal UPD of 11p13-p15.5 in BWS109 (Figure 2). We speculated the allelic imbalance statuses of these loci according to the results of MALDI-TOF MS analysis and DNA polymorphism analysis (Additional file 5: Figure S4). HB01 would harbor more paternal copies than the maternal copies in 11p13-p15.5. An abnormal allelic copy number of 20q11-q13 in HB05 would represent a higher copy number in the maternal allele than the paternal allele. HB11 would have more maternal copies of 19q13 and more paternal copies of 20q13. LOH of 7q32.2 and 11p15.5 would have occurred due to a paternal deletion and a maternal deletion, respectively. The paternal UPD was confirmed using the parents’ DNA (Additional file 5: Figure S4). The extent of the UPD was found to be at 11p11.2–pter by a SNP array analysis (data not shown).
Figure 2
Figure 2

Aberrant methylations and genetic alterations of 33 imprinted DMRs in 12 hepatoblastomas. Aberrant hypomethylation and aberrant hypermethylation found in each comparison are indicated by blue and red boxes, respectively. Aberrant methylation was identified by comparing adjacent normal liver tissue with normal control livers: AxC; tumors compared with normal control livers: TxC; and tumors compared with adjacent normal livers: TxA. The classes of these DMRs in previous reports are shown in Pat/Mat and Gametic/Somatic rows. ZDBF2-DMR was tentatively assigned as a somatic DMR based on experiments from mice [35]. MALDI-TOF MS revealed that 13 DMRs, such as TP73-DMR, SPTBN1-DMR, IGF2R-DMR2, IGF2-DMR0, IGF2-DMR2, WT1-AS-DMR, DLK1-DMR, IG-DMR-CG4, IG-DMR-CG6, TCEB3C-DMR, USP29-DMR, NNAT-DMR, and GNASXL-DMR did not show differential methylation in control livers. Most of these methylation statuses were confirmed by pyrosequencing analysis. Mat: maternally methylated DMR; Pat: paternally methylated DMR; gametic: gametic DMR; somatic: somatic DMR. ACN: abnormal copy number; LOH: loss of heterozygosity; UPD: paternal uniparental disomy; DD: difficult to decide.

The following results described are shown in Figure 2. All tumors carried aberrant methylation in at least one DMR. HB05 carried aberrant methylations at 8 DMRs, the highest number of aberrant methylations, whereas HB03 and HB08 carried aberrant methylations at only one DMR, MCTS2-DMR and KvDMR1, respectively. A total of 18 of 33 DMRs showed aberrant methylation, whereas 15 DMRs did not show such features in any tumors. There was large deviation in the incidence of aberrant methylation among the DMRs. KvDMR1 and IGF2-DMR0 were the most frequently hypomethylated DMRs in 3 of 9 tumors. The most frequently hypermethylated DMR was INPP5Fv2-DMR, which occurred in 7 of 12 samples. RB1-DMR was also hypermethylated with a high frequency, which occurred in 5 of 12 samples. In addition, GNASXL-DMR was hypermethylated in 4 of 10 samples. The following DMRs showed aberrant methylation in only one tumor: ARH1-CG1, NAP1L5-DMR, PEG10-DMR, H19 promoter, WT1-AS-DMR, MEG3-CG7-DMR, MCTS2-DMR, NESP-DMR, and GNAS1A-DMR.

Two chromosomal loci, 11p15.5 and 20q13.32, showed high frequencies of genetic and epigenetic alterations at 10 of 12 and 7 of 12, respectively. In the 11p15.5 locus, seven tumors carried the aberrant methylation and three samples carried genetic alterations. In the 20q13.32 locus, five tumors carried aberrant methylation and two carried an abnormal allelic copy number.

Comparisons of aberrantly methylated DMRs

We compared the numbers of aberrantly hypomethylated and hypermethylated DMRs in three pairs of the sample groups (Figure 3A). We excluded the DMRs harboring UPD or copy number abnormalities for the statistical analyses. Comparing adjacent normal liver tissues (A) and control livers (C), herein denoted as AxC, only hypomethylation was observed in adjacent normal liver tissue (p = 0.031). In the TxC comparison, both hypermethylation and hypomethylation were observed in tumors (no significant difference). In the TxA comparison, hypermethylation was observed more frequently than hypomethylation with statistical significance (p = 0.013). In addition, the number of hypomethylated DMRs in tumors was higher than that of adjacent normal liver tissue (p = 0.040), although Bonferroni correction did not show statistical significance of the difference with a significance level of 0.05/3 (approximately 0.0167). These results suggested a possibility that hypomethylation occurred at certain DMRs in adjacent normal liver tissue prior to tumor development, whereas hypermethylation occurred only within the tumor tissue itself.
Figure 3
Figure 3

Comparison of aberrant methylations of DMRs in hepatoblastoma. (A) Comparison of the numbers of aberrantly hypomethylated (blue bars) and hypermethylated (red bars) DMRs. Hypomethylation, but not hypermethlation, was observed in adjacent normal tissues. Hypermethylation was observed more frequently than hypomethylation in tumors compared to adjacent normal tissues. (B) Comparison of the frequencies of aberrant methylations between the paternally methylated DMRs (black bars) and the maternally methylated DMRs (white bars). Left and right panels indicate hypomethylation and hypermethylation, respectively. There was no significant difference in both hypomethylation and hypermethylation. (C) Comparison of the frequencies of aberrant methylations between gametic DMRs (yellow bars) and somatic DMRs (green bars). Left and right panels indicate hypomethylation and hypermethylation, respectively. Hypermethylation was observed more frequently at gametic DMRs with statistically significant difference in tumors compared to control livers. The frequencies were compared in (B) and (C) because of the difference in the total numbers of DMRs in the two categories compared; that is, the paternal and the maternal DMRs or the gametic and the somatic DMRs. The asterisks indicated significant difference. The p values are shown above the comparisons. b: binomial test, c: chi squared test.

We further compared the frequencies of aberrant methylation between paternally methylated DMRs and maternally methylated DMRs. As for hypomethylation, there was no significant difference between the two kinds of DMRs in each comparison (Figure 3B, left panel). In contrast, hypermethylation, which was observed only in tumors, tended to occur more frequently at maternally methylated DMRs than paternally methylated DMRs in the TxC comparison (p = 0.060) (Figure 3B, right panel). We also compared the frequencies of aberrant methylation between gametic DMRs and somatic DMRs. No significant difference in hypomethylation was found between the two kinds of DMRs in each comparison (Figure 3C, left panel). In contrast to hypomethylation, hypermethylation occurred at gametic DMRs more frequently with statistical significance (p = 0.043) (Figure 3C, right panel). This difference was mainly due to the frequent occurrence of hypermethylation at three maternally methylated and gametic DMRs, such as INPP5Fv2-DMR, RB1-DMR, and GNASXL-DMR (Figure 2).

Methylation status of LINE-1 in hepatoblastoma

We analyzed the methylation status of LINE-1 in all samples by pyrosequencing to assess the genome-wide methylation level. LINE-1 is a human repetitive element and constitutes approximately 30% of the human genome [18]. Its methylation status has been used as a surrogate marker for genome-wide DNA hypomethylation in many cancers [19, 20]. We analyzed the methylation of four CpG sites in a LINE-1 sequence that were hypomethylated in cancers [21, 22]. We compared the methylation levels of each CpG site among the three groups using Bonferroni correction with significance level of 0.0167 (Figure 4). Tumors showed slight hypomethylation only at CpG1 among four CpGs (p = 0.015 in TxA). However, other CpG sites did not show hypomethylation although bare hypermethylations (less than 2.5%) was found only in adjacent normal liver tissues at CpG2 (p = 0.001 in AxC and p = 0.010 in AxT). These results suggested that the genome-wide methylation levels were almost same among three sample groups.
Figure 4
Figure 4

Methylation status of LINE-1. Methylation levels of four successive CpG sites in LINE-1 were analyzed in hepatoblastomas (T), paired adjacent normal liver tissues (A), and three normal control livers (C) by pyrosequencing. The normal controls were analyzed in triplicate and the average values were plotted with standard deviations. The average values of methylation (%) are shown for each sample groups. The asterisks indicated significant differences (p < 0.0167). Bonfferoni correction was applied for this statistical analysis.

Discussion

In this study, we found many imprinted DMRs methylated aberrantly in hepatoblastomas and paired adjacent normal liver tissue. An important finding was that the aberrant hypomethylation occurred not only in tumor tissue but also in adjacent normal liver tissue. One possible explanation is that the occurrence of the aberrant hypomethylation at certain DMRs may be a very early and specific event prior to tumor development, although there is another possibility that the tumor may induce methylation changes in the adjacent tissues. Okamoto et al. have previously reported a similar phenomenon with respect to aberrant hypermethylation of H19-DMR that was frequently found in normal tissues adjacent to Wilms’ tumors, which carried the same aberrant methylation [23]. Based on the results, it was hypothesized that the preceding aberrant methylation may be a constitutional aberration in the onset of embryonal tumors. In contrast to the hypomethylation, the aberrant hypermethylation of the DMR occurred only in tumors. These results indicated that the hypermethylation of the DMRs, especially for INPP5Fv2-DMR, RB1-DMR, and GNASXL-DMR, was a specific event for tumor development; this suggested that the pre-cancerous cells did not carry hypermethylation at the DMRs, but acquired the aberrant methylation during tumor development.

We also analyzed the genome-wide methylation level, represented by LINE-1 methylation, and we did not find large difference among three sample groups as a whole. LINE-1 is usually hypomethylated in many adult tumors, and its methylation level correlates with clinicopathological features of the tumors [19]. The different situation concerning LINE-1 methylation between hepatoblastoma and adult tumors may reflect a different mechanism of tumorigenesis in embryonal tumors as compared to adult tumors.

Hypermethylation in tumors was frequently observed at three DMRs, INPPF5v2-DMR, RB1-DMR, and GNASXL-DMR. INPPF5v2-DMR controls the expression of INPP5F transcript variant 2, which encodes a protein of an unknown function [24, 25]. RB1-DMR, located in intron 2 of the RB1 gene, leads to maternal expression of transcript variants from exon 2B with very low expression in normal tissues [26]. The function of the variants in cell proliferation is not known. Thus, the effect of these hypermethylated DMRs on tumorigenesis would be little or unknown. GNASXL-DMR is associated with the paternal expression of GNASXL, which encodes a protein involved in signal transduction [2729]. The DMR was shown to be mostly unmethylated in control livers (Additional file 2: Figure S1). Thus, hypermethylation of GNASXL-DMR would reduce expression of GNASXL. Unfortunately, the expression of genes linked to aberrantly methylated DMRs could not be analyzed due to poor RNA quality, which was probably due to effects of chemotherapy and a limited amount of samples. Therefore, we could not assess the involvement of hypermethylation in tumorigenesis of hepatoblastoma.

Another important finding was the frequent occurrence of both genetic and epigenetic alterations at the two chromosomal loci, 11p15.5 and 20q13.3. The 11p15.5 locus is a well-known imprinted locus responsible for BWS, a tumor-predisposing imprinting disorder. The locus was found to be altered genetically and/or epigenetically in 10 of 12 tumors. Hypermethylation at H19-DMR and hypomethylation at IGF2-DMR0 associated with biallelic expression of IGF2 were reported in adult and embryological tumors, including hepatoblastoma [6, 7]. Hypermethylation at H19-DMR and the H19 promoter also reduced the expression of H19 in Wilms’ tumor [30, 31]. In addition to epigenetic alterations, genetic alterations, such as the amplification of paternal alleles leading to overexpression of IGF2 and LOH of the maternal allele leading to reduced expression of H19, were observed in sporadic Wilms’ tumors [32, 33]. In this study, in addition to the hypermethylations at H19-DMR and the H19 promoter in two tumors, hypomethylation at IGF2-DMR0 occurred in another two adjacent normal liver tissues. Further, abnormal allelic copy number, paternal UPD, and maternal LOH of 11p15.5 were observed. The overexpression of IGF2 and the reduced expression of H19 would play an important role in tumorigenesis of hepatoblastoma.

The 20q13.3 locus was also altered genetically and/or epigenetically in 7 of 12 tumors. This locus is responsible for pseudohypoparathyroidism, a condition in which pathogenesis is attributed to the tissue specific imprinting of Gsα, for example, which occurs in the proximal renal tubule. On the other hand, an extra copy of chromosome 20 has been known to be the most recurrent cytogenetic alteration in hepatoblastoma [2, 34]. We found copy number differences of the alleles in this region, suggesting the existence of non-imprinted oncogenic gene(s) in this region.

Many epigenetic and genetic alterations were found at the loci linked to the 33 imprinted DMRs in 12 hepatoblastomas. However, since sample numbers in this study were small, more hepatoblastoma samples should be analyzed to confirm the present data and to evaluate the precise role of these alterations in tumorigenesis in addition to assessing their usefulness as markers for clinical characteristics, such as stage classification, response to chemotherapy, and prognosis. Also needed are the expression analyses of the genes linked to aberrantly methylated DMRs to assess their role in tumor development, although it is very difficult to obtain hepatoblastoma samples without any chemotherapeutic history.

Conclusions

We found epigenetic and genetic characteristics of hepatoblastoma by comprehensive epigenetic and genetic analyses of 33 DMRs linked to imprinting loci in 12 hepatoblastoma samples and their adjacent normal liver tissues. These included aberrant hypomethylation in adjacent normal liver tissue, tumor-specific hypermethylation, and the frequent occurrence of both genetic and epigenetic alterations at 11p15.5 and 20q13.3 loci. Further studies using more hepatoblastoma samples are needed to confirm the present results and evaluate their roles in the tumor development.

Abbreviations

DMR: 

Differentially methylated region

LOH: 

Loss of heterozygosity

UPD: 

Uniparental disomy

LINE-1: 

Long interspersed nuclear element-1

ICR: 

Imprinting control regions

BWS: 

Beckwith-Wiedemann syndrome

MALDI-TOF MS: 

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry

SNP: 

Single nucleotide polymorphism

COBRA: 

Combined bisulfite restriction analysis.

Declarations

Acknowledgements

We thank Prof. Yutaka Kondo, Division of Epigenomics, Aichi Cancer Center Research Institute, Japan for technical advices in LINE-1 methylation analysis. This study was supported, in part, by a Grant for Research on Intractable Diseases from the Ministry of Health, Labor, and Welfare; a Grant for Child Health and Development from the National Center for Child Health and Development; and, a Grant-in-Aid for Challenging Exploratory Research and a Grant-in-Aid for Young Scientists (B) from the Japan Society for the Promotion of Science.

Authors’ Affiliations

(1)
Department of Biomolecular Sciences, Division of Molecular Genetics & Epigenetics, Faculty of Medicine, Saga University, Nabeshima 5-1-1, Saga 849-8501, Japan
(2)
Faculty of Medicine, Sam Ratulangi University, Manado, Indonesia
(3)
Department of Pediatric Surgery, Faculty of Medicine, Kyushu University, Fukuoka, Japan
(4)
Department of Pediatrics, Toho University, Omori Medical Center, Tokyo, Japan
(5)
Department of Maternal-Fetal Biology, National Research Institute for Child Health and Development, Tokyo, Japan
(6)
Department of Pathology and Microbiology, Division of Pathology, Faculty of Medicine, Saga University, Saga, Japan
(7)
Department of Anatomic Pathology, Pathological Sciences, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan
(8)
Department of Pediatric Surgery, Graduate School of Medical Science, Kyoto Prefectural University of Medicine, Kyoto, Japan

References

  1. Herzog CE, Andrassy RJ, Eftekhari F: Childhood cancers: hepatoblastoma. Oncologist. 2000, 5 (6): 445-453. 10.1634/theoncologist.5-6-445.View ArticlePubMedGoogle Scholar
  2. Tomlinson GE, Kappler R: Genetics and epigenetics of hepatoblastoma. Pediatr Blood Cancer. 2012, 59 (5): 785-792. 10.1002/pbc.24213.View ArticlePubMedGoogle Scholar
  3. Tomizawa S, Sasaki H: Genomic imprinting and its relevance to congenital disease, infertility, molar pregnancy and induced pluripotent stem cell. J Hum Genet. 2012, 57 (2): 84-91. 10.1038/jhg.2011.151.View ArticlePubMedGoogle Scholar
  4. Murrell A: Genomic imprinting and cancer: from primordial germ cells to somatic cells. ScientificWorldJournal. 2006, 6: 1888-1910.View ArticlePubMedGoogle Scholar
  5. Choufani S, Shuman C, Weksberg R: Beckwith-Wiedemann syndrome. Am J Med Genet C: Semin Med Genet. 2010, 154C (3): 343-354. 10.1002/ajmg.c.30267.View ArticleGoogle Scholar
  6. Scelfo RA, Schwienbacher C, Veronese A, Gramantieri L, Bolondi L, Querzoli P, Nenci I, Calin GA, Angioni A, Barbanti-Brodano G, et al: Loss of methylation at chromosome 11p15.5 is common in human adult tumors. Oncogene. 2002, 21 (16): 2564-2572. 10.1038/sj.onc.1205336.View ArticlePubMedGoogle Scholar
  7. Honda S, Arai Y, Haruta M, Sasaki F, Ohira M, Yamaoka H, Horie H, Nakagawara A, Hiyama E, Todo S, et al: Loss of imprinting of IGF2 correlates with hypermethylation of the H19 differentially methylated region in hepatoblastoma. Br J Cancer. 2008, 99 (11): 1891-1899. 10.1038/sj.bjc.6604754.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Chitragar S, Iyer VK, Agarwala S, Gupta SD, Sharma A, Wari MN: Loss of heterozygosity on chromosome 11p15.5 and relapse in hepatoblastomas. Eur J Pediatr Surg. 2011, 21 (1): 50-53. 10.1055/s-0030-1267208.View ArticlePubMedGoogle Scholar
  9. Holm TM, Jackson-Grusby L, Brambrink T, Yamada Y, Rideout WM, Jaenisch R: Global loss of imprinting leads to widespread tumorigenesis in adult mice. Cancer Cell. 2005, 8 (4): 275-285. 10.1016/j.ccr.2005.09.007.View ArticlePubMedGoogle Scholar
  10. Ehrich M, Nelson MR, Stanssens P, Zabeau M, Liloglou T, Xinarianos G, Cantor CR, Field JK, van den Boom D: Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc Natl Acad Sci U S A. 2005, 102 (44): 15785-15790. 10.1073/pnas.0507816102.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Woodfine K, Huddleston JE, Murrell A: Quantitative analysis of DNA methylation at all human imprinted regions reveals preservation of epigenetic stability in adult somatic tissue. Epigenetics Chromatin. 2011, 4 (1): 1-10.1186/1756-8935-4-1.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Bollati V, Baccarelli A, Hou L, Bonzini M, Fustinoni S, Cavallo D, Byun HM, Jiang J, Marinelli B, Pesatori AC, et al: Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res. 2007, 67 (3): 876-880. 10.1158/0008-5472.CAN-06-2995.View ArticlePubMedGoogle Scholar
  13. Roebuck DJ, Aronson D, Clapuyt P, Czauderna P, dJe Ville de Goyet J, Gauthier F, Mackinlay G, Maibach R, McHugh K, Olsen OE, et al: 2005 PRETEXT: a revised staging system for primary malignant liver tumours of childhood developed by the SIOPEL group. Pediatr Radiol. 2007, 37 (2): 123-132. 10.1007/s00247-006-0361-5. quiz 249–150View ArticlePubMedGoogle Scholar
  14. Higashimoto K, Nakabayashi K, Yatsuki H, Yoshinaga H, Jozaki K, Okada J, Watanabe Y, Aoki A, Shiozaki A, Saito S, et al: Aberrant methylation of H19-DMR acquired after implantation was dissimilar in soma versus placenta of patients with Beckwith-Wiedemann syndrome. Am J Med Genet A. 2012, 158A (7): 1670-1675. 10.1002/ajmg.a.35335.View ArticlePubMedGoogle Scholar
  15. Tost J, Dunker J, Gut IG: Analysis and quantification of multiple methylation variable positions in CpG islands by Pyrosequencing. Biotechniques. 2003, 35 (1): 152-156.PubMedGoogle Scholar
  16. Claus R, Wilop S, Hielscher T, Sonnet M, Dahl E, Galm O, Jost E, Plass C: A systematic comparison of quantitative high-resolution DNA methylation analysis and methylation-specific PCR. Epigenetics. 2012, 7 (7): 772-780. 10.4161/epi.20299.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Uejima H, Lee MP, Cui H, Feinberg AP: Hot-stop PCR: a simple and general assay for linear quantitation of allele ratios. Nat Genet. 2000, 25 (8): 375-376.PubMedGoogle Scholar
  18. Hancks DC, Kazazian HH: Active human retrotransposons: variation and disease. Curr Opin Genet Dev. 2012, 22 (3): 191-203. 10.1016/j.gde.2012.02.006.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Takai D, Yagi Y, Habib N, Sugimura T, Ushijima T: Hypomethylation of LINE1 retrotransposon in human hepatocellular carcinomas, but not in surrounding liver cirrhosis. Jpn J Clin Oncol. 2000, 30 (7): 306-309. 10.1093/jjco/hyd079.View ArticlePubMedGoogle Scholar
  20. Kitkumthorn N, Mutirangura A: Long interspersed nuclear element-1 hypomethylation in cancer: biology and clinical applications. Clin Epigenetics. 2011, 2 (2): 315-330. 10.1007/s13148-011-0032-8.View ArticlePubMedPubMed CentralGoogle Scholar
  21. An B, Kondo Y, Okamoto Y, Shinjo K, Kanemitsu Y, Komori K, Hirai T, Sawaki A, Tajika M, Nakamura T, et al: Characteristic methylation profile in CpG island methylator phenotype-negative distal colorectal cancers. Int J Cancer. 2010, 127 (9): 2095-2105. 10.1002/ijc.25225.View ArticlePubMedGoogle Scholar
  22. Okamoto Y, Sawaki A, Ito S, Nishida T, Takahashi T, Toyota M, Suzuki H, Shinomura Y, Takeuchi I, Shinjo K, et al: Aberrant DNA methylation associated with aggressiveness of gastrointestinal stromal tumour. Gut. 2012, 61 (3): 392-401. 10.1136/gut.2011.241034.View ArticlePubMedGoogle Scholar
  23. Okamoto K, Morison IM, Taniguchi T, Reeve AE: Epigenetic changes at the insulin-like growth factor II/H19 locus in developing kidney is an early event in Wilms tumorigenesis. Proc Natl Acad Sci U S A. 1997, 94 (10): 5367-5371. 10.1073/pnas.94.10.5367.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Choi JD, Underkoffler LA, Wood AJ, Collins JN, Williams PT, Golden JA, Schuster EF, Loomes KM, Oakey RJ: A novel variant of Inpp5f is imprinted in brain, and its expression is correlated with differential methylation of an internal CpG island. Mol Cell Biol. 2005, 25 (13): 5514-5522. 10.1128/MCB.25.13.5514-5522.2005.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Monk D, Arnaud P, Frost JM, Wood AJ, Cowley M, Martin-Trujillo A, Guillaumet-Adkins A, Iglesias Platas I, Camprubi C, Bourc’his D, et al: Human imprinted retrogenes exhibit non-canonical imprint chromatin signatures and reside in non-imprinted host genes. Nucleic Acids Res. 2011, 39 (11): 4577-4586. 10.1093/nar/gkq1230.View ArticlePubMedPubMed CentralGoogle Scholar
  26. Kanber D, Berulava T, Ammerpohl O, Mitter D, Richter J, Siebert R, Horsthemke B, Lohmann D, Buiting K: The human retinoblastoma gene is imprinted. PLoS Genet. 2009, 5 (12): e1000790-10.1371/journal.pgen.1000790.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Hayward BE, Kamiya M, Strain L, Moran V, Campbell R, Hayashizaki Y, Bonthron DT: The human GNAS1 gene is imprinted and encodes distinct paternally and biallelically expressed G proteins. Proc Natl Acad Sci USA. 1998, 95 (17): 10038-10043. 10.1073/pnas.95.17.10038.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Pasolli HA, Klemke M, Kehlenbach RH, Wang Y, Huttner WB: Characterization of the extra-large G protein alpha-subunit XLalphas. I. Tissue distribution and subcellular localization. J Biol Chem. 2000, 275 (43): 33622-33632. 10.1074/jbc.M001335200.View ArticlePubMedGoogle Scholar
  29. Klemke M, Pasolli HA, Kehlenbach RH, Offermanns S, Schultz G, Huttner WB: Characterization of the extra-large G protein alpha-subunit XLalphas. II. Signal transduction properties. J Biol Chem. 2000, 5 (43): 33633-33640.View ArticleGoogle Scholar
  30. Moulton T, Crenshaw T, Hao Y, Moosikasuwan J, Lin N, Dembitzer F, Hensle T, Weiss L, McMorrow L, Loew T: Epigenetic lesions at the H19 locus in Wilms’ tumour patients. Nat Genet. 1994, 7 (3): 440-447. 10.1038/ng0794-440.View ArticlePubMedGoogle Scholar
  31. Steenman MJ, Rainier S, Dobry CJ, Grundy P, Horon IL, Feinberg AP: Loss of imprinting of IGF2 is linked to reduced expression and abnormal methylation of H19 in Wilms’ tumour. Nat Genet. 1994, 7 (3): 433-439. 10.1038/ng0794-433.View ArticlePubMedGoogle Scholar
  32. Satoh Y, Nakadate H, Nakagawachi T, Higashimoto K, Joh K, Masaki Z, Uozumi J, Kaneko Y, Mukai T, Soejima H: Genetic and epigenetic alterations on the short arm of chromosome 11 are involved in a majority of sporadic Wilms’ tumours. Br J Cancer. 2006, 95 (4): 541-547. 10.1038/sj.bjc.6603302.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Erson A, Petty E: Kidney: Nephroblastoma (Wilms tumor). Atlas Genet Cytogenet Oncol Haematol. 2007, 11 (1): 50-53.Google Scholar
  34. Steenman M: Liver: Hepatoblastoma. Atlas Genet Cytogenet Oncol Haematol. 2002, 6 (1): 50-52.Google Scholar
  35. Kobayashi H, Sakurai T, Sato S, Nakabayashi K, Hata K, Kono T: Imprinted DNA methylation reprogramming during early mouse embryogenesis at the Gpr1-Zdbf2 locus is linked to long cis-intergenic transcription. FEBS Lett. 2012, 586 (6): 827-833. 10.1016/j.febslet.2012.01.059.View ArticlePubMedGoogle Scholar
  36. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/608/prepub

Copyright

© Rumbajan et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement