Skip to main content

Differentially methylated loci in NAFLD cirrhosis are associated with key signaling pathways


Altered DNA methylation events contribute to the pathogenesis and progression of metabolic disorders, including nonalcoholic fatty liver disease (NAFLD). Investigations of global DNA methylation patterns in liver biopsies representing severe NAFLD fibrosis have been limited. We used the HumanMethylation 450K BeadChip to analyze genome-wide methylation in patients with biopsy-proven grade 3/4 NAFLD fibrosis/cirrhosis (N = 14) and age- and sex-matched controls with normal histology (N = 15). We identified 208 CpG islands (CGIs), including 99 hypomethylated and 109 hypermethylated CGIs, showing statistically significant evidence (adjusted P value < 0.05) for differential methylation between cirrhotic and normal samples. Comparison of β values for each CGI to the read count of its corresponding gene obtained from RNA-sequencing analysis revealed negative correlation (adjusted P value < 0.05) for 34 transcripts. These findings provide supporting evidence for a role for CpG methylation in the pathogenesis of NAFLD-related cirrhosis, including confirmation of previously reported differentially methylated CGIs, and contribute new insight into the molecular mechanisms underlying the initiation and progression of liver fibrosis and cirrhosis.


Obesity and insulin resistance are commonly associated with fat accumulation in the liver, which is classified as nonalcoholic fatty liver disease (NAFLD) [1, 2]. This common disease encompasses a wide range of liver conditions ranging from hepatic steatosis to nonalcoholic steatohepatitis (NASH) in the absence of significant alcohol consumption [3]. NAFLD/NASH has become a leading cause of cryptogenic cirrhosis [4] and is currently the third leading clinical indication for liver transplantation in the USA [5]. The molecular pathogenesis underlying NAFLD appears to involve many factors, including those of genetic and environmental origins [6].

The relationship between NAFLD and obesity raises the possibility that epigenetic mechanisms associated with “metabolic memory” may be involved [7]. DNA methylation of CpG dinucleotides is a well-characterized epigenetic modification that occurs mostly within CpG islands (CGIs) in promoter regions and functions to regulate gene expression. Evidence supporting a role for DNA methylation in the development of metabolic diseases, including NAFLD, is emerging. For example, diets enriched for energy and certain macronutrients appear to induce transgenerational epigenetic changes in rodent models that predispose offspring to hepatic steatosis and NASH [8]. Changes in DNA methylation have been associated with hepatic steatosis [9, 10] and progression of fibrosis [11], while methyl-depleted diets contribute to the development of steatohepatitis, cirrhosis, and liver cancer in rodents [12, 13]. In humans, hypermethylation of NADH dehydrogenase 6 [14] and PPARGC1A has been correlated with NAFLD [15]. To our knowledge, only a few studies have investigated differential methylation of genes using a genome-wide approach in patients with NASH [16,17,18,19]. Several studies have also recently identified changes in the methylation profile of specific genes in NAFLD and NASH [20, 21], and type 2 diabetes (T2D) has been associated with epigenetic modifications in human liver [22].

To complement previous studies, we assessed DNA methylation status in liver biopsies of individuals with histologically documented NAFLD-related cirrhosis in the setting of extreme obesity. Because methylation of CpG sites in promoter regions is related to transcriptional regulation, we also examined the relationship between methylation status and gene expression. We identified over 30 CGIs that were differentially methylated in NAFLD cirrhosis and correlated with hepatic gene expression in the same liver samples. These findings contribute supporting evidence for a role for CpG methylation in the pathogenesis of NAFLD-related cirrhosis, including confirmation of almost 90 previously reported differentially methylated CpG sites, and provide new insight into the molecular mechanisms underlying the initiation and progression of liver fibrosis and cirrhosis.

Materials and methods

Study sample

Liver wedge biopsies were intraoperatively obtained from Caucasian women enrolled in the Bariatric Surgery Program at the Geisinger Clinic Center for Nutrition and Weight Management and histologically evaluated using NASH CRN criteria as described [23,24,25]. Patients with histologic or serologic evidence for other chronic liver diseases were excluded from this study. Both medical history and histological assessment excluded individuals with clinically significant alcohol intake and drug use from participation in the bariatric surgery program. Clinical data, including demographics, clinical measures, ICD-9 codes, medical history, medication use, and common lab results, were available for all study participants as described previously [26]. All study participants provided written informed consent for research, which was conducted according to The Code of Ethics of the World Medical Association (Declaration of Helsinki). The Institutional Review Boards of Geisinger Health System, the Translational Genomics Research Institute, and the Lewis Katz School of Medicine at Temple University approved the research.

Methylation analysis

Whole genome methylation profiling was performed using the 450K Infinium Methylation BeadChip Assay (Illumina; San Diego, CA, USA). Liver genomic DNA was extracted using the QIAamp DNA kit (Qiagen; Germantown, MD, USA), and DNA concentration was determined using the Quant-iT PicoGreen dsDNA Assay kit (Thermo Fisher Scientific; Waltham, MA). Bisulfite conversion of genomic DNA was performed using the EZ DNA Methylation Kit (Zymo Research; Irvine, CA, USA). Bisulfite-treated DNA was then hybridized to arrays according to the manufacturer’s protocol. Methylation levels for each CpG residue were estimated as the ratio of the methylated signal intensity over the sum of the methylated and unmethylated intensities at each locus using the minfi R package and presented as β values [27]. For quality control, we used the publicly available software CpGassoc [28] to exclude any samples with probe detection call rates < 95%, as well as those with an average intensity value of either < 50% of the experiment-wide sample mean or < 2000 arbitrary units (AU). Data points were set to missing with detection P values > 0.01. Probes overlapping with copy number variants, as well as those mapping to multiple locations with up to two mismatches, were excluded from the analysis. All samples were checked for atypical raw intensity distributions, and β value correlation among others within each respective group. Data were normalized using the SWAN algorithm provided by the minfi package (

To identify CpG sites showing differentially methylation between fibrotic and normal samples, we used the City of Hope CpG Island Analysis Pipeline (COHCAP) [29] and the minfi package. Differentially methylated site analysis encompasses testing each genomic position for association between methylation and phenotype with a F-test for categorical outcome. CpG sites were defined as methylated if they showed a percentage of methylation > 0.5 and unmethylated if they had β values < 0.5. To identify differentially methylated CGI, signal from differentially methylated CpG sites within a CGI was averaged for each sample before group comparison. We used a FDR < 0.05 [30] to define a set of differentially methylated CGIs.

RNA extraction, sequencing, and analysis

Total RNA was extracted using the RNAeasy kit (Qiagen), quantified using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific), and converted to cDNA using the Ovation RNA-Seq System V2 (NuGEN; San Carlos, CA, USA). We performed massively parallel RNA-Seq using polyA-selected RNA from biopsied liver samples using the Illumina HiSeq2000 platform, sequencing to a depth of 60 M 100 bp paired-end reads. During the generation of qseq and fastq files for alignment, low-quality reads were identified and removed, and the indexed reads identified and grouped accordingly. Filtered reads were aligned with the human genome using the Bowtie program [31]. Aligned RNA-Seq reads were imported into the Cufflinks program [32] to assemble alignments into a parsimonious set of transcripts and this set of annotated transcripts were quantified in each sample using DESeq2 [33]. To identify genes with negative expression correlations, we compared β values for each CGI to the fold-change of its corresponding gene obtained from RNA-seq analysis. The correlation was estimated using Spearman’s rho statistic, with P values computed via the asymptotic t distribution approximation.

Functional enrichment analysis

The Ingenuity Pathway Analysis (IPA) software (Qiagen) was used to identify canonical signaling pathways and network connections associated with CGIs. The significance of the association between CGI and canonical pathway was assessed using two criteria: (1) the ratio of the number of molecules mapped to the pathway and total number of molecules involved in the canonical pathway and (2) the Benjamini-Hochberg corrected P value from the right-tailed Fisher exact test.


Study sample characteristics

DNA was obtained from liver biopsies from 15 patients with NAFLD, including 4 with stage 3 bridging fibrosis (F3) and 11 with either incomplete cirrhosis or stage 4 (F4) fibrosis (cirrhosis), and 15 individuals with normal liver histology matched by age, sex, BMI at biopsy, and T2D status. The demographic information and clinical characteristics of the study participants are shown in Table 1. Individuals with NAFLD fibrosis had significantly higher serum levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT), insulin, and triglycerides compared to individuals without fibrosis. NAFLD patients with fibrosis also manifested lobular inflammation, portal inflammation, and hepatocyte ballooning.

Table 1 Patient demographics and clinical characteristics

Identification of differentially methylated CGIs

For each hepatic DNA sample, methylation levels for each CpG residue were estimated as β values, which represent the ratio of the methylated signal intensity over the sum of the methylated and unmethylated intensities. Initial hierarchical clustering analysis revealed that F3 fibrosis could be segregated as a separate group from F4, consistent with the distinct histology of stage 3 versus stage 4. To avoid the introduction of unnecessary heterogeneity, we excluded these three samples from further analysis. In addition, one NAFLD fibrosis sample failed quality control measures and was removed from further analysis. Using data from the remaining 11 cirrhosis and 15 normal samples, we identified 208 CGIs, including 99 hypomethylated and 109 hypermethylated CGIs (Additional file 1: Table S1) showing statistically significant evidence (adjusted P value < 0.05) for differential methylation. Analysis of individual sites revealed differential methylation of 4275 CpG sites, (adjusted P value < 0.05), corresponding to 1713 hypomethylated and 2562 hypermethylated CpG sites; the top 100 differentially methylated CpG sites are shown in Additional file 2: Table S2. The top CGI, prioritized by strength of statistically significant evidence for differential methylation, magnitude of difference in methylation levels (Δ beta), and biological relevance to general liver function or metabolism, are shown in Table 2 (hypermethylation) and Table 3 (hypomethylation).

Table 2 CpG islands showing the strongest evidence for hypermethylation in NAFLD cirrhosis
Table 3 Ten CpG islands showing the strongest evidence for hypomethylation in NAFLD cirrhosis

Correlation of CpG methylation with hepatic gene expression

To identify gene-methylation correlations, we compared β values for each CGI to the FPKM (fragments per kilobase of transcript per million mapped reads) data of its corresponding gene. Focusing on pairwise associations with a significant negative correlation at FDR < 0.05, we found evidence for negative correlation between 34 CGI-transcript pairs (Table 4). Box plots of mean (± standard deviation) methylation β values for CGIs associated with gene expression differences in patients with NAFLD fibrosis (n = 11) compared to individuals with normal liver histology (n = 15) are presented in Additional file 3: Figure S1.

Table 4 Integration with gene expression based on correlation statistic

Functional enrichment analysis

We used functional enrichment analyses to identify the effect of differentially methylated CGI on different canonical pathways. The analysis revealed 113 enriched canonical pathways (Additional file 4: Table S3). Pathways with relevance to NAFLD cirrhosis identified in this analysis are highlighted in Table 5. The top pathways included production of nitric oxide and reactive oxygen species, LXR/RXR activation, and FXR/RXR activation.

Table 5 Functional enrichment analysis: canonical pathways


To our knowledge, only four studies have investigated methylation profiles in NAFLD patients [16,17,18,19]. Although all four studies utilized the same array-based platform for measuring DNA methylation, they differed with respect to analytical strategy, experimental design, and study sample composition. In the first published study, levels of methylation and mRNA expression were assessed in normal-weight controls (N = 18) and obese individuals with liver histology consistent with normal (N = 18), steatosis (N = 12), or NASH (N = 15). Seventy-four differentially methylated CpG sites (FDR < 0.004) were found in comparisons of the four phenotypic groups with NAFLD-specific methylation and expression differences observed for nine genes: ACLY (ATP citrate lyase), GALNTL4 (polypeptide N-acetylgalactosaminyltransferase 18), GRID1 (glutamate ionotropic receptor delta type subunit 1), IGF1 (insulin like growth factor 1), IGFBP2 (insulin like growth factor binding protein 2), IP6K3 (inositol hexakisphosphate kinase 3), PC (pyruvate carboxylase), PLCG1 (phospholipase C gamma 1), and PRKCE (protein kinase C epsilon) [16].

In the same year, Murphy et al. [18] identified 69,247 differentially methylated CpG sites, most of which were hypomethylated, in obese patients with advanced (N = 23) versus mild (N = 33) NAFLD [18]. In that study, methylation was correlated with gene transcript abundance levels for 7% of the differentially methylated CpG sites and methylation at FGFR2, MATA1, and CASP1 was validated in a replication cohort. Despite the similar design between that study and the current work, there were notable differences between the two, which could account for the disparate findings. In the published study, (1) the mild fibrosis group was a mix of NAFLD histological types including grade 1 fibrosis; (2) the advanced NAFLD group included stage 3 fibrosis; (3) the BMI range was 32.0–33.8 across case and controls in both discovery and replication cohorts; and (4) the case and control groups were not matched for T2D. More recently, de Mello et al. [17] identified 1292 CpG sites, representing 677 genes that showed differences in DNA methylation between 26 NASH patients and 34 individuals with normal liver histology, independent of T2D status, age, sex, and BMI. In the present work, we also used a homogeneous control group without histopathological manifestations of NAFLD. The lack of a normal control group is widespread in studies of NAFLD, where liver tissue is often obtained in conjunction with clinically indicated biopsies; thus, a high pre-procedure probability of a pathological finding often results in the presence of some level of NAFLD. We also found that stage 3 fibrosis could be segregated as an epigenetically separate group using hierarchical clustering analysis; thus, we removed these samples from further analysis. This level of bioinformatic quality control is critical when analyzing large genomic datasets such as genome-wide methylation array data.

In the most recent study, Hotta et al. [19] identified a number of differentially methylated regions in a comparison of patients with either mild or advanced NAFLD fibrosis from Japan. That study differed from the current study, not only with respect to ethnic background, but also the NAFLD groups were heterogeneous with respect to fibrosis status and neither the mild nor the advanced fibrosis groups were obese [19].

Despite the differences in design, we compared our results with those reported by published studies discussed above and identified 86 common differentially methylated sites (Additional file 5: Table S4). Interestingly, four (AQP1, FGFR2, RBP5, and MGMT) overlapped with differentially methylated CpG sites in all four studies, suggesting that these genes may be a common core set with potential importance in disease pathogenesis. AQP1 (aquaporin 1) encodes a water channel that is overexpressed in fibrosis and cirrhosis and appears to promote hepatic fibrosis through mechanisms involving pathological angiogenesis [34]; association of increased AQP1 expression is associated with hypomethylation of CpG site found among these studies. FGFR2 (fibroblast growth factor receptor 2), which was found to be hypomethylated across these studies, has been linked with liver fibrosis. FGF2 levels were increased in hepatic stellate cell activation [35] and liver cirrhosis [36]. Further, FGF2 knockout in the carbon tetrachloride mouse model of hepatic injury was associated with decreased collagen expression and protection from liver fibrosis [37], while Brivanib, an inhibitor of FGFR and vascular endothelial growth factor (VEGF) was shown to inhibit activation of hepatic stellate cells in vitro and liver fibrosis in three different animal models [38]. The biological relevance of the other two genes, retinol-binding protein 5 (RBP5) and O6-methylguanine-DNA methyltransferase (MGMT), both of which were hypomethylated in these studies, to the development of NAFLD-related fibrosis remains unknown.

The most highly associated canonical pathways identified through functional enrichment analysis provide additional evidence for epigenetic regulation in the pathophysiology of NAFLD-related cirrhosis. The most highly enriched pathway, Production of Nitric Oxide and Reactive Oxygen Species in Macrophages, implicates a well-studied molecular mechanism in the progression of NAFLD to cirrhosis. In human studies and in rodent models of NAFLD/NASH, increased reactive oxygen species appear to be a central feature of lipotoxicity and inflammation [39]. Although the use of antioxidants as a therapeutic strategy has not resulted in significant effects in either humans or animal models, most studies suffer from one or more apparent weaknesses [40]. Our results suggest that epigenetic modulation of antioxidant pathways could be a promising therapeutic approach. The next two associated pathways, LXR/RXR activation and FXR/RXR activation, are mechanistically linked via heterodimerization of the nuclear receptors farnesoid X receptor (FXR) and retinoid X receptor (RXR) [41]. FXR is highly expressed in the liver and is involved in several key metabolic processes including bile acid synthesis, glucose and lipid metabolism, and regulation of inflammatory pathways [42]. Administration of the synthetic bile acid derivative and FXR agonist, obeticholic acid, has been shown to reduce the histological NAFLD fibrosis score in a multicenter, randomized, double blind, placebo-controlled study [43]. To date, the role of epigenetic regulation of this pathway in the potential clinical effectiveness of bile acid-based therapeutics has not yet been explored.

Despite the overlap in findings of differential methylation among studies, our sample size, while comparable to several other reports [16, 44, 45], is still limited. We dichotomized our cohort into extreme histologic phenotypes to increase our power to identify clinically relevant differences in DNA methylation levels and were therefore unable to stratify into phenotypic groups to assess intermediate levels of peri-sinusoidal fibrosis or portal fibrosis. Additional studies including individuals spanning the spectrum of NAFLD fibrosis will be critical to extend our findings. We also note that the cross-sectional design does not allow associations with disease progression to be drawn. In the absence of serial liver biopsies in the present cohort, we were not able to utilize a prospective design. Assessment of methylation levels in longitudinal biopsies will be necessary to determine the roles of specific candidates in disease progression. In addition, due to the small sample size in this study, we focused our investigation on females with extreme obesity, who underwent gastric bypass surgery and had liver biopsies. We utilized this design because obesity is a significant risk factor for NAFLD [46], and liver biopsy tissue was obtained without clinical indication, which effectively eliminates bias toward clinically suspect liver disease. Despite the relative phenotypic homogeneity, we were able to replicate findings of differential methylation from study samples encompassing a greater range of phenotypic diversity, indicating the influence of shared epigenetic effects. However, conclusions with regard to NAFLD-related fibrosis obtained from this cohort may still not be relevant to populations with less severe obesity or of different ethnicities without additional validation.

In summary, the results obtained in the current study not only confirm previous findings of differential methylation in NAFLD patients but also provide novel evidence of differential methylation and RNA expression profiles associated with NAFLD-related fibrosis. Future studies will be needed to determine the extent to which DNA methylation patterns in the liver are represented in other metabolically relevant tissues such as visceral and subcutaneous fat [47, 48], as well as peripheral blood leukocytes, which will be critical for the development of non-invasive markers of NAFLD stage. Additional studies, including those showing functional consequences of differentially methylated sites, i.e., disruption of transcription factor binding, will be necessary to confirm the role of specific CpG loci in liver fibrosis. These approaches are expected to yield new insights into the pathological mechanisms underlying the development of fibrosis and cirrhosis in NAFLD.



ATP citrate lyase


Alanine aminotransaminase


Aquaporin 1


Aspartate aminotransaminase


Body mass index


Caspase 1


Cytosine phosphate guanine islands


City of Hope CpG Island Analysis Pipeline


Cytosine phosphate guanine


False discovery rate


Fibroblast growth factor receptor 2


Fragments per kilobase of transcript per million mapped reads


Farnesoid X receptor


Polypeptide N-acetylgalactosaminyltransferase 18


Glutamate ionotropic receptor delta type subunit 1


International Classification of Diseases-ninth revision

IGF1 :

Insulin like growth factor 1


Insulin-like growth factor binding protein 2

IP6K3 :

Inositol hexakisphosphate kinase 3


Ingenuity Pathway Analysis


Liver X receptor


Methionine adenosyltransferase 1A


O6-Methylguanine-DNA methyltransferase


Nicotinamide adenine dinucleotide


Nonalcoholic fatty liver disease


Nonalcoholic steatohepatitis clinical research network


Nonalcoholic steatohepatitis

PC :

Pyruvate carboxylase


Phospholipase C gamma 1


Peroxisome proliferator-activated receptor gamma coactivator 1-alpha


Protein kinase C epsilon


Retinol binding protein 5


Retinoid X receptor


Type 2 diabetes


Vascular endothelial growth factor


  1. 1.

    Wang Y, Li YY, Nie YQ, Zhou YJ, Cao CY, Xu L. Association between metabolic syndrome and the development of non-alcoholic fatty liver disease. Exp Ther Med. 2013;6:77–84.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. 2.

    Asrih M, Jornayvaz FR. Metabolic syndrome and nonalcoholic fatty liver disease: is insulin resistance the link? Mol Cell Endocrinol. 2015;418(Pt 1):55–65.

    Article  PubMed  CAS  Google Scholar 

  3. 3.

    Neuschwander-Tetri BA, Caldwell SH. Nonalcoholic steatohepatitis: summary of an AASLD Single Topic Conference. Hepatology. 2003;37:1202–19.

    Article  PubMed  Google Scholar 

  4. 4.

    Starley BQ, Calcagno CJ, Harrison SA. Nonalcoholic fatty liver disease and hepatocellular carcinoma: a weighty connection. Hepatology. 2010;51:1820–32.

    Article  PubMed  Google Scholar 

  5. 5.

    Wong RJ, Cheung R, Ahmed A. Nonalcoholic steatohepatitis is the most rapidly growing indication for liver transplantation in patients with hepatocellular carcinoma in the US. Hepatology. 2014;59:2188–95.

    Article  PubMed  Google Scholar 

  6. 6.

    Oliveira CP, Stefano JT. Genetic polymorphisms and oxidative stress in non-alcoholic steatohepatitis (NASH): a mini review. Clin Res Hepatol Gastroenterol. 2015;39(Suppl 1):S35–40.

    Article  PubMed  CAS  Google Scholar 

  7. 7.

    Cencioni C, Spallotta F, Greco S, Martelli F, Zeiher AM, Gaetano C. Epigenetic mechanisms of hyperglycemic memory. Int J Biochem Cell Biol. 2014;51:155–8.

    Article  PubMed  CAS  Google Scholar 

  8. 8.

    Pruis MG, Lendvai A, Bloks VW, Zwier MV, Baller JF, de Bruin A, Groen AK, Plosch T. Maternal western diet primes non-alcoholic fatty liver disease in adult mouse offspring. Acta Physiol (Oxf). 2014;210:215–27.

    Article  CAS  Google Scholar 

  9. 9.

    Chang X, Yan H, Fei J, Jiang M, Zhu H, Lu D, Gao X. Berberine reduces methylation of the MTTP promoter and alleviates fatty liver induced by a high-fat diet in rats. J Lipid Res. 2010;51:2504–15.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. 10.

    Cordero P, Gomez-Uriz AM, Campion J, Milagro FI, Martinez JA. Dietary supplementation with methyl donors reduces fatty liver and modifies the fatty acid synthase DNA methylation profile in rats fed an obesogenic diet. Genes Nutr. 2012;

  11. 11.

    Tian W, Xu H, Fang F, Chen Q, Xu Y, Shen A. Brahma-related gene 1 bridges epigenetic regulation of proinflammatory cytokine production to steatohepatitis in mice. Hepatology. 2013;58:576–88.

    Article  PubMed  CAS  Google Scholar 

  12. 12.

    Lu SC, Alvarez L, Huang ZZ, Chen L, An W, Corrales FJ, Avila MA, Kanel G, Mato JM. Methionine adenosyltransferase 1A knockout mice are predisposed to liver injury and exhibit increased expression of genes involved in proliferation. Proc Natl Acad Sci U S A. 2001;98:5560–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. 13.

    Martinez-Chantar ML, Corrales FJ, Martinez-Cruz LA, Garcia-Trevijano ER, Huang ZZ, Chen L, Kanel G, Avila MA, Mato JM, Lu SC. Spontaneous oxidative stress and liver tumors in mice lacking methionine adenosyltransferase 1A. FASEB J. 2002;16:1292–4.

    Article  PubMed  CAS  Google Scholar 

  14. 14.

    Pirola CJ, Fernandez Gianotti T, Burgueno AL, Rey-Funes M, Loidl CF, Mallardi P, San Martino J, Castano GO, Sookoian S. Epigenetic modification of liver mitochondrial DNA is associated with histological severity of nonalcoholic fatty liver disease. Gut. 2012;

  15. 15.

    Sookoian S, Rosselli MS, Gemma C, Burgueno AL, Fernandez Gianotti T, Castano GO, Pirola CJ. Epigenetic regulation of insulin resistance in nonalcoholic fatty liver disease: impact of liver methylation of the peroxisome proliferator-activated receptor gamma coactivator 1alpha promoter. Hepatology. 2010;52:1992–2000.

    Article  PubMed  CAS  Google Scholar 

  16. 16.

    Ahrens M, Ammerpohl O, von Schonfels W, Kolarova J, Bens S, Itzel T, Teufel A, Herrmann A, Brosch M, Hinrichsen H, et al. DNA methylation analysis in nonalcoholic fatty liver disease suggests distinct disease-specific and remodeling signatures after bariatric surgery. Cell Metab. 2013;18:296–302.

    Article  PubMed  CAS  Google Scholar 

  17. 17.

    de Mello VD, Matte A, Perfilyev A, Mannisto V, Ronn T, Nilsson E, Kakela P, Ling C, Pihlajamaki J. Human liver epigenetic alterations in non-alcoholic steatohepatitis are related to insulin action. Epigenetics. 2017;12:287–95.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Murphy SK, Yang H, Moylan CA, Pang H, Dellinger A, Abdelmalek MF, Garrett ME, Ashley-Koch A, Suzuki A, Tillmann HL, et al. Relationship between methylome and transcriptome in patients with nonalcoholic fatty liver disease. Gastroenterology. 2013;145:1076–87.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. 19.

    Hotta K, Kitamoto T, Kitamoto A, Ogawa Y, Honda Y, Kessoku T, Yoneda M, Imajo K, Tomeno W, Saito S, Nakajima A. Identification of the genomic region under epigenetic regulation during non-alcoholic fatty liver disease progression. Hepatol Res. 2017;

  20. 20.

    Zeybel M, Hardy T, Robinson SM, Fox C, Anstee QM, Ness T, Masson S, Mathers JC, French J, White S, Mann J. Differential DNA methylation of genes involved in fibrosis progression in non-alcoholic fatty liver disease and alcoholic liver disease. Clin Epigenetics. 2015;7:25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. 21.

    Kitamoto T, Kitamoto A, Ogawa Y, Honda Y, Imajo K, Saito S, Yoneda M, Nakamura T, Nakajima A, Hotta K. Targeted-bisulfite sequence analysis of the methylation of CpG islands in genes encoding PNPLA3, SAMM50, and PARVB of patients with non-alcoholic fatty liver disease. J Hepatol. 2015;63:494–502.

    Article  PubMed  CAS  Google Scholar 

  22. 22.

    Nilsson E, Matte A, Perfilyev A, de Mello VD, Kakela P, Pihlajamaki J, Ling C. Epigenetic alterations in human liver from subjects with type 2 diabetes in parallel with reduced folate levels. J Clin Endocrinol Metab. 2015;100:E1491–501.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. 23.

    Leti F, Malenica I, Doshi M, Courtright A, Van Keuren-Jensen K, Legendre C, Still CD, Gerhard GS, DiStefano JK. High-throughput sequencing reveals altered expression of hepatic microRNAs in nonalcoholic fatty liver disease-related fibrosis. Transl Res. 2015;

  24. 24.

    DiStefano JK, Kingsley C, Craig Wood G, Chu X, Argyropoulos G, Still CD, Done SC, Legendre C, Tembe W, Gerhard GS. Genome-wide analysis of hepatic lipid content in extreme obesity. Acta Diabetol. 2015;52:373–82.

    Article  PubMed  CAS  Google Scholar 

  25. 25.

    Gerhard GS, Benotti P, Wood GC, Chu X, Argyropoulos G, Petrick A, Strodel WE, Gabrielsen JD, Ibele A, Still CD, et al. Identification of novel clinical factors associated with hepatic fat accumulation in extreme obesity. J Obes. 2014;2014:368210.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Wood GC, Chu X, Manney C, Strodel W, Petrick A, Gabrielsen J, Seiler J, Carey D, Argyropoulos G, Benotti P, et al. An electronic health record-enabled obesity database. BMC Med Inform Decis Mak. 2012;12:45.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Network TCGAR. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008;455:1061–8.

    Article  CAS  Google Scholar 

  28. 28.

    Barfield RT, Kilaru V, Smith AK, Conneely KN. CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics. 2012;28:1280–1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. 29.

    Warden CD, Lee H, Tompkins JD, Li X, Wang C, Riggs AD, Yu H, Jove R, Yuan YC. COHCAP: an integrative genomic pipeline for single-nucleotide resolution DNA methylation analysis. Nucleic Acids Res. 2013;41:e117.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. 30.

    Storey JD. A direct approach to false discovery rates. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2002;64:479–98.

    Article  Google Scholar 

  31. 31.

    Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. 32.

    Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. 33.

    Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. 34.

    Huebert RC, Vasdev MM, Shergill U, Das A, Huang BQ, Charlton MR, LaRusso NF, Shah VH. Aquaporin-1 facilitates angiogenic invasion in the pathological neovasculature that accompanies cirrhosis. Hepatology. 2010;52:238–48.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. 35.

    Fibbi G, Pucci M, Grappone C, Pellegrini G, Salzano R, Casini A, Milani S, Del Rosso M. Functions of the fibrinolytic system in human Ito cells and its control by basic fibroblast and platelet-derived growth factor. Hepatology. 1999;29:868–78.

    Article  PubMed  CAS  Google Scholar 

  36. 36.

    Jin-no K, Tanimizu M, Hyodo I, Kurimoto F, Yamashita T. Plasma level of basic fibroblast growth factor increases with progression of chronic liver disease. J Gastroenterol. 1997;32:119–21.

    Article  PubMed  CAS  Google Scholar 

  37. 37.

    Yu C, Wang F, Jin C, Huang X, Miller DL, Basilico C, McKeehan WL. Role of fibroblast growth factor type 1 and 2 in carbon tetrachloride-induced hepatic injury and fibrogenesis. Am J Pathol. 2003;163:1653–62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. 38.

    Nakamura I, Zakharia K, Banini BA, Mikhail DS, Kim TH, Yang JD, Moser CD, Shaleh HM, Thornburgh SR, Walters I, Roberts LR. Brivanib attenuates hepatic fibrosis in vivo and stellate cell activation in vitro by inhibition of FGF, VEGF and PDGF signaling. PLoS One. 2014;9:e92273.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. 39.

    Sunny NE, Bril F, Cusi K. Mitochondrial adaptation in nonalcoholic fatty liver disease: novel mechanisms and treatment strategies. Trends Endocrinol Metab. 2017;28:250–60.

    Article  PubMed  CAS  Google Scholar 

  40. 40.

    Torok NJ. Dysregulation of redox pathways in liver fibrosis. Am J Physiol Gastrointest Liver Physiol. 2016;311:G667–74.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Rizzo G, Renga B, Mencarelli A, Pellicciari R, Fiorucci S. Role of FXR in regulating bile acid homeostasis and relevance for human diseases. Curr Drug Targets Immune Endocr Metabol Disord. 2005;5:289–303.

    Article  PubMed  CAS  Google Scholar 

  42. 42.

    Armstrong LE, Guo GL. Role of FXR in liver inflammation during nonalcoholic steatohepatitis. Curr Pharmacol Rep. 2017;3:92–100.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. 43.

    Oseini AM, Sanyal AJ. Therapies in non-alcoholic steatohepatitis (NASH). Liver Int. 2017;37 Suppl 1:97–103.

    Article  PubMed  Google Scholar 

  44. 44.

    Cheung O, Puri P, Eicken C, Contos MJ, Mirshahi F, Maher JW, Kellum JM, Min H, Luketic VA, Sanyal AJ. Nonalcoholic steatohepatitis is associated with altered hepatic microRNA expression. Hepatology. 2008;48:1810–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. 45.

    Estep M, Armistead D, Hossain N, Elarainy H, Goodman Z, Baranova A, Chandhoke V, Younossi ZM. Differential expression of miRNAs in the visceral adipose tissue of patients with non-alcoholic fatty liver disease. Aliment Pharmacol Ther. 2010;32:487–97.

    Article  PubMed  CAS  Google Scholar 

  46. 46.

    Chiang DJ, Pritchard MT, Nagy LE. Obesity, diabetes mellitus, and liver fibrosis. Am J Physiol Gastrointest Liver Physiol. 2011;300:G697–702.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. 47.

    Duval C, Thissen U, Keshtkar S, Accart B, Stienstra R, Boekschoten MV, Roskams T, Kersten S, Muller M. Adipose tissue dysfunction signals progression of hepatic steatosis towards nonalcoholic steatohepatitis in C57BL/6 mice. Diabetes. 2010;59:3181–91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. 48.

    Lomonaco R, Ortiz-Lopez C, Orsak B, Webb A, Hardies J, Darland C, Finch J, Gastaldelli A, Harrison S, Tio F, Cusi K. Effect of adipose tissue insulin resistance on metabolic parameters and liver histology in obese patients with nonalcoholic fatty liver disease. Hepatology. 2012;55:1389–97.

    Article  PubMed  CAS  Google Scholar 

Download references


The authors acknowledge the contribution of the TGen Sequencing Core.


This work was supported by funding from the National Institutes of Health grant DK091601 (JKD and GSG) and the Translational Genomics Research Institute.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information




JKD conceived the concept and design of the study. LL, XC, ATP, and CDS contributed to the acquisition of samples and data of the study. IM and JKD analyzed and interpreted the data of the study. JKD drafted the manuscript. JKD, GSG, IM, XC, ATP, and CDS helped in critical revision of the manuscript for important intellectual content. JKD and GSG obtained funding of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Johanna K. DiStefano.

Ethics declarations

Ethics approval and consent to participate

All study participants provided written informed consent for research, which was conducted according to The Code of Ethics of the World Medical Association (Declaration of Helsinki). The Institutional Review Boards of Geisinger Health System, the Translational Genomics Research Institute, and the Lewis Katz School of Medicine at Temple University approved the research.

Consent for publication

The manuscript does not contain any personal health identifiers.

Competing interests

CDS receives grant and consulting support from Ethicon Endo-Surgery. Other authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. CpG islands differentially methylated between normal and fibrotic samples (N = 208) (DOCX 37 kb)

Additional file 2:

Table S2. Top 100 differentially methylated CpG sites (DOCX 30 kb)

Additional file 3:

Figure S1. Box plot of mean (± standard deviation) methylation β values for CGI associated with gene expression differences in patients with NAFLD fibrosis (n = 11) compared to individuals with normal liver histology (n = 15). (PDF 544 kb)

Additional file 4:

Table S3. Canonical pathways identified by pathway analysis using all significant CpG islands (N = 208). (DOCX 26 kb)

Additional file 5:

Table S4. Validation of previously reported methylation findings. (DOCX 34 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gerhard, G.S., Malenica, I., Llaci, L. et al. Differentially methylated loci in NAFLD cirrhosis are associated with key signaling pathways. Clin Epigenet 10, 93 (2018).

Download citation


  • Nonalcoholic fatty liver disease
  • Nonalcoholic steatohepatitis
  • Liver fibrosis
  • DNA methylation
  • Epigenetics
  • Severe obesity
  • Bariatric surgery