基于ChIP-seq的差异组蛋白修饰区域的筛选

组蛋白修饰是在基因组水平上起到重要调控作用的表观遗传修饰,随着ChIP-Seq的广泛使用,高通量数据的积累,为从全基因组研究组蛋白修饰模式奠定了基础。但目前缺乏在多样本间筛选疾病相关的调控区域的方法,因而本文开发了一种多细胞系的差异筛选算法来识别差异组蛋白修饰区域。本文通过窗口移动法来估计组蛋白修饰水平,并根据信息熵理论定量各个细胞系之间的差异。基于随机背景来确定差异显著性阈值。利用此算法来筛选人类全基因组9个细胞系间H3K4me3差异的区域,结果显示这些区域显著富集在基因启动子上和其他重要的染色质状态上,且与先前人们发现的活性启动子染色质状态显著重叠。通过文献挖掘进一步证实了与白血病相关的基因组标记。这些结果表明基于熵的策略可有效地挖掘多细胞系间以及与疾病相关的差异组蛋白修饰。     

Identification of Epigenetic Biomarkers of Lung Adenocarcinoma through Multi-Omics Data Analysis

Epigenetic mechanisms such as DNA methylation or histone modifications are essential for the regulation of gene expression and development of tissues. Alteration of epigenetic modifications can be used as an epigenetic biomarker for diagnosis and as promising targets for epigenetic therapy. A recent study explored cancer-cell specific epigenetic biomarkers by examining different types of epigenetic modifications simultaneously. However, it was based on microarrays and reported biomarkers that were also present in normal cells at a low frequency. Here, we first analyzed multi-omics data (including ChIP-Seq data of six types of histone modifications: H3K27ac, H3K4me1, H3K9me3, H3K36me3, H3K27me3, and H3K4me3) obtained from 26 lung adenocarcinoma cell lines and a normal cell line. We identified six genes with both H3K27ac and H3K4me3 histone modifications in their promoter regions, which were not present in the normal cell line, but present in ≥85% (22 out of 26) and ≤96% (25 out of 26) of the lung adenocarcinoma cell lines. Of these genes, NUP210 (encoding a main component of the nuclear pore complex) was the only gene in which the two modifications were not detected in another normal cell line. RNA-Seq analysis revealed that NUP210 was aberrantly overexpressed among the 26 lung adenocarcinoma cell lines, although the frequency of NUP210 overexpression was lower (19.3%) in 57 lung adenocarcinoma tissue samples studied and stored in another database. This study provides a basis to discover epigenetic biomarkers highly specific to a certain cancer, based on multi-omics data at the cell population level.     

HPRep: Quantifying Reproducibility in HiChIP and PLAC-Seq Datasets

HiChIP and PLAC-Seq are emerging technologies for studying genome-wide long-range chromatin interactions mediated by the protein of interest, enabling more sensitive and cost-efficient interrogation of protein-centric chromatin conformation. However, due to the unbalanced read distribution introduced by protein immunoprecipitation, existing reproducibility measures developed for Hi-C data are not appropriate for the analysis of HiChIP and PLAC-Seq data. Here, we present HPRep, a stratified and weighted correlation metric derived from normalized contact counts, to quantify reproducibility in HiChIP and PLAC-Seq data. We applied HPRep to multiple real datasets and demonstrate that HPRep outperforms existing reproducibility measures developed for Hi-C data. Specifically, we applied HPRep to H3K4me3 PLAC-Seq data from mouse embryonic stem cells and mouse brain tissues as well as H3K27ac HiChIP data from human lymphoblastoid cell line GM12878 and leukemia cell line K562, showing that HPRep can more clearly separate among pseudo-replicates, real replicates, and non-replicates. Furthermore, in an H3K4me3 PLAC-Seq dataset consisting of 11 samples from four human brain cell types, HPRep demonstrated the expected clustering of data that could not be achieved by existing methods developed for Hi-C data, highlighting the need for a reproducibility metric tailored to HiChIP and PLAC-Seq data.     

Genome-wide characterization of menin-dependent H3K4me3 reveals a specific role for menin in the regulation of genes implicated in MEN1-like tumors

Inactivating mutations in the MEN1 gene predisposing to the multiple endocrine neoplasia type 1 (MEN1) syndrome can also cause sporadic pancreatic endocrine tumors. MEN1 encodes menin, a subunit of MLL1/MLL2-containing histone methyltransferase complexes that trimethylate histone H3 at lysine 4 (H3K4me3). The importance of menin-dependent H3K4me3 in normal and transformed pancreatic endocrine cells is unclear. To study the role of menin-dependent H3K4me3, we performed in vitro differentiation of wild-type as well as menin-null mouse embryonic stem cells (mESCs) into pancreatic islet-like endocrine cells (PILECs). Gene expression analysis and genome-wide H3K4me3 ChIP-Seq profiling in wild-type and menin-null mESCs and PILECs revealed menin-dependent H3K4me3 at the imprinted Dlk1-Meg3 locus in mESCs, and all four Hox loci in differentiated PILECs. Specific and significant loss of H3K4me3 and gene expression was observed for genes within the imprinted Dlk1-Meg3 locus in menin-null mESCs and the Hox loci in menin-null PILECs. Given that the reduced expression of genes within the DLK1-MEG3 locus and the HOX loci is associated with MEN1-like sporadic tumors, our data suggests a possible role for menin-dependent H3K4me3 at these genes in the initiation and progression of sporadic pancreatic endocrine tumors. Furthermore, our investigation also demonstrates that menin-null mESCs can be differentiated in vitro into islet-like endocrine cells, underscoring the utility of menin-null mESC-derived specialized cell types for genome-wide high-throughput studies.     

MLL3 suppresses tumorigenesis through regulating TNS3 enhancer activity

MLL3 is a histone H3K4 methyltransferase that is frequently mutated in cancer, but the underlying molecular mechanisms remain elusive. Here, we found that MLL3 depletion by CRISPR/sgRNA significantly enhanced cell migration, but did not elevate the proliferation rate of cancer cells. Through RNA-Seq and ChIP-Seq approaches, we identified TNS3 as the potential target gene for MLL3. MLL3 depletion caused downregulation of H3K4me1 and H3K27ac on an enhancer ~ 7 kb ahead of TNS3. 3C assay indicated the identified enhancer interacts with TNS3 promoter and repression of enhancer activity by dCas9-KRAB system impaired TNS3 expression. Exogenous expression of TNS3 in MLL3 deficient cells completely blocked the enhanced cell migration phenotype. Taken together, our study revealed a novel mechanism for MLL3 in suppressing cancer, which may provide novel targets for diagnosis or drug development.Subject terms: Tumour-suppressor proteins, Chromosomes     

Role of hPHF1 in H3K27 methylation and Hox gene silencing

Polycomb group (PcG) proteins are required for maintaining the silent state of the homeotic genes and other important developmental regulators. The silencing function of the PcG proteins has been linked to their intrinsic histone modifying enzymatic activities. The EED-EZH2 complex, containing the core subunits EZH2, EED, SUZ12, and RbAp48, functions as a histone H3K27-specific methyltransferase. Here we describe the identification and characterization of a related EED-EZH2 protein complex which is distinguished from the previous complex by the presence of another PcG protein, hPHF1. Consistent with the ability of hPHF1 to stimulate the enzymatic activity of the core EED-EZH2 complex in vitro, manipulation of mPcl1, the mouse counterpart of hPHF1, in NIH 3T3 cells and cells of the mouse male germ cell line GC1spg results in global alteration of H3K27me2 and H3K27me3 levels and Hox gene expression. Small interfering RNA-mediated knockdown of mPcl1 affects association of the Eed-Ezh2 complex with certain Hox genes, such as HoxA10, as well as Hox gene expression concomitant with an alteration on the H3K27me2 levels of the corresponding promoters. Therefore, our results reveal hPHF1 as a component of a novel EED-EZH2 complex and demonstrate its important role in H3K27 methylation and Hox gene silencing.     

Reprogramming of H3K27me3 Is Critical for Acquisition of Pluripotency from Cultured Arabidopsis Tissues

In plants, multiple detached tissues are capable of forming a pluripotent cell mass, termed callus, when cultured on media containing appropriate plant hormones. Recent studies demonstrated that callus resembles the root-tip meristem, even if it is derived from aerial organs. This finding improves our understanding of the regeneration process of plant cells; however, the molecular mechanism that guides cells of different tissue types to form a callus still remains elusive. Here, we show that genome-wide reprogramming of histone H3 lysine 27 trimethylation (H3K27me3) is a critical step in the leaf-to-callus transition. The Polycomb Repressive Complex 2 (PRC2) is known to function in establishing H3K27me3. By analyzing callus formation of mutants corresponding to different histone modification pathways, we found that leaf blades and/or cotyledons of the PRC2 mutants curly leaf swinger (clf swn) and embryonic flower2 (emf2) were defective in callus formation. We identified the H3K27me3-covered loci in leaves and calli by a ChIP-chip assay, and we found that in the callus H3K27me3 levels decreased first at certain auxin-pathway genes. The levels were then increased at specific leaf genes but decreased at a number of root-regulatory genes. Changes in H3K27me3 levels were negatively correlated with expression levels of the corresponding genes. One possible role of PRC2-mediated H3K27me3 in the leaf-to-callus transition might relate to elimination of leaf features by silencing leaf-regulatory genes, as most leaf-preferentially expressed regulatory genes could not be silenced in the leaf explants of clf swn. In contrast to the leaf explants, the root explants of both clf swn and emf2 formed calli normally, possibly because the root-to-callus transition bypasses the leaf gene silencing process. Furthermore, our data show that PRC2-mediated H3K27me3 and H3K27 demethylation act in parallel in the reprogramming of H3K27me3 during the leaf-to-callus transition, suggesting a general mechanism for cell fate transition in plants.     

A chromatin-dependent mechanism regulates gene expression at the core of the Arabidopsis circadian clock

The mechanisms of circadian clock function in Arabidopsis rely on the complex relationships among core clock components. The current model of the Arabidopsis oscillator comprises a myriad of repressors but the mechanisms responsible for activation remain largely unknown. In our recent studies, we have demonstrated that the rhythms in H3 acetylation (H3ac) and H3K4 trimethylation (H3K4me3) are a key mechanism at the positive arm of the oscillator. H3K4me3 rhythmic accumulation is delayed compared to that of H3ac, which opens the possibility for separate roles for each mark. Indeed, the use of inhibitors that block H3K4me3 accumulation was concomitant with increased clock repressor binding, suggesting that H3K4me3 might control the timing from activation to repression. Plants mis-expressing the histone methyltransferase SET DOMAIN GROUP 2 (SDG2/ATXR3) displayed altered H3K4me3 accumulation, oscillator gene expression and clock repressor binding, suggesting that SDG2/ATXR3 is a key component contributing to proper circadian expression.     

Proteome-wide analysis in Saccharomyces cerevisiae identifies several PHD fingers as novel direct and selective binding modules of histone H3 methylated at either lysine 4 or lysine 36

The PHD finger motif is a signature chromatin-associated motif that is found throughout eukaryotic proteomes. Here we have determined the histone methyl-lysine binding activity of the PHD fingers present within the Saccharomyces cerevisiae proteome. We provide evidence on the genomic scale that PHD fingers constitute a general class of effector modules for histone H3 trimethylated at lysine 4 (H3K4me3) and histone H3 trimethylated at lysine 36 (H3K36me3). Structural modeling of PHD fingers demonstrates a conserved mechanism for recognizing the trimethyl moiety and provides insight into the molecular basis of affinity for the different methyl-histone ligands. Together, our study suggests that a common function for PHD fingers is to transduce methyl-lysine events and sheds light on how a single histone modification can be linked to multiple biological outcomes.     

Dynamic distribution of linker histone h1.5 in cellular differentiation

Linker histones are essential components of chromatin, but the distributions and functions of many during cellular differentiation are not well understood. Here, we show that H1.5 binds to genic and intergenic regions, forming blocks of enrichment, in differentiated human cells from all three embryonic germ layers but not in embryonic stem cells. In differentiated cells, H1.5, but not H1.3, binds preferentially to genes that encode membrane and membrane-related proteins. Strikingly, 37% of H1.5 target genes belong to gene family clusters, groups of homologous genes that are located in proximity to each other on chromosomes. H1.5 binding is associated with gene repression and is required for SIRT1 binding, H3K9me2 enrichment, and chromatin compaction. Depletion of H1.5 results in loss of SIRT1 and H3K9me2, increased chromatin accessibility, deregulation of gene expression, and decreased cell growth. Our data reveal for the first time a specific and novel function for linker histone subtype H1.5 in maintenance of condensed chromatin at defined gene families in differentiated human cells.     

Histone H3 lysine 27 and 9 hypermethylation within the Bad promoter region mediates 5-Aza-2′-deoxycytidine-induced Leydig cell apoptosis: implications of 5-Aza-2′-deoxycytidine toxicity to male reproduction

5-Aza-2′-deoxycitidine (5-Aza), an anticancer agent, results in substantial toxicity to male reproduction, causing a decline in sperm quality associated with reduced testosterone. Here, we report that 5-Aza increased the apoptotic protein Bad epigenetically in the testosterone-producing mouse TM3 Leydig cell line. 5-Aza decreased cell viability in a dose- and time-dependent manner with concomitant increase in Bad protein. This increase is accompanied by increased cleavages of both poly ADP ribose polymerase and caspase-3. Flow cytometric analysis further supported 5-Aza-derived apoptosis in TM3 cells. Bisulfite sequencing analysis failed to identify putative methylcytosine site(s) in CpG islands of the Bad promoter. A chromatin immunoprecipitation assay revealed decreased levels of trimethylation at lysine 27 of histone H3 (H3K27-3me) and H3K9-3me in the Bad promoter region in response to 5-Aza treatment. Knock-down by siRNA of enhancer of zeste homologue 2 (EZH2), a histone methyltransferase responsible for H3K27-3me, or demethylation of H3K9-3me by BIX-01294 showed significantly increased levels in Bad expression and consequent Leydig cell apoptosis. In conclusion, our results demonstrate for the first time that Bad expression is regulated at least by EZH2-mediated H3K27-3me or G9a-like protein/euchromatic histone methyltransferase 1 (GLP/Eu-HMTase1)-mediated H3K9-3me in mouse TM3 Leydig cells, which may be implicated in 5-Aza-derived toxicity to male reproduction.