Methylation Profile of Single Hepatocytes Derived from Hepatitis B Virus-Related Hepatocellular Carcinoma

Background With the development of high-throughput screening, a variety of genetic alterations has been found in hepatocellular carcinoma (HCC). Although previous studies on HCC methylation profiles have focused on liver tissue, studies using isolated hepatocytes are rare. The heterogeneity of liver composition may impact the genuine methylation status of HCC; therefore, it is important to clarify the methylation profile of hepatocytes to aid in understanding the process of tumorigenesis. Methods and Findings The global methylation profile of single hepatocytes isolated from liver tissue of hepatitis B virus (HBV) related HCC (HBHC) was analyzed using Illumina Infinium Human Methylation27 BeadChips, and combined bisulfite restriction analysis (COBRA) and bisulfite sequencing were used to validate the 20 significant hypermethylated genes identified. In this study, we found many noteworthy differences in the genome-wide methylation profiles of single hepatocytes of HBHC. Unsupervised hierarchical clustering analysis showed that hepatocyte methylation profiles could be classified according to three cell types: hepatocytes of HCC, adjacent hepatocytes and normal hepatocytes. Among the 20 most hypermethylated genes in the hepatocytes of HBHC, 7 novel genes (WNK2, EMILIN2, TLX3, TM6SF1, TRIM58, HIST1H4Fand GRASP) were found to be hypermethylated in HBHC and hypomethylated in paired adjacent liver tissues; these findings have not been reported in previous studies on tissue samples. Conclusion The genome-wide methylation profile of purified single hepatocytes of HBHC was aided in understanding the process of tumorigenesis, and a series of novel methylated genes found in this study have the potential to be biomarkers for the diagnosis and prognosis of HBHC.


Introduction
Hepatocellular carcinoma (HCC) is one of the most common cancers in the world and is commonly developed from liver cirrhosis that are secondary to viral infection [1]. China is one of the highest prevalent areas of HCC as chronic hepatitis B carriers account for more than 10% of its population. The prognosis of HCC patients is poor even after partial hepatectomy due to the high rates of metastasis and relapse [2,3], and there is an urgent need to elucidate the mechanism of hepatocarcinogenesis.
Epigenetic changes have been reported recently as one of the mechanisms of tumorigenesis, and DNA methylation, which is one of the major epigenetic alterations, has also been reported in HCC. Aberrant DNA methylation of promoter CpG islands has been associated with global hypomethylation and specific loci hypermethylation, which are believed to have potential as diagnostic markers in the progression of malignant tumors. In HCC, DNA hypermethylation of candidate genes, including GSTP1 [4], RASSF1A [5], RIZ1, APC [6], SOCS1 [7], APC [8], and E-cadherin [9], has been examined using methylation sensitive polymerase chain reaction, combined bisulfite restriction analysis (COBRA), and bisulfite sequencing techniques in previous studies. With the development of genome-wide screening, the global methylation profiling of cancer has become clearer than ever before. A group of potential tumor biomarker genes has been reported in HCC [10]. Different gene methylation profiles have been found in HCCs with different etiological backgrounds [11,12,13], pathophysiological processes [10], and histological characteristics.
However, most of the large-scale screening profiles have been performed on tissue samples. Tissues often consist of a heterogeneous population of cells. The same is also true in cancerous tissues. In addition to tumor cells, HCC tissues can have fibroblasts, stromal cells, vasculature, inflammatory cells and other non-parenchymal cells [14]. Although these non-parenchymal cells have some effect on the progression of cancer, this mixed population may distort or even hide the genuine methylation status of cancer cells, making diagnosis difficult to determine and experimental results hard to interpret.
Several studies from independent labs have reported gene expression and proteomic profile changes in isolated cells compared to primary tissues. Harrell et al. compared the gene expression of whole organs to isolated cells from the organ and showed a distinctly different gene expression profile, with only minor overlap between the two samples [15]. Waanders et al. reported that the proteomic analysis of single pancreatic islets was different from that of pancreatic tissue, suggesting that single-cell analysis is more readily accessible to proteomic measurements [16]. Even the tumor epithelium and tumor-associated stroma have shown different methylation profiles, as GSTP1 and RARbeta2 were found to be hypermethylated in the tumor epithelium and less so in tumor-associated stroma [17]. In the context of the liver, hepatocytes from HCC have shown reduced heterogeneity and contamination in tissue samples and present less protein spots than HCC tissue [18]. Considering the difference between isolated hepatocytes and liver tissue, the gene expression and proteomic profiles of hepatocytes isolated from cirrhotic nodules [19], fetal livers [20], HCC lymphatic metastases [21] and neoplasms [18] have been studied. The analysis of these hepatocytes has revealed not only specific patterns that may reflect the prognosis and drug sensitivity of cancer cell [22,23] but also the identity of genes involved in malignant transformation, progression [24], and tumor metastasis [21]. The homogeneity of isolated hepatocytes may allow the elucidation of the mechanism of tumorigenesis in HCC, with the global genome-wide methylation profile of HCC cells being of particular interest.
In the present study, we used a four-step collagenase perfusion method to separate hepatocytes from liver tissues and analyzed the aberrant promoter methylation of three paired liver cancer samples and adjacent liver samples using Illumina Infinium Human Methylation27 BeadChips that contained 27,000 CpGs. And to eliminate the effect of different clinical backgrounds, we focused on HBV related HCC (HBHC). The CpG island methylation profile of hepatocytes isolated from primary tissues and distinct pathways in single hepatocytes from HBHC were obtained and revealed a group of potential novel biomarkers in HBHC.

Isolation and characterization of separated cells from liver tissues
For clarification, in this study, the hepatocytes separated from live cancer tissues were called HHCs, while HAs were those that were collected from adjacent non-cancerous liver tissues and HNs were those that were isolated from normal liver tissues. The isolated hepatocyte motility rate was similar in the three groups (95%64%) using a four-step collagenase perfusion method; therefore, a cell motility rate .90% was chosen for further analysis. The total number of separated hepatocytes was 1.3860.386106 per gram of hepatic tissue, and the contamination rate of non-parenchymal cells in the HHC population was ,2%, with a small cellular fraction as determined using microscopic analysis.
Genome-wide DNA methylation profiling identifies methylated genes in freshly isolated hepatocytes and liver tissues Global DNA methylation profiles were measured using Illumina Infinium Human Methylation27 BeadChips, which target 14,475 total refseq genes, 12,833 well-annotated genes described in the NCBI CCDS database, 144 methylation hotspots in cancer genes, 982 cancer-related targets and 110 miRNA promoters. Three pairs of HHCs and HAs were included for comparison using the HNs as the control. Using this approach, the average methylation rate of HHCs was higher than that of matched HAs (Figure 1) as reported in previous study [25].
To distinguish the genes that were differentially methylated in HHCs, methylated genes were filtered according to P,0.05 and a delta b.0.15. HHCs showed 792 genes that were hypermethylated and 533 genes that were hypomethylated compared with HNs, and 345 and 796 genes were hypermethylated and hypomethylated, when compared with HAs respectively. When compared to both HAs and HNs combined, HHCs showed 315 and 460 genes that were hypermethylated and hypomethylated, respectively. KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses were performed to identify the pathways with which the aberrant genes in HHCs were associated. Hypermethylated genes in HHCs were enriched in the KEGG pathway, most of which have been reported to be involved in the following cancer pathways: gap junction (p = 9.1610-5), calcium signaling (p = 5.8610-4), natural killer cell-mediated cytotoxicity (p = 3. We further restricted inclusion criteria as P,0.001 and the delta b.0.15, which resulted in the identification of 63 CpG sites corresponding to 62 genes (Table S1); the average distance from the 63 CpG sites to the transcription start site (TSS) was 256.58 bp. A total of 50 CpG sites showed significantly higher methylation in HHCs relative to the other two groups of hepatocytes. While 10 CpG sites outside of the CpG islands showed less methylation in HHCs, there were only 2 CpG sites inside CpG islands that showed hypomethylation. A total of 45 CpG sites exhibited a significant monotonic decreasing trend in the methylation rate (HHCs.HAs.HNs), while 8 CpG sites showed increasing trends (HHCs,HAs,HNs). Among the 63 genes showing the most significant differences in methylation, 38 genes have been associated with cancer in previous studies. Among them, 10 genes have been reported in HCC, and ZMYND10 [26], SYK [27], DAB2IP [28], and DLEC1 [29] have been reported to be tumor-suppressor genes in HCC. Unsupervised hierarchical clustering analysis with a P,0.001 showed that in the single hepatocytes group, the cells could be classified into three types: HN, HA, HHC ( Figure 3).
The methylation profile in HAs was also compared with that of HNs, and 124 and 115 genes were identified as differently methylated (p,0.05) and hypermethylated, respectively. A number of these genes were associated with MAPK (FGF2, BDNF, TGFB3, FGF23, and FGF17) and colorectal cancer (dcc, cycs, and tgfb3) signaling pathway.

Methylation analysis in paired normal and liver cancer tissues
To further investigate whether the methylation profiles of isolated hepatocytes were consistent with tissue samples, we chose the 20 most significantly hypermethylated genes in five pairs of HBHC tissues and adjacent liver tissues for further validation using COBRA. Two previously reported genes (ZMYND10 and SYK) and one unknown gene (DKFZp434I1020) were excluded from the analysis. The PCR products of these seventeen genes amplified from bisulfite-treated cancer and adjacent tissues were digested using the methylation-sensitive restriction enzyme BstUI. Two of the genes (KIF17 and ATP8A2) showed no methylation in either cancer or matched tissues, and eight genes (EPHA4, OVOL1, UTF1, ADAM8, BOLL, FOXE3, GU-CY1A2, and ADCY5) were equally methylated in the two groups. These ten genes were removed from further analysis, and the remaining seven genes were analyzed in ten paired liver tissues. These genes were frequently methylated in HBHC, while those in the matched non-cancerous liver tissues were  Neuroactive ligandreceptor interaction   15  GRP83, CRHR2, ADRA2B, DAB2IP, ADRB3, DRD4,  EDNRB, F2R, GALR1, NPBWR1, HTR7, NMBR,   either methylated to a significantly lower extent or not methylated at all. A summary of the promoter methylation analysis is shown in Tables 3 and 4 and the cloning and sequencing analysis results are shown in Figure 4. These results indicate that cancer tissues have much higher methylation rates than the related non-cancerous tissues, which confirm the results from COBRA.

Expression analysis in cell lines with 5-aza-dc
To confirm the role of methylation in gene expression silencing, we treated four hepatoma cell lines with the demethylating agent 5-aza-dc. Seven genes that were validated in ten paired HBHC and adjacent liver tissues using COBRA were selected for this experiment, and their expression was analyzed using real-time quantitative polymerase chain reaction (RT-PCR). All of the   Table 5.

Discussion
Carcinogenesis is a multi-stage, multi-mechanistic process and the details of the progression of a normal cell as it transforms into a cancer that exhibits excessive proliferation and reduced apoptotic potential are still unclear. The global methylation profile helps us to clarify the progression of hepatocarcinogenesis. In this study, we analyzed the methylation profiles of single hepatocyte isolated from HBHC using Illumina Infinium Human Methylation27 BeadChips to better understand the molecular mechanisms underlying tumorigenesis.
The collagenase perfusion method is a classic method to isolate hepatocytes from liver tissues. This method were used to isolated hepatocyte from normal liver [30,31], liver cirrhosis [32] and HCC [20,22,23,33] to study molecular mechanisms. Fibrous liver connective tissue was removed through successive perfusion with EDTA, dispase, collagenase type IV, and DNase I. Centrifuged at low speed to eliminate non-parenchymal cells. Single isolated hepatocyte is not only different from liver tissue without the latter's heterogeneity which might influence the methylation profile [13], but is also different from cell lines. Because cell lines were cultured in vitro and did not have the background of chronic viral infection, hypermethylation in cell lines could be a background event caused by repeated passages and the culture environment, which could be significantly different from primary malignancies [34]. Therefore, we focused our study on freshly isolated hepatocytes, which not only exhibit the pathology background but are also homogeneous.
Using the isolated single hepatocytes from HBHC, we detected novel genes that were specifically methylated. However, among these genes, only a few overlapped with previous studies. Beside the differences between single hepatocytes and liver tissue, one explanation for these differences is that the studies varied in pathogen, pathology, tumor grade, race and other clinical data. In this study, we selected all samples according to the same criteria to eliminate the effect of different clinical backgrounds. In addition, all HCC specimens were moderately or poor differentiated, and all patients suffered from hepatitis B virus (HBV) and exhibited the same underlying cirrhosis etiology. Another explanation for the variability of the findings from the multiple studies could be that the methylation microarrays from each company detected methylated genes in the human genome using different detection methods with variations in the sensitivity of detection and filter range.
In this study, the HHC methylation profile detected genes that were linked with several cancer-associated pathways. Gap junction genes (ADCY5, DAB2IP, ADCY4 CUCY1A2, CX36, and TUBB6) were significantly hypermethylated in HHCs (p = 0.000091) and have shown significantly reduced expression in HCC compared with normal livers in previous studies [35] [31]. Cancer was one of the first pathologies to be associated with gap junction defects, and tumor formation, invasiveness and cancer cell dissemination have all been closely correlated with gap junctions [36]. The inhibition of gap junctional intercellular communication (GJIC) results in a block in cell communication and the loss of the control of endocrine regulation between normal cells to regulate growth and terminally differentiate, thus resulting in immortalization and escape from apoptosis [37,38]. GJIC is postulated to be an important mechanism to control growth and differentiation and eliminate the interference from non-parenchymal liver cells, and it was significantly associated with isolated hepatocyte methylation profiles.
Calcium signaling, cell adhesion molecules and apoptosis are all closely related to gap junctions. Calcium signaling represents a general mechanism of intercellular signaling via the gap junctionmediated diffusion of Ca2+/inositol 1,4,5-triphosphate [37]; cancer cell transformation, cell proliferation and cell death are associated with a major rearrangement of Ca2+ pumps, Na/Ca exchangers and Ca2+ channels, which leads to enhanced proliferation and an impaired ability to die that ultimately determines cell fate [38]. Cell adhesion molecules play an important role in regulating cell survival, proliferation and motility   [39], and apoptosis is an important link between inflammation and liver cancer. Inflammatory hepatitis might activate NF-kB by the concerted action of cytokines that are resistant to apoptosis and promote cell survival of pre-cancerous hepatocytes [27]. Another pathway that is associated with cancer is natural killer cellmediated cytotoxicity. Nature killer cell(NK cell) play a role in mediating innate host immune responses, which are regulated by a balance between activating and inhibitory receptors. The dysregulation of natural killer cell-mediated cytotoxicity breaks this balance because upregulation of the ligands for inhibitory receptors and/or the loss of ligands for activating receptors aids in tumor cell escape from NK cell recognition [40]. Except the distinct pathways identified in HHCs, a group of aberrantly methylated genes was selected. EPHA4, KIF17, OVOL1, UTF1, ADAM8 BOLL, FOXE3, GUCY1A2, and ADCY5 exhibited the same promoter methylation status in Figure 4. DNA methylation status of seven selected genes in HBHC and adjacent liver tissue. COBRA analysis were performed as amplified PCR of sodium bisulfite-converted genomic DNAs, PCR products incubated with BstUI to measure the methylation rate and samples without incubation with BstUI were listed as control. Selected PCR products were cloned into T-A clones and sequencing. Five samples of each alleles listed were shown by black and white circles represent methylated and unmethylated CpG dinucleotides respectively. "MI" was the abbreviation of methylation index, calculated as the percentage of the methylated dinucleotides versus total. Detail methylated sites were marked in red, the position of COBRA primers is also shown in blue, covering the region was respect to the transcription start site. AJ was the abbreviation of adjacent liver tissue. doi:10.1371/journal.pone.0019862.g004 HBHC and adjacent liver tissue, suggesting that the HHC methylation is acquired during cirrhosis and chronic inflammation with persistent viral infection during the pre-cancerous state before the formation of the tumor. SYK and ZMYND10 (also known as BLU) were hypermethylated in HHCs in our studies. These genes have been reported to be methylated in HCC, indicative of tumorsuppressor genes. SYK has been highly correlated with lower overall survival rates. [26,27]. WNK2, emilin2, tlx3, tm6sf1, trim58, hist1h4f, and grasp were significantly hypermethylated in HBHC. Wnk2 belongs to a serine/threonine kinase subfamily that has the catalytic lysine in a non-canonical position. WNK2 point mutations have been found in lung cancer, adult glial tumors and other solid tumor types, and it has also been reported to be a tumor suppressor [41]. EMILIN2 belongs to a family of extracellular matrix glycoproteins that contain an EMI domain at the N-terminus and a C1q domain at the C-terminus. EMILIN2 has been shown to suppress the growth of breast and other cancer cells and trigger apoptosis in cancer cells via the extrinsic apoptotic pathway following EMILIN2 binding to the trail receptors DR4 and DR5 [42]. TLX3 has been shown to be hypermethylated in prostate cancer [43]. Tm6sf1, trim58, hist1h4f, and grasp were identified as novel hypermethylated genes in HBHC; these genes were methylated in HHCs and HBHC tissues, suggesting that the methylation of these genes might have critical role in tumorigenesis and cancer progression. These novel genes were upregulated in response to 5-aza-dc in at least one of the tested cell lines, and based on the different invasive behaviors of these cell lines, the low expression was regulated by hypermethylation.
Another observation is the distinct methylation profile of isolated hepatocytes from non-cancerous livers, which was different from that of normal hepatocytes, as reported in previous studies with tissue samples [6]; some pathways associated with cancer were also found in adjacent cells. According to a recent theory about tumorigenesis by Slaughter [44], field cancerization results from the expansion of pre-neoplastic cells in early tumorigenesis. Aberrantly methylated genes in adjacent cells would affect the progress of hepatocarcinogenesis and would be useful in early detection and risk assessment.
This study is the first report on the global methylation profiling of hepatocytes isolated from tissues. Even with the small sample size, the aberrant methylation profile of HHCs is distinct from that of normal hepatocytes and reveals potential biomarkers in the diagnosis and prognosis of HBHC. These novel methylated genes should be further analyzed in larger clinical samples to elucidate their clinical implication, and more samples with different etiological backgrounds, pathophysiological processes, and histological characteristics should be used in further studies.

Patient specimens
All samples in this study were taken from a surgical specimen, obtained through hepatectomy or liver transplantation. All the patients have signed an informed-consent. No prisoners or organs from prisoners were used in this study. All patients included in this Figure 5. Gene re-expression after treated with 5-aza-dc. Gene expression level was analyzed by quantitative reverse transcriptionpolymerase chain reaction. All selected genes were unregulated after treated with 10 uM 5-aza-dcC at least in one cell lines, the expression of cell lines before 5-aza-dc treatment was analyzed as control. doi:10.1371/journal.pone.0019862.g005 A total of 13 patients with HCC and 2 patients with hepatic hemangioma were selected for analysis. All patients were male to avoid sex-based variation, and the mean age was 46.84 years old (5865.99). All HCC patients were hepatitis B virus (HBV) surface antigen-positive without hepatitis C virus (HCV) infection and exhibited the same underlying cirrhosis etiology. Clinical information for these patients were summarized in Tables 6 and 7. HBV surface antigen and anti-HCV antibody were measured serologically. The collection of tumor tissues was strictly limited inside the boundary of tumor, and adjacent non-cancerous liver tissues was collected at least 3 cm from the tumor margins, all the collections were under the supervision of a same pathologist. Each sample was devided into two parts, one was used for methylation profiling analysis, the other was histologically examined by two experienced pathologists, confirming that under the microscope there was no detectable normal liver cells inside the tumor samples and no tumor cells in the surrounding liver tissues.
Three paired HBHC and adjacent liver tissues were perfused to harvest single hepatocytes from tissues, and were then used for the Illumina Infinium BeadChips analysis. Ten paired HBHC and matching tissues were available for further validation. Two normal adult liver tissues from the patients with hepatic hemangioma who underwent partial hepatectomy without HBV or HCV infection were used for comparison purposes. Tissue samples were freshly frozen at 280uC or perfused immediately to harvest single cells.

Isolation of hepatocytes from liver tissues
Three paired tumors and adjacent liver tissues and normal livers were isolated using a modified four-step collagenase perfusion method as described previously [30,31]. The liver segments were perfused by the extracorporeal circulating perfusion apparatus in a sterile cell room. Four different buffer solutions supplemented, respectively, with 0.37 mg/ml EDTA, 0.5% dispase, 0.05% collagenase type IV (Gibco BRL, Grand Island, NY), and 40 ug/ml DNase I (Sigma, St. Louis, MO) were perfused successively into small hepatic vessels for 40 minutes. Solid tumor remnants were dispersed mechanically with a scalpel blade, and then incubated in collagenase buffer for 20 min at 37uC. The pooled suspensions were filtered by 100 um nylon mesh. Isolated cells were harvested and centrifuged at 4uC at 50 g for 3 minutes, and then the liver cells were resuspended in Dulbecco's modified Eagle's medium (DMEM; Invitrogen, Carlsbad, CA) and centrifuged twice at 4uC at 50 g for 1 min.
The hepatocytes were then stained using 0.4% trypan blue, and cell motility and the number of isolated cells were calculated. DNA was extracted from isolated hepatocytes and frozen liver tissues for methylation analysis using the QIAamp DNA Mini Kit (Qiagen, Valencia, Ca) according to the manufacturer's instructions.

Illumina Infinium methylation assay
The Illumina Infinium 27k Human DNA methylation Bead-Chip v1.2 was used to obtain genome-wide DNA methylation profiles across approximately 27,000 CpGs, which were focused on the promoter regions of 14,495 genes, in purified cell samples.
Bisulfite conversion of genomic DNA was performed using the EZ DNA methylation Kit (Zymo Research, D5002, USA) according to the protocol provided by the Illumina Infinium Methylation Assay. Bisulfite-converted and -unconverted (i.e., methylated) sites were simultaneously evaluated using DNA hybridization to site-specific probes attached to beads (one set for unmethylated and the other for methylated sites), followed by allele-specific base extension that included a fluorescent label. Two different labels were used, and fluorescent signals were specific for either the T (unmethylated) or C (methylated) alleles. Methylation scores represented as b values were generated for each site using  IlluminaHGenome Studio software v1.0 (Illumina Inc., USA) and were computed based on the ratio of methylated to methylated plus unmethylated signal outputs. Thus, the b values ranged from 0 (unmethylated) to 1 (fully methylated) on a continuous scale. Initial array results were visualized using IlluminaHGenome Studio software v1.0 (Illumina Inc., USA). All computations and statistical analyses were performed using the R package (R v2.10.1; http://www.r-project.org) and Bioconductor [45]. Significant differences between the sample groups were identified using the limma (Linear Models for Microarray Analysis) package of the Bioconductor suite, and an empirical Bayesian-moderated t-statistic hypothesis test between the two specified phenotypic groups was performed [46] using significance at P value ,0.05 with the absolute delta b value cutoff set at 0.15.
Heat maps were created using Mev software [47]. The Euclidean distance between two groups of samples was calculated by the average linkage measure (the mean of all pair-wise distances (linkages) between members of the two groups concerned). Gene annotation and enrichment analyses were performed using KEGG databases using the DAVID Bioinformatics Resources (2008; http://david.abcc.ncifcrf.gov/) [48] interfaces and WebGestalt (http://bioinfo.vanderbilt.edu/webgestalt/) [49], respectively.

Combined bisulfite restriction analysis (COBRA) and bisulfite sequencing
Genomic DNA (1 mg) from freshly purified hepatocytes and liver tissues was bisulfite-converted using EpiTect Bisulfite Kits (Qiagen, Valencia, CA, USA) and used for PCR amplification. PCR primers were designed using Methyl Primer Express to amplify regions of CpG islands that overlapped or were close to the transcription start sites. All primer sequences and the mythylated sites are listed in Table 8. The PCR products were then digested with the methylation-sensitive endonuclease BstUI at 60uC for 2 hours. Digested PCR fragments were visualized using a 2% agarose gel, and the PCR products for bisulfite sequencing were cloned into the pGEM-T easy vector (Promage, USA) following the manufacturer's instructions. After being transfected into Escherichia coli, up to 10 clones were selected and sequenced using M13 reverse primers. The methylation status of each allele was determined using the methylation index (MI), which was calculated as the number of methylated CpG dinucleotides divided by the total number of CpG dinucleotides sequenced multiplied by 100.

Cell line culture conditions and 5-aza-dc treatment
Four hepatoma cell lines, HepG2, 7402, 97L and 97H, which were all routinely cultured in our institution, were grown in DMEM supplemented with 10% FBS at 37uC with humidified 5% CO2. Cells were seeded in six-well plates 24 hours before treatment with 10 mM 5-aza-29-deoxycytidine (5-aza-dc; Sigma, USA) or phosphate-buffered saline as a control. Media and drugs were replaced every 24 h and the cells were harvested at 72 h post-treatment.

Real-time quantitative polymerase chain reaction (RT-PCR) Analysis
Total RNA from the four hepatoma cell lines was extracted before and after treatment with 5-aza-dc using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's protocols. RNA (1 mg) was reverse-transcribed using the AMV Reverse Transcription System(Invitrogen, Carlsbad, CA, USA). Primer sequences are listed in Table 9. RT-PCR was performed with SYBR Premix Ex TaqTMH (Takara, Japan) on an ABI prism 7500 sequencer. The relative expression levels were analyzed using the DDCt method, and GAPDH was used as an internal control.

Statistical analysis
Data were evaluated using Student's t-test for each pair using SPSS software v16.0. The level of significance for all statistical analyses was P,0.05.

Supporting Information
Table S1 List of differentially methylated genes in single hepatocyte derived from HBHC (P,0.001).