Epigenetic basis of hepatocellular carcinoma: A network-based integrative meta-analysis

Abstract

AIM

To identify the key epigenetically modulated genes and pathways in HCC by performing an integrative meta-analysis of all major, well-annotated and publicly available methylation datasets using tools of network analysis.

METHODS

PubMed and Gene Expression Omnibus were searched for genome-wide DNA methylation datasets. Patient clinical and demographic characteristics were obtained. DNA methylation data were integrated using the Ingenuity Pathway Analysis, a software package for visualizing and analyzing biological networks. Pathway enrichment analysis was performed using IPA, which also provides literature-driven and computationally-predicted annotations for significant association of genes to curated molecular pathways.

RESULTS

From an initial 928 potential abstracts, we identified and analyzed 11 eligible high-throughput methylation datasets representing 354 patients. A significant proportion of studies did not provide concomitant clinical data. In the promoter region, HIST1H2AJ and SPDYA were the most commonly methylated, whereas HRNBP3 gene was the most commonly hypomethylated. ESR1 and ERK were central genes in the principal networks. The pathways most associated with the frequently methylated genes were G-protein coupled receptor and cAMP-mediated signalling.

CONCLUSION

Using an integrative network-based analysis approach of genome-wide DNA methylation data of both the promoter and body of genes, we identified G-protein coupled receptor signalling as the most highly associated with HCC. This encompasses a diverse range of cancer pathways, such as the PI3K/Akt/mTOR and Ras/Raf/MAPK pathways, and is therefore supportive of previous literature on gene expression in HCC. However, there are novel targetable genes such as HIST1H2AJ that are epigenetically modified, suggesting their potential as biomarkers and for therapeutic targeting of the HCC epigenome.

Key Words: Network analysis; Hepatocellular carcinoma; Methylation

Core tip: Hepatocellular carcinoma (HCC) is a high-fatality cancer with limited screening biomarkers and therapeutic options. It arises in the context of chronic liver disease, having accumulated epigenetic changes over time. The goal of this study was to perform an integrative network-based meta-analysis of all genome-wide DNA methylation data in HCC. Using bioinformatics tools, we identified the most important aberrantly methylated genes and associated pathways. G-protein receptor signaling was the most significantly associated with HCC based on differential methylation of involved genes, which is consistent with the implication of the Ras/Raf/MAPK and mTOR pathways. The identification of novel epigenetically modified genes such as HIST1H2AJ within known pathways suggests targeting of the epigenome as a potential therapeutic avenue for HCC.

INTRODUCTION

Hepatocellular carcinoma (HCC) arises in the context of chronic liver disease, where there is ongoing injury over decades. HCC incidence in North America has been increasing in recent years, in the setting of a higher prevalence of cirrhosis secondary to hepatitis C and fatty liver disease[1]. It is the fifth most common cancer worldwide, and five-year survival is the second worst worldwide among all cancers at 8.9%. HCC is often diagnosed at later stages, and there is an inability to tolerate chemotherapy in patients with cirrhosis[1]. Curative treatment with resection, radiofrequency ablation or transplantation is only possible in early stage disease[2]. When diagnosed at a later stage, the first-line chemotherapeutic agent is sorafenib, which extends survival only by 3 mo[2]. Several other trials of chemotherapeutic regimens have been developed based on studies of genomic and transcriptomic data, with no further improvement in overall survival[3]. There is therefore a dire need to better understand HCC pathogenesis, elucidate screening biomarkers in patients at risk, and develop more optimal therapeutic agents.

Epigenetic changes are of significant interest to this malignancy, given that it results from mutations accumulating over time with exposure to various insults such as viral hepatitis, alcohol or fatty liver. Epigenetic modifications are heritable states of gene expression without altering DNA sequences[4]. They encompass processes such as DNA methylation, histone modifications, non-coding RNAs and nucleosome positioning. These changes are passed along faithfully to daughter cells during cell division[4]. Among these, DNA methylation has been the most studied, regulating gene expression through a stable silencing mechanism[5]. Covalent modification by DNA methyltransferases of cytosine residues with methyl groups in CpG dinucleotides occurs preferentially at the 5’ end in promoter regions. Transcriptional gene silencing results from this through two mechanisms: steric hindrance of transcription factors being able to access their cognate binding sites on gene promoters[5], and direct binding of methyl CpG binding domain containing proteins to the methylated DNA causing transcriptional repression[6]. The recent advent of genome-wide methylation analysis has enabled an appreciation of methylation status in genes of interest to cancer: Hypermethylation of tumor suppressor genes, hypomethylation of oncogenes, and methylation of repetitive elements[7]. The extensive reprogramming of the epigenome in cancer has led to a growing interest in epigenetic therapy. Specifically in HCC, aberrant DNA methylation of tumor suppressor gene promoters has been documented[8]. These epigenetic changes have been closely correlated with disease stage and clinical outcome[9]. There has been significant variability in the reported frequency of hypermethylated loci in HCC[10-12]. CDKN2A is methylated in 30%-70% of HCCs[10,12,13], RASSF1A in up to 85%[11,12], GSTP1 in 50-90%[14] and MGMT in 40%[15]. DNA methylation loci have also been reported as significantly enriched in the signaling networks of cellular development, gene expression, cell death, and cancer[16].

Our study represents the first comprehensive network-based attempt to integrate all relevant, publicly available, high-throughput genome-wide DNA methylation data to better understand the epigenetic landscape in HCC. Network and pathway analysis tools can enable identification of the most commonly methylated genes across studies and associated pathways, and propose novel treatment options using network-based analysis[17,18].

The goal of this study was to identify key epigenetically modulated genes and pathways in HCC by integrating all major, well-annotated and publicly available methylation datasets datasets using tools of network analysis.

MATERIALS AND METHODS

Data collection, analysis and database compiling

Genome-wide methylation profiles related to HCC samples were downloaded from published datasets (PubMed, http://www.ncbi.nlm.nih.gov/PubMed) using the following MeSH terms: “{[“methylation”(MeSH Terms) or “methylation” (all fields)] and [“carcinoma, hepatocellular” (MeSH terms) or [“carcinoma” (all fields) and “hepatocellular” (all fields)] or “hepatocellular carcinoma” (All Fields) or ["hepatocellular" (all fields) and "carcinoma" (all fields)]} and ["humans" (MeSH Terms) and English (lang)]. All entries on PubMed since 2002, which represents the advent of high-throughput profiling, were considered for inclusion. A second search was performed using Gene Expression Omnibus (GEO), a public functional genomics data repository containing genome-wide methylation profile array data (https://http://www.ncbi.nlm.nih.gov/geo). This search was performed using the following MeSH terms {[“methylation” (MeSH Terms) or methylation (all fields)] and [“carcinoma, hepatocellular” (MeSH terms) or hepatocellular carcinoma (all fields) and “Homo sapiens” (porgn) and “Homo sapiens” (porgn) and “Homo sapiens” (porgn)] and “Homo sapiens” (porgn)}, covering all HCC high-throughput methylation profiling datasets comparing HCC to adjacent non-tumoral tissue.

The study workflow is illustrated in Figure 1A. Results were retrieved from both databases: GEO and PubMed. The exclusion criteria listed in Figure 1A were applied to identify papers reporting quantitative results of methylation profile performed on HCC patients and the relative adjacent tissue as control.

Figure 1 Flow chart. A: Workflow of Data collection, analysis and database compiling; B: Bioinformatics flow for selecting differentially methylated CpG sites from HCC meta-analysis; C: Bioinformatics flow for selecting differentially methylated genes from HCC meta-analysis. HCC: Hepatocellular carcinoma.

Available patient data, including etiology of liver disease (HCV, HBV, alcohol, fatty liver disease) on the basis of which the HCC tumors developed, presence of cirrhosis, the Model for End-stage Liver Disease score (MELD score, an assessment of the severity of liver dysfunction), tumor histology, stage of cancer, alpha-fetoprotein (AFP) level, overall and recurrence-free survival following treatment were also documented.

We identified 928 abstracts retrieved by the search on PubMed and 233 results were obtained from GEO. The flow chart outlining the selection process is detailed in Figure 1A. Details regarding the 11 included studies[8,19-30], together with the information on number of samples per group, per study are provided in Table 1.

Table 1 Datasets used for the meta-analysis.

No	GEO Accession	PubMed ID	HCC samples	Adjacent tissue samples
1		21500188	13	12
2		243066621	45	45
3		253762921	22	22
4	GSE592602	25945129	8	8
5	GSE297202	21747116	12	12
6	GSE180812	20165882	20	20
7	GSE37988	22234943	62	62
8	GSE44970	24012984	20	8
9	GSE54503	23208076	66	66
10	GSE57956	25093504	59	59
11	GSE607532	25294808	27	27

¹Only genes are reported, without information on CpG sites;

²Differential methylation data analysis performed in-house.

Demographic and clinical characteristics

Demographic and clinical patient information pertaining to each dataset are presented in Tables 2 and 3.

Table 2 Availability of clinicopathological information on the 11 datasets used for the integrative analysis of genome-wide DNA methylation.

Clinical-pathological information	n = 11	%
Cirrhosis status in HCC samples	6	55
Child-Pugh/MELD Score	3	27
HCC Etiology	10	91
Alphafetoprotein level	5	45
Tumor grade	5	45
Tumor stage	5	45
Survival data	3	27

Table 3 Clinicopathological information of the individual datasets and methodology used for methylation analysis.

Dataset	Year	PMID	GEO dataset	HCC (n)	Controls (n)	Liver cirrhosis in HCC samples	Etiology of liver disease (n)	Method
1	2011	21500188		13	12	Y (12)	HBV (3), HCV (4), alcoholic (6)	Human methylation 27 DNA analysis bead-chip
2	2014	24306662		45	45	Y (120), N (34)	HBV (149), HCV (1), nonviral (4)	Illumina GoldenGate Methylation Beadarray Cancer Panel I
3	2014	25376292		22	22	N/A	HBV (1), HCV (9), alcohol (4), other (8)	Infinium Humanmethylation 27 Beadchip
4	2015	25945129	GSE59260	8	8	N/A	HBV-HCV-(8)	Nimblegen Human DNA Methylation 3 x 720K CpG Island PI
5	2011	21747116	GSE29720	12	12	N/A	N/A	Agilent-017075 human hg 18 promoter 800-200
6	2010	20165882	GSE18081	20	20	Y (20)	HCV (20)	Illumina Golden Gate Methylation Beadarray Cancer Panel I
7	2012	22234943	GSE37988	62	62	N/A	HBV-HCV-(7), HBV + HCV-(36), HBV-HCV+(6), HBV + HCV+(13)	Illumina Human Methylation27 Beadchip
8	2013	24012984	GSE44970	20	8	N/A	HCV (8)	Human Methylation27 Beadchip
9	2013	23208076	GSE54503	66	66	Y (48), N (17), missing (1)	HBV-HCV-(19), HCV (19), HBV (13), HBV + HCV (4), missing (11)	Infinium Human Methylation 450K Beadchip
10	2014	25093504	GSE57956	59	59	Y (37), N (21)	HBV+(36), HBV-(23)	Intinium Ilumanmethylation27 Beadchip
11	2014	25294808	GSE60753	27	27	Y (26), N (1)	HBV (1), HCV (7), alcohol (9), other (10)	Iinfinium -450K Human Methy lation Beadchip

Only 6 out of 11 papers (55%) included details regarding presence/absence of cirrhosis, and 10 out 11 papers (91%) provided details regarding the etiology of liver disease. MELD or Child-Pugh score were provided in 27% of papers. Less than half of the selected publications, 5 out of 11 (45%), included details regarding the stage of cancer, although all studies were performed using hepatectomy patient samples. The same trend is identified for information regarding the histologic grade of the tumor (well-, moderately-, or poorly-differentiated tumors). Only 5 out of 11 papers (45%) had alpha-fetoprotein levels available. Overall survival and HCC recurrence statistics as follow-up data were available in 3/11 studies (27%).

Genomic region selection

For the final network-based integrative analysis, we selected 11 datasets (Table 1). Out of these, raw data from eight studies were available on the GEO website (https://http://www.ncbi.nlm.nih.gov/geo/). Except for the GSE60753[25] dataset, for which we have performed our own analysis with R[31] (due to comparison between sample groups required being different from the main paper), we selected the CpG sites or genes reported to be hyper-or hypo- methylated in the corresponding publications. In 6/11 datasets, CpGs and the mapped genes were selected, whereas in the remaining 5/11 datasets, only the genes without CpG sites were found. Therefore, we separated our analysis into two parts: (1) Taking into account approximately 13500 CpG sites provided or obtained with R analysis (Figure 1B); and (2) considering only the 4122 differentially methylated genes without information on the corresponding CpG sites (Figure 1C). In both cases, the genomic region is considered differentially methylated between HCC tissue and the adjacent non-tumoral sample, if the FDR[32] corrected P-value < 0.01. Furthermore, we filtered out everything that did not satisfy the criteria: ∆β ≥ 0.20 or ∆β ≤ -0.20, where ∆β = β_HCC - β_adjacent was the difference in methylation between above specified groups. When the CpG sites were considered, the Illumina HumanMethylation450K and 27K platforms were used for mapping to the genes. When multiple sites or genes were found having the same sense of differential methylation, the mean value of ∆β was calculated. The CpGs in the 5’UTR, 1^st Exon, TSS200, TSS1500 or in CpG islands were considered in the promoter and all other CpGs were considered to be in the body of the gene.

Pathway and network analysis

Two final lists of differentially methylated genes corresponding to CpGs in the promoter (n = 765, Supplementary Table 1) or to the body of the gene (n = 411, Supplementary Table 2) and their corresponding mean (∆β) were uploaded into IPA (Ingenuity Systems^®, http://www.ingenuity.com). Based on the manually-curated Ingenuity Knowledge Base derived from experiments and findings published in top peer-reviewed journals, IPA identifies a series of canonical pathways, diseases and functions or networks associated with the molecules in the input list. For each of these, a P-value is calculated with the right-tailed Fisher’s exact test[33], which takes into account the number of focus molecules (input genes) in the network and the total number of molecules in the IPA database that could be included in the corresponding networks.

RESULTS

Based on datasets with known CpG sites, the most frequently hypermethylated genes in the promoter region included HIST1H2AJ, which is a histone protein and SPDYA, a cell cycle regulator known to trigger transition from G1 to S phase. The HRNBP3 gene, an RNA-binding protein, was the most commonly hypomethylated in the promoter region. Further details are provided in Supplementary Tables 1 and 2.

Using the differentially methylated CpG sites in the promoter

Canonical pathways: The most significantly associated pathways with our list of differentially methylated genes are given in Table 4. G-protein coupled receptor signaling, Transcriptional Regulatory Network in Embryonic Cells, cAMP-mediated signaling were the top hits.

Table 4 Top canonical pathways identified by IPA for the genes corresponding to CpG sites in promoter.

Ingenuity canonical pathways	-log (P-value)	Molecules
G-protein coupled receptor signaling	3.84E+00	DRD5, GNA11, VIPR2, ADCY5, ADRB1, CNR1, PIK3R5, NPY1R, FPR1, FFAR3, NFKBID, MC2R, PDPK1, GRM4, MC3R, CXCR2, PRKAR1B, DRD4, PDE6B, HCAR2, DUSP4, PTGDR
Transcriptional regulatory network in embryonic stem cells	3.42E+00	MYF5, SIX3, PAX6, GBX2, CDYL, FOXD3, ONECUT1, FOXC1
cAMP-mediated signaling	3.22E+00	DRD5, VIPR2, ADCY5, ADRB1, CNR1, NPY1R, FPR1, FFAR3, MC2R, GRM4, MC3R, CXCR2, PRKAR1B, DRD4, PDE6B, HCAR2, DUSP4, PTGDR

Diseases and functions: Not surprisingly Cancer, Organismal Injury and Abnormalities were identified as some of the top diseases and functions. Among different types of cancer, abdominal cancer (P = 1.7E-210), digestive system cancer (P = 7.88E-14), abdominal carcinoma and digestive organ tumor (P = 8.58E-13–1.76E-12) were listed. Cellular development (P-value=3.22E-03–2.66E-08), growth and proliferation (P = 3.22E-03–6.99E-06) were the most important molecular and cellular functions associated with HCC methylated data.

Networks: Among the most statistically and biologically significant networks associated with the genes differentially methylated in the promoter region, three of them captured our attention: Organismal Development, Organismal Injury and Abnormalities, Cellular Development (Figure 2A), Lipid Metabolism, Small Molecule Biochemistry, Cell Death and Survival (Figure 2B) and Cell-to-Cell Signaling and Interaction, Drug Metabolism, Small Molecule Biochemistry (Figure 2C). We noted that networks 1 and 3 were mainly formed by the hyper-methylated genes, whereas network 2 is constituted by both hyper- and hypo-methylated genes approximately equally.

Figure 2 Networks associated with differentially methylated CpG sites in the promoter regions of genes. A: Organismal development, organismal injury and abnormalities, cellular development; B: Lipid metabolism, small molecule biochemistry, cell death and survival; C: Cell-to-cell signaling and interaction, drug metabolism, small molecule biochemistry.

Using the differentially methylated CpG sites in the body of the gene

Canonical pathways: G-protein coupled receptor signaling and cAMP-mediated signaling were among the top 20 hits (Table 5), which was similar to the pathways generated by the genes with methylated CpG sites in the promoter region.

Table 5 Canonical pathways identified by IPA for the genes with methylation differences in the body of the gene in hepatocellular carcinoma.

Ingenuity canonical pathways	-log (P-value)	Molecules
Aryl hydrocarbon receptor signaling	2.48E+00	GSTM1, CCND2, TFF1, ALDH1L2, GSTM2, TP73, GSTP1, ALDH3A1
G-protein coupled receptor signaling	2.03E+00	CAMK2B, RGS7, GABBR1, PDE4D, ADCY2, PDE1C, ADRA1D, NPR3, PDE10A, GRM6, PRKCG
cAMP-mediated signaling	1.73E+00	CAMK2B, RGS7, GABBR1, PDE4D, ADCY2, PDE1C, NPR3, PDE10A, GRM6

Diseases and functions: Our study shows that DNA methylation in HCC patients leads to the same diseases and functions, regardless of the CpG site position (promoter or body), with cancer and, in particular, abdominal/digestive system cancer among the top listed by IPA.

Networks: From the top detected networks, the Cell-To-Cell Signaling and Interaction, Cellular Assembly and Organization, Cellular Function and Maintenance (Figure 3A) and Drug Metabolism, Glutathione Depletion in Liver, Small Molecule Biochemistry (Figure 3B) are closely related to HCC.

Figure 3 Networks associated with differentially methylated CpG sites in body of the gene. A: Cellular assembly and organization, cellular function and maintenance; B: Drug metabolism, glutathione depletion in liver, small molecule biochemistry.

Validation using only the reported differentially methylated genes

Having identified the differentially methylated genes corresponding to the CpG sites in the promoter or the body of the genes through this integrative meta-analysis, we wanted to verify how many of these overlapped with the reported differentially methylated genes in those studies that did not include information on CpG sites. The Venn diagram in Figure 4 shows 165 genes reported genes in common with CpG sites in the promoter and 82 genes in common with CpG sites in the body (Supplementary Table 3).

Figure 4 Venn diagram intersecting three lists: (1) reported differentially methylated genes in HCC from studies not providing information on the corresponding CpG sites; (2) identified differentially methylated genes in HCC corresponding to CpG sites in promoter; and (3) identified differentially methylated genes in HCC corresponding to CpG sites in body of the gene.

DISCUSSION

The literature on the epigenome in HCC has grown since the advent of tools permitting genome-wide methylation analysis. Epigenetic changes in HCC arise in the context of various etiologies of chronic liver disease, and have been revealed to contribute to tumorigenesis and cancer progression. Therapeutic targeting of HCC has not been as successful as in other malignancies, and requires exploration of a different approach[34,35]. Given that this cancer is driven by various known environmental factors, targeting epigenetic changes in HCC represents a potentially promising therapeutic avenue[36].

The current study is the largest network-based integrative meta-analysis of all publicly available genome-wide DNA methylation data in HCC, with 354 HCC samples represented. These HCCs had arisen mainly in the context of viral hepatitis B and C, with only a few occurring in patients with alcoholic cirrhosis. Therefore, the literature on methylation in HCC is heavily weighted towards epigenetic changes from viral infection, and the aberrantly methylated genes in our analysis will likewise be influenced by the greater proportion of hepatitis B and C. Clinical information regarding tumor grade, disease stage, and survival were only available in around half of the datasets, thereby limiting the ability to correlate with histopathological characteristics and disease outcome. The genome-wide DNA methylation datasets included in our integrative analysis were published from 2010 to 2015.

Network-based tools offer a different and unique perspective into the key genes and pathways implicated in disease pathogenesis and progression[37]. Network-based medicine is critical to a broader understanding of HCC, whose pathogenesis has been difficult to elucidate given the multiplicity of underlying liver disease etiologies[30]. Epigenetic changes impact genetic networks, and a network-based integrative meta-analysis is ideally suited to integrating and exploring effects of networks on disease pathogenesis[38]. Using IPA, we performed this integrative analysis in order to identify the most commonly aberrantly methylated genes and associated pathways. The most commonly hyper- and hypomethylated genes were identified. These included HIST1H2AJ, which is a histone-coding cell cycle gene previously also identified as hypermethylated in patient lung adenocarcinoma samples[37] and head and neck squamous cell carcinomas[39]. In a study investigating the genetic-and-epigenetic cell cycle network in HeLa cancer cells, methylation of HIST1H2AJ (among other genes) was found to result in cell proliferation and anti-apoptosis through NFκB, TGF-β, and PI3K/Akt/mTOR pathways[40]. HIST1H2AJ has not previously been highlighted as a gene of interest in HCC, which is a novel finding of our integrative analysis of methylation datasets. This illustrates the power of integrating all available high-throughput data to better understand important genes in cancer. SPDYA, a cell cycle regulator known to trigger transition from G1 to S phase, was differentially hypermethylated. The HRNBP3 gene, an RNA-binding protein, was the most commonly hypomethylated, as had been reported in the integrative analysis of epigenetic data by Song et al[16] in 2012. One would thereby anticipate increased gene expression of HRNBP3 in HCC. These genes with differential methylation have not previously been highlighted in the HCC literature, and serve as potential new biomarkers and therapeutic targets[41].

Using IPA, we then determined the most commonly affected networks in HCC.

G-protein coupled receptor signaling, Transcriptional Regulatory Network in Embryonic Cells, cAMP-mediated signaling were the top hits, which was in perfect agreement with the work of Song et al[16], wherein they used IPA to analyze methylation profiling for a set of 27 HCC tumors compared with 20 normal patients. G-protein coupled receptor signalling is common to various principal pathways known to be implicated in HCC, including the PI3K/Akt/mTOR and Ras/Raf/MAPK pathways based on genomic and gene expression analyses. Therefore, our results reinforce the biological rationale of targeting these pathways. We also elucidated the crosstalk between proteins within the networks of interest to HCC. This analysis revealed ESR1 and ERK to be proteins central to they key networks.

A unique aspect of our study was the analysis of methylation at CpG sites in both the promoter and body of genes. Whereas methylation in the gene promoter is known to cause transcriptional repression, methylation in the body has the opposite effect, promoting gene expression. A novel finding was that the genes with the greatest differential methylation in the promoter were the same as those with the greatest differential methylation in the body, further confirming the importance of these genes. We were also able to validate the identity of several genes with data on CpG sites within datasets without such data available.

Limitations of our study include the lack of methylation data on individual HCC samples. Given the relatively recent advent of genome-wide methylation analysis methods, with the earliest dataset in HCC being released in 2010, this analysis was representative of only 354 samples in comparison to a similar number of non-cancerous liver samples. Nonetheless, our study is the largest integrative network-based analysis of DNA methylation in HCC. Clinicopathological characteristics such as grade, stage and survival were available only for half of the datasets, thereby limiting the ability to correlate these data points with the most aberrantly methylated genes. Finally, these data were most representative of hepatitis B and C, as described above.

In conclusion, our integrative analysis of genome-wide DNA methylation represents the largest such study in HCC. By integrating all genome-wide DNA methylation data with network-based tools, we have systematically elucidated the landscape of epigenetic DNA modifications in HCC and identified novel potential biomarkers and targetable genes within known pathways of interest to HCC. Therapeutic targeting of the epigenome in HCC is a potential avenue to address this malignancy that arises in the context of various etiologies of chronic liver disease.

ARTICLE HIGHLIGHTS

Research background

The advent of high-throughput technologies in epigenetics has led to improved characterization of methylation status and its impact on development of Hepatocellular carcinoma (HCC).

Research motivation

HCC is a malignancy that arises in the context of ongoing liver injury from various causes, such as hepatitis B, hepatitis C, alcoholic and non-alcoholic liver disease. Therefore, epigenetic changes are very likely to contribute to the pathogenesis of this malignancy.

Research objectives

We aimed to identify the key epigenetically modulated genes and pathways in HCC by performing an integrative meta-analysis of all major, well-annotated and publicly available methylation datasets using tools of network analysis.

Research methods

PubMed and Gene Expression Omnibus were searched for genome-wide DNA methylation datasets. Patient clinical and demographic characteristics were obtained. DNA methylation data were integrated using the Ingenuity Pathway Analysis, a software package for visualizing and analyzing biological networks. Pathway enrichment analysis was performed using IPA, which also provides literature-driven and computationally-predicted annotations for significant association of genes to curated molecular pathways.

Research results

From an initial 928 potential abstracts, we identified and analyzed 11 eligible high-throughput methylation datasets representing 354 patients. A significant proportion of studies did not provide concomitant clinical data. In the promoter region, HIST1H2AJ and SPDYA were the most commonly methylated, whereas HRNBP3 gene was the most commonly hypomethylated. ESR1 and ERK were central genes in the principal networks. The pathways most associated with the frequently methylated genes were G-protein coupled receptor and cAMP-mediated signalling.

Research conclusions

Using an integrative network-based analysis approach of genome-wide DNA methylation data of both the promoter and body of genes, we identified G-protein coupled receptor signalling as the most highly associated with HCC. This encompasses a diverse range of cancer pathways, such as the PI3K/Akt/mTOR and Ras/Raf/MAPK pathways, and is therefore supportive of previous literature on gene expression in HCC. However, there are novel targetable genes such as HIST1H2AJ that are epigenetically modified, suggesting their potential as biomarkers and for therapeutic targeting of the HCC epigenome.

Research perspectives

Our integrative analysis of genome-wide DNA methylation represents the largest such study in HCC. By integrating all genome-wide DNA methylation data with network-based tools, we have systematically elucidated the landscape of epigenetic DNA modifications in HCC and identified novel potential biomarkers and targetable genes within known pathways of interest to HCC. Therapeutic targeting of the epigenome in HCC is a potential avenue to address this malignancy that arises in the context of various etiologies of chronic liver disease.