Epigenetic memory of environmental exposures as a mediator of liver disease.

Epigenetic changes are a common feature of human disease, including liver disease and its progression to liver cancer. The most frequent form of liver cancer, HCC, is unusual because most of its causes, or etiologic drivers, are known and are dominated by environmental exposures, including viral infection, alcohol abuse, and overnutrition/metabolic syndrome. The epigenome is a regulatory system overlayed on the genetic material that regulates when, where, and to what extent genes are expressed in developmental, cell type, and disease-associated contexts. Deregulation of the epigenome has emerged as a major player in the pathologic effects of liver disease driving exposures, particularly during their early phases when genetic changes are uncommon. Although it is inherent in the definition of an epigenetic process to be reversible, emerging evidence indicates that epigenetic changes persist after the removal of the exposure and contribute to long-term risk of disease progression. In other systems, environmental exposures lead to beneficial adaptive changes in expression that facilitate processes such as wound healing, and these too are driven by epigenetic changes. What remains unclear, however, is what drives the transition from a beneficial epigenetic memory to a maladaptive scar, the epigenetic processes involved in forming these memories, and whether this process can be modulated for therapeutic purposes. In this review, we discuss these concepts in relation to liver disease and more broadly using examples from other tissue types and diseases, and finally consider how epigenetic therapies could be used to reprogram maladaptive epigenetic memories to delay and/or prevent hepatocarcinogenesis.

Titration-based normalization of antibody amount improves consistency of ChIP-seq experiments.

Chromatin immunoprecipitation (ChIP) is an antibody-based approach that is frequently utilized in chromatin biology and epigenetics. The challenge in experimental variability by unpredictable nature of usable input amounts from samples and undefined antibody titer in ChIP reaction still remains to be addressed. Here, we introduce a simple and quick method to quantify chromatin inputs and demonstrate its utility for normalizing antibody amounts to the optimal titer in individual ChIP reactions. For a proof of concept, we utilized ChIP-seq validated antibodies against the key enhancer mark, acetylation of histone H3 on lysine 27 (H3K27ac), in the experiments. The results indicate that the titration-based normalization of antibody amounts improves assay outcomes including the consistency among samples both within and across experiments for a broad range of input amounts.

Interferon drives HCV scarring of the epigenome and creates targetable vulnerabilities following viral clearance.

BACKGROUND AND AIMS: Chronic HCV infection is a leading etiologic driver of cirrhosis and ultimately HCC. Of the approximately 71 million individuals chronically infected with HCV, 10%-20% are expected to develop severe liver complications in their lifetime. Epigenetic mechanisms including DNA methylation and histone modifications become profoundly disrupted in disease processes including liver disease. APPROACH AND RESULTS: To understand how HCV infection influences the epigenome and whether these events remain as "scars" following cure of chronic HCV infection, we mapped genome-wide DNA methylation, four key regulatory histone modifications (H3K4me3, H3K4me1, H3K27ac, and H3K27me3), and open chromatin in parental and HCV-infected immortalized hepatocytes and the Huh7.5 HCC cell line, along with DNA methylation and gene-expression analyses following elimination of HCV in these models through treatment with interferon-α (IFN-α) or a direct-acting antiviral (DAA). Our data demonstrate that HCV infection profoundly affects the epigenome (particularly enhancers); HCV shares epigenetic targets with interferon-α targets; and an overwhelming majority of epigenetic changes induced by HCV remain as "scars" on the epigenome following viral cure. Similar findings are observed in primary human patient samples cured of chronic HCV infection. Supplementation of IFN-α/DAA antiviral regimens with DNA methyltransferase inhibitor 5-aza-2'-deoxycytidine synergizes in reverting aberrant DNA methylation induced by HCV. Finally, both HCV-infected and cured cells displayed a blunted immune response, demonstrating a functional effect of epigenetic scarring. CONCLUSIONS: Integration of epigenetic and transcriptional data elucidate key gene deregulation events driven by HCV infection and how this may underpin the long-term elevated risk for HCC in patients cured of HCV due to epigenome scarring.

CpGtools: a python package for DNA methylation analysis.

MOTIVATION: DNA methylation can be measured at the single CpG level using sodium bisulfite conversion of genomic DNA followed by sequencing or array hybridization. Many analytic tools have been developed, yet there is still a high demand for a comprehensive and multifaceted tool suite to analyze, annotate, QC and visualize the DNA methylation data. RESULTS: We developed the CpGtools package to analyze DNA methylation data generated from bisulfite sequencing or Illumina methylation arrays. The CpGtools package consists of three types of modules: (i) 'CpG position modules' focus on analyzing the genomic positions of CpGs, including associating other genomic and epigenomic features to a given list of CpGs and generating the DNA motif logo enriched in the genomic contexts of a given list of CpGs; (ii) 'CpG signal modules' are designed to analyze DNA methylation values, such as performing the PCA or t-SNE analyses, using Bayesian Gaussian mixture modeling to classify CpG sites into fully methylated, partially methylated and unmethylated groups, profiling the average DNA methylation level over user-specified genomics regions and generating the bean/violin plots and (iii) 'differential CpG analysis modules' focus on identifying differentially methylated CpGs between groups using different statistical methods including Fisher's Exact Test, Student's t-test, ANOVA, non-parametric tests, linear regression, logistic regression, beta-binomial regression and Bayesian estimation. AVAILABILITY AND IMPLEMENTATION: CpGtools is written in Python under the open-source GPL license. The source code and documentation are freely available at https://github.com/liguowang/cpgtools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integrating the Epigenome to Identify Drivers of Hepatocellular Carcinoma.

Disruption of epigenetic mechanisms has been intimately linked to the etiology of human cancer. Understanding how these epigenetic mechanisms (including DNA methylation [5mC], hydroxymethylation [5hmC], and histone post-translational modifications) work in concert to drive cancer initiation and progression remains unknown. Hepatocellular carcinoma (HCC) is increasing in frequency in Western countries but lacks efficacious treatments. The epigenome of HCC remains understudied. To better understand the epigenetic underpinnings of HCC, we performed a genome-wide assessment of 5mC, 5hmC, four histone modifications linked to promoter/enhancer function (H3K4me1, H3K27ac, H3K4me3, and H3K27me3), and transcription across normal, cirrhotic, and HCC liver tissue. Implementation of bioinformatic strategies integrated these epigenetic marks with each other and with transcription to provide a comprehensive epigenetic profile of how and when the liver epigenome is perturbed during progression to HCC. Our data demonstrate significant deregulation of epigenetic regulators combined with disruptions in the epigenome hallmarked by profound loss of 5hmC, locus-specific gains in 5mC and 5hmC, and markedly altered histone modification profiles, particularly remodeling of enhancers. Data integration demonstrates that these marks collaborate to influence transcription (e.g., hyper-5hmC in HCC-gained active enhancers is linked to elevated expression) of genes regulating HCC proliferation. Two such putative epigenetic driver loci identified through our integrative approach, COMT and FMO3, increase apoptosis and decrease cell viability in liver-derived cancer cell lines when ectopically re-expressed. Conclusion: Altogether, integration of multiple epigenetic parameters is a powerful tool for identifying epigenetically regulated drivers of HCC and elucidating how epigenome deregulation contributes to liver disease and HCC.

Alpha-1 Antitrypsin Deficiency Liver Disease, Mutational Homogeneity Modulated by Epigenetic Heterogeneity With Links to Obesity.

Alpha-1 antitrypsin deficiency (AATD) liver disease is characterized by marked heterogeneity in presentation and progression, despite a common underlying gene mutation, strongly suggesting the involvement of other genetic and/or epigenetic modifiers. Variation in clinical phenotype has added to the challenge of detection, diagnosis, and testing of new therapies in patients with AATD. We examined the contribution of DNA methylation (5-methylcytosine [5mC]) to AATD liver disease heterogeneity because 5mC responds to environmental and genetic cues and its deregulation is a major driver of liver disease. Using liver biopsies from adults with early-stage AATD and the ZZ genotype, genome-wide 5mC patterns were interrogated. We compared DNA methylation among patients with early AATD, and among patients with normal liver, cirrhosis, and hepatocellular carcinoma derived from multiple etiologic exposures, and linked patient clinical/demographic features. Global analysis revealed significant genomic hypomethylation in AATD liver-impacting genes related to liver cancer, cell cycle, and fibrosis, as well as key regulatory molecules influencing growth, migration, and immune function. Further analysis indicated that 5mC changes are localized, with hypermethylation occurring within a background of genome-wide 5mC loss and with patients with AATD manifesting distinct epigenetic landscapes despite their mutational homogeneity. By integrating clinical data with 5mC landscapes, we observed that CpGs differentially methylated among patients with AATD disease are linked to hallmark clinical features of AATD (e.g., hepatocyte degeneration and polymer accumulation) and further reveal links to well-known sex-specific effects of liver disease progression. Conclusion: Our data reveal molecular epigenetic signatures within this mutationally homogeneous group that point to ways to stratify patients for liver disease risk.

ZBTB24 is a transcriptional regulator that coordinates with DNMT3B to control DNA methylation.

The interplay between transcription factors and epigenetic writers like the DNA methyltransferases (DNMTs), and the role of this interplay in gene expression, is being increasingly appreciated. ZBTB24, a poorly characterized zinc-finger protein, or the de novo methyltransferase DNMT3B, when mutated, cause Immunodeficiency, Centromere Instability, and Facial anomalies (ICF) syndrome, suggesting an underlying mechanistic link. Chromatin immunoprecipitation coupled with loss-of-function approaches in model systems revealed common loci bound by ZBTB24 and DNMT3B, where they function to regulate gene body methylation. Genes coordinately regulated by ZBTB24 and DNMT3B are enriched for molecular mechanisms essential for cellular homeostasis, highlighting the importance of the ZBTB24-DNMT3B interplay in maintaining epigenetic patterns required for normal cellular function. We identify a ZBTB24 DNA binding motif, which is contained within the promoters of most of its transcriptional targets, including CDCA7, AXIN2, and OSTC. Direct binding of ZBTB24 at the promoters of these genes targets them for transcriptional activation. ZBTB24 binding at the promoters of RNF169 and CAMKMT, however, targets them for transcriptional repression. The involvement of ZBTB24 targets in diverse cellular programs, including the VDR/RXR and interferon regulatory pathways, suggest that ZBTB24's role as a transcriptional regulator is not restricted to immune cells.

Loss of SETD2 Induces a Metabolic Switch in Renal Cell Carcinoma Cell Lines toward Enhanced Oxidative Phosphorylation.

SETD2, a histone H3 lysine trimethyltransferase, is frequently inactivated and associated with recurrence of clear cell renal cell carcinoma (ccRCC). However, the impact of SETD2 loss on metabolic alterations in ccRCC is still unclear. In this study, SETD2 null isogenic 38E/38F clones derived from 786-O cells were generated by zinc finger nucleases, and subsequent metabolic, genomic, and cellular phenotypic changes were analyzed by targeted metabolomics, RNA sequencing, and biological methods, respectively. Our results showed that compared with parental 786-O cells, 38E/38F cells had elevated levels of MTT/Alamar blue levels, ATP, glycolytic/mitochondrial respiratory capacity, citrate synthase (CS) activity, and TCA metabolites such as aspartate, malate, succinate, fumarate, and α-ketoglutarate. The 38E/38F cells also utilized alternative sources beyond pyruvate to generate acetyl-CoA for the TCA cycle. Moreover, 38E/38F cells showed disturbed gene networks mainly related to mitochondrial metabolism and the oxidation of fatty acids and glucose, which was associated with increased PGC1α, mitochondrial mass, and cellular size/complexity. Our results indicate that SETD2 deficiency induces a metabolic switch toward enhanced oxidative phosphorylation in ccRCC, which can be related to PGC1α-mediated metabolic networks. Therefore, this current study lays the foundation for the further development of a global metabolic analysis of cancer cells in individual patients, which ultimately will have significant potential for the discovery of novel therapeutics and precision medicine in SETD2-inactivated ccRCC.

Enhanced and controlled chromatin extraction from FFPE tissues and the application to ChIP-seq.

BACKGROUND: Epigenetic dysregulation is involved in the etiology and progression of various human diseases. Formalin-fixed paraffin-embedded (FFPE) samples represent the gold standard for archiving pathology samples, and thus FFPE samples are a major resource of samples in clinical research. However, chromatin-based epigenetic assays in the clinical settings are limited to fresh or frozen samples, and are hampered by low chromatin yield in FFPE samples due to the lack of a reliable and efficient chromatin preparation method. Here, we introduce a new chromatin extraction method from FFPE tissues (Chrom-EX PE) for chromatin-based epigenetic assays. RESULTS: During rehydration of FFPE tissues, applying a tissue-level cross-link reversal into the deparaffinized tissue at 65 °C dramatically increased chromatin yield in the soluble fraction. The resulting chromatin is compatible with targeted ChIP-qPCR and genome-wide ChIP-seq approaches. The chromatin prepared by Chrom-EX PE showed a gradual fragmentation pattern with varying incubation temperature. At temperatures below 37 °C, the majority of soluble chromatin is over 1 kb. The soluble chromatin prepared in the range of 45-60 °C showed a typical nucleosomal pattern. And the majority of chromatin prepared at 65 °C is close to mononucleosomal size. These observations indicate that chromatin preparation from FFPE samples can be controlled for downstream chromatin-based epigenetic assays. CONCLUSIONS: This study provided a new method that achieves efficient extraction of high-quality chromatin suitable for chromatin-based epigenetic assays with less damage on chromatin. This approach may provide a way to circumvent the over-fixed nature of FFPE tissues for future technology development.

Genetic and Epigenetic Heterogeneity in Normal Liver Homeostasis and Its Implications for Liver Disease and Hepatocellular Cancer.

Hepatocellular carcinoma (HCC) is the most prevalent primary tumor of the liver, and is steadily becoming one of the most lethal cancers worldwide. Liver resection, which is the recommended procedure for early localized HCC, results in frequent recurrence (50-70%), while the standard of care for late-stage HCC, multikinase inhibitors, only improves survival by a few months. The lack of success for these treatment modalities is attributable, at least in part, to marked phenotypic heterogeneity within the tumor. Intratumoral heterogeneity (ITH) has emerged as a defining characteristic of human tumors, with individual cancer cells displaying distinct differences in properties including growth rate, metastatic capacity, and response to treatment. This heterogeneity, which is unlikely to be captured from a biopsy, impacts outcome because a single treatment targeting one cancer-specific pathway would spare tumor cells having distinct characteristics. Development of effective biomarkers remains a major challenge for similar reasons. Understanding, interpreting, and circumventing the impact of ITH is therefore paramount for developing reliable biomarkers and designing effective individualized treatment strategies for HCC.

CAME: identification of chromatin accessibility from nucleosome occupancy and methylome sequencing.

Motivation: Chromatin accessibility plays a key role in epigenetic regulation of gene activation and silencing. Open chromatin regions allow regulatory elements such as transcription factors and polymerases to bind for gene expression while closed chromatin regions prevent the activity of transcriptional machinery. Recently, Methyltransferase Accessibility Protocol for individual templates-Bisulfite Genome Sequencing (MAPit-BGS) and nucleosome occupancy and methylome sequencing (NOMe-seq) have been developed for simultaneously profiling chromatin accessibility and DNA methylation on single molecules. Therefore, there is a great demand in developing computational methods to identify chromatin accessibility from MAPit-BGS and NOMe-seq. Results: In this article, we present CAME (Chromatin Accessibility and Methylation), a seed-extension based approach that identifies chromatin accessibility from NOMe-seq. The efficiency and effectiveness of CAME were demonstrated through comparisons with other existing techniques on both simulated and real data, and the results show that our method not only can precisely identify chromatin accessibility but also outperforms other methods. Availability and Implementation: CAME is implemented in java and the program is freely available online at http://sourceforge.net/projects/came/. Contacts: jechoi@gru.edu or khryu@dblab.chungbuk.ac.kr. Supplementary information: Supplementary data are available at Bioinformatics online.

Purification of nanogram-range immunoprecipitated DNA in ChIP-seq application.

BACKGROUND: Chromatin immunoprecipitation-sequencing (ChIP-seq) is a widely used epigenetic approach for investigating genome-wide protein-DNA interactions in cells and tissues. The approach has been relatively well established but several key steps still require further improvement. As a part of the procedure, immnoprecipitated DNA must undergo purification and library preparation for subsequent high-throughput sequencing. Current ChIP protocols typically yield nanogram quantities of immunoprecipitated DNA mainly depending on the target of interest and starting chromatin input amount. However, little information exists on the performance of reagents used for the purification of such minute amounts of immunoprecipitated DNA in ChIP elution buffer and their effects on ChIP-seq data. Here, we compared DNA recovery, library preparation efficiency, and ChIP-seq results obtained with several commercial DNA purification reagents applied to 1 ng ChIP DNA and also investigated the impact of conditions under which ChIP DNA is stored. RESULTS: We compared DNA recovery of ten commercial DNA purification reagents and phenol/chloroform extraction from 1 to 50 ng of immunopreciptated DNA in ChIP elution buffer. The recovery yield was significantly different with 1 ng of DNA while similar in higher DNA amounts. We also observed that the low nanogram range of purified DNA is prone to loss during storage depending on the type of polypropylene tube used. The immunoprecipitated DNA equivalent to 1 ng of purified DNA was subject to DNA purification and library preparation to evaluate the performance of four better performing purification reagents in ChIP-seq applications. Quantification of library DNAs indicated the selected purification kits have a negligible impact on the efficiency of library preparation. The resulting ChIP-seq data were comparable with the dataset generated by ENCODE consortium and were highly correlated between the data from different purification reagents. CONCLUSIONS: This study provides comparative data on commercial DNA purification reagents applied to nanogram-range immunopreciptated ChIP DNA and evidence for the importance of storage conditions of low nanogram-range purified DNA. We verified consistent high performance of a subset of the tested reagents. These results will facilitate the improvement of ChIP-seq methodology for low-input applications.

iTagPlot: an accurate computation and interactive drawing tool for tag density plot.

MOTIVATION: Tag density plots are very important to intuitively reveal biological phenomena from capture-based sequencing data by visualizing the normalized read depth in a region. RESULTS: We have developed iTagPlot to compute tag density across functional features in parallel using multicores and a grid engine and to interactively explore it in a graphical user interface. It allows us to stratify features by defining groups based on biological function and measurement, summary statistics and unsupervised clustering. AVAILABILITY AND IMPLEMENTATION: http://sourceforge.net/projects/itagplot/.

SETD2 loss in renal epithelial cells drives epithelial-to-mesenchymal transition in a TGF-ß-independent manner.

Histone-lysine N-methyltransferase SETD2 (SETD2), the sole histone methyltransferase that catalyzes trimethylation of lysine 36 on histone H3 (H3K36me3), is often mutated in clear cell renal cell carcinoma (ccRCC). SETD2 mutation and/or loss of H3K36me3 is linked to metastasis and poor outcome in ccRCC patients. Epithelial-to-mesenchymal transition (EMT) is a major pathway that drives invasion and metastasis in various cancer types. Here, using novel kidney epithelial cell lines isogenic for SETD2, we discovered that SETD2 inactivation drives EMT and promotes migration, invasion, and stemness in a transforming growth factor-beta-independent manner. This newly identified EMT program is triggered in part through secreted factors, including cytokines and growth factors, and through transcriptional reprogramming. RNA-seq and assay for transposase-accessible chromatin sequencing uncovered key transcription factors upregulated upon SETD2 loss, including SOX2, POU2F2 (OCT2), and PRRX1, that could individually drive EMT and stemness phenotypes in SETD2 wild-type (WT) cells. Public expression data from SETD2 WT/mutant ccRCC support the EMT transcriptional signatures derived from cell line models. In summary, our studies reveal that SETD2 is a key regulator of EMT phenotypes through cell-intrinsic and cell-extrinsic mechanisms that help explain the association between SETD2 loss and ccRCC metastasis.

Histone Modifications and miRNA Perturbations Contribute to Transcriptional Dysregulation of Hypertrophy in Obstructive Hypertrophic Cardiomyopathy.

Background: Recently, we demonstrated transcriptional downregulation of hypertrophy pathways in myectomy tissue derived from patients with obstructive hypertrophic cardiomyopathy (HCM) despite translational activation of hypertrophy pathways. The mechanisms and modifiers of this transcriptional dysregulation in HCM remain unexplored. We hypothesized that miRNA and post-translational modifications of histones contribute to transcriptional dysregulation in HCM. Methods: First, miRNA-sequencing and chromatin immunoprecipitation sequencing (ChIP-seq) were performed on HCM myectomy tissue and control donor hearts to characterize miRNA and differential histone marks across the genome. Next, the differential miRNA and histone marks were integrated with RNA-sequencing (RNA-seq) data. Finally, the effects of miRNA and histones were removed in silico to determine their necessity for transcriptional dysregulation of pathways. Results: miRNA-analysis identified 19 differentially expressed miRNA. ChIP-seq analysis identified 2,912 (7%) differential H3K4me3 peaks, 23,339 (21%) differential H3K9ac peaks, 33 (0.05%) differential H3K9me3 peaks, 58,837 (42%) differential H3K27ac peaks, and 853 (3%) differential H3K27me3 peaks. Univariate analysis of concordance between H3K9ac with RNA-seq data showed activation of cardiac hypertrophy signaling, while H3K27me showed downregulation of cardiac hypertrophy signaling. Similarly, miRNAs were predicted to result in downregulation of cardiac hypertrophy signaling. In silico knock-out that effects either miRNA or histones attenuated transcriptional downregulation while knocking out both abolished downregulation of hypertrophy pathways completely. Conclusion: Myectomy tissue from patients with obstructive HCM shows transcriptional dysregulation, including transcriptional downregulation of hypertrophy pathways mediated by miRNA and post-translational modifications of histones. Cardiac hypertrophy loci showed activation via changes in H3K9ac and a mix of activation and repression via H3K27ac.

DNA methylation at DLGAP2 and risk for relapse in alcohol dependence during acamprosate treatment.

BACKGROUND: Alcohol use disorders are prevalent mental disorders with significant health implications. Epigenetic alterations may play a role in their pathogenesis, as DNA methylation at several genes has been associated with these disorders. We have previously shown that methylation in the DLGAP2 gene, coding for a synaptic density protein, is associated with alcohol dependence. In this study, we aimed to examine the association between DLGAP2 methylation and treatment response among patients undergoing acamprosate treatment. METHODS: 102 patients under acamprosate treatment were included. DNA methylation analysis at DLGAP2 was performed by bisulfite pyrosequencing at the start and after 3-month treatment. Treatment outcomes were having a relapse during the treatment and severity of craving at the end of three months. Cox proportional hazard and linear regression models were performed. RESULTS: Patients whose methylation levels were decreased during the treatment showed an increased risk for relapse within three months in comparison to the ones without methylation change (hazard ratio [HR]=2.44; 95% confidence interval [CI]=1.04, 5.73; p=0.04). For the same group, a positive association for the severity of craving was observed, yet statistical significance was not reached (ß=2.97; 95% CI=-0.41, 6.34; p=0.08). CONCLUSION: We demonstrate that patients whose DLGAP2 methylation levels decrease during acamprosate treatment are more likely to relapse compared to the ones without changes. This is in line with our previous findings showing that DLGAP2 methylation is lower in alcohol dependent subjects compared to controls, and might suggest a role for changes in DLGAP2 methylation in treatment response.

Rapid and transient recruitment of DNMT1 to DNA double-strand breaks is mediated by its interaction with multiple components of the DNA damage response machinery.

DNA methylation is an epigenetic mark critical for regulating transcription, chromatin structure and genome stability. Although many studies have shed light on how methylation impacts transcription and interfaces with the histone code, far less is known about how it regulates genome stability. We and others have shown that DNA methyltransferase 1 (DNMT1), the maintenance methyltransferase, contributes to the cellular response to DNA damage, yet DNMT1's exact role in this process remains unclear. DNA damage, particularly in the form of double-strand breaks (DSBs), poses a major threat to genome integrity. Cells therefore possess a potent system to respond to and repair DSBs, or to initiate cell death. In the current study, we used a near-infrared laser microirradiation system to directly study the link between DNMT1 and DSBs. Our results demonstrate that DNMT1 is rapidly but transiently recruited to DSBs. DNMT1 recruitment is dependent on its ability to interact with both PCNA and the ATR effector kinase CHK1, but is independent of its catalytic activity. In addition, we show for the first time that DNMT1 interacts with the 9-1-1 PCNA-like sliding clamp and that this interaction also contributes to DNMT1 localization to DNA DSBs. Finally, we demonstrate that DNMT1 modulates the rate of DSB repair and is essential for suppressing abnormal activation of the DNA damage response in the absence of exogenous damage. Taken together, our studies provide compelling additional evidence for DNMT1 acting as a regulator of genome integrity and as an early responder to DNA DSBs.