H3K27me3 forms BLOCs over silent genes and intergenic regions and specifies a histone banding pattern on a mouse autosomal chromosome

Advancing personalized medicine in neurodegenerative diseases: The role of epigenetics and pharmacoepigenomics in pharmacotherapy

About 80 % of brain disorders have a genetic basis. The pathogenesis of most neurodegenerative diseases is associated with a myriad of genetic defects, epigenetic alterations (DNA methylation, histone/chromatin remodeling, miRNA dysregulation), and environmental factors. The emergence of new sequencing technologies and tools to study the epigenome has led to identifying predictive biomarkers for earlier diagnosis, opening up the possibility of prophylactical interventions. As a result, advances in pharmacogenetics and pharmacoepigenomics now allow for personalized treatments based on the profile of each patient and the specific genetic and epigenetic mechanisms involved. This Review highlights the complexity of neurodegenerative diseases and the variability in patient responses to pharmacotherapy, emphasizing the influence of genetic polymorphisms on the pharmacokinetics and pharmacodynamics of drugs used to treat those conditions. We specifically discuss the potential modulatory effect of several genetic polymorphisms associated with an increased risk of developing different neurodegenerative diseases. We explore genetic and genomic technologies and the potential of analyzing individual-specific drug metabolism to predict and influence drug response and associated clinical outcomes. We also provide insights into the mechanism of action of the drugs under investigation and their potential impact on disease-modifying pathways. Finally, the Review underscores the great potential of this field to enhance the effectiveness and safety of drug treatments through personalized medicine.

Increased EZH2 function in regulatory T cells promotes their capacity to suppress autoimmunity by driving effector differentiation prior to activation

The immunosuppressive function of regulatory T (Treg) cells is essential for maintaining immune homeostasis. Enhancer of zeste homolog 2 (EZH2), a histone H3 lysine 27 (H3K27) methyltransferase, plays a key role in maintaining Treg cell function upon CD28 co-stimulation, and Ezh2 deletion in Treg cells causes autoimmunity. Here we assessed whether increased EZH2 activity in Treg cells would improve Treg cell function. Using an Ezh2 gain-of-function mutation, Ezh2 Y641F , we found that Treg cells expressing Ezh2 Y641F displayed an increased effector Treg phenotype and were poised for improved homing to organ tissues. Expression of Ezh2 Y641F in Treg cells led to more rapid remission from autoimmunity. H3K27me3 profiling and transcriptomic analysis revealed a redistribution of H3K27me3, which prompted a gene expression profile in naïve Ezh2 Y641F Treg cells that recapitulated aspects of CD28-activated Ezh2 WT Treg cells. Altogether, increased EZH2 activity promotes the differentiation of effector Treg cells that can better suppress autoimmunity.

Highlights

EZH2 function promotes effector differentiation of Treg cells.EZH2 function promotes Treg cell migration to organ tissues.EZH2 function in Treg cells improves remission from autoimmunity.EZH2 function poises naïve Treg cells to adopt a CD28-activated phenotype.

Analyzing histone ChIP-seq data with a bin-based probability of being signal

Histone ChIP-seq is one of the primary methods for charting the cellular epigenomic landscape, the components of which play a critical regulatory role in gene expression. Analyzing the activity of regulatory elements across datasets and cell types can be challenging due to shifting peak positions and normalization artifacts resulting from, for example, differing read depths, ChIP efficiencies, and target sizes. Moreover, broad regions of enrichment seen in repressive histone marks often evade detection by commonly used peak callers. Here, we present a simple and versatile method for identifying enriched regions in ChIP-seq data that relies on estimating a gamma distribution fit to non-overlapping 5kB genomic bins to establish a global background. We use this distribution to assign a probability of being signal (PBS) between zero and one to each 5 kB bin. This approach, while lower in resolution than typical peak-calling methods, provides a straightforward way to identify enriched regions and compare enrichments among multiple datasets, by transforming the data to values that are universally normalized and can be readily visualized and integrated with downstream analysis methods. We demonstrate applications of PBS for both broad and narrow histone marks, and provide several illustrations of biological insights which can be gleaned by integrating PBS scores with downstream data types.

H4K20me3 is important for Ash1-mediated H3K36me3 and transcriptional silencing in facultative heterochromatin in a fungal pathogen

Facultative heterochromatin controls development and differentiation in many eukaryotes. In metazoans, plants, and many filamentous fungi, facultative heterochromatin is characterized by transcriptional repression and enrichment with nucleosomes that are trimethylated at histone H3 lysine 27 (H3K27me3). While loss of H3K27me3 results in derepression of transcriptional gene silencing in many species, additional up- and downstream layers of regulation are necessary to mediate control of transcription in chromosome regions enriched with H3K27me3. Here, we investigated the effects of one histone mark on histone H4, namely H4K20me3, in the fungus Zymoseptoria tritici, a globally important pathogen of wheat. Deletion of kmt5, the gene encoding the sole methyltransferase responsible for H4K20 methylation, resulted in global derepression of transcription, especially in regions of facultative heterochromatin. Derepression in the absence of H4K20me3 not only affected known genes but also a large number of novel, previously undetected transcripts generated from regions of facultative heterochromatin on accessory chromosomes. Transcriptional activation in kmt5 deletion strains was accompanied by a complete loss of Ash1-mediated H3K36me3 and chromatin reorganization affecting H3K27me3 and H3K4me2 distribution in regions of facultative heterochromatin. Strains with H4K20L, M or Q mutations in the single histone H4 gene of Z. tritici recapitulated these chromatin changes, suggesting that H4K20me3 is important for Ash1-mediated H3K36me3. The ∆kmt5 mutants we obtained were more sensitive to genotoxic stressors than wild type and both, ∆kmt5 and ∆ash1, showed greatly increased rates of accessory chromosome loss. Taken together, our results provide insights into an unsuspected mechanism involved in the assembly and maintenance of facultative heterochromatin.

The role of heterochromatin in 3D genome organization during preimplantation development

During the early stages of mammalian development, the epigenetic state of the parental genome is completely reprogrammed to give rise to the totipotent embryo. An important aspect of this remodeling concerns the heterochromatin and the spatial organization of the genome. While heterochromatin and genome organization are intricately linked in pluripotent and somatic systems, little is known about their relationship in the totipotent embryo. In this review, we summarize the current knowledge on the reprogramming of both regulatory layers. In addition, we discuss available evidence on their relationship and put this in the context of findings in other systems.

Single-cell sortChIC identifies hierarchical chromatin dynamics during hematopoiesis

Post-translational histone modifications modulate chromatin activity to affect gene expression. How chromatin states underlie lineage choice in single cells is relatively unexplored. We develop sort-assisted single-cell chromatin immunocleavage (sortChIC) and map active (H3K4me1 and H3K4me3) and repressive (H3K27me3 and H3K9me3) histone modifications in the mouse bone marrow. During differentiation, hematopoietic stem and progenitor cells (HSPCs) acquire active chromatin states mediated by cell-type-specifying transcription factors, which are unique for each lineage. By contrast, most alterations in repressive marks during differentiation occur independent of the final cell type. Chromatin trajectory analysis shows that lineage choice at the chromatin level occurs at the progenitor stage. Joint profiling of H3K4me1 and H3K9me3 demonstrates that cell types within the myeloid lineage have distinct active chromatin but share similar myeloid-specific heterochromatin states. This implies a hierarchical regulation of chromatin during hematopoiesis: heterochromatin dynamics distinguish differentiation trajectories and lineages, while euchromatin dynamics reflect cell types within lineages.

scChIX-seq infers dynamic relationships between histone modifications in single cells

Regulation of chromatin states involves the dynamic interplay between different histone modifications to control gene expression. Recent advances have enabled mapping of histone marks in single cells, but most methods are constrained to profile only one histone mark per cell. Here, we present an integrated experimental and computational framework, scChIX-seq (single-cell chromatin immunocleavage and unmixing sequencing), to map several histone marks in single cells. scChIX-seq multiplexes two histone marks together in single cells, then computationally deconvolves the signal using training data from respective histone mark profiles. This framework learns the cell-type-specific correlation structure between histone marks, and therefore does not require a priori assumptions of their genomic distributions. Using scChIX-seq, we demonstrate multimodal analysis of histone marks in single cells across a range of mark combinations. Modeling dynamics of in vitro macrophage differentiation enables integrated analysis of chromatin velocity. Overall, scChIX-seq unlocks systematic interrogation of the interplay between histone modifications in single cells.

DMRscaler: a scale-aware method to identify regions of differential DNA methylation spanning basepair to multi-megabase features

BACKGROUND

Pathogenic mutations in genes that control chromatin function have been implicated in rare genetic syndromes. These chromatin modifiers exhibit extraordinary diversity in the scale of the epigenetic changes they affect, from single basepair modifications by DNMT1 to whole genome structural changes by PRM1/2. Patterns of DNA methylation are related to a diverse set of epigenetic features across this full range of epigenetic scale, making DNA methylation valuable for mapping regions of general epigenetic dysregulation. However, existing methods are unable to accurately identify regions of differential methylation across this full range of epigenetic scale directly from DNA methylation data.

RESULTS

To address this, we developed DMRscaler, a novel method that uses an iterative windowing procedure to capture regions of differential DNA methylation (DMRs) ranging in size from single basepairs to whole chromosomes. We benchmarked DMRscaler against several DMR callers in simulated and natural data comparing XX and XY peripheral blood samples. DMRscaler was the only method that accurately called DMRs ranging in size from 100 bp to 1 Mb (pearson's r = 0.94) and up to 152 Mb on the X-chromosome. We then analyzed methylation data from rare-disease cohorts that harbor chromatin modifier gene mutations in NSD1, EZH2, and KAT6A where DMRscaler identified novel DMRs spanning gene clusters involved in development.

CONCLUSION

Taken together, our results show DMRscaler is uniquely able to capture the size of DMR features across the full range of epigenetic scale and identify novel, co-regulated regions that drive epigenetic dysregulation in human disease.

Epigenetics and Early Development

The epigenome controls all aspect of eukaryotic development as the packaging of DNA greatly affects gene expression. Epigenetic changes are reversible and do not affect the DNA sequence itself but rather control levels of gene expression. As a result, the science of epigenetics focuses on the physical configuration of chromatin in the proximity of gene promoters rather than on the mechanistic effects of gene sequences on transcription and translation. In the present review we discuss three prominent epigenetic modifications, DNA methylation, histone methylation/acetylation, and the effects of chromatin remodeling complexes. Specifically, we introduce changes to the methylated state of DNA through DNA methyltransferases and DNA demethylases, discuss the effects of histone tail modifications such as histone acetylation and methylation on gene expression and present the functions of major ATPase subunit containing chromatin remodeling complexes. We also introduce examples of how changes in these epigenetic factors affect early development in humans and mice. In summary, this review provides an overview over the most important epigenetic mechanisms and provides examples of the dramatic effects of epigenetic changes in early mammalian development.

Time makes histone H3 modifications drift in mouse liver

To detect the epigenetic drift of time passing, we determined the genome-wide distributions of mono- and tri-methylated lysine 4 and acetylated and tri-methylated lysine 27 of histone H3 in the livers of healthy 3, 6 and 12 months old C57BL/6 mice. The comparison of different age profiles of histone H3 marks revealed global redistribution of histone H3 modifications with time, in particular in intergenic regions and near transcription start sites, as well as altered correlation between the profiles of different histone modifications. Moreover, feeding mice with caloric restriction diet, a treatment known to retard aging, reduced the extent of changes occurring during the first year of life in these genomic regions.

Mitotically heritable, RNA polymerase II-independent H3K4 dimethylation stimulates INO1 transcriptional memory

For some inducible genes, the rate and molecular mechanism of transcriptional activation depend on the prior experiences of the cell. This phenomenon, called epigenetic transcriptional memory, accelerates reactivation, and requires both changes in chromatin structure and recruitment of poised RNA polymerase II (RNAPII) to the promoter. Memory of inositol starvation in budding yeast involves a positive feedback loop between transcription factor-dependent interaction with the nuclear pore complex and histone H3 lysine 4 dimethylation (H3K4me2). While H3K4me2 is essential for recruitment of RNAPII and faster reactivation, RNAPII is not required for H3K4me2. Unlike RNAPII-dependent H3K4me2 associated with transcription, RNAPII-independent H3K4me2 requires Nup100, SET3C, the Leo1 subunit of the Paf1 complex and, upon degradation of an essential transcription factor, is inherited through multiple cell cycles. The writer of this mark (COMPASS) physically interacts with the potential reader (SET3C), suggesting a molecular mechanism for the spreading and re-incorporation of H3K4me2 following DNA replication.

Distinctive aspects of the placental epigenome and theories as to how they arise

The placenta has a methylome dramatically unlike that of any somatic cell type. Among other distinctions, it features low global DNA methylation, extensive “partially methylated domains” packed in dense heterochromatin and methylation of hundreds of CpG islands important in somatic development. These features attract interest in part because a substantial fraction of human cancers feature the exact same phenomena, suggesting parallels between epigenome formation in placentation and cancer. Placenta also features an expanded set of imprinted genes, some of which come about by distinctive developmental pathways. Recent discoveries, some from far outside the placental field, shed new light on how the unusual placental epigenetic state may arise. Nonetheless, key questions remain unresolved.

Mapping cis-regulatory elements in the midgestation mouse placenta

The placenta is a temporary organ that provides the developing fetus with nutrients, oxygen, and protection in utero. Defects in its development, which may be caused by misregulated gene expression, can lead to devastating outcomes for the mother and fetus. In mouse, placental defects during midgestation commonly lead to embryonic lethality. However, the regulatory mechanisms controlling expression of genes during this period have not been thoroughly investigated. Therefore, we generated and analyzed ChIP-seq data for multiple histone modifications known to mark cis-regulatory regions. We annotated active and poised promoters and enhancers, as well as regions generally associated with repressed gene expression. We found that poised promoters were associated with neuronal development genes, while active promoters were largely associated with housekeeping genes. Active and poised enhancers were associated with placental development genes, though only active enhancers were associated with genes that have placenta-specific expression. Motif analysis within active enhancers identified a large network of transcription factors, including those that have not been previously studied in the placenta and are candidates for future studies. The data generated and genomic regions annotated provide researchers with a foundation for future studies, aimed at understanding how specific genes in the midgestation mouse placenta are regulated.

Bivalent Regulation and Related Mechanisms of H3K4/27/9me3 in Stem Cells

The "bivalent domain" is a unique histone modification region consisting of two histone tri-methylation modifications. Over the years, it has been revealed that the maintenance and dynamic changes of the bivalent domains play a vital regulatory role in the differentiation of various stem cell systems, as well as in other cells, such as immunomodulation. Tri-methylation modifications involved in the formation of the bivalent domains are interrelated and mutually regulated, thus regulating many life processes of cells. Tri-methylation of histone H3 at lysine 4 (H3K4me3), tri-methylation of histone H3 at lysine 9 (H3K9me3) and tri-methylation of histone H3 at lysine 27 (H3K27me3) are the main tri-methylation modifications involved in the formation of bivalent domains. The three form different bivalent domains in pairs. Furthermore, it is equally clear that H3K4me3 is a positive regulator of transcription and that H3K9me3/H3K27me3 are negative regulators. Enzymes related to the regulation of histone methylation play a significant role in the "homeostasis" and "breaking homeostasis" of the bivalent domains. Bivalent domains regulate target genes, upstream transcription, downstream targeting regulation and related cytokines during the establishment and breakdown of homeostasis, and exert the specific regulation of stem cells. Indeed, a unified mechanism to explain the bivalent modification in all stem cells has been difficult to define, and whether the bivalent modification is antagonistic in inducing the differentiation of homologous stem cells is controversial. In this review, we focus on the different bivalent modifications in several key stem cells and explore the main mechanisms and effects of these modifications involved. Finally, we discussed the close relationship between bivalent domains and immune cells, and put forward the prospect of the application of bivalent domains in the field of stem cells.

Characterizing Genetic Regulatory Elements in Ovine Tissues

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

Sequence deeper without sequencing more: Bayesian resolution of ambiguously mapped reads

Next-generation sequencing (NGS) has transformed molecular biology and contributed to many seminal insights into genomic regulation and function. Apart from whole-genome sequencing, an NGS workflow involves alignment of the sequencing reads to the genome of study, after which the resulting alignments can be used for downstream analyses. However, alignment is complicated by the repetitive sequences; many reads align to more than one genomic locus, with 15-30% of the genome not being uniquely mappable by short-read NGS. This problem is typically addressed by discarding reads that do not uniquely map to the genome, but this practice can lead to systematic distortion of the data. Previous studies that developed methods for handling ambiguously mapped reads were often of limited applicability or were computationally intensive, hindering their broader usage. In this work, we present SmartMap: an algorithm that augments industry-standard aligners to enable usage of ambiguously mapped reads by assigning weights to each alignment with Bayesian analysis of the read distribution and alignment quality. SmartMap is computationally efficient, utilizing far fewer weighting iterations than previously thought necessary to process alignments and, as such, analyzing more than a billion alignments of NGS reads in approximately one hour on a desktop PC. By applying SmartMap to peak-type NGS data, including MNase-seq, ChIP-seq, and ATAC-seq in three organisms, we can increase read depth by up to 53% and increase the mapped proportion of the genome by up to 18% compared to analyses utilizing only uniquely mapped reads. We further show that SmartMap enables the analysis of more than 140,000 repetitive elements that could not be analyzed by traditional ChIP-seq workflows, and we utilize this method to gain insight into the epigenetic regulation of different classes of repetitive elements. These data emphasize both the dangers of discarding ambiguously mapped reads and their power for driving biological discovery.

Genome Architecture Facilitates Phenotypic Plasticity in the Honeybee (Apis mellifera)

Phenotypic plasticity, the ability of an organism to alter its phenotype in response to an environmental cue, facilitates rapid adaptation to changing environments. Plastic changes in morphology and behavior are underpinned by widespread gene expression changes. However, it is unknown if, or how, genomes are structured to ensure these robust responses. Here, we use repression of honeybee worker ovaries as a model of plasticity. We show that the honeybee genome is structured with respect to plasticity; genes that respond to an environmental trigger are colocated in the honeybee genome in a series of gene clusters, many of which have been assembled in the last 80 My during the evolution of the Apidae. These clusters are marked by histone modifications that prefigure the gene expression changes that occur as the ovary activates, suggesting that these genomic regions are poised to respond plastically. That the linear sequence of the honeybee genome is organized to coordinate widespread gene expression changes in response to environmental influences and that the chromatin organization in these regions is prefigured to respond to these influences is perhaps unexpected and has implications for other examples of plasticity in physiology, evolution, and human disease.

The deubiquitinase Usp9x regulates PRC2-mediated chromatin reprogramming during mouse development

Pluripotent cells of the mammalian embryo undergo extensive chromatin rewiring to prepare for lineage commitment after implantation. Repressive H3K27me3, deposited by Polycomb Repressive Complex 2 (PRC2), is reallocated from large blankets in pre-implantation embryos to mark promoters of developmental genes. The regulation of this global redistribution of H3K27me3 is poorly understood. Here we report a post-translational mechanism that destabilizes PRC2 to constrict H3K27me3 during lineage commitment. Using an auxin-inducible degron system, we show that the deubiquitinase Usp9x is required for mouse embryonic stem (ES) cell self-renewal. Usp9x-high ES cells have high PRC2 levels and bear a chromatin and transcriptional signature of the pre-implantation embryo, whereas Usp9x-low ES cells resemble the post-implantation, gastrulating epiblast. We show that Usp9x interacts with, deubiquitinates and stabilizes PRC2. Deletion of Usp9x in post-implantation embryos results in the derepression of genes that normally gain H3K27me3 after gastrulation, followed by the appearance of morphological abnormalities at E9.5, pointing to a recurrent link between Usp9x and PRC2 during development. Usp9x is a marker of "stemness" and is mutated in various neurological disorders and cancers. Our results unveil a Usp9x-PRC2 regulatory axis that is critical at peri-implantation and may be redeployed in other stem cell fate transitions and disease states.

Histone modifications form a cell-type-specific chromosomal bar code that persists through the cell cycle

Chromatin configuration influences gene expression in eukaryotes at multiple levels, from individual nucleosomes to chromatin domains several Mb long. Post-translational modifications (PTM) of core histones seem to be involved in chromatin structural transitions, but how remains unclear. To explore this, we used ChIP-seq and two cell types, HeLa and lymphoblastoid (LCL), to define how changes in chromatin packaging through the cell cycle influence the distributions of three transcription-associated histone modifications, H3K9ac, H3K4me3 and H3K27me3. We show that chromosome regions (bands) of 10-50 Mb, detectable by immunofluorescence microscopy of metaphase (M) chromosomes, are also present in G1 and G2. They comprise 1-5 Mb sub-bands that differ between HeLa and LCL but remain consistent through the cell cycle. The same sub-bands are defined by H3K9ac and H3K4me3, while H3K27me3 spreads more widely. We found little change between cell cycle phases, whether compared by 5 Kb rolling windows or when analysis was restricted to functional elements such as transcription start sites and topologically associating domains. Only a small number of genes showed cell-cycle related changes: at genes encoding proteins involved in mitosis, H3K9 became highly acetylated in G2M, possibly because of ongoing transcription. In conclusion, modified histone isoforms H3K9ac, H3K4me3 and H3K27me3 exhibit a characteristic genomic distribution at resolutions of 1 Mb and below that differs between HeLa and lymphoblastoid cells but remains remarkably consistent through the cell cycle. We suggest that this cell-type-specific chromosomal bar-code is part of a homeostatic mechanism by which cells retain their characteristic gene expression patterns, and hence their identity, through multiple mitoses.

H3K27me3-rich genomic regions can function as silencers to repress gene expression via chromatin interactions

The mechanisms underlying gene repression and silencers are poorly understood. Here we investigate the hypothesis that H3K27me3-rich regions of the genome, defined from clusters of H3K27me3 peaks, may be used to identify silencers that can regulate gene expression via proximity or looping. We find that H3K27me3-rich regions are associated with chromatin interactions and interact preferentially with each other. H3K27me3-rich regions component removal at interaction anchors by CRISPR leads to upregulation of interacting target genes, altered H3K27me3 and H3K27ac levels at interacting regions, and altered chromatin interactions. Chromatin interactions did not change at regions with high H3K27me3, but regions with low H3K27me3 and high H3K27ac levels showed changes in chromatin interactions. Cells with H3K27me3-rich regions knockout also show changes in phenotype associated with cell identity, and altered xenograft tumor growth. Finally, we observe that H3K27me3-rich regions-associated genes and long-range chromatin interactions are susceptible to H3K27me3 depletion. Our results characterize H3K27me3-rich regions and their mechanisms of functioning via looping.

RECOGNICER: A coarse-graining approach for identifying broad domains from ChIP-seq data

Background

Histone modifications are major factors that define chromatin states and have functions in regulating gene expression in eukaryotic cells. Chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) technique has been widely used for profiling the genome-wide distribution of chromatin-associating protein factors. Some histone modifications, such as H3K27me3 and H3K9me3, usually mark broad domains in the genome ranging from kilobases (kb) to megabases (Mb) long, resulting in diffuse patterns in the ChIP-seq data that are challenging for signal separation. While most existing ChIP-seq peak-calling algorithms are based on local statistical models without account of multi-scale features, a principled method to identify scale-free board domains has been lacking.

Methods

Here we present RECOGNICER (Recursive coarse-graining identification for ChIP-seq enriched regions), a computational method for identifying ChIP-seq enriched domains on a large range of scales. The algorithm is based on a coarse-graining approach, which uses recursive block transformations to determine spatial clustering of local enriched elements across multiple length scales.

Results

We apply RECOGNICER to call H3K27me3 domains from ChIP-seq data, and validate the results based on H3K27me3's association with repressive gene expression. We show that RECOGNICER outperforms existing ChIP-seq broad domain calling tools in identifying more whole domains than separated pieces.

Conclusion

RECOGNICER can be a useful bioinformatics tool for next-generation sequencing data analysis in epigenomics research.

Morphine leads to global genome changes in H3K27me3 levels via a Polycomb Repressive Complex 2 (PRC2) self-regulatory mechanism in mESCs

Background

Environmentally induced epigenetic changes can lead to health problems or disease, but the mechanisms involved remain unclear. Morphine can pass through the placental barrier leading to abnormal embryo development. However, the mechanism by which morphine causes these effects and how they sometimes persist into adulthood is not well known. To unravel the morphine-induced chromatin alterations involved in aberrant embryo development, we explored the role of the H3K27me3/PRC2 repressive complex in gene expression and its transmission across cellular generations in response to morphine.

Results

Using mouse embryonic stem cells as a model system, we found that chronic morphine treatment induces a global downregulation of the histone modification H3K27me3. Conversely, ChIP-Seq showed a remarkable increase in H3K27me3 levels at specific genomic sites, particularly promoters, disrupting selective target genes related to embryo development, cell cycle and metabolism. Through a self-regulatory mechanism, morphine downregulated the transcription of PRC2 components responsible for H3K27me3 by enriching high H3K27me3 levels at the promoter region. Downregulation of PRC2 components persisted for at least 48 h (4 cell cycles) following morphine removal, though promoter H3K27me3 levels returned to control levels.

Conclusions

Morphine induces targeting of the PRC2 complex to selected promoters, including those of PRC2 components, leading to characteristic changes in gene expression and a global reduction in H3K27me3. Following morphine removal, enhanced promoter H3K27me3 levels revert to normal sooner than global H3K27me3 or PRC2 component transcript levels. We suggest that H3K27me3 is involved in initiating morphine-induced changes in gene expression, but not in their maintenance.

Graphic abstract

Model of Polycomb repressive complex 2 (PRC2) and H3K27me3 alterations induced by chronic morphine exposure. Morphine induces H3K27me3 enrichment at promoters of genes encoding core members of the PRC2 complex and is associated with their transcriptional downregulation.

Broad genic repression domains signify enhanced silencing of oncogenes

Cancers result from a set of genetic and epigenetic alterations. Most known oncogenes were identified by gain-of-function mutations in cancer, yet little is known about their epigenetic features. Through integrative analysis of 11,596 epigenomic profiles and mutations from >8200 tumor-normal pairs, we discover broad genic repression domains (BGRD) on chromatin as an epigenetic signature for oncogenes. A BGRD is a widespread enrichment domain of the repressive histone modification H3K27me3 and is further enriched with multiple other repressive marks including H3K9me3, H3K9me2, and H3K27me2. Further, BGRD displays widespread enrichment of repressed cis-regulatory elements. Shortening of BGRDs is linked to derepression of transcription. BGRDs at oncogenes tend to be conserved across normal cell types. Putative tumor-promoting genes and lncRNAs defined using BGRDs are experimentally verified as required for cancer phenotypes. Therefore, BGRDs play key roles in epigenetic regulation of cancer and provide a direction for mutation-independent discovery of oncogenes.

Epigenetically altered genes can have a key role in cancer pathobiology but epigenetic signatures that distinguish oncogenes are not yet known. Here, the authors identify broad genic repression domains, defined by widespread H3K27me3 modification, as an epigenetic signature to provide mutation-independent information for discovery of potential oncogenes.

A physical basis for quantitative ChIP-sequencing

ChIP followed by next-generation sequencing (ChIP-Seq) is a key technique for mapping the distribution of histone posttranslational modifications (PTMs) and chromatin-associated factors across genomes. There is a perceived challenge to define a quantitative scale for ChIP-Seq data, and as such, several approaches making use of exogenous additives, or "spike-ins," have recently been developed. Herein, we report on the development of a quantitative, physical model defining ChIP-Seq. The quantitative scale on which ChIP-Seq results should be compared emerges from the model. To test the model and demonstrate the quantitative scale, we examine the impacts of an EZH2 inhibitor through the lens of ChIP-Seq. We report a significant increase in immunoprecipitation of presumed off-target histone PTMs after inhibitor treatment, a trend predicted by the model but contrary to spike-in-based indications. Our work also identifies a sensitivity issue in spike-in normalization that has not been considered in the literature, placing limitations on its utility and trustworthiness. We call our new approach the sans-spike-in method for quantitative ChIP-sequencing (siQ-ChIP). A number of changes in community practice of ChIP-Seq, data reporting, and analysis are motivated by this work.

A Chromodomain-Helicase-DNA-Binding Factor Functions in Chromatin Modification and Gene Regulation

Proteins in the Chromodomain-Helicase/ATPase-DNA-binding domain (CHD) family are divided into three groups. The function of group I CHD proteins in nucleosome positioning is well established, while that of group II members (represented by CHD3/Mi2) remains unclear. Using high-throughput approaches, we investigated the function of the group II rice (Oryza sativa) CHD protein CHR729 in nucleosome positioning, gene expression, histone methylation, and binding. Our data revealed that the chr729 mutation led to increased nucleosome occupancy in the rice genome and altered the expression and histone H3K4me3 modification of many, mainly underexpressed, genes. Further analysis showed that the mutation affected both the deposition and depletion of H3K4me3 in distinct chromatin regions, with concomitant changes in H3K27me3 modification. Genetic and genomic analyses revealed that CHR729 and JMJ703, an H3K4 demethylase, had agonistic, antagonistic, and independent functions in modulating H3K4me3 and the expression of subsets of genes. In addition, CHR729 binding was enriched in H3K4me3-marked genic and H3K27me3-marked intergenic regions. The results indicate that CHR729 has distinct functions in regulating H3K4me3 and H3K27me3 modifications and gene expression at different chromatin domains and provide insight into chromatin regulation of bivalent genes marked by both H3K4me3 and H3K27me3.

Systematic characterization of chromatin modifying enzymes identifies KDM3B as a critical regulator in castration resistant prostate cancer

Androgen deprivation therapy (ADT) is the standard care for prostate cancer (PCa) patients who fail surgery or radiotherapy. While initially effective, the cancer almost always recurs as a more aggressive castration resistant prostate cancer (CRPC). Previous studies have demonstrated that chromatin modifying enzymes can play a critical role in the conversion to CRPC. However, only a handful of these potential pharmacological targets have been tested. Therefore, in this study, we conducted a focused shRNA screen of chromatin modifying enzymes previously shown to be involved in cellular differentiation. We found that altering the balance between histone methylation and demethylation impacted growth and proliferation. Of all genes tested, KDM3B, a histone H3K9 demethylase, was found to have the most antiproliferative effect. These results were phenocopied with a KDM3B CRISPR/Cas9 knockout. When tested in several PCa cell lines, the decrease in proliferation was remarkably specific to androgen-independent cells. Genetic rescue experiments showed that only the enzymatically active KDM3B could recover the phenotype. Surprisingly, despite the decreased proliferation of androgen-independent cell no alterations in the cell cycle distribution were observed following KDM3B knockdown. Whole transcriptome analyses revealed changes in the gene expression profile following loss of KDM3B, including downregulation of metabolic enzymes such as ARG2 and RDH11. Metabolomic analysis of KDM3B knockout showed a decrease in several critical amino acids. Overall, our work reveals, for the first time, the specificity and the dependence of KDM3B in CRPC proliferation.

Broad Heterochromatic Domains Open in Gonocyte Development Prior to De Novo DNA Methylation

Facultative heterochromatin forms and reorganizes in response to external stimuli. However, how the initial establishment of such a chromatin state is regulated in cell-cycle-arrested cells remains unexplored. Mouse gonocytes are arrested male germ cells, at which stage the genome-wide DNA methylome forms. Here, we discovered transiently accessible heterochromatin domains of several megabases in size in gonocytes and named them differentially accessible domains (DADs). Open DADs formed in gene desert and gene cluster regions, primarily at transposons, with the reprogramming of histone marks, suggesting DADs as facultative heterochromatin. De novo DNA methylation took place with two waves in gonocytes: the first region specific and the second genome-wide. DADs were resistant to the first wave and their opening preceded the second wave. In addition, the higher-order chromosome architecture was reorganized with less defined chromosome compartments in gonocytes. These findings suggest that multiple layers of chromatin reprogramming facilitate de novo DNA methylation.

A unified encyclopedia of human functional DNA elements through fully automated annotation of 164 human cell types

Semi-automated genome annotation methods such as Segway take as input a set of genome-wide measurements such as of histone modification or DNA accessibility and output an annotation of genomic activity in the target cell type. Here we present annotations of 164 human cell types using 1615 data sets. To produce these annotations, we automated the label interpretation step to produce a fully automated annotation strategy. Using these annotations, we developed a measure of the importance of each genomic position called the “conservation-associated activity score.” We further combined all annotations into a single, cell type-agnostic encyclopedia that catalogs all human regulatory elements.

H3K9me3, H3K36me3, and H4K20me3 Expression Correlates with Patient Outcome in Esophageal Squamous Cell Carcinoma as Epigenetic Markers

BACKGROUND

Histone methylation, as an essential pattern of posttranslational modifications, contributes to multiple cancer-related biological processes. Dysregulation of histone methylation is now considered a biomarker for cancer prognosis.

AIMS

This study investigated and evaluated the potential role of four histone lysine trimethylation markers as biomarkers for esophageal squamous cell carcinoma (ESCC) prognosis.

METHODS

Tissue arrays were made from 135 paraffin-embedded ESCC samples and examined for histone markers by immunohistochemistry, and 10 pairs of cancer and noncancerous mucosa tissues from ESCC patients were investigated with Western blot. Chi-squared test, Kaplan-Meier analysis with log-rank test, and Cox proportional hazard trend analyses were performed to assess the prognostic values of the markers.

RESULTS

Histone 3 lysine 4 trimethylation (H3K4me3), histone 3 lysine 9 trimethylation (H3K9me3), and histone 4 lysine 20 trimethylation (H4K20me3), but not histone 3 lysine 36 trimethylation (H3K36me3), showed stronger immunostaining signals in tumor tissues than in the corresponding adjacent non-neoplastic mucosa tissues. The expression patterns of H3K36me3, H3K9me3, and H4K20me3 correlated with tumor infiltrating depth, lymph node involvement, and pTNM stage. Low-scoring H3K9me3 and H4K20me3 predicted better prognosis, while H3K36me3 manifested the opposite trend. Poor prognosis occurred in ESCC patients with expression patterns of high levels of H3K9me3, high levels of H4K20me3, and low levels of H3K36me3 expression.

CONCLUSIONS

H3K9me3, H4K20me3, and H3K36me3 showed a close relationship with clinical features and were considered independent risk factors for survival of ESCC patients. The combination of H3K9me3, H4K20me3, and H3K36me3 expression, rather than the expression of a single histone marker, is believed to further enhance evaluations of ESCC prognosis and management.

The Airn lncRNA does not require any DNA elements within its locus to silence distant imprinted genes

Long non-coding (lnc) RNAs are numerous and found throughout the mammalian genome, and many are thought to be involved in the regulation of gene expression. However, the majority remain relatively uncharacterised and of uncertain function making the use of model systems to uncover their mode of action valuable. Imprinted lncRNAs target and recruit epigenetic silencing factors to a cluster of imprinted genes on the same chromosome, making them one of the best characterized lncRNAs for silencing distant genes in cis. In this study we examined silencing of the distant imprinted gene Slc22a3 by the lncRNA Airn in the Igf2r imprinted cluster in mouse. Previously we proposed that imprinted lncRNAs may silence distant imprinted genes by disrupting promoter-enhancer interactions by being transcribed through the enhancer, which we called the enhancer interference hypothesis. Here we tested this hypothesis by first using allele-specific chromosome conformation capture (3C) to detect interactions between the Slc22a3 promoter and the locus of the Airn lncRNA that silences it on the paternal chromosome. In agreement with the model, we found interactions enriched on the maternal allele across the entire Airn gene consistent with multiple enhancer-promoter interactions. Therefore, to test the enhancer interference hypothesis we devised an approach to delete the entire Airn gene. However, the deletion showed that there are no essential enhancers for Slc22a2, Pde10a and Slc22a3 within the Airn gene, strongly indicating that the Airn RNA rather than its transcription is responsible for silencing distant imprinted genes. Furthermore, we found that silent imprinted genes were covered with large blocks of H3K27me3 on the repressed paternal allele. Therefore we propose an alternative hypothesis whereby the chromosome interactions may initially guide the lncRNA to target imprinted promoters and recruit repressive chromatin, and that these interactions are lost once silencing is established.

Long non-coding (lnc) RNAs are numerous in the mammalian genome and many have been implicated in gene regulation. However, the vast majority are uncharacterised and of uncertain function making known functional lncRNAs valuable models for understanding their mechanism of action. One mode of lncRNA action is to recruit epigenetic silencing to target distant genes on the same chromosome. A well-characterized group of lncRNAs that act in this way to silence genes are imprinted lncRNAs. In this study we examined how the imprinted lncRNA Airn silences genes in the Igf2r imprinted cluster, focusing primarily on silencing of the distant imprinted gene Slc22a3. We found that Airn expression blocks chromosome interactions between the Slc22a3 promoter and the Airn gene locus. By making a large genomic deletion including the Airn gene we showed that these interactions are not essential enhancer/promoter interactions, but may help to guide the Airn RNA to target genes to recruit epigenetic silencing. Our study adds to the understanding of how lncRNAs may act to silence distant genes.

Does the Presence of Transposable Elements Impact the Epigenetic Environment of Human Duplicated Genes?

Epigenetic modifications have an important role to explain part of the intra- and inter-species variation in gene expression. They also have a role in the control of transposable elements (TEs) whose activity may have a significant impact on genome evolution by promoting various mutations, which are expected to be mostly deleterious. A change in the local epigenetic landscape associated with the presence of TEs is expected to affect the expression of neighboring genes since these modifications occurring at TE sequences can spread to neighboring sequences. In this work, we have studied how the epigenetic modifications of genes are conserved and what the role of TEs is in this conservation. For that, we have compared the conservation of the epigenome associated with human duplicated genes and the differential presence of TEs near these genes. Our results show higher epigenome conservation of duplicated genes from the same family when they share similar TE environment, suggesting a role for the differential presence of TEs in the evolutionary divergence of duplicates through variation in the epigenetic landscape.

Misregulation of the IgH Locus in Thymocytes

Functional antigen receptor genes are assembled by somatic rearrangements that are largely lymphocyte lineage specific. The immunoglobulin heavy chain (IgH) gene locus is unique amongst the seven antigen receptor loci in undergoing partial gene rearrangements in the wrong lineage. Here we demonstrate that breakdown of lineage-specificity is associated with inappropriate activation of the Eμ enhancer during T cell development by a different constellation of transcription factors than those used in developing B cells. This is reflected in reduced enhancer-induced epigenetic changes, eRNAs, formation of the RAG1/2-rich recombination center, attenuated chromatin looping and markedly different utilization of D_H gene segments in CD4⁺CD8⁺ (DP) thymocytes. Additionally, CTCF-dependent V_H locus compaction is disrupted in DP cells despite comparable transcription factor binding in both lineages. These observations identify multiple mechanisms that contribute to lineage-specific antigen receptor gene assembly.

N-terminal acetylation and methylation differentially affect the function of MYL9

Deciphering the histone code has illustrated that acetylation or methylation on the same residue can have analogous or opposing roles. However, little is known about the interplay between these post-translational modifications (PTMs) on the same nonhistone residues. We have recently discovered that N-terminal acetyltransferases (NATs) and N-terminal methyltransferases (NRMTs) can have overlapping substrates and identified myosin regulatory light chain 9 (MYL9) as the first confirmed protein to occur in either α-amino-methylated (Nα-methyl) or α-amino-acetylated (Nα-acetyl) states in vivo Here we aim to determine if these PTMs function similarly or create different MYL9 proteoforms with distinct roles. We use enzymatic assays to directly verify MYL9 is a substrate of both NRMT1 and NatA and generate mutants of MYL9 that are exclusive for Nα-acetylation or Nα-methylation. We then employ eukaryotic cell models to probe the regulatory functions of these Nα-PTMs on MYL9. Our results show that, contrary to prevailing dogma, neither of these modifications regulate the stability of MYL9. Rather, exclusive Nα-acetylation promotes cytoplasmic roles of MYL9, while exclusive Nα-methylation promotes the nuclear role of MYL9 as a transcription factor. The increased cytoplasmic activity of Nα-acetylated MYL9 corresponds with increased phosphorylation at serine 19, a key MYL9 activating PTM. Increased nuclear activity of Nα-methylated MYL9 corresponds with increased DNA binding. Nα-methylation also results in a decrease of interactions between the N-terminus of MYL9 and a host of cytoskeletal proteins. These results confirm that Nα-acetylation and Nα-methylation differentially affect MYL9 function by creating distinct proteoforms with different internal PTM patterns and binding properties.

Global analysis of H3K4me3/H3K27me3 in Brachypodium distachyon reveals VRN3 as critical epigenetic regulation point in vernalization and provides insights into epigenetic memory

Vernalization, the requirement of plants for long-term exposure to low environmental temperature for flowering, is an epigenetic phenomenon. Histone modification regulation has been revealed in vernalization, but is limited to key genes. Now, we know that VRN1 is epigenetically critical for monocots. Genome-wide analysis is still unavailable, however. We performed chromatin immunoprecipitation-sequencing for H3K4me3/H3K27me3 in Brachypodium distachyon to obtain a global view of histone modifications in vernalization on a genome-wide scale and for different pathways/genes. Our data showed that H3K4me3 and H3K27me3 play distinct roles in vernalization. Unlike H3K4me3, H3K27me3 exhibited regional regulation, showed main regulation targets in vernalization and contributed to epigenetic memory. For genes in four flowering regulation pathways, only FT2 (functional ortholog of VRN3 in B. distachyon) and VRN1 showed coordinated changes in H3K4me3/H3K27me3. The epigenetic response at VRN3 was weaker under short-day than under long-day conditions. VRN3 was revealed as an epigenetic regulation point integrating vernalization and day length signals. We globally identified genes maintaining vernalization-induced epigenetic changes. Most of these genes showed dose-dependent vernalization responses, revealing a quantitative 'recording system' for vernalization. Our studies shed light on the epigenetic role of VRN3 and H3K4me3/H3K27me3 in vernalization and reveal genes underlying epigenetic memory, laying the foundation for further study.

Epigenetic analyses of planarian stem cells demonstrate conservation of bivalent histone modifications in animal stem cells

Planarian flatworms have an indefinite capacity to regenerate missing or damaged body parts owing to a population of pluripotent adult stems cells called neoblasts (NBs). Currently, little is known about the importance of the epigenetic status of NBs and how histone modifications regulate homeostasis and cellular differentiation. We have developed an improved and optimized ChIP-seq protocol for NBs in Schmidtea mediterranea and have generated genome-wide profiles for the active marks H3K4me3 and H3K36me3, and suppressive marks H3K4me1 and H3K27me3. The genome-wide profiles of these marks were found to correlate well with NB gene expression profiles. We found that genes with little transcriptional activity in the NB compartment but which switch on in post-mitotic progeny during differentiation are bivalent, being marked by both H3K4me3 and H3K27me3 at promoter regions. In further support of this hypothesis, bivalent genes also have a high level of paused RNA Polymerase II at the promoter-proximal region. Overall, this study confirms that epigenetic control is important for the maintenance of a NB transcriptional program and makes a case for bivalent promoters as a conserved feature of animal stem cells and not a vertebrate-specific innovation. By establishing a robust ChIP-seq protocol and analysis methodology, we further promote planarians as a promising model system to investigate histone modification–mediated regulation of stem cell function and differentiation.

Histone methylation changes are required for life cycle progression in the human parasite Schistosoma mansoni

Epigenetic mechanisms and chromatin structure play an important role in development. Their impact is therefore expected to be strong in parasites with complex life cycles and multiple, strikingly different, developmental stages, i.e. developmental plasticity. Some studies have already described how the chromatin structure, through histone modifications, varies from a developmental stage to another in a few unicellular parasites. While H3K4me3 profiles remain relatively constant, H3K27 trimethylation and bivalent methylation show strong variation. Inhibitors (A366 and GSK343) of H3K27 histone methyltransferase activity in S. mansoni efficiently blocked miracidium to sporocyst transition indicating that H3K27 trimethylation is required for life cycle progression. As S. mansoni is a multicellular parasite that significantly affects both the health and economy of endemic areas, a better understanding of fluke developmental processes within the definitive host will likely highlight novel disease control strategies. Towards this goal, we also studied H4K20me1 in female cercariae and adults. In particular, we found that bivalent trimethylation of H3K4 and H3K27 at the transcription start site of genes is a landmark of the cercarial stage. In cercariae, H3K27me3 presence and strong enrichment in H4K20me1 over long regions (10–100 kb) is associated with development related genes. Here, we provide a broad overview of the chromatin structure of a metazoan parasite throughout its most important lifecycle stages. The five developmental stages studied here present distinct chromatin structures, indicating that histone methylation plays an important role during development. Hence, components of the histone methylation (and demethylation) machinery may provide suitable Schistosomiasis control targets.

Schistosoma mansoni is a parasitic flatworm and causative agent of intestinal schistosomiasis, a neglected tropical disease affecting 67 million people worldwide. The parasite has a complex life cycle involving two consecutive obligate hosts (a poikilotherm snail and a homeotherm mammal) and two transitions between these hosts as free-swimming larvae. Here, we show that the chromatin structure of five different developmental stages is characterized by specific changes in chemical modifications of histones, basic proteins that are closely associated with DNA (trimethylation of lysines 4 and 27 and histone H3, and monomethylation of lysine 20 on histone H4). These modifications occur around protein coding genes as well as within repetitive genomic elements. A functional role for histone methylation during schistosome development was elucidated by the use of epi-drugs targeting G9a/GLP and EZH2 histone methyltransferase orthologs in S. mansoni. Our results indicate that histone methylation plays an important role during schistosome development and suggest that the enzymes responsible for maintaining these chromatin modifications are suitable targets for anti-schistosomal drugs.

Genes Predisposed to DNA Hypermethylation during Acquired Resistance to Chemotherapy Are Identified in Ovarian Tumors by Bivalent Chromatin Domains at Initial Diagnosis

Bivalent chromatin domains containing both active H3K4me3 and repressive H3K27me3 histone marks define gene sets poised for expression or silencing in differentiating embryonic stem (ES) cells. In cancer cells, aberrantly poised genes may facilitate changes in transcriptional states after exposure to anticancer drugs. In this study, we used ChIP-seq to characterize genome-wide positioning of H3K4me3- and H3K27me3-associated chromatin in primary high-grade serous ovarian carcinomas and in normal ovarian surface and fallopian tube tissue. Gene sets with proximal bivalent marks defined in this manner were evaluated subsequently as signatures of systematic change in DNA methylation and gene expression, comparing pairs of tissue samples taken from patients at primary presentation and relapse following chemotherapy. We found that gene sets harboring bivalent chromatin domains at their promoters in tumor tissue, but not normal epithelia, overlapped with Polycomb-repressive complex target genes as well as transcriptionally silenced genes in normal ovarian and tubal stem cells. The bivalently marked genes we identified in tumors before chemotherapy displayed increased promoter CpG methylation and reduced gene expression at relapse after chemotherapy of ovarian cancer. Overall, our results support the hypothesis that preexisting histone modifications at genes in a poised chromatin state may lead to epigenetic silencing during acquired drug resistance.Significance: These results suggest epigenetic targets for intervention to prevent the emergence of cancer drug resistance. Cancer Res; 78(6); 1383-91. ©2018 AACR.

Transcriptional decomposition reveals active chromatin architectures and cell specific regulatory interactions

Transcriptional regulation is tightly coupled with chromosomal positioning and three-dimensional chromatin architecture. However, it is unclear what proportion of transcriptional activity is reflecting such organisation, how much can be informed by RNA expression alone and how this impacts disease. Here, we develop a computational transcriptional decomposition approach separating the proportion of expression associated with genome organisation from independent effects not directly related to genomic positioning. We show that positionally attributable expression accounts for a considerable proportion of total levels and is highly informative of topological associating domain activities and organisation, revealing boundaries and chromatin compartments. Furthermore, expression data alone accurately predict individual enhancer-promoter interactions, drawing features from expression strength, stabilities, insulation and distance. We characterise predictions in 76 human cell types, observing extensive sharing of domains, yet highly cell-type-specific enhancer-promoter interactions and strong enrichments in relevant trait-associated variants. Overall, our work demonstrates a close relationship between transcription and chromatin architecture.

Roles of H3K27me2 and H3K27me3 Examined during Fate Specification of Embryonic Stem Cells

The polycomb repressive complex 2 (PRC2) methylates lysine 27 of histone H3 (H3K27) through its catalytic subunit Ezh2. PRC2-mediated di- and tri-methylation (H3K27me2/H3K27me3) have been interchangeably associated with gene repression. However, it remains unclear whether these two degrees of H3K27 methylation have different functions. In this study, we have generated isogenic mouse embryonic stem cells (ESCs) with a modified H3K27me2/H3K27me3 ratio. Our findings document dynamic developmental control in the genomic distribution of H3K27me2 and H3K27me3 at regulatory regions in ESCs. They also reveal that modifying the ratio of H3K27me2 and H3K27me3 is sufficient for the acquisition and repression of defined cell lineage transcriptional programs and phenotypes and influences induction of the ESC ground state.

Gene bivalency at Polycomb domains regulates cranial neural crest positional identity

The cranial neural crest cells are multipotent cells that provide head skeletogenic mesenchyme and are crucial for craniofacial patterning. We analyzed the chromatin landscapes of mouse cranial neural crest subpopulations in vivo. Early postmigratory subpopulations contributing to distinct mouse craniofacial structures displayed similar chromatin accessibility patterns yet differed transcriptionally. Accessible promoters and enhancers of differentially silenced genes carried H3K27me3/H3K4me2 bivalent chromatin marks embedded in large enhancer of zeste homolog 2-dependent Polycomb domains, indicating transcriptional poising. These postmigratory bivalent chromatin regions were already present in premigratory progenitors. At Polycomb domains, H3K27me3 antagonized H3K4me2 deposition, which was restricted to accessible sites. Thus, bivalent Polycomb domains provide a chromatin template for the regulation of cranial neural crest cell positional identity in vivo, contributing insights into the epigenetic regulation of face morphogenesis.

Broad Chromatin Domains: An Important Facet of Genome Regulation

Chromatin composition differs across the genome, with distinct compositions characterizing regions associated with different properties and functions. Whereas many histone modifications show local enrichment over genes or regulatory elements, marking can also span large genomic intervals defining broad chromatin domains. Here we highlight structural and functional features of chromatin domains marked by histone modifications, with a particular emphasis on the potential roles of H3K27 methylation domains in the organization and regulation of genome activity in metazoans.

c-Myc is Required for BRAFV600E-Induced Epigenetic Silencing by H3K27me3 in Tumorigenesis

BRAF^V600E mutation is frequently found in human cancers particularly thyroid cancer and melanoma, and is involved in the regulation of gene expression through activating MAPK/Erk signaling. Trimethylation of histone 3 lysine 27 (H3K27me3) is a critical epigenetic mark for the maintenance of gene silencing in tumorigenesis. However, molecular mechanism underlying the complex interplay between these two molecular events remains to be explored. In the present study, we conducted chromatin immunoprecipitation combined with next-generation sequencing (ChIP-Seq) and expression microarray analysis in NIH3T3 cells to explore the relationship between H3K27me3 and transcriptional regulation by BRAF^V600E mutation. Our results showed that activated MAPK/Erk signaling by BRAF^V600E mutation was a trigger of this epigenetic processing at many downstream target genes in cancer cell lines and Braf^V600E-induced thyroid cancer of transgenetic mice. By integrating ChIP-Seq and gene expression microarray data, we identified 150 down-regulated loci with increased levels of H3K27me3 in BRAF-mutant cells relative to BRAF wild-type cells. Our data also demonstrated that c-Myc, a downstream key effector of BRAF^V600E signaling, was required for BRAF^V600E-induced changes in H3K27me3 through regulating the components of the polycomb repressive complex 2 (PRC2) genes Ezh2, Suz12 and Jarid2 at both transcriptional levels via direct binding to their regulatory elements and post-transcriptional levels via repressing the miR-26a, miR-200b and miR-155. In addition, BRAF^V600E also caused gene silencing through Erk1/2-induced RNA polymerase II (RNAPII) poising and chromatin architecture. Collectively, our data uncover a previously unknown epigenetic mechanism in the tumorigenesis of BRAF^V600E-driven cancers.

Regulation of hepatic microRNA expression by hepatocyte nuclear factor 4 alpha

AIM

To uncover the role of hepatocyte nuclear factor 4 alpha (HNF4α) in regulating hepatic expression of microRNAs.

METHODS

Microarray and real-time PCR were used to determine hepatic expression of microRNAs in young-adult mice lacking Hnf4α expression in liver (Hnf4α-LivKO). Integrative genomics viewer software was used to analyze the public chromatin immunoprecipitation-sequencing datasets for DNA-binding of HNF4α, RNA polymerase-II, and histone modifications to loci of microRNAs in mouse liver and human hepatoma cells. Dual-luciferase reporter assay was conducted to determine effects of HNF4α on the promoters of mouse and human microRNAs as well as effects of microRNAs on the untranslated regions (3’UTR) of two genes in human hepatoma cells.

RESULTS

Microarray data indicated that most microRNAs remained unaltered by Hnf4α deficiency in Hnf4α-LivKO mice. However, certain liver-predominant microRNAs were down-regulated similarly in young-adult male and female Hnf4α-LivKO mice. The down-regulation of miR-101, miR-192, miR-193a, miR-194, miR-215, miR-802, and miR-122 as well as induction of miR-34 and miR-29 in male Hnf4α-LivKO mice were confirmed by real-time PCR. Analysis of public chromatin immunoprecipitation-sequencing data indicates that HNF4α directly binds to the promoters of miR-101, miR-122, miR-194-2/miR-192 and miR-193, which is associated with histone marks of active transcription. Luciferase reporter assay showed that HNF4α markedly activated the promoters of mouse and human miR-101b/miR-101-2 and the miR-194/miR-192 cluster. Additionally, miR-192 and miR-194 significantly decreased activities of luciferase reporters for the 3’UTR of histone H3F3 and chromodomain helicase DNA binding protein 1 (CHD1), respectively, suggesting that miR-192 and miR-194 might be important in chromosome remodeling through directly targeting H3F3 and CHD1.

CONCLUSION

HNF4α is essential for hepatic basal expression of a group of liver-enriched microRNAs, including miR-101, miR-192, miR-193a, miR-194 and miR-802, through which HNF4α may play a major role in the post-transcriptional regulation of gene expression and maintenance of the epigenome in liver.

Retrieving Chromatin Patterns from Deep Sequencing Data Using Correlation Functions

Epigenetic modifications and other chromatin features partition the genome on multiple length scales. They define chromatin domains with distinct biological functions that come in sizes ranging from single modified DNA bases to several megabases in the case of heterochromatic histone modifications. Due to chromatin folding, domains that are well separated along the linear nucleosome chain can form long-range interactions in three-dimensional space. It has now become a routine task to map epigenetic marks and chromatin structure by deep sequencing methods. However, assessing and comparing the properties of chromatin domains and their positional relationships across data sets without a priori assumptions remains challenging. Here, we introduce multiscale correlation evaluation (MCORE), which uses the fluctuation spectrum of mapped sequencing reads to quantify and compare chromatin patterns over a broad range of length scales in a model-independent manner. We applied MCORE to map the chromatin landscape in mouse embryonic stem cells and differentiated neural cells. We integrated sequencing data from chromatin immunoprecipitation, RNA expression, DNA methylation, and chromosome conformation capture experiments into network models that reflect the positional relationships among these features on different genomic scales. Furthermore, we used MCORE to compare our experimental data to models for heterochromatin reorganization during differentiation. The application of correlation functions to deep sequencing data complements current evaluation schemes and will support the development of quantitative descriptions of chromatin networks.

The H3K9 dimethyltransferases EHMT1/2 protect against pathological cardiac hypertrophy

Cardiac hypertrophic growth in response to pathological cues is associated with reexpression of fetal genes and decreased cardiac function and is often a precursor to heart failure. In contrast, physiologically induced hypertrophy is adaptive, resulting in improved cardiac function. The processes that selectively induce these hypertrophic states are poorly understood. Here, we have profiled 2 repressive epigenetic marks, H3K9me2 and H3K27me3, which are involved in stable cellular differentiation, specifically in cardiomyocytes from physiologically and pathologically hypertrophied rat hearts, and correlated these marks with their associated transcriptomes. This analysis revealed the pervasive loss of euchromatic H3K9me2 as a conserved feature of pathological hypertrophy that was associated with reexpression of fetal genes. In hypertrophy, H3K9me2 was reduced following a miR-217-mediated decrease in expression of the H3K9 dimethyltransferases EHMT1 and EHMT2 (EHMT1/2). miR-217-mediated, genetic, or pharmacological inactivation of EHMT1/2 was sufficient to promote pathological hypertrophy and fetal gene reexpression, while suppression of this pathway protected against pathological hypertrophy both in vitro and in mice. Thus, we have established a conserved mechanism involving a departure of the cardiomyocyte epigenome from its adult cellular identity to a reprogrammed state that is accompanied by reexpression of fetal genes and pathological hypertrophy. These results suggest that targeting miR-217 and EHMT1/2 to prevent H3K9 methylation loss is a viable therapeutic approach for the treatment of heart disease.

Hypoxia increases genome-wide bivalent epigenetic marking by specific gain of H3K27me3

Background

Trimethylation at histone H3 lysine 4 (H3K4me3) and lysine 27 (H3K27me3) controls gene activity during development and differentiation. Whether H3K4me3 and H3K27me3 changes dynamically in response to altered microenvironmental conditions, including low-oxygen conditions commonly present in solid tumors, is relatively unknown. Demethylation of H3K4me3 and H3K27me3 is mediated by oxygen and 2-oxoglutarate dioxygenases enzymes, suggesting that oxygen deprivation (hypoxia) may influence histone trimethylation. Using the MCF7 breast epithelial adenocarcinoma cell model, we have determined the relationship between epigenomic and transcriptomic reprogramming as a function of fluctuating oxygen tension.

Results

We find that in MCF7, H3K4me3 and H3K27me3 marks rapidly increase at specific locations throughout the genome and are largely reversed upon reoxygenation. Whereas dynamic changes are relatively highest for H3K27me3 marking under hypoxic conditions, H3K4me3 occupation is identified as the defining epigenetic marker of transcriptional control. In agreement with the global increase of H3K27 trimethylation, we provide direct evidence that the histone H3K27me3 demethylase KDM6B/JMJD3 is inactivated by limited oxygen. In situ immunohistochemical analysis confirms a marked rise of histone trimethylation in hypoxic tumor areas. Acquisition of H3K27me3 at H3K4me3-marked loci results in a striking increase in “bivalent” epigenetic marking. Hypoxia-induced bivalency substantially overlaps with embryonal stem cell-associated genic bivalency and is retained at numerous loci upon reoxygenation. Transcriptional activity is selectively and progressively dampened at bivalently marked loci upon repeated exposure to hypoxia, indicating that this subset of genes uniquely maintains the potential for epigenetic regulation by KDM activity.

Conclusions

These data suggest that dynamic regulation of the epigenetic state within the tumor environment may have important consequences for tumor plasticity and biology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0086-0) contains supplementary material, which is available to authorized users.

Stable Caenorhabditis elegans chromatin domains separate broadly expressed and developmentally regulated genes

Eukaryotic genomes are organized into domains of differing structure and activity. There is evidence that the domain organization of the genome regulates its activity, yet our understanding of domain properties and the factors that influence their formation is poor. Here, we use chromatin state analyses in early embryos and third-larval stage (L3) animals to investigate genome domain organization and its regulation in Caenorhabditis elegans At both stages we find that the genome is organized into extended chromatin domains of high or low gene activity defined by different subsets of states, and enriched for H3K36me3 or H3K27me3, respectively. The border regions between domains contain large intergenic regions and a high density of transcription factor binding, suggesting a role for transcription regulation in separating chromatin domains. Despite the differences in cell types, overall domain organization is remarkably similar in early embryos and L3 larvae, with conservation of 85% of domain border positions. Most genes in high-activity domains are expressed in the germ line and broadly across cell types, whereas low-activity domains are enriched for genes that are developmentally regulated. We find that domains are regulated by the germ-line H3K36 methyltransferase MES-4 and that border regions show striking remodeling of H3K27me1, supporting roles for H3K36 and H3K27 methylation in regulating domain structure. Our analyses of C. elegans chromatin domain structure show that genes are organized by type into domains that have differing modes of regulation.