Nuclease hypersensitive sites in chromatin

The histone variant H2A.Z is an important regulator of enhancer activity

Gene regulatory programs in different cell types are largely defined through cell-specific enhancers activity. The histone variant H2A.Z has been shown to play important roles in transcription mainly by controlling proximal promoters, but its effect on enhancer functions remains unclear. Here, we demonstrate by genome-wide approaches that H2A.Z is present at a subset of active enhancers bound by the estrogen receptor alpha (ERα). We also determine that H2A.Z does not influence the local nucleosome positioning around ERα enhancers using ChIP sequencing at nucleosomal resolution and unsupervised pattern discovery. We further highlight that H2A.Z-enriched enhancers are associated with chromatin accessibility, H3K122ac enrichment and hypomethylated DNA. Moreover, upon estrogen stimulation, the enhancers occupied by H2A.Z produce enhancer RNAs (eRNAs), and recruit RNA polymerase II as well as RAD21, a member of the cohesin complex involved in chromatin interactions between enhancers and promoters. Importantly, their recruitment and eRNAs production are abolished by H2A.Z depletion, thereby revealing a novel functional link between H2A.Z occupancy and enhancer activity. Taken together, our findings suggest that H2A.Z acts as an important player for enhancer functions by establishing and maintaining a chromatin environment required for RNA polymerase II recruitment, eRNAs transcription and enhancer-promoters interactions, all essential attributes of enhancer activity.

Mining the coding and non-coding genome for cancer drivers

Progress in next-generation sequencing provides unprecedented opportunities to fully characterize the spectrum of somatic mutations of cancer genomes. Given the large number of somatic mutations identified by such technologies, the prioritization of cancer-driving events is a consistent bottleneck. Most bioinformatics tools concentrate on driver mutations in the coding fraction of the genome, those causing changes in protein products. As more non-coding pathogenic variants are identified and characterized, the development of computational approaches to effectively prioritize cancer-driving variants within the non-coding fraction of human genome is becoming critical. After a short summary of methods for coding variant prioritization, we here review the highly diverse non-coding elements that may act as cancer drivers and describe recent methods that attempt to evaluate the deleteriousness of sequence variation in these elements. With such tools, the prioritization and identification of cancer-implicated regulatory elements and non-coding RNAs is becoming a reality.

Tissue-Specific Enrichment of Lymphoma Risk Loci in Regulatory Elements

Though numerous polymorphisms have been associated with risk of developing lymphoma, how these variants function to promote tumorigenesis is poorly understood. Here, we report that lymphoma risk SNPs, especially in the non-Hodgkin's lymphoma subtype chronic lymphocytic leukemia, are significantly enriched for co-localization with epigenetic marks of active gene regulation. These enrichments were seen in a lymphoid-specific manner for numerous ENCODE datasets, including DNase-hypersensitivity as well as multiple segmentation-defined enhancer regions. Furthermore, we identify putatively functional SNPs that are both in regulatory elements in lymphocytes and are associated with gene expression changes in blood. We developed an algorithm, UES, that uses a Monte Carlo simulation approach to calculate the enrichment of previously identified risk SNPs in various functional elements. This multiscale approach integrating multiple datasets helps disentangle the underlying biology of lymphoma, and more broadly, is generally applicable to GWAS results from other diseases as well.

Epigenomic mapping and effect sizes of noncoding variants associated with psychotropic drug response

AIM

To provide insight into potential regulatory mechanisms of gene expression underlying addiction, analgesia, psychotropic drug response and adverse drug events, genome-wide association studies searching for variants associated with these phenotypes has been undertaken with limited success. We undertook analysis of these results with the aim of applying epigenetic knowledge to aid variant discovery and interpretation.

METHODS

We applied conditional imputation to results from 26 genome-wide association studies and three candidate gene-association studies. The analysis workflow included data from chromatin conformation capture, chromatin state annotation, DNase I hypersensitivity, hypomethylation, anatomical localization and biochronicity. We also made use of chromatin state data from the epigenome roadmap, transcription factor-binding data, spatial maps from published Hi-C datasets and 'guilt by association' methods.

RESULTS

We identified 31 pharmacoepigenomic SNPs from a total of 2024 variants in linkage disequilibrium with lead SNPs, of which only 6% were coding variants. Interrogation of chromatin state using our workflow and the epigenome roadmap showed agreement on 34 of 35 tissue assignments to regulatory elements including enhancers and promoters. Loop boundary domains were inferred by association with CTCF (CCCTC-binding factor) and cohesin, suggesting proximity to topologically associating domain boundaries and enhancer clusters. Spatial interactions between enhancer-promoter pairs detected both known and previously unknown mechanisms. Addiction and analgesia SNPs were common in relevant populations and exhibited large effect sizes, whereas a SNP located in the promoter of the SLC1A2 gene exhibited a moderate effect size for lithium response in bipolar disorder in patients of European ancestry. SNPs associated with drug-induced organ injury were rare but exhibited the largest effect sizes, consistent with the published literature.

CONCLUSION

This work demonstrates that an in silico bioinformatics-based approach using integrative analysis of a diversity of molecular and morphological data types can discover pharmacoepigenomic variants that are suitable candidates for further validation in cell lines, animal models and human clinical trials.

Genome-Wide Profiling of PARP1 Reveals an Interplay with Gene Regulatory Regions and DNA Methylation

Poly (ADP-ribose) polymerase-1 (PARP1) is a nuclear enzyme involved in DNA repair, chromatin remodeling and gene expression. PARP1 interactions with chromatin architectural multi-protein complexes (i.e. nucleosomes) alter chromatin structure resulting in changes in gene expression. Chromatin structure impacts gene regulatory processes including transcription, splicing, DNA repair, replication and recombination. It is important to delineate whether PARP1 randomly associates with nucleosomes or is present at specific nucleosome regions throughout the cell genome. We performed genome-wide association studies in breast cancer cell lines to address these questions. Our studies show that PARP1 associates with epigenetic regulatory elements genome-wide, such as active histone marks, CTCF and DNase hypersensitive sites. Additionally, the binding of PARP1 to chromatin genome-wide is mutually exclusive with DNA methylation pattern suggesting a functional interplay between PARP1 and DNA methylation. Indeed, inhibition of PARylation results in genome-wide changes in DNA methylation patterns. Our results suggest that PARP1 controls the fidelity of gene transcription and marks actively transcribed gene regions by selectively binding to transcriptionally active chromatin. These studies provide a platform for developing our understanding of PARP1’s role in gene regulation.

Tandem repeat variation in human and great ape populations and its impact on gene expression divergence

Tandem repeats (TRs) are stretches of DNA that are highly variable in length and mutate rapidly. They are thus an important source of genetic variation. This variation is highly informative for population and conservation genetics. It has also been associated with several pathological conditions and with gene expression regulation. However, genome-wide surveys of TR variation in humans and closely related species have been scarce due to technical difficulties derived from short-read technology. Here we explored the genome-wide diversity of TRs in a panel of 83 human and nonhuman great ape genomes, in a total of six different species, and studied their impact on gene expression evolution. We found that population diversity patterns can be efficiently captured with short TRs (repeat unit length, 1–5 bp). We examined the potential evolutionary role of TRs in gene expression differences between humans and primates by using 30,275 larger TRs (repeat unit length, 2–50 bp). Genes that contained TRs in the promoters, in their 3′ untranslated region, in introns, and in exons had higher expression divergence than genes without repeats in the regions. Polymorphic small repeats (1–5 bp) had also higher expression divergence compared with genes with fixed or no TRs in the gene promoters. Our findings highlight the potential contribution of TRs to human evolution through gene regulation.

Taking Stock of Regulatory Variation

Three recent studies measure individual variation in regulatory DNA accessibility. What do they tell us about the prospects of assessing variation in single cells and across populations?

Chromatin Dynamics in the Regulation of CFTR Expression

The contribution of chromatin dynamics to the regulation of human disease-associated loci such as the cystic fibrosis transmembrane conductance regulator (CFTR) gene has been the focus of intensive experimentation for many years. Recent technological advances in the analysis of transcriptional mechanisms across the entire human genome have greatly facilitated these studies. In this review we describe the complex machinery of tissue-specific regulation of CFTR expression, and put earlier observations in context by incorporating them into datasets generated by the most recent genomics methods. Though the gene promoter is required for CFTR expression, cell-type specific regulatory elements are located elsewhere in the gene and in flanking intergenic regions. Probably within its own topological domain established by the architectural proteins CTCF and cohesin, the CFTR locus utilizes chromatin dynamics to remodel nucleosomes, recruit cell-selective transcription factors, and activate intronic enhancers. These cis-acting elements are then brought to the gene promoter by chromatin looping mechanisms, which establish long-range interactions across the locus. Despite its complexity, the CFTR locus provides a paradigm for elucidating the critical role of chromatin dynamics in the transcription of individual human genes.

Functional compensation among HMGN variants modulates the DNase I hypersensitive sites at enhancers

DNase I hypersensitive sites (DHSs) are a hallmark of chromatin regions containing regulatory DNA such as enhancers and promoters; however, the factors affecting the establishment and maintenance of these sites are not fully understood. We now show that HMGN1 and HMGN2, nucleosome-binding proteins that are ubiquitously expressed in vertebrate cells, maintain the DHS landscape of mouse embryonic fibroblasts (MEFs) synergistically. Loss of one of these HMGN variants led to a compensatory increase of binding of the remaining variant. Genome-wide mapping of the DHSs in Hmgn1(-/-), Hmgn2(-/-), and Hmgn1(-/-)n2(-/-) MEFs reveals that loss of both, but not a single HMGN variant, leads to significant remodeling of the DHS landscape, especially at enhancer regions marked by H3K4me1 and H3K27ac. Loss of HMGN variants affects the induced expression of stress-responsive genes in MEFs, the transcription profiles of several mouse tissues, and leads to altered phenotypes that are not seen in mice lacking only one variant. We conclude that the compensatory binding of HMGN variants to chromatin maintains the DHS landscape, and the transcription fidelity and is necessary to retain wild-type phenotypes. Our study provides insight into mechanisms that maintain regulatory sites in chromatin and into functional compensation among nucleosome binding architectural proteins.

From Structural Variation of Gene Molecules to Chromatin Dynamics and Transcriptional Bursting

Transcriptional activation of eukaryotic genes is accompanied, in general, by a change in the sensitivity of promoter chromatin to endonucleases. The structural basis of this alteration has remained elusive for decades; but the change has been viewed as a transformation of one structure into another, from "closed" to "open" chromatin. In contradistinction to this static and deterministic view of the problem, a dynamical and probabilistic theory of promoter chromatin has emerged as its solution. This theory, which we review here, explains observed variation in promoter chromatin structure at the level of single gene molecules and provides a molecular basis for random bursting in transcription-the conjecture that promoters stochastically transition between transcriptionally conducive and inconducive states. The mechanism of transcriptional regulation may be understood only in probabilistic terms.

High throughput technologies for the functional discovery of mammalian enhancers: new approaches for understanding transcriptional regulatory network dynamics

Completion of the human and mouse genomes has inspired new initiatives to obtain a global understanding of the functional regulatory networks governing gene expression. Enhancers are primary regulatory DNA elements determining precise spatio- and temporal gene expression patterns, but the observation that they can function at any distance from the gene(s) they regulate has made their genome-wide characterization challenging. Since traditional, single reporter approaches would be unable to accomplish this enormous task, high throughput technologies for mapping chromatin features associated with enhancers have emerged as an effective surrogate for enhancer discovery. However, the last few years have witnessed the development of several new innovative approaches that can effectively screen for and discover enhancers based on their functional activation of transcription using massively parallel reporter systems. In addition to their application for genome annotation, these new high throughput functional approaches open new and exciting avenues for modeling gene regulatory networks.

Determinants of nucleosome positioning and their influence on plant gene expression

Nucleosome positioning influences the access of transcription factors (TFs) to their binding sites and gene expression. Studies in plant, animal, and fungal models demonstrate similar nucleosome positioning patterns along genes and correlations between occupancy and expression. However, the relationships among nucleosome positioning, cis-regulatory element accessibility, and gene expression in plants remain undefined. Here we showed that plant nucleosome depletion occurs on specific 6-mer motifs and this sequence-specific nucleosome depletion is predictive of expression levels. Nucleosome-depleted regions in Arabidopsis thaliana tend to have higher G/C content, unlike yeast, and are centered on specific G/C-rich 6-mers, suggesting that intrinsic sequence properties, such as G/C content, cannot fully explain plant nucleosome positioning. These 6-mer motif sites showed higher DNase I hypersensitivity and are flanked by strongly phased nucleosomes, consistent with known TF binding sites. Intriguingly, this 6-mer-specific nucleosome depletion pattern occurs not only in promoter but also in genic regions and is significantly correlated with higher gene expression level, a phenomenon also found in rice but not in yeast. Among the 6-mer motifs enriched in genes responsive to treatment with the defense hormone jasmonate, there are no significant changes in nucleosome occupancy, suggesting that these sites are potentially preconditioned to enable rapid response without changing chromatin state significantly. Our study provides a global assessment of the joint contribution of nucleosome occupancy and motif sequences that are likely cis-elements to the control of gene expression in plants. Our findings pave the way for further understanding the impact of chromatin state on plant transcriptional regulatory circuits.

Dynamics of chromatin accessibility and long-range interactions in response to glucocorticoid pulsing

Although physiological steroid levels are often pulsatile (ultradian), the genomic effects of this pulsatility are poorly understood. By utilizing glucocorticoid receptor (GR) signaling as a model system, we uncovered striking spatiotemporal relationships between receptor loading, lifetimes of the DNase I hypersensitivity sites (DHSs), long-range interactions, and gene regulation. We found that hormone-induced DHSs were enriched within ± 50 kb of GR-responsive genes and displayed a broad spectrum of lifetimes upon hormone withdrawal. These lifetimes dictate the strength of the DHS interactions with gene targets and contribute to gene regulation from a distance. Our results demonstrate that pulsatile and constant hormone stimulations induce unique, treatment-specific patterns of gene and regulatory element activation. These modes of activation have implications for corticosteroid function in vivo and for steroid therapies in various clinical settings.

The Proportion of Chromatin Graded between Closed and Open States Determines the Level of Transcripts Derived from Distinct Promoters in the CYP19 Gene

The human CYP19 gene encodes aromatase, which converts androgens to estrogens. CYP19 mRNA variants are transcribed mainly from three promoters. Quantitative RT-PCR was used to measure the relative amounts of each of the three transcripts and determine the on/off state of the promoters. While some of the promoters were silent, CYP19 mRNA production differed among the other promoters, whose estimated transcription levels were 0.001% to 0.1% of that of the TUBB control gene. To investigate the structural aspects of chromatin that were responsible for this wide range of activity of the CYP19 promoters, we used a fractionation protocol, designated SEVENS, which sequentially separates densely packed nucleosomes from dispersed nucleosomes. The fractional distribution of each inactive promoter showed a similar pattern to that of the repressed reference loci; the inactive regions were distributed toward lower fractions, in which closed chromatin comprising packed nucleosomes was enriched. In contrast, active CYP19 promoters were raised toward upper fractions, including dispersed nucleosomes in open chromatin. Importantly, these active promoters were moderately enriched in the upper fractions as compared to active reference loci, such as the TUBB promoter; the proportion of open chromatin appeared to be positively correlated to the promoter strength. These results, together with ectopic transcription accompanied by an increase in the proportion of open chromatin in cells treated with an H3K27me inhibitor, indicate that CYP19 mRNA could be transcribed from a promoter in which chromatin is shifted toward an open state in the equilibrium between closed and open chromatin.

Survey of protein-DNA interactions in Aspergillus oryzae on a genomic scale

The genome-scale delineation of in vivo protein–DNA interactions is key to understanding genome function. Only ∼5% of transcription factors (TFs) in the Aspergillus genus have been identified using traditional methods. Although the Aspergillus oryzae genome contains >600 TFs, knowledge of the in vivo genome-wide TF-binding sites (TFBSs) in aspergilli remains limited because of the lack of high-quality antibodies. We investigated the landscape of in vivo protein–DNA interactions across the A. oryzae genome through coupling the DNase I digestion of intact nuclei with massively parallel sequencing and the analysis of cleavage patterns in protein–DNA interactions at single-nucleotide resolution. The resulting map identified overrepresented de novo TF-binding motifs from genomic footprints, and provided the detailed chromatin remodeling patterns and the distribution of digital footprints near transcription start sites. The TFBSs of 19 known Aspergillus TFs were also identified based on DNase I digestion data surrounding potential binding sites in conjunction with TF binding specificity information. We observed that the cleavage patterns of TFBSs were dependent on the orientation of TF motifs and independent of strand orientation, consistent with the DNA shape features of binding motifs with flanking sequences.

Unraveling the 3D genome: genomics tools for multiscale exploration

A decade of rapid method development has begun to yield exciting insights into the 3D architecture of the metazoan genome and the roles it may play in regulating transcription. Here we review core methods and new tools in the modern genomicist's toolbox at three length scales, ranging from single base pairs to megabase-scale chromosomal domains, and discuss the emerging picture of the 3D genome that these tools have revealed. Blind spots remain, especially at intermediate length scales spanning a few nucleosomes, but thanks in part to new technologies that permit targeted alteration of chromatin states and time-resolved studies, the next decade holds great promise for hypothesis-driven research into the mechanisms that drive genome architecture and transcriptional regulation.

The 'dark matter' in the plant genomes: non-coding and unannotated DNA sequences associated with open chromatin

Sequencing of complete plant genomes has become increasingly more routine since the advent of the next-generation sequencing technology. Identification and annotation of large amounts of noncoding but functional DNA sequences, including cis-regulatory DNA elements (CREs), have become a new frontier in plant genome research. Genomic regions containing active CREs bound to regulatory proteins are hypersensitive to DNase I digestion and are called DNase I hypersensitive sites (DHSs). Several recent DHS studies in plants illustrate that DHS datasets produced by DNase I digestion followed by next-generation sequencing (DNase-seq) are highly valuable for the identification and characterization of CREs associated with plant development and responses to environmental cues. DHS-based genomic profiling has opened a door to identify and annotate the 'dark matter' in sequenced plant genomes.

Enhancer-associated RNAs as therapeutic targets

INTRODUCTION

Regulation of gene expression involves a variety of mechanisms driven by a complex regulatory network of factors. Control of transcription is an important step in gene expression regulation, which integrates the function of cis-acting and trans-acting elements. Among cis-regulatory elements, enhancer RNA (eRNA)-producing domains recently emerged as widespread and potent regulators of transcription and cell fate decision. Thus, manipulation of eRNA levels becomes a novel and appealing avenue for the design of new therapeutic treatments.

AREAS COVERED

In this review, we focus on eRNA-producing domains. We describe mechanisms involved in their cell-type specific selection and activation as well as their epigenetic features. In addition, we present their function and the growing evidences of their deregulation in human diseases. Finally, we discuss eRNAs as potential therapeutic targets.

EXPERT OPINION

As key factors in the control of transcription, eRNAs appear to possess a great potential for the establishment of new therapy options. However, thorough testing as well as providing the genetic toolbox to target eRNAs will be needed to fully assess the practical and clinical possibilities.

Using the ENCODE Resource for Functional Annotation of Genetic Variants

This article illustrates the use of the Encyclopedia of DNA Elements (ENCODE) resource to generate or refine hypotheses from genomic data on disease and other phenotypic traits. First, the goals and history of ENCODE and related epigenomics projects are reviewed. Second, the rationale for ENCODE and the major data types used by ENCODE are briefly described, as are some standard heuristics for their interpretation. Third, the use of the ENCODE resource is examined. Standard use cases for ENCODE, accessing the ENCODE resource, and accessing data from related projects are discussed. Although the focus of this article is the use of ENCODE data, some of the same approaches can be used with data from other projects.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

An integrative analysis of TFBS-clustered regions reveals new transcriptional regulation models on the accessible chromatin landscape

DNase I hypersensitive sites (DHSs) define the accessible chromatin landscape and have revolutionised the discovery of distinct cis-regulatory elements in diverse organisms. Here, we report the first comprehensive map of human transcription factor binding site (TFBS)-clustered regions using Gaussian kernel density estimation based on genome-wide mapping of the TFBSs in 133 human cell and tissue types. Approximately 1.6 million distinct TFBS-clustered regions, collectively spanning 27.7% of the human genome, were discovered. The TFBS complexity assigned to each TFBS-clustered region was highly correlated with genomic location, cell selectivity, evolutionary conservation, sequence features, and functional roles. An integrative analysis of these regions using ENCODE data revealed transcription factor occupancy, transcriptional activity, histone modification, DNA methylation, and chromatin structures that varied based on TFBS complexity. Furthermore, we found that we could recreate lineage-branching relationships by simple clustering of the TFBS-clustered regions from terminally differentiated cells. Based on these findings, a model of transcriptional regulation determined by TFBS complexity is proposed.

Tissue-resident macrophage enhancer landscapes are shaped by the local microenvironment

Macrophages are critical for innate immune defense and also control organ homeostasis in a tissue-specific manner. They provide a fitting model to study the impact of ontogeny and microenvironment on chromatin state and whether chromatin modifications contribute to macrophage identity. Here, we profile the dynamics of four histone modifications across seven tissue-resident macrophage populations. We identify 12,743 macrophage-specific enhancers and establish that tissue-resident macrophages have distinct enhancer landscapes beyond what can be explained by developmental origin. Combining our enhancer catalog with gene expression profiles and open chromatin regions, we show that a combination of tissue- and lineage-specific transcription factors form the regulatory networks controlling chromatin specification in tissue-resident macrophages. The environment is capable of shaping the chromatin landscape of transplanted bone marrow precursors, and even differentiated macrophages can be reprogrammed when transferred into a new microenvironment. These results provide a comprehensive view of macrophage regulatory landscape and highlight the importance of the microenvironment, along with pioneer factors in orchestrating identity and plasticity.

DNA topoisomerases beyond the standard role

Chromatin is dynamically changing its structure to accommodate and control DNA-dependent processes inside of eukaryotic cells. These changes are necessarily linked to changes of DNA topology, which might itself serve as a regulatory signal to be detected by proteins. Thus, DNA Topoisomerases may contribute to the regulation of many events occurring during the transcription cycle. In this review we will focus on DNA Topoisomerase functions in transcription, with particular emphasis on the multiplicity of tasks beyond their widely appreciated role in solving topological problems associated with transcription elongation.

Principles of regulatory information conservation between mouse and human

To broaden our understanding of the evolution of gene regulation mechanisms, we generated occupancy profiles for 34 orthologous transcription factors (TFs) in human-mouse erythroid progenitor, lymphoblast and embryonic stem-cell lines. By combining the genome-wide transcription factor occupancy repertoires, associated epigenetic signals, and co-association patterns, here we deduce several evolutionary principles of gene regulatory features operating since the mouse and human lineages diverged. The genomic distribution profiles, primary binding motifs, chromatin states, and DNA methylation preferences are well conserved for TF-occupied sequences. However, the extent to which orthologous DNA segments are bound by orthologous TFs varies both among TFs and with genomic location: binding at promoters is more highly conserved than binding at distal elements. Notably, occupancy-conserved TF-occupied sequences tend to be pleiotropic; they function in several tissues and also co-associate with many TFs. Single nucleotide variants at sites with potential regulatory functions are enriched in occupancy-conserved TF-occupied sequences.

ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide

This unit describes Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq), a method for mapping chromatin accessibility genome-wide. This method probes DNA accessibility with hyperactive Tn5 transposase, which inserts sequencing adapters into accessible regions of chromatin. Sequencing reads can then be used to infer regions of increased accessibility, as well as to map regions of transcription-factor binding and nucleosome position. The method is a fast and sensitive alternative to DNase-seq for assaying chromatin accessibility genome-wide, or to MNase-seq for assaying nucleosome positions in accessible regions of the genome.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

A novel virus-inducible enhancer of the interferon-β gene with tightly linked promoter and enhancer activities

Long-range enhancers of transcription are a key component of the genomic regulatory architecture. Recent studies have identified bi-directionally transcribed RNAs emanating from these enhancers known as eRNAs. However, it remains unclear how tightly coupled eRNA production is with enhancer activity. Through our systematic search for long-range elements that interact with the interferon-β gene, a model system for studying inducible transcription, we have identified a novel enhancer, which we have named L2 that regulates the expression of interferon-β. We have demonstrated its virus-inducible enhancer activity by analyzing epigenomic profiles, transcription factor association, nascent RNA production and activity in reporter assays. This enhancer exhibits intimately linked virus-inducible enhancer and bidirectional promoter activity that is largely dependent on a conserved Interferon Stimulated Response Element and robustly generates virus inducible eRNAs. Notably, its enhancer and promoter activities are fully retained in reporter assays even upon a complete elimination of its associated eRNA sequences. Finally, we show that L2 regulates IFNB1 expression by siRNA knockdown of eRNAs, and the deletion of L2 in a BAC transfection assay. Thus, L2 is a novel enhancer that regulates IFNB1 and whose eRNAs exert significant activity in vivo that is distinct from those activities recapitulated in the luciferase reporter assays.

Expression QTL-based analyses reveal candidate causal genes and loci across five tumor types

The majority of trait-associated loci discovered through genome-wide association studies are located outside of known protein coding regions. Consequently, it is difficult to ascertain the mechanism underlying these variants and to pinpoint the causal alleles. Expression quantitative trait loci (eQTLs) provide an organizing principle to address both of these issues. eQTLs are genetic loci that correlate with RNA transcript levels. Large-scale data sets such as the Cancer Genome Atlas (TCGA) provide an ideal opportunity to systematically evaluate eQTLs as they have generated multiple data types on hundreds of samples. We evaluated the determinants of gene expression (germline variants and somatic copy number and methylation) and performed cis-eQTL analyses for mRNA expression and miRNA expression in five tumor types (breast, colon, kidney, lung and prostate). We next tested 149 known cancer risk loci for eQTL effects, and observed that 42 (28.2%) were significantly associated with at least one transcript. Lastly, we described a fine-mapping strategy for these 42 eQTL target-gene associations based on an integrated strategy that combines the eQTL level of significance and the regulatory potential as measured by DNaseI hypersensitivity. For each of the risk loci, our analyses suggested 1 to 81 candidate causal variants that may be prioritized for downstream functional analysis. In summary, our study provided a comprehensive landscape of the genetic determinants of gene expression in different tumor types and ranked the genes and loci for further functional assessment of known cancer risk loci.

Enhancer alterations in cancer: a source for a cell identity crisis

Enhancers are selectively utilized to orchestrate gene expression programs that first govern pluripotency and then proceed to highly specialized programs required for the process of cellular differentiation. Whereas gene-proximal promoters are typically active across numerous cell types, distal enhancer activation is cell-type-specific and central to cell fate determination, thereby accounting for cell identity. Recent studies have highlighted the diversity of enhancer usage, cataloguing millions of such elements in the human genome. The disruption of enhancer activity, through genetic or epigenetic alterations, can impact cell-type-specific functions, resulting in a wide range of pathologies. In cancer, these alterations can promote a 'cell identity crisis', in which enhancers associated with oncogenes and multipotentiality are activated, while those promoting cell fate commitment are inactivated. Overall, these alterations favor an undifferentiated cellular phenotype. Here, we review the current knowledge regarding the role of enhancers in normal cell function, and discuss how genetic and epigenetic changes in enhancer elements potentiate oncogenesis. In addition, we discuss how understanding the mechanisms regulating enhancer activity can inform therapeutic opportunities in cancer cells and highlight key challenges that remain in understanding enhancer biology as it relates to oncology.

Prediction of DNase I hypersensitive sites by using pseudo nucleotide compositions

DNase I hypersensitive sites (DHS) associated with a wide variety of regulatory DNA elements. Knowledge about the locations of DHS is helpful for deciphering the function of noncoding genomic regions. With the acceleration of genome sequences in the postgenomic age, it is highly desired to develop cost-effective computational methods to identify DHS. In the present work, a support vector machine based model was proposed to identify DHS by using the pseudo dinucleotide composition. In the jackknife test, the proposed model obtained an accuracy of 83%, which is competitive with that of the existing method. This result suggests that the proposed model may become a useful tool for DHS identifications.

Identification of novel transcription factors in osteoclast differentiation using genome-wide analysis of open chromatin determined by DNase-seq

Clarification of the mechanisms underlying osteoclast differentiation enables us to understand the physiology of bone metabolism as well as the pathophysiology of bone diseases such as osteoporosis. Recently, it has been reported that epigenetics can determine cell fate and regulate cell type-specific gene expression. However, little is known about epigenetics during osteoclastogenesis. To reveal a part of epigenetics, especially focused on chromatin dynamics, during early osteoclastogenesis and to identify novel transcription factors involved in osteoclastogenesis, we performed a genome-wide analysis of open chromatin during receptor activator of NF-κB ligand (RANKL)-induced osteoclastogenesis using DNase I hypersensitive sites sequencing (DNase-seq). DNase-seq was performed using the extracted nuclei from RAW264 cells treated with or without RANKL for 24 hours, followed by several bioinformatic analyses. DNase I hypersensitive sites (DHSs) were dynamically changed during RANKL-induced osteoclastogenesis and they accumulated in promoter regions. The distributions of DHSs among cis-regulatory DNA regions were identical regardless of RANKL stimulation. Motif discovery analysis successfully identified well-known osteoclastogenic transcription factors including Jun, CREB1, FOS, ATF2, and ATF4, but also novel transcription factors for osteoclastogenesis such as Zscan10, Atf1, Nrf1, and Srebf2. siRNA knockdown of these identified novel transcription factors impaired osteoclastogenesis. Taken together, DNase-seq is a useful tool for comprehension of epigenetics, especially chromatin dynamics during osteoclastogenesis and for identification of novel transcription factors involved in osteoclastogenesis. This study may reveal underlying mechanisms that determine cell type-specific differentiation of bone cells and may lead to investigation of novel therapeutic targets for osteoporosis.

Effective automated feature construction and selection for classification of biological sequences

BACKGROUND

Many open problems in bioinformatics involve elucidating underlying functional signals in biological sequences. DNA sequences, in particular, are characterized by rich architectures in which functional signals are increasingly found to combine local and distal interactions at the nucleotide level. Problems of interest include detection of regulatory regions, splice sites, exons, hypersensitive sites, and more. These problems naturally lend themselves to formulation as classification problems in machine learning. When classification is based on features extracted from the sequences under investigation, success is critically dependent on the chosen set of features.

METHODOLOGY

We present an algorithmic framework (EFFECT) for automated detection of functional signals in biological sequences. We focus here on classification problems involving DNA sequences which state-of-the-art work in machine learning shows to be challenging and involve complex combinations of local and distal features. EFFECT uses a two-stage process to first construct a set of candidate sequence-based features and then select a most effective subset for the classification task at hand. Both stages make heavy use of evolutionary algorithms to efficiently guide the search towards informative features capable of discriminating between sequences that contain a particular functional signal and those that do not.

RESULTS

To demonstrate its generality, EFFECT is applied to three separate problems of importance in DNA research: the recognition of hypersensitive sites, splice sites, and ALU sites. Comparisons with state-of-the-art algorithms show that the framework is both general and powerful. In addition, a detailed analysis of the constructed features shows that they contain valuable biological information about DNA architecture, allowing biologists and other researchers to directly inspect the features and potentially use the insights obtained to assist wet-laboratory studies on retainment or modification of a specific signal. Code, documentation, and all data for the applications presented here are provided for the community at http://www.cs.gmu.edu/~ashehu/?q=OurTools.

Open chromatin profiling in mice livers reveals unique chromatin variations induced by high fat diet

Metabolic diseases result from multiple genetic and environmental factors. We report here that one manner in which environmental factors can contribute to metabolic disease progression is through modification to chromatin. We demonstrate that high fat diet leads to chromatin remodeling in the livers of C57BL/6J mice, as compared with mice fed a control diet, and that these chromatin changes are associated with changes in gene expression. We further show that the regions of greatest variation in chromatin accessibility are targeted by liver transcription factors, including HNF4α, CCAAT/enhancer-binding protein α (CEBP/α), and FOXA1. Repeating the chromatin and gene expression profiling in another mouse strain, DBA/2J, revealed that the regions of greatest chromatin change are largely strain-specific and that integration of chromatin, gene expression, and genetic data can be used to characterize regulatory regions. Our data indicate dramatic changes in the epigenome due to diet and demonstrate strain-specific dynamics in chromatin remodeling.

An enhancer polymorphism at the cardiomyocyte intercalated disc protein NOS1AP locus is a major regulator of the QT interval

QT interval variation is assumed to arise from variation in repolarization as evidenced from rare Na- and K-channel mutations in Mendelian QT prolongation syndromes. However, in the general population, common noncoding variants at a chromosome 1q locus are the most common genetic regulators of QT interval variation. In this study, we use multiple human genetic, molecular genetic, and cellular assays to identify a functional variant underlying trait association: a noncoding polymorphism (rs7539120) that maps within an enhancer of NOS1AP and affects cardiac function by increasing NOS1AP transcript expression. We further localized NOS1AP to cardiomyocyte intercalated discs (IDs) and demonstrate that overexpression of NOS1AP in cardiomyocytes leads to altered cellular electrophysiology. We advance the hypothesis that NOS1AP affects cardiac electrical conductance and coupling and thereby regulates the QT interval through propagation defects. As further evidence of an important role for propagation variation affecting QT interval in humans, we show that common polymorphisms mapping near a specific set of 170 genes encoding ID proteins are significantly enriched for association with the QT interval, as compared to genome-wide markers. These results suggest that focused studies of proteins within the cardiomyocyte ID are likely to provide insights into QT prolongation and its associated disorders.

Generation of GFP Reporter Human Induced Pluripotent Stem Cells Using AAVS1 Safe Harbor Transcription Activator-Like Effector Nuclease

Generation of a fluorescent GFP reporter line in human induced pluripotent stem cells (hiPSCs) provides enormous potentials in both basic stem cell research and regenerative medicine. A protocol for efficiently generating such an engineered reporter line by gene targeting is highly desired. Transcription activator-like effector nucleases (TALENs) are a new class of artificial restriction enzymes that have been shown to significantly promote homologous recombination by >1000-fold. The AAVS1 (adeno-associated virus integration site 1) locus is a "safe harbor" and has an open chromatin structure that allows insertion and stable expression of transgene. Here, we describe a step-by-step protocol from determination of TALENs activity, hiPSC culture, and delivery of a donor into AAVS1 targeting site, to validation of targeted integration by PCR and Southern blot analysis using hiPSC line, and a pair of open-source AAVS1 TALENs.

The chromatin "landscape" of a murine adult β-globin gene is unaffected by deletion of either the gene promoter or a downstream enhancer

In mammals, the complex tissue- and developmental-specific expression of genes within the β-globin cluster is known to be subject to control by the gene promoters, by a locus control region (LCR) located upstream of the cluster, and by sequence elements located across the intergenic regions. Despite extensive investigation, however, the complement of sequences that is required for normal regulation of chromatin structure and gene expression within the cluster is not fully defined. To further elucidate regulation of the adult β-globin genes, we investigate the effects of two deletions engineered within the endogenous murine β-globin locus. First, we find that deletion of the β2-globin gene promoter, while eliminating β2-globin gene expression, results in no additional effects on chromatin structure or gene expression within the cluster. Notably, our observations are not consistent with competition among the β-globin genes for LCR activity. Second, we characterize a novel enhancer located 3′ of the β2-globin gene, but find that deletion of this sequence has no effect whatsoever on gene expression or chromatin structure. This observation highlights the difficulty in assigning function to enhancer sequences identified by the chromatin “landscape” or even by functional assays.

Defining functional DNA elements in the human genome

With the completion of the human genome sequence, attention turned to identifying and annotating its functional DNA elements. As a complement to genetic and comparative genomics approaches, the Encyclopedia of DNA Elements Project was launched to contribute maps of RNA transcripts, transcriptional regulator binding sites, and chromatin states in many cell types. The resulting genome-wide data reveal sites of biochemical activity with high positional resolution and cell type specificity that facilitate studies of gene regulation and interpretation of noncoding variants associated with human disease. However, the biochemically active regions cover a much larger fraction of the genome than do evolutionarily conserved regions, raising the question of whether nonconserved but biochemically active regions are truly functional. Here, we review the strengths and limitations of biochemical, evolutionary, and genetic approaches for defining functional DNA segments, potential sources for the observed differences in estimated genomic coverage, and the biological implications of these discrepancies. We also analyze the relationship between signal intensity, genomic coverage, and evolutionary conservation. Our results reinforce the principle that each approach provides complementary information and that we need to use combinations of all three to elucidate genome function in human biology and disease.

Uncoupling transcription from covalent histone modification

It is widely accepted that transcriptional regulation of eukaryotic genes is intimately coupled to covalent modifications of the underlying chromatin template, and in certain cases the functional consequences of these modifications have been characterized. Here we present evidence that gene activation in the silent heterochromatin of the yeast Saccharomyces cerevisiae can occur in the context of little, if any, covalent histone modification. Using a SIR-regulated heat shock-inducible transgene, hsp82-2001, and a natural drug-inducible subtelomeric gene, YFR057w, as models we demonstrate that substantial transcriptional induction (>200-fold) can occur in the context of restricted histone loss and negligible levels of H3K4 trimethylation, H3K36 trimethylation and H3K79 dimethylation, modifications commonly linked to transcription initiation and elongation. Heterochromatic gene activation can also occur with minimal H3 and H4 lysine acetylation and without replacement of H2A with the transcription-linked variant H2A.Z. Importantly, absence of histone modification does not stem from reduced transcriptional output, since hsp82-ΔTATA, a euchromatic promoter mutant lacking a TATA box and with threefold lower induced transcription than heterochromatic hsp82-2001, is strongly hyperacetylated in response to heat shock. Consistent with negligible H3K79 dimethylation, dot1Δ cells lacking H3K79 methylase activity show unimpeded occupancy of RNA polymerase II within activated heterochromatic promoter and coding regions. Our results indicate that large increases in transcription can be observed in the virtual absence of histone modifications often thought necessary for gene activation.

The proper regulation of gene expression is of fundamental importance in the maintenance of normal growth and development. Misregulation of genes can lead to such outcomes as cancer, diabetes and neurodegenerative disease. A key step in gene regulation occurs during the transcription of the chromosomal DNA into messenger RNA by the enzyme, RNA polymerase II. Histones are small, positively charged proteins that package genomic DNA into arrays of bead-like particles termed nucleosomes, the principal components of chromatin. Increasing evidence suggests that nucleosomal histones play an active role in regulating transcription, and that this is derived in part from reversible chemical (“covalent”) modifications that take place on their amino acids. These histone modifications create novel surfaces on nucleosomes that can serve as docking sites for other proteins that control a gene's expression state. In this study we present evidence that contrary to the general case, covalent modifications typically associated with transcription are minimally used by genes embedded in a specialized, condensed chromatin structure termed heterochromatin in the model organism baker's yeast. Our observations are significant, for they suggest that gene transcription can occur in a living cell in the virtual absence of covalent modification of the chromatin template.

STAT3, STAT4, NFATc1, and CTCF regulate PD-1 through multiple novel regulatory regions in murine T cells

Programmed death-1 (PD-1) is a crucial negative regulator of CD8 T cell development and function, yet the mechanisms that control its expression are not fully understood. Through a nonbiased DNase I hypersensitivity assay, four novel regulatory regions within the Pdcd1 locus were identified. Two of these elements flanked the locus, bound the transcriptional insulator protein CCCTC-binding factor, and interacted with each other, creating a potential regulatory compartmentalization of the locus. In response to T cell activation signaling, NFATc1 bound to two of the novel regions that function as independent regulatory elements. STAT binding sites were identified in these elements as well. In splenic CD8 T cells, TCR-induced PD-1 expression was augmented by IL-6 and IL-12, inducers of STAT3 and STAT4 activity, respectively. IL-6 or IL-12 on its own did not induce PD-1. Importantly, STAT3/4 and distinct chromatin modifications were associated with the novel regulatory regions following cytokine stimulation. The NFATc1/STAT regulatory regions were found to interact with the promoter region of the Pdcd1 gene, providing a mechanism for their action. Together these data add multiple novel distal regulatory regions and pathways to the control of PD-1 expression and provide a molecular mechanism by which proinflammatory cytokines, such as IL-6 or IL-12, can augment PD-1 expression.

Enhancer RNAs and regulated transcriptional programs

A large portion of the human genome is transcribed into RNAs without known protein-coding functions, far outnumbering coding transcription units. Extensive studies of long noncoding RNAs (lncRNAs) have clearly demonstrated that they can play critical roles in regulating gene expression, development, and diseases, acting both as transcriptional activators and repressors. More recently, enhancers have been found to be broadly transcribed, resulting in the production of enhancer-derived RNAs, or eRNAs. Here, we review emerging evidence suggesting that at least some eRNAs contribute to enhancer function. We discuss these findings with respect to potential mechanisms of action of eRNAs and other ncRNAs in regulated gene expression.

The chromatin architectural proteins HMGD1 and H1 bind reciprocally and have opposite effects on chromatin structure and gene regulation

BACKGROUND

Chromatin architectural proteins interact with nucleosomes to modulate chromatin accessibility and higher-order chromatin structure. While these proteins are almost certainly important for gene regulation they have been studied far less than the core histone proteins.

RESULTS

Here we describe the genomic distributions and functional roles of two chromatin architectural proteins: histone H1 and the high mobility group protein HMGD1 in Drosophila S2 cells. Using ChIP-seq, biochemical and gene specific approaches, we find that HMGD1 binds to highly accessible regulatory chromatin and active promoters. In contrast, H1 is primarily associated with heterochromatic regions marked with repressive histone marks. We find that the ratio of HMGD1 to H1 binding is a better predictor of gene activity than either protein by itself, which suggests that reciprocal binding between these proteins is important for gene regulation. Using knockdown experiments, we show that HMGD1 and H1 affect the occupancy of the other protein, change nucleosome repeat length and modulate gene expression.

CONCLUSION

Collectively, our data suggest that dynamic and mutually exclusive binding of H1 and HMGD1 to nucleosomes and their linker sequences may control the fluid chromatin structure that is required for transcriptional regulation. This study provides a framework to further study the interplay between chromatin architectural proteins and epigenetics in gene regulation.

Rapid and unambiguous detection of DNase I hypersensitive site in rare population of cells

DNase I hypersensitive (DHS) sites are important for understanding cis regulation of gene expression. However, existing methods for detecting DHS sites in small numbers of cells can lead to ambiguous results. Here we describe a simple new method, in which DNA fragments with ends generated by DNase I digestion are isolated and used as templates for two PCR reactions. In the first PCR, primers are derived from sequences up- and down-stream of the DHS site. If the DHS site exists in the cells, the first PCR will not produce PCR products due to the cuts of the templates by DNase I between the primer sequences. In the second PCR, one primer is derived from sequence outside the DHS site and the other from the adaptor. This will produce a smear of PCR products of different sizes due to cuts by DNase I at different positions at the DHS site. With this design, we detected a DHS site at the CD4 gene in two CD4 T cell populations using as few as 2×10(4) cells. We further validated this method by detecting a DHS site of the IL-4 gene that is specifically present in type 2 but not type 1 T helper cells. Overall, this method overcomes the interference by genomic DNA not cut by DNase I at the DHS site, thereby offering unambiguous detection of DHS sites in the cells.

Genome-scale mapping of DNase I hypersensitivity

DNase I-seq is a global and high-resolution method that uses the nonspecific endonuclease DNase I to map chromatin accessibility. These accessible regions, designated as DNase I hypersensitive sites (DHSs), define the regulatory features, (e.g., promoters, enhancers, insulators, and locus control regions) of complex genomes. In this unit, methods are described for nuclei isolation, digestion of nuclei with limiting concentrations of DNase I, and the biochemical fractionation of DNase I hypersensitive sites in preparation for high-throughput sequencing. DNase I-seq is an unbiased and robust method that is not predicated on an a priori understanding of regulatory patterns or chromatin features.

Enhancer-targeted genome editing selectively blocks innate resistance to oncokinase inhibition

Thousands of putative enhancers are characterized in the human genome, yet few have been shown to have a functional role in cancer progression. Inhibiting oncokinases, such as EGFR, ALK, ERBB2, and BRAF, is a mainstay of current cancer therapy but is hindered by innate drug resistance mediated by up-regulation of the HGF receptor, MET. The mechanisms mediating such genomic responses to targeted therapy are unknown. Here, we identify lineage-specific enhancers at the MET locus for multiple common tumor types, including a melanoma lineage-specific enhancer 63 kb downstream from the MET TSS. This enhancer displays inducible chromatin looping with the MET promoter to up-regulate MET expression upon BRAF inhibition. Epigenomic analysis demonstrated that the melanocyte-specific transcription factor, MITF, mediates this enhancer function. Targeted genomic deletion (<7 bp) of the MITF motif within the MET enhancer suppressed inducible chromatin looping and innate drug resistance, while maintaining MITF-dependent, inhibitor-induced melanoma cell differentiation. Epigenomic analysis can thus guide functional disruption of regulatory DNA to decouple pro- and anti-oncogenic functions of a dominant transcription factor and block innate resistance to oncokinase therapy.