library_construction_protocol: Liver was flushed with physiological salt solution, removed and snap-frozen in liquid nitrogen and stored at −80 °C. Total DNA from liver tissue samples was extracted using Machery Nagel Nucleospin tissue kit and used for whole genome bisulfite sequencing (WGBS) by NXTGNT (www.nxtgnt.com). Concentration of the extracted DNA was measured using Quant-iT PicoGreen kit (Invitrogen). DNA quality was checked on a 1% agarose gel (E-gel EX Invitrogen). 500ng DNA was used from each sample for further analysis. Fragmentation was performed using Covaris model S2 to obtain fragments with a size of approximately 400 bp, followed by bisulfite conversion with the EZ DNA methylation Gold kit (Zymo Research, Irvine, CA, USA) according to manufacturer's protocol. For library preparation, the NEBNext Ultra II DNA library prep kit (New England Biolabs, Ipswich, MA, USA) was used according to protocol, followed by sequencing on three Hiseq3000 lanes (Illumina), generating PE2x150bp reads.