Wang, Liguo

Significance The multiple-domain protein ZMYND8 is well known as a chromatin reader that recognizes methylated, acetylated, and/or mutated histones. Here we report that ZMYND8 is overexpressed in human ccRCC, and that it preferentially binds T487-phosphorylated EZH2. We show that in hypoxia-stimulated or von Hippel–Lindau–deficient cells, ZMYND8 binding of EZH2 is important for maintenance of EZH2 phosphorylation level, inhibition of PRC2 complex formation, and enhanced interaction between EZH2 and FOXM1. Our results reveal a reader function of ZMYND8 in recognizing a phosphorylation “code” on a nonhistone protein. Our findings also identify ZMYND8 binding as a facet of regulation and function of EZH2 phosphorylation in the Polycomb-dependent to -independent switch of EZH2/PRC2 activities.

Abstract The tumor-suppressor protein RB acts as a transcription repressor via interaction of its pocket domain with an LXCXE motif in histone deacetylase (HDAC) proteins such as HDAC1. Here, we demonstrate that HDAC5 deficient for the LXCXE motif interacts with both RB-N (via an FXXXV motif) and RB-C segments, and such interactions are diminished by phosphorylation of RB serine-249/threonine-252 and threonine-821. HDAC5 was frequently downregulated or deleted in human cancers such as prostate cancer. Loss of HDAC5 increased histone H3 lysine 27 acetylation (H3K27-ac) and circumvented RB-mediated repression of cell-cycle–related pro-oncogenic genes. HDAC5 loss also conferred resistance to CDK4/6 inhibitors such as palbociclib in prostate and breast cancer cells in vitro and prostate tumors in vivo, but this effect was overcome by the BET-CBP/p300 dual inhibitor NEO2734. Our findings reveal an unknown role of HDAC5 in RB-mediated histone deacetylation and gene repression and define a new mechanism modulating CDK4/6 inhibitor therapeutic sensitivity in cancer cells. Significance: This study defines a previously uncharacterized role of HDAC5 in tumor suppression and provides a viable strategy to overcome CDK4/6 inhibitor resistance in HDAC5-deficent cancer.

Abstract 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.

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.

Abstract Motivation The significance of long non-coding RNAs (lncRNAs) in many biological processes and diseases has gained intense interests over the past several years. However, computational identification of lncRNAs in a wide range of species remains challenging; it requires prior knowledge of well-established sequences and annotations or species-specific training data, but the reality is that only a limited number of species have high-quality sequences and annotations. Results Here we first characterize lncRNAs in contrast to protein-coding RNAs based on feature relationship and find that the feature relationship between open reading frame length and guanine-cytosine (GC) content presents universally substantial divergence in lncRNAs and protein-coding RNAs, as observed in a broad variety of species. Based on the feature relationship, accordingly, we further present LGC, a novel algorithm for identifying lncRNAs that is able to accurately distinguish lncRNAs from protein-coding RNAs in a cross-species manner without any prior knowledge. As validated on large-scale empirical datasets, comparative results show that LGC outperforms existing algorithms by achieving higher accuracy, well-balanced sensitivity and specificity, and is robustly effective (>90% accuracy) in discriminating lncRNAs from protein-coding RNAs across diverse species that range from plants to mammals. To our knowledge, this study, for the first time, differentially characterizes lncRNAs and protein-coding RNAs based on feature relationship, which is further applied in computational identification of lncRNAs. Taken together, our study represents a significant advance in characterization and identification of lncRNAs and LGC thus bears broad potential utility for computational analysis of lncRNAs in a wide range of species. Availability and implementation LGC web server is publicly available at http://bigd.big.ac.cn/lgc/calculator. The scripts and data can be downloaded at http://bigd.big.ac.cn/biocode/tools/BT000004. Supplementary information Supplementary data are available at Bioinformatics online.

A microfluidic device was designed to efficiently fragment chromatin for downstream epigenetic assays.

Background: Transcription factors are DNA-binding proteins that play key roles in many fundamental biological processes. Unraveling their interactions with DNA is essential to identify their target genes and understand the regulatory network. Genome-wide identification of their binding sites became feasible thanks to recent progress in experimental and computational approaches. ChIP-chip, ChIP-seq, and ChIP-exo are three widely used techniques to demarcate genome-wide transcription factor binding sites. Objective: This review aims to provide an overview of these three techniques including their experiment procedures, computational approaches, and popular analytic tools. Conclusion: ChIP-chip, ChIP-seq, and ChIP-exo have been the major techniques to study genome- wide in vivo protein-DNA interaction. Due to the rapid development of next-generation sequencing technology, array-based ChIP-chip is deprecated and ChIP-seq has become the most widely used technique to identify transcription factor binding sites in genome-wide. The newly developed ChIP-exo further improves the spatial resolution to single nucleotide. Numerous tools have been developed to analyze ChIP-chip, ChIP-seq and ChIP-exo data. However, different programs may employ different mechanisms or underlying algorithms thus each will inherently include its own set of statistical assumption and bias. So choosing the most appropriate analytic program for a given experiment needs careful considerations. Moreover, most programs only have command line interface so their installation and usage will require basic computation expertise in Unix/Linux.

Abstract Purpose: Homozygous deletions play important roles in carcinogenesis. The genome-wide screening for homozygously deleted genes in many different cancer types with a large number of patient specimens representing the tumor heterogeneity has not been done. Experimental Design: We performed integrative analyses of the copy-number profiles of 10,759 patients across 31 cancer types from The Cancer Genome Atlas project. Results: We found that the type-I interferon, α-, and β-defensin genes were homozygously deleted in 19 cancer types with high frequencies (7%–31%, median = 12%; interquartile range = 10%–16.5%). Patients with homozygous deletion of interferons exhibited significantly shortened overall or disease-free survival time in a number of cancer types, whereas patients with homozygous deletion of defensins did not significantly associate with worse overall or disease-free survival. Gene expression analyses suggested that homozygous deletion of interferon and defensin genes could activate genes involved in oncogenic and cell-cycle pathways but repress other genes involved in immune response pathways, suggesting their roles in promoting tumorigenesis and helping cancer cells evade immune surveillance. Further analysis of the whole exomes of 109 patients with melanoma demonstrated that the homozygous deletion of interferon (P = 0.0029, OR = 11.8) and defensin (P = 0.06, OR = 2.79) genes are significantly associated with resistance to anti-CTLA4 immunotherapy. Conclusions: Our analysis reveals that the homozygous deletion of interferon and defensin genes is prevalent in human cancers, and importantly this feature can be used as a novel prognostic biomarker for immunotherapy resistance. Clin Cancer Res; 24(14); 3299–308. ©2018 AACR.

Tumor tissues are chronically exposed to hypoxia owing to aberrant vascularity. Lipid droplet (LD) accumulation is a hallmark of hypoxic cancer cells, yet how LDs form and function during hypoxia remains poorly understood. Herein, we report that in various cancer cells upon oxygen deprivation, HIF-1 activation down-modulates LD catabolism mediated by adipose triglyceride lipase (ATGL), the key enzyme for intracellular lipolysis. Proteomics and functional analyses identified hypoxia-inducible gene 2 (HIG2), a HIF-1 target, as a new inhibitor of ATGL. Knockout of HIG2 enhanced LD breakdown and fatty acid (FA) oxidation, leading to increased ROS production and apoptosis in hypoxic cancer cells as well as impaired growth of tumor xenografts. All of these effects were reversed by co-ablation of ATGL. Thus, by inhibiting ATGL, HIG2 acts downstream of HIF-1 to sequester FAs in LDs away from the mitochondrial pathways for oxidation and ROS generation, thereby sustaining cancer cell survival in hypoxia.