Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity

Genome-wide identification of regulatory elements in Sertoli cells

A current goal of molecular biology is to identify transcriptional networks that regulate cell differentiation. However, identifying functional gene regulatory elements has been challenging in the context of developing tissues where material is limited and cell types are mixed. To identify regulatory sites during sex determination, we subjected Sertoli cells from mouse fetal testes to DNaseI-seq and ChIP-seq for H3K27ac. DNaseI-seq identified putative regulatory sites around genes enriched in Sertoli and pregranulosa cells; however, active enhancers marked by H3K27ac were enriched proximal to only Sertoli-enriched genes. Sequence analysis identified putative binding sites of known and novel transcription factors likely controlling Sertoli cell differentiation. As a validation of this approach, we identified a novel Sertoli cell enhancer upstream of Wt1, and used it to drive expression of a transgenic reporter in Sertoli cells. This work furthers our understanding of the complex genetic network that underlies sex determination and identifies regions that potentially harbor non-coding mutations underlying disorders of sexual development.

Computational Approaches for Mining GRO-Seq Data to Identify and Characterize Active Enhancers

Transcriptional enhancers are DNA regulatory elements that are bound by transcription factors and act to positively regulate the expression of nearby or distally located target genes. Enhancers have many features that have been discovered using genomic analyses. Recent studies have shown that active enhancers recruit RNA polymerase II (Pol II) and are transcribed, producing enhancer RNAs (eRNAs). GRO-seq, a method for identifying the location and orientation of all actively transcribing RNA polymerases across the genome, is a powerful approach for monitoring nascent enhancer transcription. Furthermore, the unique pattern of enhancer transcription can be used to identify enhancers in the absence of any information about the underlying transcription factors. Here, we describe the computational approaches required to identify and analyze active enhancers using GRO-seq data, including data pre-processing, alignment, and transcript calling. In addition, we describe protocols and computational pipelines for mining GRO-seq data to identify active enhancers, as well as known transcription factor binding sites that are transcribed. Furthermore, we discuss approaches for integrating GRO-seq-based enhancer data with other genomic data, including target gene expression and function. Finally, we describe molecular biology assays that can be used to confirm and explore further the function of enhancers that have been identified using genomic assays. Together, these approaches should allow the user to identify and explore the features and biological functions of new cell type-specific enhancers.

Epigenetic regulation of bone cells

Bone is a major organ in the skeletal system that supports and protects muscles and other organs, facilitates movement and hematopoiesis, and forms a reservoir of minerals including calcium. The cells in the bone, such as osteoblasts, osteoclasts, and osteocytes, orchestrate sequential and balanced regulatory mechanisms to maintain bone and are capable of differentiating in bones. Bone development and remodeling require a precise regulation of gene expressions in bone cells, a process governed by epigenetic mechanisms such as histone modification, DNA methylation, and chromatin structure. Importantly, lineage-specific transcription factors can determine the epigenetic regulation of bone cells. Emerging data suggest that perturbation of epigenetic programs can affect the function and activity of bone cells and contributes to pathogenesis of bone diseases, including osteoporosis. Thus, understanding epigenetic regulations in bone cells would be important for early diagnosis and future therapeutic approaches.

Editing the Neuronal Genome: a CRISPR View of Chromatin Regulation in Neuronal Development, Function, and Plasticity

The dynamic orchestration of gene expression is crucial for the proper differentiation, function, and adaptation of cells. In the brain, transcriptional regulation underlies the incredible diversity of neuronal cell types and contributes to the ability of neurons to adapt their function to the environment. Recently, novel methods for genome and epigenome editing have begun to revolutionize our understanding of gene regulatory mechanisms. In particular, the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system has proven to be a particularly accessible and adaptable technique for genome engineering. Here, we review the use of CRISPR/Cas9 in neurobiology and discuss how these studies have advanced understanding of nervous system development and plasticity. We cover four especially salient applications of CRISPR/Cas9: testing the consequences of enhancer mutations, tagging genes and gene products for visualization in live cells, directly activating or repressing enhancers in vivo, and manipulating the epigenome. In each case, we summarize findings from recent studies and discuss evolving adaptations of the method.

DNA Hypomethylation in Intragenic and Intergenic Enhancer Chromatin of Muscle-Specific Genes Usually Correlates with their Expression

Tissue-specific enhancers are critical for gene regulation. In this study, we help elucidate the contribution of muscle-associated differential DNA methylation to the enhancer activity of highly muscle-specific genes. By bioinformatic analysis of 44 muscle-associated genes, we show that preferential gene expression in skeletal muscle (SkM) correlates with SkM-specific intragenic and intergenic enhancer chromatin and overlapping foci of DNA hypomethylation. Some genes, e.g., CASQ1 and FBXO32, displayed broad regions of both SkM- and heart-specific enhancer chromatin but exhibited focal SkM-specific DNA hypomethylation. Half of the genes had SkM-specific super-enhancers. In contrast to simple enhancer/gene-expression correlations, a super-enhancer was associated with the myogenic MYOD1 gene in both SkM and myoblasts even though SkM has < 1 percent as much MYOD1 expression. Local chromatin differences in this super-enhancer probably contribute to the SkM/myoblast differential expression. Transfection assays confirmed the tissue-specificity of the 0.3-kb core enhancer within MYOD1's super-enhancer and demonstrated its repression by methylation of its three CG dinucleotides. Our study suggests that DNA hypomethylation increases enhancer tissue-specificity and that SkM super-enhancers sometimes are poised for physiologically important, rapid up-regulation.

Identification of the Functional Variant(s) that Explain the Low-Density Lipoprotein Receptor (LDLR) GWAS SNP rs6511720 Association with Lower LDL-C and Risk of CHD

Background

The Low-Density Lipoprotein Receptor (LDLR) SNP rs6511720 (G>T), located in intron-1 of the gene, has been identified in genome-wide association studies (GWAS) as being associated with lower plasma levels of LDL-C and a lower risk of coronary heart disease (CHD). Whether or not rs6511720 is itself functional or a marker for a functional variant elsewhere in the gene is not known.

Methods

The association of LDLR SNP rs6511720 with incidence of CHD and levels of LDL-C was determined by reference to CARDIoGRAM, C4D and Global lipids genetics consortium (GLGC) data. SNP annotation databases were used to identify possible SNP function and prioritization. Luciferase reporter assays in the liver cell line Huh7 were used to measure the effect of variant genotype on gene expression. Electrophoretic Mobility Shift Assays (EMSAs) were used to identify the Transcription Factors (TFs) involved in gene expression regulation.

Results

The phenotype-genotype analysis showed that the rs6511720 minor allele is associated with lower level of LDL-C [beta = -0.2209, p = 3.85 x10^-262], and lower risk of CHD [log (OR) = 0.1155, p = 1.04 x10^-7]. Rs6511720 is in complete linkage. Rs6511720 is in complete linkage disequilibrium (LD) with three intron-1 SNPs (rs141787760, rs60173709, rs57217136). Luciferase reporter assays in Huh7 cells showed that the rare alleles of both rs6511720 and rs57217136 caused a significant increase in LDLR expression compared to the common alleles (+29% and +24%, respectively). Multiplex Competitor-EMSAs (MC-EMSA) identified that the transcription factor serum response element (SRE) binds to rs6511720, while retinoic acid receptor (RAR) and signal transducer and activator of transcription 1 (STAT1) bind to rs57217136.

Conclusion

Both LDLR rs6511720 and rs57217136 are functional variants. Both these minor alleles create enhancer-binding protein sites for TFs and may contribute to increased LDLR expression, which is consequently associated with reduced LDL-C levels and 12% lower CHD risk.

RGmatch: matching genomic regions to proximal genes in omics data integration

Background

The integrative analysis of multiple genomics data often requires that genome coordinates-based signals have to be associated with proximal genes. The relative location of a genomic region with respect to the gene (gene area) is important for functional data interpretation; hence algorithms that match regions to genes should be able to deliver insight into this information.

Results

In this work we review the tools that are publicly available for making region-to-gene associations. We also present a novel method, RGmatch, a flexible and easy-to-use Python tool that computes associations either at the gene, transcript, or exon level, applying a set of rules to annotate each region-gene association with the region location within the gene. RGmatch can be applied to any organism as long as genome annotation is available. Furthermore, we qualitatively and quantitatively compare RGmatch to other tools.

Conclusions

RGmatch simplifies the association of a genomic region with its closest gene. At the same time, it is a powerful tool because the rules used to annotate these associations are very easy to modify according to the researcher’s specific interests. Some important differences between RGmatch and other similar tools already in existence are RGmatch’s flexibility, its wide range of user options, compatibility with any annotatable organism, and its comprehensive and user-friendly output.

Electronic supplementary material

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

Tet2 and Tet3 cooperate with B-lineage transcription factors to regulate DNA modification and chromatin accessibility

Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine, facilitating DNA demethylation and generating new epigenetic marks. Here we show that concomitant loss of Tet2 and Tet3 in mice at early B cell stage blocked the pro- to pre-B cell transition in the bone marrow, decreased Irf4 expression and impaired the germline transcription and rearrangement of the Igκ locus. Tet2/3-deficient pro-B cells showed increased CpG methylation at the Igκ 3' and distal enhancers that was mimicked by depletion of E2A or PU.1, as well as a global decrease in chromatin accessibility at enhancers. Importantly, re-expression of the Tet2 catalytic domain in Tet2/3-deficient B cells resulted in demethylation of the Igκ enhancers and restored their chromatin accessibility. Our data suggest that TET proteins and lineage-specific transcription factors cooperate to influence chromatin accessibility and Igκ enhancer function by modulating the modification status of DNA.

Genome-scale high-resolution mapping of activating and repressive nucleotides in regulatory regions

Massively parallel reporter assays (MPRAs) enable nucleotide-resolution dissection of transcriptional regulatory regions, such as enhancers, but only few regions at a time. Here we present a combined experimental and computational approach, Systematic high-resolution activation and repression profiling with reporter tiling using MPRA (Sharpr-MPRA), that allows high-resolution analysis of thousands of regions simultaneously. Sharpr-MPRA combines dense tiling of overlapping MPRA constructs with a probabilistic graphical model to recognize functional regulatory nucleotides, and to distinguish activating and repressive nucleotides, using their inferred contribution to reporter gene expression. We used Sharpr-MPRA to test 4.6 million nucleotides spanning 15,000 putative regulatory regions tiled at 5-nucleotide resolution in two human cell types. Our results recovered known cell-type-specific regulatory motifs and evolutionarily conserved nucleotides, and distinguished known activating and repressive motifs. Our results also showed that endogenous chromatin state and DNA accessibility are both predictive of regulatory function in reporter assays, identified retroviral elements with activating roles, and uncovered 'attenuator' motifs with repressive roles in active chromatin.

Plant Enhancers: A Call for Discovery

Higher eukaryotes typically contain many different cell types, displaying different cellular functions that are influenced by biotic and abiotic cues. The different functions are characterized by specific gene expression patterns mediated by regulatory sequences such as transcriptional enhancers. Recent genome-wide approaches have identified thousands of enhancers in animals, reviving interest in enhancers in gene regulation. Although the regulatory roles of plant enhancers are as crucial as those in animals, genome-wide approaches have only very recently been applied to plants. Here we review characteristics of enhancers at the DNA and chromatin level in plants and other species, their similarities and differences, and techniques widely used for genome-wide discovery of enhancers in animal systems that can be implemented in plants.

Widespread Chromatin Accessibility at Repetitive Elements Links Stem Cells with Human Cancer

Chromatin regulation is critical for differentiation and disease. However, features linking the chromatin environment of stem cells with disease remain largely unknown. We explored chromatin accessibility in embryonic and multipotent stem cells and unexpectedly identified widespread chromatin accessibility at repetitive elements. Integrating genomic and biochemical approaches, we demonstrate that these sites of increased accessibility are associated with well-positioned nucleosomes marked by distinct histone modifications. Differentiation is accompanied by chromatin remodeling at repetitive elements associated with altered expression of genes in relevant developmental pathways. Remarkably, we found that the chromatin environment of Ewing sarcoma, a mesenchymally derived tumor, is shared with primary mesenchymal stem cells (MSCs). Accessibility at repetitive elements in MSCs offers a permissive environment that is exploited by the critical oncogene responsible for this cancer. Our data demonstrate that stem cells harbor a unique chromatin landscape characterized by accessibility at repetitive elements, a feature associated with differentiation and oncogenesis.

Understanding the genetic liability to schizophrenia through the neuroepigenome

The Psychiatric Genomics Consortium-Schizophrenia Workgroup (PGC-SCZ) recently identified 108 loci associated with increased risk for schizophrenia (SCZ). The vast majority of these variants reside within non-coding sequences of the genome and are predicted to exert their effects by affecting the mechanism of action of cis regulatory elements (CREs), such as promoters and enhancers. Although a number of large-scale collaborative efforts (e.g. ENCODE) have achieved a comprehensive mapping of CREs in human cell lines or tissue homogenates, it is becoming increasingly evident that many risk-associated variants are enriched for expression Quantitative Trait Loci (eQTLs) and CREs in specific tissues or cells. As such, data derived from previous research endeavors may not capture fully cell-type and/or region specific changes associated with brain diseases. Coupling recent technological advances in genomics with cell-type specific methodologies, we are presented with an unprecedented opportunity to better understand the genetics of normal brain development and function and, in turn, the molecular basis of neuropsychiatric disorders. In this review, we will outline ongoing efforts towards this goal and will discuss approaches with the potential to shed light on the mechanism(s) of action of cell-type specific cis regulatory elements and their putative roles in disease, with particular emphasis on understanding the manner in which the epigenome and CREs influence the etiology of SCZ.

Nucleosome assembly and disassembly activity of GRWD1, a novel Cdt1-binding protein that promotes pre-replication complex formation

GRWD1 was previously identified as a novel Cdt1-binding protein that possesses histone-binding and nucleosome assembly activities and promotes MCM loading, probably by maintaining chromatin openness at replication origins. However, the molecular mechanisms underlying these activities remain unknown. We prepared reconstituted mononucleosomes from recombinant histones and a DNA fragment containing a nucleosome positioning sequence, and investigated the effects of GRWD1 on them. GRWD1 could disassemble these preformed mononucleosomes in vitro in an ATP-independent manner. Thus, our data suggest that GRWD1 facilitates removal of H2A-H2B dimers from nucleosomes, resulting in formation of hexasomes. The activity was compromised by deletion of the acidic domain, which is required for efficient histone binding. In contrast, nucleosome assembly activity of GRWD1 was not affected by deletion of the acidic domain. In HeLa cells, the acidic domain of GRWD1 was necessary to maintain chromatin openness and promote MCM loading at replication origins. Taken together, our results suggest that GRWD1 promotes chromatin fluidity by influencing nucleosome structures, e.g., by transient eviction of H2A-H2B, and thereby promotes efficient MCM loading at replication origins.

Punctuated chromatin states regulate Plasmodium falciparum antigenic variation at the intron and 2 kb upstream regions

Background

Understanding the regulation mechanism of var gene expression is crucial for explaining antigenic variation in Plasmodium falciparum. Recent work observed that while all var genes produce transcripts, only a few var genes exhibit high expression levels. However, the global regulation of var expression and the relationship between epigenetic and genetic control remains to be established.

Result

We have systematically reanalyzed the existing genomic data including chromatin configurations and gene expressions; and for the first time used robust statistical methods to show that the intron and 2 kb upstream regions of each endogenous var gene always maintain high chromatin accessibility, with high potential to bind transcription factors (TFs). The levels of transcripts for different var gene family members are associated with this chromatin accessibility. Any given var gene thus shows punctuated chromatin states throughout the asexual life cycle. This is demonstrated by three independent transcript datasets. Chromatin accessibility in the var intron and 2 kb upstream regions are also positively correlated with their GC content, suggesting the level of var genes silencing might be encoded in their intron sequences. Interestingly, both var intron and 2 kb upstream regions exhibit higher chromatin accessibility when the genes have relatively lower transcription levels, suggesting a punctuated repressive function for these regulatory elements.

Conclusion

By integrating and analyzing epigenomic, genomic and transcriptomic data, our work reveals a novel distal element in var control. We found dynamic modulations of specific epigenetic marks around the var intron and distal upstream regions are involved in the general var gene expression patterns in malarial antigenic variation.

Electronic supplementary material

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

Phenotypic Plasticity through Transcriptional Regulation of the Evolutionary Hotspot Gene tan in Drosophila melanogaster

Phenotypic plasticity is the ability of a given genotype to produce different phenotypes in response to distinct environmental conditions. Phenotypic plasticity can be adaptive. Furthermore, it is thought to facilitate evolution. Although phenotypic plasticity is a widespread phenomenon, its molecular mechanisms are only beginning to be unravelled. Environmental conditions can affect gene expression through modification of chromatin structure, mainly via histone modifications, nucleosome remodelling or DNA methylation, suggesting that phenotypic plasticity might partly be due to chromatin plasticity. As a model of phenotypic plasticity, we study abdominal pigmentation of Drosophila melanogaster females, which is temperature sensitive. Abdominal pigmentation is indeed darker in females grown at 18°C than at 29°C. This phenomenon is thought to be adaptive as the dark pigmentation produced at lower temperature increases body temperature. We show here that temperature modulates the expression of tan (t), a pigmentation gene involved in melanin production. t is expressed 7 times more at 18°C than at 29°C in female abdominal epidermis. Genetic experiments show that modulation of t expression by temperature is essential for female abdominal pigmentation plasticity. Temperature modulates the activity of an enhancer of t without modifying compaction of its chromatin or level of the active histone mark H3K27ac. By contrast, the active mark H3K4me3 on the t promoter is strongly modulated by temperature. The H3K4 methyl-transferase involved in this process is likely Trithorax, as we show that it regulates t expression and the H3K4me3 level on the t promoter and also participates in female pigmentation and its plasticity. Interestingly, t was previously shown to be involved in inter-individual variation of female abdominal pigmentation in Drosophila melanogaster, and in abdominal pigmentation divergence between Drosophila species. Sensitivity of t expression to environmental conditions might therefore give more substrate for selection, explaining why this gene has frequently been involved in evolution of pigmentation.

Environmental conditions can strongly modulate the phenotype produced by a particular genotype. This process, called phenotypic plasticity, has major implications in medicine and agricultural sciences, and is thought to facilitate evolution. Phenotypic plasticity is observed in many animals and plants but its mechanisms are only partially understood. As a model of phenotypic plasticity, we study the effect of temperature on female abdominal pigmentation in the fruit fly Drosophila melanogaster. Here we show that temperature affects female abdominal pigmentation by modulating the expression of tan (t), a gene involved in melanin production, in female abdominal epidermis. This effect is mediated at least partly by a particular regulatory sequence of t, the t_MSE enhancer. However we detected no modulation of chromatin structure of t_MSE by temperature. By contrast, the level of the active chromatin mark H3K4me3 on the t promoter is strongly increased at lower temperature. We show that the H3K4 methyl-transferase Trithorax is involved in female abdominal pigmentation and its plasticity and regulates t expression and H3K4me3 level on the t promoter. Several studies have linked t to pigmentation evolution within and between Drosophila species. Our results suggest that sensitivity of t expression to temperature might facilitate its role in pigmentation evolution.

Identification of Regulatory DNA Elements Using Genome-wide Mapping of DNase I Hypersensitive Sites during Tomato Fruit Development

Development and ripening of tomato fruit are precisely controlled by transcriptional regulation, which depends on the orchestrated accessibility of regulatory proteins to promoters and other cis-regulatory DNA elements. This accessibility and its effect on gene expression play a major role in defining the developmental process. To understand the regulatory mechanism and functional elements modulating morphological and anatomical changes during fruit development, we generated genome-wide high-resolution maps of DNase I hypersensitive sites (DHSs) from the fruit tissues of the tomato cultivar "Moneymaker" at 20 days post anthesis as well as break stage. By exploring variation of DHSs across fruit development stages, we pinpointed the most likely hypersensitive sites related to development-specific genes. By detecting binding motifs on DHSs of these development-specific genes or genes in the ascorbic acid biosynthetic pathway, we revealed the common regulatory elements contributing to coordinating gene transcription of plant ripening and specialized metabolic pathways. Our results contribute to a better understanding of the regulatory dynamics of genes involved in tomato fruit development and ripening.

The role of enhancers in cancer

Enhancer elements function as the logic gates of the genetic regulatory circuitry. One of their most important functions is the integration of extracellular signals with intracellular cell fate information to generate cell type-specific transcriptional responses. Mutations occurring in cancer often misregulate enhancers that normally control the signal-dependent expression of growth-related genes. This misregulation can result from trans-acting mechanisms, such as activation of the transcription factors or epigenetic regulators that control enhancer activity, or can be caused in cis by direct mutations that alter the activity of the enhancer or its target gene specificity. These processes can generate tumour type-specific super-enhancers and establish a 'locked' gene regulatory state that drives the uncontrolled proliferation of cancer cells. Here, we review the role of enhancers in cancer, and their potential as therapeutic targets.

A synergistic DNA logic predicts genome-wide chromatin accessibility

Enhancers and promoters commonly occur in accessible chromatin characterized by depleted nucleosome contact; however, it is unclear how chromatin accessibility is governed. We show that log-additive cis-acting DNA sequence features can predict chromatin accessibility at high spatial resolution. We develop a new type of high-dimensional machine learning model, the Synergistic Chromatin Model (SCM), which when trained with DNase-seq data for a cell type is capable of predicting expected read counts of genome-wide chromatin accessibility at every base from DNA sequence alone, with the highest accuracy at hypersensitive sites shared across cell types. We confirm that a SCM accurately predicts chromatin accessibility for thousands of synthetic DNA sequences using a novel CRISPR-based method of highly efficient site-specific DNA library integration. SCMs are directly interpretable and reveal that a logic based on local, nonspecific synergistic effects, largely among pioneer TFs, is sufficient to predict a large fraction of cellular chromatin accessibility in a wide variety of cell types.

Chromatin variation associated with liver metabolism is mediated by transposable elements

Background

Functional regulatory regions in eukaryotic genomes are characterized by the disruption of nucleosomes leading to accessible chromatin. The modulation of chromatin accessibility is one of the key mediators of transcriptional regulation, and variation in chromatin accessibility across individuals has been linked to complex traits and disease susceptibility. While mechanisms responsible for chromatin variation across individuals have been investigated, the overwhelming majority of chromatin variation remains unexplained. Furthermore, the processes through which the variation of chromatin accessibility contributes to phenotypic diversity remain poorly understood.

Results

We profiled chromatin accessibility in liver from seven strains of mice with phenotypic diversity in response to a high-fat/high-sucrose (HF/HS) diet and identified reproducible chromatin variation across the individuals. We found that sites of variable chromatin accessibility were more likely to coincide with particular classes of transposable elements (TEs) than sites with common chromatin signatures. Evolutionarily younger long interspersed nuclear elements (LINEs) are particularly likely to harbor variable chromatin sites. These younger LINEs are enriched for binding sites of immune-associated transcription factors, whereas older LINEs are enriched for liver-specific transcription factors. Genomic region enrichment analysis indicates that variable chromatin sites at TEs may function to regulate liver metabolic pathways. CRISPR-Cas9 deletion of a number of variable chromatin sites at TEs altered expression of nearby metabolic genes. Finally, we show that polymorphism of TEs and differential DNA methylation at TEs can both influence chromatin variation.

Conclusions

Our results demonstrate that specific classes of TEs show variable chromatin accessibility across strains of mice that display phenotypic diversity in response to a HF/HS diet. These results indicate that chromatin variation at TEs is an important contributor to phenotypic variation among populations.

Electronic supplementary material

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

Defining the identity of mouse embryonic dermal fibroblasts

Embryonic dermal fibroblasts in the skin have the exceptional ability to initiate hair follicle morphogenesis and contribute to scarless wound healing. Activation of the Wnt signaling pathway is critical for dermal fibroblast fate selection and hair follicle induction. In humans, mutations in Wnt pathway components and target genes lead to congenital focal dermal hypoplasias with diminished hair. The gene expression signature of embryonic dermal fibroblasts during differentiation and its dependence on Wnt signaling is unknown. Here we applied Shannon entropy analysis to identify the gene expression signature of mouse embryonic dermal fibroblasts. We used available human DNase-seq and histone modification ChiP-seq data on various cell-types to demonstrate that genes in the fibroblast cell identity signature can be epigenetically repressed in other cell-types. We found a subset of the signature genes whose expression is dependent on Wnt/β-catenin activity in vivo. With our approach, we have defined and validated a statistically derived gene expression signature that may mediate dermal fibroblast identity and function in development and disease. genesis 54:415-430, 2016. © 2016 Wiley Periodicals, Inc.

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Grabbing the genome by the NADs

The regions of the genome that interact frequently with the nucleolus have been termed nucleolar-associated domains (NADs). Deep sequencing and DNA-fluorescence in situ hybridization (FISH) experiments have revealed that these domains are enriched for repetitive elements, regions of the inactive X chromosome (Xi), and several RNA polymerase III-transcribed genes. NADs are often marked by chromatin modifications characteristic of heterochromatin, including H3K27me3, H3K9me3, and H4K20me3, and artificial targeting of genes to this area is correlated with reduced expression. It has therefore been hypothesized that NAD localization to the nucleolar periphery contributes to the establishment and/or maintenance of heterochromatic silencing. Recently published studies from several multicellular eukaryotes have begun to reveal the trans-acting factors involved in NAD localization, including the insulator protein CCCTC-binding factor (CTCF), chromatin assembly factor (CAF)-1 subunit p150, several nucleolar proteins, and two long non-coding RNAs (lncRNAs). The mechanisms by which these factors coordinate with one another in regulating NAD localization and/or silencing are still unknown. This review will summarize recently published studies, discuss where additional research is required, and speculate about the mechanistic and functional implications of genome organization around the nucleolus.

High-throughput cis-regulatory element discovery in the vector mosquito Aedes aegypti

Background

Despite substantial progress in mosquito genomic and genetic research, few cis-regulatory elements (CREs), DNA sequences that control gene expression, have been identified in mosquitoes or other non-model insects. Formaldehyde-assisted isolation of regulatory elements paired with DNA sequencing, FAIRE-seq, is emerging as a powerful new high-throughput tool for global CRE discovery. FAIRE results in the preferential recovery of open chromatin DNA fragments that are not bound by nucleosomes, an evolutionarily conserved indicator of regulatory activity, which are then sequenced. Despite the power of the approach, FAIRE-seq has not yet been applied to the study of non-model insects. In this investigation, we utilized FAIRE-seq to profile open chromatin and identify likely regulatory elements throughout the genome of the human disease vector mosquito Aedes aegypti. We then assessed genetic variation in the regulatory elements of dengue virus susceptible (Moyo-S) and refractory (Moyo-R) mosquito strains.

Results

Analysis of sequence data obtained through next generation sequencing of FAIRE DNA isolated from A. aegypti embryos revealed >121,000 FAIRE peaks (FPs), many of which clustered in the 1 kb 5’ upstream flanking regions of genes known to be expressed at this stage. As expected, known transcription factor consensus binding sites were enriched in the FPs, and of these FoxA1, Hunchback, Gfi, Klf4, MYB/ph3 and Sox9 are most predominant. All of the elements tested in vivo were confirmed to drive gene expression in transgenic Drosophila reporter assays. Of the >13,000 single nucleotide polymorphisms (SNPs) recently identified in dengue virus-susceptible and refractory mosquito strains, 3365 were found to map to FPs.

Conclusion

FAIRE-seq analysis of open chromatin in A. aegypti permitted genome-wide discovery of CREs. The results of this investigation indicate that FAIRE-seq is a powerful tool for identification of regulatory DNA in the genomes of non-model organisms, including human disease vector mosquitoes.

Electronic supplementary material

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

Distinct transcription factor complexes act on a permissive chromatin landscape to establish regionalized gene expression in CNS stem cells

Spatially distinct gene expression profiles in neural stem cells (NSCs) are a prerequisite to the formation of neuronal diversity, but how these arise from the regulatory interactions between chromatin accessibility and transcription factor activity has remained unclear. Here, we demonstrate that, despite their distinct gene expression profiles, NSCs of the mouse cortex and spinal cord share the majority of their DNase I hypersensitive sites (DHSs). Regardless of this similarity, domain-specific gene expression is highly correlated with the relative accessibility of associated DHSs, as determined by sequence read density. Notably, the binding pattern of the general NSC transcription factor SOX2 is also largely cell type specific and coincides with an enrichment of LHX2 motifs in the cortex and HOXA9 motifs in the spinal cord. Interestingly, in a zebrafish reporter gene system, these motifs were critical determinants of patterned gene expression along the rostral-caudal axis. Our findings establish a predictive model for patterned NSC gene expression, whereby domain-specific expression of LHX2 and HOX proteins act on their target motifs within commonly accessible cis-regulatory regions to specify SOX2 binding. In turn, this binding correlates strongly with these DHSs relative accessibility—a robust predictor of neighboring gene expression.

Network-Based Interpretation of Diverse High-Throughput Datasets through the Omics Integrator Software Package

High-throughput, ‘omic’ methods provide sensitive measures of biological responses to perturbations. However, inherent biases in high-throughput assays make it difficult to interpret experiments in which more than one type of data is collected. In this work, we introduce Omics Integrator, a software package that takes a variety of ‘omic’ data as input and identifies putative underlying molecular pathways. The approach applies advanced network optimization algorithms to a network of thousands of molecular interactions to find high-confidence, interpretable subnetworks that best explain the data. These subnetworks connect changes observed in gene expression, protein abundance or other global assays to proteins that may not have been measured in the screens due to inherent bias or noise in measurement. This approach reveals unannotated molecular pathways that would not be detectable by searching pathway databases. Omics Integrator also provides an elegant framework to incorporate not only positive data, but also negative evidence. Incorporating negative evidence allows Omics Integrator to avoid unexpressed genes and avoid being biased toward highly-studied hub proteins, except when they are strongly implicated by the data. The software is comprised of two individual tools, Garnet and Forest, that can be run together or independently to allow a user to perform advanced integration of multiple types of high-throughput data as well as create condition-specific subnetworks of protein interactions that best connect the observed changes in various datasets. It is available at http://fraenkel.mit.edu/omicsintegrator and on GitHub at https://github.com/fraenkel-lab/OmicsIntegrator.

HNF1 regulates critical processes in the human epididymis epithelium

The luminal environment of the epididymis participates in sperm maturation and impacts male fertility. It is dependent on the coordinated expression of many genes encoding proteins with a role in epithelial transport. We identified cis-regulatory elements for critical genes in epididymis function, by mapping open chromatin genome-wide in human epididymis epithelial (HEE) cells. Bioinformatic predictions of transcription factors binding to the regulatory elements suggested an important role for hepatocyte nuclear factor 1 (HNF1) in the transcriptional program of these cells. Chromatin immunoprecipitation and deep sequencing (ChIP-seq) revealed HNF1 target genes in HEE cells. In parallel, the contribution of HNF1 to the transcriptome of HEE cells was determined by RNA-seq, following siRNA-mediated depletion of both HNF1α and HNF1β transcription factors. Repression of these factors caused differential expression of 1892 transcripts (902 were downregulated and 990 upregulated) in comparison to non-targeting siRNAs. Differentially expressed genes with HNF1 ChIP-seq peaks within 20 kb were subject to gene ontology process enrichment analysis. Among the most significant processes associated with down-regulated genes were epithelial transport of water, phosphate and bicarbonate, all critical processes in epididymis epithelial function. Measurements of intracellular pH (pHi) confirmed a role for HNF1 in regulating the epididymis luminal environment.

Genome Wide Mapping of NR4A Binding Reveals Cooperativity with ETS Factors to Promote Epigenetic Activation of Distal Enhancers in Acute Myeloid Leukemia Cells

Members of the NR4A subfamily of orphan nuclear receptors regulate cell fate decisions via both genomic and non-genomic mechanisms in a cell and tissue selective manner. NR4As play a key role in maintenance of hematopoietic stem cell homeostasis and are critical tumor suppressors of acute myeloid leukemia (AML). Expression of NR4As is broadly silenced in leukemia initiating cell enriched populations from human patients relative to normal hematopoietic stem/progenitor cells. Rescue of NR4A expression in human AML cells inhibits proliferation and reprograms AML gene signatures via transcriptional mechanisms that remain to be elucidated. By intersecting an acutely regulated NR4A1 dependent transcriptional profile with genome wide NR4A binding distribution, we now identify an NR4A targetome of 685 genes that are directly regulated by NR4A1. We show that NR4As regulate gene transcription primarily through interaction with distal enhancers that are co-enriched for NR4A1 and ETS transcription factor motifs. Using a subset of NR4A activated genes, we demonstrate that the ETS factors ERG and FLI-1 are required for activation of NR4A bound enhancers and NR4A target gene induction. NR4A1 dependent recruitment of ERG and FLI-1 promotes binding of p300 histone acetyltransferase to epigenetically activate NR4A bound enhancers via acetylation at histone H3K27. These findings disclose novel epigenetic mechanisms by which NR4As and ETS factors cooperate to drive NR4A dependent gene transcription in human AML cells.

High-throughput small molecule screen identifies inhibitors of aberrant chromatin accessibility

Mutations in chromatin-modifying proteins and transcription factors are commonly associated with a wide variety of cancers. Through gain- or loss-of-function, these mutations may result in characteristic alterations of accessible chromatin, indicative of shifts in the landscape of regulatory elements genome-wide. The identification of compounds that reverse a specific chromatin signature could lead to chemical probes or potential therapies. To explore whether chromatin accessibility could serve as a platform for small molecule screening, we adapted formaldehyde-assisted isolation of regulatory elements (FAIRE), a chemical method to enrich for nucleosome-depleted genomic regions, as a high-throughput, automated assay. After demonstrating the validity and robustness of this approach, we applied this method to screen an epigenetically targeted small molecule library by evaluating regions of aberrant nucleosome depletion mediated by EWSR1-FLI1, the chimeric transcription factor critical for the bone and soft tissue tumor Ewing sarcoma. As a class, histone deacetylase inhibitors were greatly overrepresented among active compounds. These compounds resulted in diminished accessibility at targeted sites by disrupting transcription of EWSR1-FLI1. Capitalizing on precise differences in chromatin accessibility for drug discovery efforts offers significant advantages because it does not depend on the a priori selection of a single molecular target and may detect novel biologically relevant pathways.

Characterization of chromatin accessibility with a transposome hypersensitive sites sequencing (THS-seq) assay

Chromatin accessibility captures in vivo protein-chromosome binding status, and is considered an informative proxy for protein-DNA interactions. DNase I and Tn5 transposase assays require thousands to millions of fresh cells for comprehensive chromatin mapping. Applying Tn5 tagmentation to hundreds of cells results in sparse chromatin maps. We present a transposome hypersensitive sites sequencing assay for highly sensitive characterization of chromatin accessibility. Linear amplification of accessible DNA ends with in vitro transcription, coupled with an engineered Tn5 super-mutant, demonstrates improved sensitivity on limited input materials, and accessibility of small regions near distal enhancers, compared with ATAC-seq.

Small Genetic Circuits and MicroRNAs: Big Players in Polymerase II Transcriptional Control in Plants

RNA Polymerase II (Pol II) regulatory cascades involving transcription factors (TFs) and their targets orchestrate the genetic circuitry of every eukaryotic organism. In order to understand how these cascades function, they can be dissected into small genetic networks, each containing just a few Pol II transcribed genes, that generate specific signal-processing outcomes. Small RNA regulatory circuits involve direct regulation of a small RNA by a TF and/or direct regulation of a TF by a small RNA and have been shown to play unique roles in many organisms. Here, we will focus on small RNA regulatory circuits containing Pol II transcribed microRNAs (miRNAs). While the role of miRNA-containing regulatory circuits as modular building blocks for the function of complex networks has long been on the forefront of studies in the animal kingdom, plant studies are poised to take a lead role in this area because of their advantages in probing transcriptional and posttranscriptional control of Pol II genes. The relative simplicity of tissue- and cell-type organization, miRNA targeting, and genomic structure make the Arabidopsis thaliana plant model uniquely amenable for small RNA regulatory circuit studies in a multicellular organism. In this Review, we cover analysis, tools, and validation methods for probing the component interactions in miRNA-containing regulatory circuits. We then review the important roles that plant miRNAs are playing in these circuits and summarize methods for the identification of small genetic circuits that strongly influence plant function. We conclude by noting areas of opportunity where new plant studies are imminently needed.

Core pathway mutations induce de-differentiation of murine astrocytes into glioblastoma stem cells that are sensitive to radiation but resistant to temozolomide

BACKGROUND

Glioma stem cells (GSCs) from human glioblastomas (GBMs) are resistant to radiation and chemotherapy and may drive recurrence. Treatment efficacy may depend on GSCs, expression of DNA repair enzymes such as methylguanine methyltransferase (MGMT), or transcriptome subtype.

METHODS

To model genetic alterations in human GBM core signaling pathways, we induced Rb knockout, Kras activation, and Pten deletion mutations in cortical murine astrocytes. Neurosphere culture, differentiation, and orthotopic transplantation assays were used to assess whether these mutations induced de-differentiation into GSCs. Genome-wide chromatin landscape alterations and expression profiles were examined by formaldehyde-assisted isolation of regulatory elements (FAIRE) seq and RNA-seq. Radiation and temozolomide efficacy were examined in vitro and in an allograft model in vivo. Effects of radiation on transcriptome subtype were examined by microarray expression profiling.

RESULTS

Cultured triple mutant astrocytes gained unlimited self-renewal and multilineage differentiation capacity. These cells harbored significantly altered chromatin landscapes that were associated with downregulation of astrocyte- and upregulation of stem cell-associated genes, particularly the Hoxa locus of embryonic transcription factors. Triple-mutant astrocytes formed serially transplantable glioblastoma allografts that were sensitive to radiation but expressed MGMT and were resistant to temozolomide. Radiation induced a shift in transcriptome subtype of GBM allografts from proneural to mesenchymal.

CONCLUSION

A defined set of core signaling pathway mutations induces de-differentiation of cortical murine astrocytes into GSCs with altered chromatin landscapes and transcriptomes. This non-germline genetically engineered mouse model mimics human proneural GBM on histopathological, molecular, and treatment response levels. It may be useful for dissecting the mechanisms of treatment resistance and developing more effective therapies.

Nucleosome architecture throughout the cell cycle

Nucleosomes provide additional regulatory mechanisms to transcription and DNA replication by mediating the access of proteins to DNA. During the cell cycle chromatin undergoes several conformational changes, however the functional significance of these changes to cellular processes are largely unexplored. Here, we present the first comprehensive genome-wide study of nucleosome plasticity at single base-pair resolution along the cell cycle in Saccharomyces cerevisiae. We determined nucleosome organization with a specific focus on two regulatory regions: transcription start sites (TSSs) and replication origins (ORIs). During the cell cycle, nucleosomes around TSSs display rearrangements in a cyclic manner. In contrast to gap (G1 and G2) phases, nucleosomes have a fuzzier organization during S and M phases, Moreover, the choreography of nucleosome rearrangements correlate with changes in gene expression during the cell cycle, indicating a strong association between nucleosomes and cell cycle-dependent gene functionality. On the other hand, nucleosomes are more dynamic around ORIs along the cell cycle, albeit with tighter regulation in early firing origins, implying the functional role of nucleosomes on replication origins. Our study provides a dynamic picture of nucleosome organization throughout the cell cycle and highlights the subsequent impact on transcription and replication activity.

Genome-Wide Epigenetic Studies in Human Disease: A Primer on -Omic Technologies

Epigenetic information encoded in covalent modifications of DNA and histone proteins regulates fundamental biological processes through the action of chromatin regulators, transcription factors, and noncoding RNA species. Epigenetic plasticity enables an organism to respond to developmental and environmental signals without genetic changes. However, aberrant epigenetic control plays a key role in pathogenesis of disease. Normal epigenetic states could be disrupted by detrimental mutations and expression alteration of chromatin regulators or by environmental factors. In this primer, we briefly review the epigenetic basis of human disease and discuss how recent discoveries in this field could be translated into clinical diagnosis, prevention, and treatment. We introduce platforms for mapping genome-wide chromatin accessibility, nucleosome occupancy, DNA-binding proteins, and DNA methylation, primarily focusing on the integration of DNA methylation and chromatin immunoprecipitation-sequencing technologies into disease association studies. We highlight practical considerations in applying high-throughput epigenetic assays and formulating analytical strategies. Finally, we summarize current challenges in sample acquisition, experimental procedures, data analysis, and interpretation and make recommendations on further refinement in these areas. Incorporating epigenomic testing into the clinical research arsenal will greatly facilitate our understanding of the epigenetic basis of disease and help identify novel therapeutic targets.

Mapping nucleosome positions using DNase-seq

Although deoxyribonuclease I (DNase I) was used to probe the structure of the nucleosome in the 1960s and 1970s, in the current high-throughput sequencing era, DNase I has mainly been used to study genomic regions devoid of nucleosomes. Here, we reveal for the first time that DNase I can be used to precisely map the (translational) positions of in vivo nucleosomes genome-wide. Specifically, exploiting a distinctive DNase I cleavage profile within nucleosome-associated DNA—including a signature 10.3 base pair oscillation that corresponds to accessibility of the minor groove as DNA winds around the nucleosome—we develop a Bayes-factor–based method that can be used to map nucleosome positions along the genome. Compared to methods that require genetically modified histones, our DNase-based approach is easily applied in any organism, which we demonstrate by producing maps in yeast and human. Compared to micrococcal nuclease (MNase)-based methods that map nucleosomes based on cuts in linker regions, we utilize DNase I cuts both outside and within nucleosomal DNA; the oscillatory nature of the DNase I cleavage profile within nucleosomal DNA enables us to identify translational positioning details not apparent in MNase digestion of linker DNA. Because the oscillatory pattern corresponds to nucleosome rotational positioning, it also reveals the rotational context of transcription factor (TF) binding sites. We show that potential binding sites within nucleosome-associated DNA are often centered preferentially on an exposed major or minor groove. This preferential localization may modulate TF interaction with nucleosome-associated DNA as TFs search for binding sites.

Epigenetics of the myotonic dystrophy-associated DMPK gene neighborhood

Aim:

Identify epigenetic marks in the vicinity of DMPK (linked to myotonic dystrophy, DM1) that help explain tissue-specific differences in its expression.

Materials & methods:

At DMPK and its flanking genes (DMWD, SIX5, BHMG1 and RSPH6A), we analyzed many epigenetic and transcription profiles from myoblasts, myotubes, skeletal muscle, heart and 30 nonmuscle samples.

Results:

In the DMPK gene neighborhood, muscle-associated DNA hypermethylation and hypomethylation, enhancer chromatin, and CTCF binding were seen. Myogenic DMPK hypermethylation correlated with high expression and decreased alternative promoter usage. Testis/sperm hypomethylation of BHMG1 and RSPH6A was associated with testis-specific expression. G-quadruplex (G4) motifs and sperm-specific hypomethylation were found near the DM1-linked CTG repeats within DMPK.

Conclusion:

Tissue-specific epigenetic features in DMPK and neighboring genes help regulate its expression. G4 motifs in DMPK DNA and RNA might contribute to DM1 pathology.

A STATISTICAL MODEL TO ASSESS (ALLELE-SPECIFIC) ASSOCIATIONS BETWEEN GENE EXPRESSION AND EPIGENETIC FEATURES USING SEQUENCING DATA

Sequencing techniques have been widely used to assess gene expression (i.e., RNA-seq) or the presence of epigenetic features (e.g., DNase-seq to identify open chromatin regions). In contrast to traditional microarray platforms, sequencing data are typically summarized in the form of discrete counts, and they are able to delineate allele-specific signals, which are not available from microarrays. The presence of epigenetic features are often associated with gene expression, both of which have been shown to be affected by DNA polymorphisms. However, joint models with the flexibility to assess interactions between gene expression, epigenetic features and DNA polymorphisms are currently lacking. In this paper, we develop a statistical model to assess the associations between gene expression and epigenetic features using sequencing data, while explicitly modeling the effects of DNA polymorphisms in either an allele-specific or nonallele-specific manner. We show that in doing so we provide the flexibility to detect associations between gene expression and epigenetic features, as well as conditional associations given DNA polymorphisms. We evaluate the performance of our method using simulations and apply our method to study the association between gene expression and the presence of DNase I Hypersensitive sites (DHSs) in HapMap individuals. Our model can be generalized to exploring the relationships between DNA polymorphisms and any two types of sequencing experiments, a useful feature as the variety of sequencing experiments continue to expand.

Genome-Scale Analysis of Cell-Specific Regulatory Codes Using Nuclear Enzymes

High-throughput sequencing technologies have made it possible for biologists to generate genome-wide profiles of chromatin features at the nucleotide resolution. Enzymes such as nucleases or transposes have been instrumental as a chromatin-probing agent due to their ability to target accessible chromatin for cleavage or insertion. On the scale of a few hundred base pairs, preferential action of the nuclear enzymes on accessible chromatin allows mapping of cell state-specific accessibility in vivo. Such accessible regions contain functionally important regulatory sites, including promoters and enhancers, which undergo active remodeling for cells adapting in a dynamic environment. DNase-seq and the more recent ATAC-seq are two assays that are gaining popularity. Deep sequencing of DNA libraries from these assays, termed genomic footprinting, has been proposed to enable the comprehensive construction of protein occupancy profiles over the genome at the nucleotide level. Recent studies have discovered limitations of genomic footprinting which reduce the scope of detectable proteins. In addition, the identification of putative factors that bind to the observed footprints remains challenging. Despite these caveats, the methodology still presents significant advantages over alternative techniques such as ChIP-seq or FAIRE-seq. Here we describe computational approaches and tools for analysis of chromatin accessibility and genomic footprinting. Proper experimental design and assay-specific data analysis ensure the detection sensitivity and maximize retrievable information. The enzyme-based chromatin profiling approaches represent a powerful and evolving methodology which facilitates our understanding of how the genome is regulated.

Genome-wide identification of hypoxia-induced enhancer regions

Here we present a genome-wide method for de novo identification of enhancer regions. This approach enables massively parallel empirical investigation of DNA sequences that mediate transcriptional activation and provides a platform for discovery of regulatory modules capable of driving context-specific gene expression. The method links fragmented genomic DNA to the transcription of randomer molecule identifiers and measures the functional enhancer activity of the library by massively parallel sequencing. We transfected a Drosophila melanogaster library into S2 cells in normoxia and hypoxia, and assayed 4,599,881 genomic DNA fragments in parallel. The locations of the enhancer regions strongly correlate with genes up-regulated after hypoxia and previously described enhancers. Novel enhancer regions were identified and integrated with RNAseq data and transcription factor motifs to describe the hypoxic response on a genome-wide basis as a complex regulatory network involving multiple stress-response pathways. This work provides a novel method for high-throughput assay of enhancer activity and the genome-scale identification of 31 hypoxia-activated enhancers in Drosophila.

A novel ATAC-seq approach reveals lineage-specific reinforcement of the open chromatin landscape via cooperation between BAF and p63

Background

Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation.

Results

Using a novel “on-plate” ATAC-sequencing approach for profiling open chromatin landscapes with a low number of adherent cells, we demonstrate that the BAF complex is essential for maintaining 11.6 % of open chromatin regions in epidermal differentiation. These BAF-dependent open chromatin regions are highly cell-type-specific and are strongly enriched for binding sites for p63, a master epidermal transcription factor. The DNA sequences of p63 binding sites intrinsically favor nucleosome formation and are inaccessible in other cell types without p63 to prevent ectopic activation. In epidermal cells, BAF and p63 mutually recruit each other to maintain 14,853 open chromatin regions. We further demonstrate that BAF and p63 cooperatively position nucleosomes away from p63 binding sites and recruit transcriptional machinery to control tissue differentiation.

Conclusions

BAF displays high specificity in controlling the open chromatin landscape during epidermal differentiation by cooperating with the master transcription factor p63 to maintain lineage-specific open chromatin regions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0840-9) contains supplementary material, which is available to authorized users.

Functionally distinct patterns of nucleosome remodeling at enhancers in glucocorticoid-treated acute lymphoblastic leukemia

Background

Precise nucleosome positioning is an increasingly recognized feature of promoters and enhancers, reflecting complex contributions of DNA sequence, nucleosome positioning, histone modification and transcription factor binding to enhancer activity and regulation of gene expression. Changes in nucleosome position and occupancy, histone variants and modifications, and chromatin remodeling are also critical elements of dynamic transcriptional regulation, but poorly understood at enhancers. We investigated glucocorticoid receptor-associated (GR) nucleosome dynamics at enhancers in acute lymphoblastic leukemia.

Results

For the first time, we demonstrate functionally distinct modes of nucleosome remodeling upon chromatin binding by GR, which we term central, non-central, phased, and minimal. Central and non-central remodeling reflect nucleosome eviction by GR and cofactors, respectively. Phased remodeling involves nucleosome repositioning and is associated with rapidly activated enhancers and induction of gene expression. Minimal remodeling sites initially have low levels of enhancer-associated histone modification, but the majority of these regions gain H3K4me2 or H3K27Ac to become de novo enhancers. Minimal remodeling regions are associated with gene ontologies specific to decreased B cell number and mTOR inhibition and may make unique contributions to glucocorticoid-induced leukemia cell death.

Conclusions

Our findings form a novel framework for understanding the dynamic interplay between transcription factor binding, nucleosome remodeling, enhancer function, and gene expression in the leukemia response to glucocorticoids.

Electronic supplementary material

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

Genome-wide predictors of NF-κB recruitment and transcriptional activity

BACKGROUND

Inducible transcription factors (TFs) mediate transcriptional responses to environmental cues. In response to multiple inflammatory signals active NF-κB dimers enter the nucleus and trigger cell-type-, and stimulus-specific transcriptional programs. Although much is known about NF-κB inducing pathways and about locus-specific mechanisms of transcriptional control, it is poorly understood how the pre-existing chromatin landscape determines NF-κB target selection and activation. Specifically, it is not known which epigenetic marks and pre-bound TFs serve genome-wide as positive (negative) cues for active NF-κB.

RESULTS

We applied multivariate and combinatorial data mining techniques on a comprehensive dataset of DNA methylation, DNase I hypersensitivity, eight epigenetic marks, and 34 TFs to arrive at genome-wide patterns that predict NF-κB binding. Strikingly, we observed NF-κB recruitment to accessible and nucleosome-bound sites. Within nucleosomal DNA NF-κB binding was primed by H3K4me1 and H2A.Z, but also hyper-methylated DNA outside of promoters and CpG-islands. Many of these predictors showed combinatorial cooperativity and statistically significant interactions. Recruitment to pre-accessible sites was more frequent and influenced by chromatin-associated TFs. We observed that specific TF-combinations are greatly enriched for (or depleted of) NF-κB binding events.

CONCLUSIONS

We provide evidence of NF-κB binding within genomic regions that lack classical marks of activity. These pioneer binding events are relatively often associated with transcriptional regulation. Further, our predictive models indicate that specific combinations of epigenetic marks and transcription factors predetermine the NF-κB cistrome, supporting the feasibility of using statistical approaches to identify "histone codes".

Super Enhancers in Cancers, Complex Disease, and Developmental Disorders

Recently, unique areas of transcriptional regulation termed super-enhancers have been identified and implicated in human disease. Defined by their magnitude of size, transcription factor density, and binding of transcriptional machinery, super-enhancers have been associated with genes driving cell differentiation. While their functions are not completely understood, it is clear that these regions driving high-level transcription are susceptible to perturbation, and trait-associated single nucleotide polymorphisms (SNPs) occur within super-enhancers of disease-relevant cell types. Here we review evidence for super-enhancer involvement in cancers, complex diseases, and developmental disorders and discuss interactions between super-enhancers and cofactors/chromatin regulators.

Intersection of genetics and epigenetics in monozygotic twin genomes

As a final function of various epigenetic mechanisms, chromatin regulation is a transcription control process that especially demonstrates active interaction with genetic elements. Thus, chromatin structure has become a principal focus in recent genomics researches that strive to characterize regulatory functions of DNA variants related to diseases or other traits. Although researchers have been focusing on DNA methylation when studying monozygotic (MZ) twins, a great model in epigenetics research, interactions between genetics and epigenetics in chromatin level are expected to be an imperative research trend in the future. In this review, we discuss how the genome, epigenome, and transcriptome of MZ twins can be studied in an integrative manner from this perspective.

Identification of Predictive Cis-Regulatory Elements Using a Discriminative Objective Function and a Dynamic Search Space

The generation of genomic binding or accessibility data from massively parallel sequencing technologies such as ChIP-seq and DNase-seq continues to accelerate. Yet state-of-the-art computational approaches for the identification of DNA binding motifs often yield motifs of weak predictive power. Here we present a novel computational algorithm called MotifSpec, designed to find predictive motifs, in contrast to over-represented sequence elements. The key distinguishing feature of this algorithm is that it uses a dynamic search space and a learned threshold to find discriminative motifs in combination with the modeling of motifs using a full PWM (position weight matrix) rather than k-mer words or regular expressions. We demonstrate that our approach finds motifs corresponding to known binding specificities in several mammalian ChIP-seq datasets, and that our PWMs classify the ChIP-seq signals with accuracy comparable to, or marginally better than motifs from the best existing algorithms. In other datasets, our algorithm identifies novel motifs where other methods fail. Finally, we apply this algorithm to detect motifs from expression datasets in C. elegans using a dynamic expression similarity metric rather than fixed expression clusters, and find novel predictive motifs.

Transcriptional enhancers: functional insights and role in human disease

In recent years, studies of cis-regulatory mechanisms have evolved from a predominant focus on promoter regions to the realization that spatial and temporal gene regulation is frequently driven by long-range enhancer clusters that operate within chromosomal compartments. This increased understanding of genome function, together with the emergence of technologies that enable whole-genome sequencing of patients’ DNAs, open the prospect of dissecting the role of cis-regulatory defects in human disease. In this review we discuss how recent epigenomic studies have provided insights into the function of transcriptional enhancers. We then present examples that illustrate how integrative genomics can help uncover enhancer sequence variants underlying Mendelian and common polygenic human disease.

Comparative FAIRE-seq analysis reveals distinguishing features of the chromatin structure of ground state- and primed-pluripotent cells

Both pluripotent embryonic stem cells (ESCs), established from preimplantation murine blastocysts, and epiblast stem cells (EpiSCs), established from postimplantation embryos, can self-renew in culture or differentiate into each of the primary germ layers. While the core transcription factors (TFs) OCT4, SOX2, and NANOG are expressed in both cell types, the gene expression profiles and other features suggest that ESCs and EpiSCs reflect distinct developmental maturation stages of the epiblast in vivo. Accordingly, "naïve" or "ground state" ESCs resemble cells of the inner cell mass, whereas "primed" EpiSCs resemble cells of the postimplantation egg cylinder. To gain insight into the relationship between naïve and primed pluripotent cells, and of each of these pluripotent states to that of nonpluripotent cells, we have used FAIRE-seq to generate a comparative atlas of the accessible chromatin regions within ESCs, EpiSCs, multipotent neural stem cells, and mouse embryonic fibroblasts. We find a distinction between the accessible chromatin patterns of pluripotent and somatic cells that is consistent with the highly related phenotype of ESCs and EpiSCs. However, by defining cell-specific and shared regions of open chromatin, and integrating these data with published gene expression and ChIP analyses, we also illustrate unique features of the chromatin of naïve and primed cells. Functional studies suggest that multiple stage-specific enhancers regulate ESC- or EpiSC-specific gene expression, and implicate auxiliary TFs as important modulators for stage-specific activation by the core TFs. Together these observations provide insights into the chromatin structure dynamics accompanying transitions between these pluripotent states.