Nuclease hypersensitive sites in chromatin

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.

Reduced local mutation density in regulatory DNA of cancer genomes is linked to DNA repair

Carcinogenesis and neoplastic progression are mediated by the accumulation of somatic mutations. Here we report that the local density of somatic mutations in cancer genomes is highly reduced specifically in accessible regulatory DNA defined by DNase I hypersensitive sites. This reduction is independent of any known factors influencing somatic mutation density and is observed in diverse cancer types, suggesting a general mechanism. By analyzing individual cancer genomes¹, we show that the reduced local mutation density within regulatory DNA is linked to intact global genome repair machinery, with nearly complete abrogation of the hypomutation phenomenon in individual cancers that possess mutations in multiple nucleotide excision repair components. Together, our results connect chromatin structure, gene regulation and cancer-associated somatic mutation.

Analysis of the SWI/SNF chromatin-remodeling complex during early heart development and BAF250a repression cardiac gene transcription during P19 cell differentiation

The regulatory networks of differentiation programs and the molecular mechanisms of lineage-specific gene regulation in mammalian embryos remain only partially defined. We document differential expression and temporal switching of BRG1-associated factor (BAF) subunits, core pluripotency factors and cardiac-specific genes during post-implantation development and subsequent early organogenesis. Using affinity purification of BRG1 ATPase coupled to mass spectrometry, we characterized the cardiac-enriched remodeling complexes present in E8.5 mouse embryos. The relative abundance and combinatorial assembly of the BAF subunits provides functional specificity to Switch/Sucrose NonFermentable (SWI/SNF) complexes resulting in a unique gene expression profile in the developing heart. Remarkably, the specific depletion of the BAF250a subunit demonstrated differential effects on cardiac-specific gene expression and resulted in arrhythmic contracting cardiomyocytes in vitro. Indeed, the BAF250a physically interacts and functionally cooperates with Nucleosome Remodeling and Histone Deacetylase (NURD) complex subunits to repressively regulate chromatin structure of the cardiac genes by switching open and poised chromatin marks associated with active and repressed gene expression. Finally, BAF250a expression modulates BRG1 occupancy at the loci of cardiac genes regulatory regions in P19 cell differentiation. These findings reveal specialized and novel cardiac-enriched SWI/SNF chromatin-remodeling complexes, which are required for heart formation and critical for cardiac gene expression regulation at the early stages of heart development.

The 5p12 breast cancer susceptibility locus affects MRPS30 expression in estrogen-receptor positive tumors

Genome-wide association studies have identified numerous loci linked to breast cancer susceptibility, but the mechanism by which variations at these loci influence susceptibility is usually unknown. Some variants are only associated with particular clinical subtypes of breast cancer. Understanding how and why these variants influence subtype-specific cancer risk contributes to our understanding of cancer etiology. We conducted a genome-wide expression Quantitative Trait Locus (eQTL) study in a discovery set of 287 breast tumors and 97 normal mammary tissue samples and a replication set of 235 breast tumors. We found that the risk-associated allele of rs7716600 in the 5p12 estrogen receptor-positive (ER-positive) susceptibility locus was associated with elevated expression of the nearby gene MRPS30 exclusively in ER-positive tumors. We replicated this finding in 235 independent tumors. Further, we showed the rs7716600 risk genotype was associated with decreased MRPS30 promoter methylation exclusively in ER-positive breast tumors. In vitro studies in MCF-7 cells carrying the protective genotype showed that estrogen stimulation decreased MRPS30 promoter chromatin availability and mRNA levels. In contrast, in 600MPE cells carrying the risk genotype, estrogen increased MRPS30 expression and did not affect promoter availability. Our data suggest the 5p12 risk allele affects MRPS30 expression in estrogen-responsive tumor cells after tumor initiation by a mechanism affecting chromatin availability. These studies emphasize that the genetic architecture of breast cancer is context-specific, and integrated analysis of gene expression and chromatin remodeling in normal and tumor tissues will be required to explain the mechanisms of risk alleles.

Insights into chromatin structure and dynamics in plants

The packaging of chromatin into the nucleus of a eukaryotic cell requires an extraordinary degree of compaction and physical organization. In recent years, it has been shown that this organization is dynamically orchestrated to regulate responses to exogenous stimuli as well as to guide complex cell-type-specific developmental programs. Gene expression is regulated by the compartmentalization of functional domains within the nucleus, by distinct nucleosome compositions accomplished via differential modifications on the histone tails and through the replacement of core histones by histone variants. In this review, we focus on these aspects of chromatin organization and discuss novel approaches such as live cell imaging and photobleaching as important tools likely to give significant insights into our understanding of the very dynamic nature of chromatin and chromatin regulatory processes. We highlight the contribution plant studies have made in this area showing the potential advantages of plants as models in understanding this fundamental aspect of biology.

Predicting tissue specific transcription factor binding sites

BACKGROUND

Studies of gene regulation often utilize genome-wide predictions of transcription factor (TF) binding sites. Most existing prediction methods are based on sequence information alone, ignoring biological contexts such as developmental stages and tissue types. Experimental methods to study in vivo binding, including ChIP-chip and ChIP-seq, can only study one transcription factor in a single cell type and under a specific condition in each experiment, and therefore cannot scale to determine the full set of regulatory interactions in mammalian transcriptional regulatory networks.

RESULTS

We developed a new computational approach, PIPES, for predicting tissue-specific TF binding. PIPES integrates in vitro protein binding microarrays (PBMs), sequence conservation and tissue-specific epigenetic (DNase I hypersensitivity) information. We demonstrate that PIPES improves over existing methods on distinguishing between in vivo bound and unbound sequences using ChIP-seq data for 11 mouse TFs. In addition, our predictions are in good agreement with current knowledge of tissue-specific TF regulation.

CONCLUSIONS

We provide a systematic map of computationally predicted tissue-specific binding targets for 284 mouse TFs across 55 tissue/cell types. Such comprehensive resource is useful for researchers studying gene regulation.

Developmental fate and cellular maturity encoded in human regulatory DNA landscapes

Cellular-state information between generations of developing cells may be propagated via regulatory regions. We report consistent patterns of gain and loss of DNase I-hypersensitive sites (DHSs) as cells progress from embryonic stem cells (ESCs) to terminal fates. DHS patterns alone convey rich information about cell fate and lineage relationships distinct from information conveyed by gene expression. Developing cells share a proportion of their DHS landscapes with ESCs; that proportion decreases continuously in each cell type as differentiation progresses, providing a quantitative benchmark of developmental maturity. Developmentally stable DHSs densely encode binding sites for transcription factors involved in autoregulatory feedback circuits. In contrast to normal cells, cancer cells extensively reactivate silenced ESC DHSs and those from developmental programs external to the cell lineage from which the malignancy derives. Our results point to changes in regulatory DNA landscapes as quantitative indicators of cell-fate transitions, lineage relationships, and dysfunction.

Decoding a signature-based model of transcription cofactor recruitment dictated by cardinal cis-regulatory elements in proximal promoter regions

Genome-wide maps of DNase I hypersensitive sites (DHSs) reveal that most human promoters contain perpetually active cis-regulatory elements between −150 bp and +50 bp (−150/+50 bp) relative to the transcription start site (TSS). Transcription factors (TFs) recruit cofactors (chromatin remodelers, histone/protein-modifying enzymes, and scaffold proteins) to these elements in order to organize the local chromatin structure and coordinate the balance of post-translational modifications nearby, contributing to the overall regulation of transcription. However, the rules of TF-mediated cofactor recruitment to the −150/+50 bp promoter regions remain poorly understood. Here, we provide evidence for a general model in which a series of cis-regulatory elements (here termed ‘cardinal’ motifs) prefer acting individually, rather than in fixed combinations, within the −150/+50 bp regions to recruit TFs that dictate cofactor signatures distinctive of specific promoter subsets. Subsequently, human promoters can be subclassified based on the presence of cardinal elements and their associated cofactor signatures. In this study, furthermore, we have focused on promoters containing the nuclear respiratory factor 1 (NRF1) motif as the cardinal cis-regulatory element and have identified the pervasive association of NRF1 with the cofactor lysine-specific demethylase 1 (LSD1/KDM1A). This signature might be distinctive of promoters regulating nuclear-encoded mitochondrial and other particular genes in at least some cells. Together, we propose that decoding a signature-based, expanded model of control at proximal promoter regions should lead to a better understanding of coordinated regulation of gene transcription.

Human cells exploit different mechanisms to coordinate the expression of both protein-coding and non-coding RNAs. Elucidating these mechanisms is essential to understanding normal physiology and disease. In our attempt to identify new regulatory layers acting particularly at proximal promoters, we have computationally analyzed the genomic sequences located from −150 bp to +50 bp relative to the transcriptional start site (TSS), which are often at the center of ‘open’ chromatin regions in human promoters. We have confirmed the presence of a series of cis-regulatory elements (here referred to as ‘cardinal’ motifs) that show a strong preference for these short regions. Interestingly, these elements tend to act independently rather than in fixed combinations. Therefore, we propose that they confer unique regulatory features to the human promoter subsets that contain each of these particular elements. In agreement with this model, we have identified a large repertoire of preferential partnerships between transcription factors recognizing cardinal motifs and their associated proteins (cofactors), thus decoding a signature-based model that distinguishes distinctive regulatory types of promoters based on cardinal motifs. These signatures may underlie a new layer of transcriptional regulation to orchestrate coordinated gene expression in human promoters.

Coupling transcription factor occupancy to nucleosome architecture with DNase-FLASH

It is currently not possible to resolve the genome-wide relationship of transcription factors (TFs) and nucleosomes at the level of individual chromatin templates despite rapidly increasing data on TF and nucleosome occupancy in the human genome. Here we describe DNase I-released fragment-length analysis of hypersensitivity (DNase-FLASH), an approach that directly couples mapping of TF occupancy, via quantification of DNA microfragments released from individual TF recognition sites in regulatory DNA, to the surrounding nucleosome architecture, via analysis of larger DNA fragments, in a single assay. DNase-FLASH enables coupling of individual TF footprints to nucleosome occupancy, identifying TFs that precisely demarcate the regulatory DNA-nucleosome interface.

Gene transcription in the zebrafish embryo: regulators and networks

The precise spatial and temporal control of gene expression is a key process in the development, maintenance and regeneration of the vertebrate body. A substantial proportion of vertebrate genomes encode genes that control the transcription of the genetic information into mRNA. The zebrafish is particularly well suited to investigate gene regulatory networks underlying the control of gene expression during development due to the external development of its transparent embryos and the increasingly sophisticated tools for genetic manipulation available for this model system. We review here recent data on the analysis of cis-regulatory modules, transcriptional regulators and their integration into gene regulatory networks in the zebrafish, using the developing spinal cord as example.

A quantitative analysis of the impact on chromatin accessibility by histone modifications and binding of transcription factors in DNase I hypersensitive sites

It is known that chromatin features such as histone modifications and the binding of transcription factors exert a significant impact on the "openness" of chromatin. In this study, we present a quantitative analysis of the genome-wide relationship between chromatin features and chromatin accessibility in DNase I hypersensitive sites. We found that these features show distinct preference to localize in open chromatin. In order to elucidate the exact impact, we derived quantitative models to directly predict the "openness" of chromatin using histone modification features and transcription factor binding features, respectively. We show that these two types of features are highly predictive for chromatin accessibility in a statistical viewpoint. Moreover, our results indicate that these features are highly redundant and only a small number of features are needed to achieve a very high predictive power. Our study provides new insights into the true biological phenomena and the combinatorial effects of chromatin features to differential DNase I hypersensitivity.

Combining Hi-C data with phylogenetic correlation to predict the target genes of distal regulatory elements in human genome

Defining the target genes of distal regulatory elements (DREs), such as enhancer, repressors and insulators, is a challenging task. The recently developed Hi-C technology is designed to capture chromosome conformation structure by high-throughput sequencing, and can be potentially used to determine the target genes of DREs. However, Hi-C data are noisy, making it difficult to directly use Hi-C data to identify DRE-target gene relationships. In this study, we show that DREs-gene pairs that are confirmed by Hi-C data are strongly phylogenetic correlated, and have thus developed a method that combines Hi-C read counts with phylogenetic correlation to predict long-range DRE-target gene relationships. Analysis of predicted DRE-target gene pairs shows that genes regulated by large number of DREs tend to have essential functions, and genes regulated by the same DREs tend to be functionally related and co-expressed. In addition, we show with a couple of examples that the predicted target genes of DREs can help explain the causal roles of disease-associated single-nucleotide polymorphisms located in the DREs. As such, these predictions will be of importance not only for our understanding of the function of DREs but also for elucidating the causal roles of disease-associated noncoding single-nucleotide polymorphisms.

Comprehensive characterization of erythroid-specific enhancers in the genomic regions of human Krüppel-like factors

Background

Mapping of DNase I hypersensitive sites (DHSs) is a powerful tool to experimentally identify cis-regulatory elements (CREs). Among CREs, enhancers are abundant and predominantly act in driving cell-specific gene expression. Krüppel-like factors (KLFs) are a family of eukaryotic transcription factors. Several KLFs have been demonstrated to play important roles in hematopoiesis. However, transcriptional regulation of KLFs via CREs, particularly enhancers, in erythroid cells has been poorly understood.

Results

In this study, 23 erythroid-specific or putative erythroid-specific DHSs were identified by DNase-seq in the genomic regions of 17 human KLFs, and their enhancer activities were evaluated using dual-luciferase reporter (DLR) assay. Of the 23 erythroid-specific DHSs, the enhancer activities of 15 DHSs were comparable to that of the classical enhancer HS2 in driving minimal promoter (minP). Fifteen DHSs, some overlapping those that increased minP activities, acted as enhancers when driving the corresponding KLF promoters (KLF-Ps) in erythroid cells; of these, 10 DHSs were finally characterized as erythroid-specific KLF enhancers. These 10 erythroid-specific KLF enhancers were further confirmed using chromatin immunoprecipitation coupled to sequencing (ChIP-seq) data-based bioinformatic and biochemical analyses.

Conclusion

Our present findings provide a feasible strategy to extensively identify gene- and cell-specific enhancers from DHSs obtained by high-throughput sequencing, which will help reveal the transcriptional regulation and biological functions of genes in some specific cells.

The open chromatin landscape of Kaposi's sarcoma-associated herpesvirus

Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus which establishes latent infection in endothelial and B cells, as well as in primary effusion lymphoma (PEL). During latency, the viral genome exists as a circular DNA minichromosome (episome) and is packaged into chromatin analogous to human chromosomes. Only a small subset of promoters, those which drive latent RNAs, are active in latent episomes. In general, nucleosome depletion ("open chromatin") is a hallmark of eukaryotic regulatory elements such as promoters and transcriptional enhancers or insulators. We applied formaldehyde-assisted isolation of regulatory elements (FAIRE) followed by next-generation sequencing to identify regulatory elements in the KSHV genome and integrated these data with previously identified locations of histone modifications, RNA polymerase II occupancy, and CTCF binding sites. We found that (i) regions of open chromatin were not restricted to the transcriptionally defined latent loci; (ii) open chromatin was adjacent to regions harboring activating histone modifications, even at transcriptionally inactive loci; and (iii) CTCF binding sites fell within regions of open chromatin with few exceptions, including the constitutive LANA promoter and the vIL6 promoter. FAIRE-identified nucleosome depletion was similar among B and endothelial cell lineages, suggesting a common viral genome architecture in all forms of latency.

Migraine genetics--a review: Part I

Significant progress in molecular genetics has advanced our understanding of the genetic basis of migraine. The fundamentals of molecular genetics, and the recent advances in the field, are important for clinicians to understand, as they provide a foundation for critical appraisal of the literature, unprecedented insights into the pathogenesis of the disorder, and reveal promising treatment targets for future drug development. This paper provides a primer of molecular genetics and will be followed by a companion paper on the genetic advances in migraine, the methodology of genome wide association studies, and the potential clinical implications.

Causes and consequences of chromatin variation between inbred mice

Variation at regulatory elements, identified through hypersensitivity to digestion by DNase I, is believed to contribute to variation in complex traits, but the extent and consequences of this variation are poorly characterized. Analysis of terminally differentiated erythroblasts in eight inbred strains of mice identified reproducible variation at approximately 6% of DNase I hypersensitive sites (DHS). Only 30% of such variable DHS contain a sequence variant predictive of site variation. Nevertheless, sequence variants within variable DHS are more likely to be associated with complex traits than those in non-variant DHS, and variants associated with complex traits preferentially occur in variable DHS. Changes at a small proportion (less than 10%) of variable DHS are associated with changes in nearby transcriptional activity. Our results show that whilst DNA sequence variation is not the major determinant of variation in open chromatin, where such variants exist they are likely to be causal for complex traits.

Regulatory sites of the genome affect gene expression and complex traits, including disease susceptibility. Variable regulatory sites are potentially interesting because they are a likely cause of phenotypic variation, providing a bridge between sequence and transcriptional variation. In this paper we identify regions of the genome where DNA is not wrapped up in chromatin (hence potentially regulatory) in eight inbred strains of mice. We compare sites that vary among strains and compare them to non-variable sites. We show that more than half of variable sites cannot be attributed to local sequence variation. Functional consequences (in terms of readily detectable changes in gene expression) are associated with less than 10% of variable DNase I hypersensitive sites. We show that variable sites are enriched for sequence variants contributing to complex traits in mice.

A detailed protocol for formaldehyde-assisted isolation of regulatory elements (FAIRE)

Nucleosome displacement is a key event in the regulation of gene expression in the eukaryotic genome. This unit details an approach called Formaldehyde-Assisted Isolation of Regulatory Elements (FAIRE) for isolating nucleosome-depleted regions. FAIRE does not rely on the use of antibodies or enzymes, and has proven successful in most eukaryotic cells and tissues. The set of regulatory elements enriched by FAIRE is similar to those identified through DNase hypersensitivity. The enriched fragments can be detected by quantitative PCR, tiling DNA microarrays, or next-generation sequencing. Although the signal-to-noise ratio is typically lower than that observed for DNase assays, FAIRE has high sample-to-sample reproducibility, requires very low amounts of input material, is inexpensive, is amenable to high-throughput adaptations, and is a relatively simple procedure with a high rate of success, even for those without extensive experience in molecular biology protocols.

Conserved non-coding elements and cis regulation: actions speak louder than words

It is a truth (almost) universally acknowledged that conserved non-coding genomic sequences function in the cis regulation of neighbouring genes. But is this a misconception? The literature is strewn with examples of conserved non-coding sequences being able to drive reporter expression, but the extent to which such sequences are actually used endogenously in vivo is only now being rigorously explored using unbiased genome-scale approaches. Here, we review the emerging picture, examining the extent to which conserved non-coding sequences equivalently regulate gene expression in different species, or at different developmental stages, and how genomics approaches are revealing the relationship between sequence conservation and functional use of cis-regulatory elements.

DNase-seq predicts regions of rotational nucleosome stability across diverse human cell types

DNase-seq is primarily used to identify nucleosome-depleted DNase I hypersensitive (DHS) sites genome-wide that correspond to active regulatory elements. However, ∼40 yr ago it was demonstrated that DNase I also digests with a ∼10-bp periodicity around nucleosomes matching the exposure of the DNA minor groove as it wraps around histones. Here, we use DNase-seq data from 49 samples representing diverse cell types to reveal this digestion pattern at individual loci and predict genomic locations where nucleosome rotational positioning, the orientation of DNA with respect to the histone surface, is stably maintained. We call these regions DNase I annotated regions of nucleosome stability (DARNS). Compared to MNase-seq experiments, we show DARNS correspond well to annotated nucleosomes. Interestingly, many DARNS are positioned over only one side of annotated nucleosomes, suggesting that the periodic digestion pattern attenuates over the nucleosome dyad. DARNS reproduce the arrangement of nucleosomes around transcription start sites and are depleted at ubiquitous DHS sites. We also generated DARNS from multiple lymphoblast cell line (LCL) samples. We found that LCL DARNS were enriched at DHS sites present in most of the original 49 samples but absent in LCLs, while multi-cell-type DARNS were enriched at LCL-specific DHS sites. This indicates that variably open DHS sites are often occupied by rotationally stable nucleosomes in cell types where the DHS site is closed. DARNS provide additional information about precise DNA orientation within individual nucleosomes not available from other nucleosome positioning assays and contribute to understanding the role of chromatin in gene regulation.

Identification of a secondary promoter within the human B cell receptor component gene hCD79b

The human B cell-specific protein, CD79b (also known as Igβ and B29) constitutes an essential signal transduction component of the B cell receptor. Although its function is central to the triggering of B cell terminal differentiation in response to antigen stimulation, the transcriptional determinants that control CD79b gene expression remain poorly defined. In the present study, we explored these determinants using a series of hCD79b transgenic mouse models. Remarkably, we observed that the previously described hCD79b promoter along with its associated enhancer elements and first exon could be deleted without appreciable loss of hCD79b transcriptional activity or tissue specificity. In this deletion setting, a secondary promoter located within exon 2 maintained full levels and specificity of hCD79b transcription. Of note, this secondary promoter was also active, albeit at lower levels, in the wild-type hCD79b locus. The activity of the secondary promoter was dependent on the action(s) of a conserved sequence element mapping to a chromatin DNase I hypersensitive site located within intron 1. mRNA generated from this secondary promoter is predicted to encode an Igβ protein lacking a signal sequence and thus unable to serve normal B cell receptor function. Although the physiologic role of the hCD79b secondary promoter and its encoded protein remain unclear, the current data suggest that it has the capacity to play a role in normal as well as pathologic states in B cell proliferation and function.

From remote enhancers to gene regulation: charting the genome's regulatory landscapes

Vertebrate genes are characterized by the presence of cis-regulatory elements located at great distances from the genes they control. Alterations of these elements have been implicated in human diseases and evolution, yet little is known about how these elements interact with their surrounding sequences. A recent survey of the mouse genome with a regulatory sensor showed that the regulatory activities of these elements are not organized in a gene-centric manner, but instead are broadly distributed along chromosomes, forming large regulatory landscapes with distinct tissue-specific activities. A large genome-wide collection of expression data from this regulatory sensor revealed some basic principles of this complex genome regulatory architecture, including a substantial interplay between enhancers and other types of activities to modulate gene expression. We discuss the implications of these findings for our understanding of non-coding transcription, and of the possible consequences of structural genomic variations in disease and evolution.

High-throughput identification of long-range regulatory elements and their target promoters in the human genome

Enhancer elements are essential for tissue-specific gene regulation during mammalian development. Although these regulatory elements are often distant from their target genes, they affect gene expression by recruiting transcription factors to specific promoter regions. Because of this long-range action, the annotation of enhancer element-target promoter pairs remains elusive. Here, we developed a novel analysis methodology that takes advantage of Hi-C data to comprehensively identify these interactions throughout the human genome. To do this, we used a geometric distribution-based model to identify DNA-DNA interaction hotspots that contact gene promoters with high confidence. We observed that these promoter-interacting hotspots significantly overlap with known enhancer-associated histone modifications and DNase I hypersensitive sites. Thus, we defined thousands of candidate enhancer elements by incorporating these features, and found that they have a significant propensity to be bound by p300, an enhancer binding transcription factor. Furthermore, we revealed that their target genes are significantly bound by RNA Polymerase II and demonstrate tissue-specific expression. Finally, we uncovered that these elements are generally found within 1 Mb of their targets, and often regulate multiple genes. In total, our study presents a novel high-throughput workflow for confident, genome-wide discovery of enhancer-target promoter pairs, which will significantly improve our understanding of these regulatory interactions.

Insulator protein Su(Hw) recruits SAGA and Brahma complexes and constitutes part of Origin Recognition Complex-binding sites in the Drosophila genome

Despite increasing data on the properties of replication origins, molecular mechanisms underlying origin recognition complex (ORC) positioning in the genome are still poorly understood. The Su(Hw) protein accounts for the activity of best-studied Drosophila insulators. Here, we show that Su(Hw) recruits the histone acetyltransferase complex SAGA and chromatin remodeler Brahma to Su(Hw)-dependent insulators, which gives rise to regions with low nucleosome density and creates conditions for ORC binding. Depletion in Su(Hw) leads to a dramatic drop in the levels of SAGA, Brahma and ORC subunits and a significant increase in nucleosome density on Su(Hw)-dependent insulators, whereas artificial Su(Hw) recruitment itself is sufficient for subsequent SAGA, Brahma and ORC binding. In contrast to the majority of replication origins that associate with promoters of active genes, Su(Hw)-binding sites constitute a small proportion (6%) of ORC-binding sites that are localized preferentially in transcriptionally inactive chromatin regions termed BLACK and BLUE chromatin. We suggest that the key determinants of ORC positioning in the genome are DNA-binding proteins that constitute different DNA regulatory elements, including insulators, promoters and enhancers. Su(Hw) is the first example of such a protein.

Long-range transcriptional regulation by the p110 CUX1 homeodomain protein on the ENCODE array

Background

Overexpression of the Cut homeobox 1 gene, CUX1, inversely correlates with patient survival in breast cancers. Cell-based assays and molecular studies have revealed that transcriptional regulation by CUX1 involves mostly the proteolytically processed p110 isoform. As there is no antibody specific to p110 CUX1 only, an alternate strategy must be employed to identify its targets.

Results

We expressed physiological levels of a tagged-p110 CUX1 protein and performed chromatin affinity purification followed by hybridization on ENCODE and promoter arrays. Targets were validated by chromatin immunoprecipitation and transcriptional regulation by CUX1 was analyzed in expression profiling and RT-qPCR assays following CUX1 knockdown or p110 CUX1 overexpression. Approximately 47% and 14% of CUX1 binding sites were respectively mapped less than 4 Kbp, or more than 40 Kbp, away from a transcription start site. More genes exhibited changes in expression following CUX1 knockdown than p110 CUX1 overexpression. CUX1 directly activated or repressed 7.4% and 8.4% of putative targets identified on the ENCODE and promoter arrays respectively. This proportion increased to 11.2% for targets with 2 binding sites or more. Transcriptional repression was observed in a slightly higher proportion of target genes. The CUX1 consensus binding motif, ATCRAT, was found at 47.2% of the CUX1 binding sites, yet only 8.3% of the CUX1 consensus motifs present on the array were bound in vivo. The presence of a consensus binding motif did not have an impact on whether a target gene was repressed or activated. Interestingly, the distance between a binding site and a transcription start site did not significantly reduced the ability of CUX1 to regulate a target gene. Moreover, CUX1 not only was able to regulate the next adjacent gene, but also regulated the gene located beyond this one as well as the gene located further away in the opposite direction.

Conclusion

Our results demonstrate that p110 CUX1 can activate or repress transcription when bound at a distance and can regulate more than one gene on certain genomic loci.

Non-coding transcription at cis-regulatory elements: computational and experimental approaches

Mammalian genomes are pervasively transcribed, generating mostly RNAs with no coding potential that display different size, structure and interspecies sequence conservation. A prominent contribution to the ncRNA pool comes from the transcription of cis-regulatory elements, namely promoters, enhancers and locus control regions. While this phenomenon has been extensively documented, possible roles of such ncRNAs in gene regulation are still unclear. Addressing this issue will require experimental strategies dealing with the low abundance of enhancer-templated ncRNAs and aimed at specifically dissecting the relative role of transcription per se vs. RNA products. In this review, we first focus on the identification and characterization of cis-regulatory elements, highlighting the differences between emerging classes of ncRNAs associated to specific chromatin signatures. We then discuss current experimental strategies to dissect the function of nc transcription and computational approaches to the analysis and classification of regulatory sequences identified in next-generation sequencing experiments.

A NF-κB-dependent dual promoter-enhancer initiates the lipopolysaccharide-mediated transcriptional activation of the chicken lysozyme in macrophages

The transcriptional activation of the chicken lysozyme gene (cLys) by lipopolysaccharide (LPS) in macrophages is dependent on transcription of a LPS-Inducible Non-Coding RNA (LINoCR) triggering eviction of the CCCTC-binding factor (CTCF) from a negative regulatory element upstream of the lysozyme transcription start site. LINoCR is transcribed from a promoter originally characterized as a hormone response enhancer in the oviduct. Herein, we report the characterization of this cis-regulatory element (CRE). In activated macrophages, a 60 bp region bound by NF-κB, AP1 and C/EBPβ controls this CRE, which is strictly dependent on NF-κB binding for its activity in luciferase assays. Moreover, the serine/threonine kinase IKKα, known to be recruited by NF-κB to NF-κB-dependent genes is found at the CRE and within the transcribing regions of both cLys and LINoCR. Such repartition suggests a simultaneous promoter and enhancer activity of this CRE, initiating cLys transcriptional activation and driving CTCF eviction. This recruitment was transient despite persistence of both cLys transcription and NF-κB binding to the CRE. Finally, comparing cLys with other LPS-inducible genes indicates that IKKα detection within transcribing regions can be correlated with the presence of the elongating form of RNA polymerase II or concentrated in the 3′ end of the gene.

Latent enhancers activated by stimulation in differentiated cells

According to current models, once the cell has reached terminal differentiation, the enhancer repertoire is completely established and maintained by cooperatively acting lineage-specific transcription factors (TFs). TFs activated by extracellular stimuli operate within this predetermined repertoire, landing close to where master regulators are constitutively bound. Here, we describe latent enhancers, defined as regions of the genome that in terminally differentiated cells are unbound by TFs and lack the histone marks characteristic of enhancers but acquire these features in response to stimulation. Macrophage stimulation caused sequential binding of stimulus-activated and lineage-determining TFs to these regions, enabling deposition of enhancer marks. Once unveiled, many of these enhancers did not return to a latent state when stimulation ceased; instead, they persisted and mediated a faster and stronger response upon restimulation. We suggest that stimulus-specific expansion of the cis-regulatory repertoire provides an epigenomic memory of the exposure to environmental agents.

Modification of enhancer chromatin: what, how, and why?

Emergence of form and function during embryogenesis arises in large part through cell-type- and cell-state-specific variation in gene expression patterns, mediated by specialized cis-regulatory elements called enhancers. Recent large-scale epigenomic mapping revealed unexpected complexity and dynamics of enhancer utilization patterns, with 400,000 putative human enhancers annotated by the ENCODE project alone. These large-scale efforts were largely enabled through the understanding that enhancers share certain stereotypical chromatin features. However, an important question still lingers: what is the functional significance of enhancer chromatin modification? Here we give an overview of enhancer-associated modifications of histones and DNA and discuss enzymatic activities involved in their dynamic deposition and removal. We describe potential downstream effectors of these marks and propose models for exploring functions of chromatin modification in regulating enhancer activity during development.

Helper T-cell identity and evolution of differential transcriptomes and epigenomes

CD4(+) T cells are critical for the elimination of an immense array of microbial pathogens. Among the ways they accomplish this task is to generate progeny with specialized, characteristic patterns of gene expression. From this perspective, helper cells can be viewed as pluripotent precursors that adopt distinct cell fates. Although there are aspects of helper cell differentiation that can be modeled as a classic cell fate commitment, CD4(+) T cells also maintain considerable flexibility in their transcriptional program. This makes sense in terms of host defense, but raises the question of how these remarkable cells balance both these requirements, a high degree of specific gene expression and the capacity for plasticity. In this review, we discuss recent advances in our understanding of CD4(+) T-cell specification, focusing on how genomic perspectives have influenced our views of these processes. The relative contributions of sensors of the cytokine milieu, especially the signal transducer and activator of transcription family transcription factors, 'master regulators', and other transcription factors are considered as they relate to the helper cell transcriptome and epigenome.

Neurokinin B causes acute GnRH secretion and repression of GnRH transcription in GT1-7 GnRH neurons

Genetic studies in human patients with idiopathic hypogonadotropic hypogonadism (IHH) identified mutations in the genes that encode neurokinin B (NKB) and the neurokinin 3 receptor (NK3R). However, determining the mechanism whereby NKB regulates gonadotropin secretion has been difficult because of conflicting results from in vivo studies investigating the luteinizing hormone (LH) response to senktide, a NK3R agonist. NK3R is expressed in a subset of GnRH neurons and in kisspeptin neurons that are known to regulate GnRH secretion. Thus, one potential source of inconsistency is that NKB could produce opposing direct and indirect effects on GnRH secretion. Here, we employ the GT1-7 cell model to elucidate the direct effects of NKB on GnRH neuron function. We find that GT1-7 cells express NK3R and respond to acute senktide treatment with c-Fos induction and increased GnRH secretion. In contrast, long-term senktide treatment decreased GnRH secretion. Next, we focus on the examination of the mechanism underlying the long-term decrease in secretion and determine that senktide treatment represses transcription of GnRH. We further show that this repression of GnRH transcription may involve enhanced c-Fos protein binding at novel activator protein-1 (AP-1) half-sites identified in enhancer 1 and the promoter, as well as chromatin remodeling at the promoter of the GnRH gene. These data indicate that NKB could directly regulate secretion from NK3R-expressing GnRH neurons. Furthermore, whether the response is inhibitory or stimulatory toward GnRH secretion could depend on the history or length of exposure to NKB because of a repressive effect on GnRH transcription.

Personal and population genomics of human regulatory variation

The characteristics and evolutionary forces acting on regulatory variation in humans remains elusive because of the difficulty in defining functionally important noncoding DNA. Here, we combine genome-scale maps of regulatory DNA marked by DNase I hypersensitive sites (DHSs) from 138 cell and tissue types with whole-genome sequences of 53 geographically diverse individuals in order to better delimit the patterns of regulatory variation in humans. We estimate that individuals likely harbor many more functionally important variants in regulatory DNA compared with protein-coding regions, although they are likely to have, on average, smaller effect sizes. Moreover, we demonstrate that there is significant heterogeneity in the level of functional constraint in regulatory DNA among different cell types. We also find marked variability in functional constraint among transcription factor motifs in regulatory DNA, with sequence motifs for major developmental regulators, such as HOX proteins, exhibiting levels of constraint comparable to protein-coding regions. Finally, we perform a genome-wide scan of recent positive selection and identify hundreds of novel substrates of adaptive regulatory evolution that are enriched for biologically interesting pathways such as melanogenesis and adipocytokine signaling. These data and results provide new insights into patterns of regulatory variation in individuals and populations and demonstrate that a large proportion of functionally important variation lies beyond the exome.

Linking disease associations with regulatory information in the human genome

Genome-wide association studies have been successful in identifying single nucleotide polymorphisms (SNPs) associated with a large number of phenotypes. However, an associated SNP is likely part of a larger region of linkage disequilibrium. This makes it difficult to precisely identify the SNPs that have a biological link with the phenotype. We have systematically investigated the association of multiple types of ENCODE data with disease-associated SNPs and show that there is significant enrichment for functional SNPs among the currently identified associations. This enrichment is strongest when integrating multiple sources of functional information and when highest confidence disease-associated SNPs are used. We propose an approach that integrates multiple types of functional data generated by the ENCODE Consortium to help identify “functional SNPs” that may be associated with the disease phenotype. Our approach generates putative functional annotations for up to 80% of all previously reported associations. We show that for most associations, the functional SNP most strongly supported by experimental evidence is a SNP in linkage disequilibrium with the reported association rather than the reported SNP itself. Our results show that the experimental data sets generated by the ENCODE Consortium can be successfully used to suggest functional hypotheses for variants associated with diseases and other phenotypes.

Spermiogenic Germ Cell Phase-Specific DNA Damage Following Cyclophosphamide Exposure

The production of genetically competent spermatozoa is essential for normal embryo development. The chemotherapeutic drug cyclophosphamide creates cross-links and DNA strand breaks in many cell types, including germ cells. This study assessed the phase specificity of the susceptibility of spermiogenic germ cells to genetic damage induced by cyclophosphamide. Adult male rats were given cyclophosphamide using one of four schedules: 1) high dose/acute- day 1, 100 mg/kg; 2) low dose/subchronic, 4 days-days 1-4, 6.0 mg/kg/d; 3) high dose/subchronic, 4 days-day 1, 100 mg/kg, and days 2-4, 50 mg/kg/d; and 4) low dose/chronic-daily, 6.0 mg/kg/d for 14-28 days. To capture cauda epididymal spermatozoa exposed to cyclophosphamide during late, mid-, and early spermiogenesis, animals were sacrificed on days 14, 21, and 28, respectively. Spermatozoa were analyzed for DNA strand breaks using the comet assay. No dramatic increases in damage were seen after high-dose/acute exposure to cyclophosphamide. Subchronic exposure showed a dose-related increase in DNA damage; maximal damage, as demonstrated by comet tail parameters, was seen after 21 days, reflecting an increased susceptibility of step 9-14 spermatids. Low-dose chronic exposure to cyclophosphamide induced DNA damage, which reached a plateau by day 21. The magnitude of damage at all time points after low-dose chronic exposure was much greater than that following low-dose exposure for 4 days, indicating an accumulation of damage over time. Thus, the DNA damage induced by cyclophosphamide is germ cell phase-specific. The most damaging effects of cyclophosphamide occurred during a key point of sperm chromatin remodeling (histone hyperacetylation and transition protein deposition). We speculate that strand breaks disrupt chromatin remodeling, hence affecting chromatin structure and embryo development.

Isolation of nuclei for use in genome-wide DNase hypersensitivity assays to probe chromatin structure

DNase hypersensitivity (DHS) analysis coupled with high-throughput DNA sequencing (DNase-seq) has emerged as a powerful tool to analyze chromatin accessibility and identify regulatory sequences in genomic DNA on a global scale. In this method, intact nuclei are isolated from fresh tissue or cultured cells and then subjected to limited digestion using DNase I. The resulting short DNA fragments released by DNase digestion, which correspond to regions of open chromatin structure, are subsequently purified and identified by high throughput next generation DNA sequencing. This chapter describes methods used to isolate intact nuclei from mouse liver suitable for DNase-seq studies. The following chapter presents a detailed protocol for DNase I digestion of liver nuclei followed by the isolation of DNase-released fragments for sequencing and genome-wide mapping of DHS sites.

Gene regulation

This review concisely highlights the complexity of regulatory events. It provides examples of how interconnectivity of regulatory hubs could maintain transcriptional synergy and orchestrate the proper spatial and temporal patterns of gene expression.

Upstream distal regulatory elements contact the Lmo2 promoter in mouse erythroid cells

The Lim domain only 2 (Lmo2) gene encodes a transcriptional cofactor critical for the development of hematopoietic stem cells. Several distal regulatory elements have been identified upstream of the Lmo2 gene in the human and mouse genomes that are capable of enhancing reporter gene expression in erythroid cells and may be responsible for the high level transcription of Lmo2 in the erythroid lineage. In this study we investigate how these elements regulate transcription of Lmo2 and whether or not they function cooperatively in the endogenous context. Chromosome conformation capture (3C) experiments show that chromatin-chromatin interactions exist between upstream regulatory elements and the Lmo2 promoter in erythroid cells but that these interactions are absent from kidney where Lmo2 is transcribed at twelve fold lower levels. Specifically, long range chromatin-chromatin interactions occur between the Lmo2 proximal promoter and two broad regions, 3–31 and 66–105 kb upstream of Lmo2, which we term the proximal and distal control regions for Lmo2 (pCR and dCR respectively). Each of these regions is bound by several transcription factors suggesting that multiple regulatory elements cooperate in regulating high level transcription of Lmo2 in erythroid cells. Binding of CTCF and cohesin which support chromatin loops at other loci were also found within the dCR and at the Lmo2 proximal promoter. Intergenic transcription occurs throughout the dCR in erythroid cells but not in kidney suggesting a role for these intergenic transcripts in regulating Lmo2, similar to the broad domain of intergenic transcription observed at the human β-globin locus control region. Our data supports a model in which the dCR functions through a chromatin looping mechanism to contact and enhance Lmo2 transcription specifically in erythroid cells. Furthermore, these chromatin loops are supported by the cohesin complex recruited to both CTCF and transcription factor bound regions.

Chromatin remodeling in cancer: a gateway to regulate gene transcription

Cancer cells are remarkably adaptive to diverse survival strategies, probably due to its ability to interpret signaling cues differently than the normal cells. It appears as if cancer cells are constantly sampling, selecting and adapting signaling pathways to favor its proliferation. This process of successful adaptive evolution eventually renders a retractile nature to therapeutic regimens, fueling to the process of cancer progression. Based on plethora of available information, it is now evident that multiple signaling pathways eventually converge, perhaps, in a tempo-spatial manner, onto DNA template-dependent dynamic processes. Considering the complexity and packaging of eukaryotic genome, this process involves energy-dependent sub-events mediated by chromatin remodelers. Chromatin remodeler proteins function as gatekeepers and constitute a major determinant of accessibility of accessory factors to nucleosome DNA, allowing a wide repertoire of biological functions. And thus, aberrant expression or epigenetic modulation of remodeler proteins confers a unique ability to cancer cells to reprogram its genome for the maintenance of oncogenic phenotypes. Cancer cells can uniquely select a multi-subunit remodeler proteome for oncogenic advantage. This review summarizes our current understanding and importance of remodeler and chromatin proteins in cancer biology and also highlights the paradoxical role of proteins with or without dual-regulator functions. It is our hope that an in-depth understanding of these events is likely to provide a next set of opportunities for novel strategies for targeted cancer therapeutics.

NCBI Epigenomics: what's new for 2013

The Epigenomics resource at the National Center for Biotechnology Information (NCBI) has been created to serve as a comprehensive public repository for whole-genome epigenetic data sets (www.ncbi.nlm.nih.gov/epigenomics). We have constructed this resource by selecting the subset of epigenetics-specific data from the Gene Expression Omnibus (GEO) database and then subjecting them to further review and annotation. Associated data tracks can be viewed using popular genome browsers or downloaded for local analysis. We have performed extensive user testing throughout the development of this resource, and new features and improvements are continuously being implemented based on the results. We have made substantial usability improvements to user interfaces, enhanced functionality, made identification of data tracks of interest easier and created new tools for preliminary data analyses. Additionally, we have made efforts to enhance the integration between the Epigenomics resource and other NCBI databases, including the Gene database and PubMed. Data holdings have also increased dramatically since the initial publication describing the NCBI Epigenomics resource and currently consist of >3700 viewable and downloadable data tracks from 955 biological sources encompassing five well-studied species. This updated manuscript highlights these changes and improvements.

Long distance relationships: enhancer-promoter communication and dynamic gene transcription

The three-dimensional regulation of gene transcription involves loop formation between enhancer and promoter elements, controlling spatiotemporal gene expression in multicellular organisms. Enhancers are usually located in non-coding DNA and can activate gene transcription by recruiting transcription factors, chromatin remodeling factors and RNA Polymerase II. Research over the last few years has revealed that enhancers have tell-tale characteristics that facilitate their detection by several approaches, although the hallmarks of enhancers are not always uniform. Enhancers likely play an important role in the activation of genes by functioning as a primary point of contact for transcriptional activators, and by making physical contact with gene promoters often by means of a chromatin loop. Although numerous transcriptional regulators participate in the formation of chromatin loops that bring enhancers into proximity with promoters, the mechanism(s) of enhancer-promoter connectivity remain enigmatic. Here we discuss enhancer function, review some of the many proteins shown to be involved in establishing enhancer-promoter loops, and describe the dynamics of enhancer-promoter contacts during development, differentiation and in specific cell types.

ChIP-seq and beyond: new and improved methodologies to detect and characterize protein-DNA interactions

Chromatin immunoprecipitation experiments followed by sequencing (ChIP-seq) detect protein-DNA binding events and chemical modifications of histone proteins. Challenges in the standard ChIP-seq protocol have motivated recent enhancements in this approach, such as reducing the number of cells that are required and increasing the resolution. Complementary experimental approaches - for example, DNaseI hypersensitive site mapping and analysis of chromatin interactions that are mediated by particular proteins - provide additional information about DNA-binding proteins and their function. These data are now being used to identify variability in the functions of DNA-binding proteins across genomes and individuals. In this Review, I describe the latest advances in methods to detect and functionally characterize DNA-bound proteins.

Identifying and characterizing regulatory sequences in the human genome with chromatin accessibility assays

After finishing a human genome reference sequence in 2002, the genomics community has turned to the task of interpreting it. A primary focus is to identify and characterize not only protein-coding genes, but all functional elements in the genome. The effort includes both individual investigators and large-scale projects like the Encyclopedia of DNA Elements (ENCODE) project. As part of the ENCODE project, several groups have identified millions of regulatory elements in hundreds of human cell-types using DNase-seq and FAIRE-seq experiments that detect regions of nucleosome-free open chromatin. ChIP-seq experiments have also been used to discover transcription factor binding sites and map histone modifications. Nearly all identified elements are found in non-coding DNA, hypothesizing a function for previously unannotated sequence. In this review, we provide an overview of the ENCODE effort to define regulatory elements, summarize the main results, and discuss implications of the millions of regulatory elements distributed throughout the genome.

Chromatin accessibility reveals insights into androgen receptor activation and transcriptional specificity

BACKGROUND

Epigenetic mechanisms such as chromatin accessibility impact transcription factor binding to DNA and transcriptional specificity. The androgen receptor (AR), a master regulator of the male phenotype and prostate cancer pathogenesis, acts primarily through ligand-activated transcription of target genes. Although several determinants of AR transcriptional specificity have been elucidated, our understanding of the interplay between chromatin accessibility and AR function remains incomplete.

RESULTS

We used deep sequencing to assess chromatin structure via DNase I hypersensitivity and mRNA abundance, and paired these datasets with three independent AR ChIP-seq datasets. Our analysis revealed qualitative and quantitative differences in chromatin accessibility that corresponded to both AR binding and an enrichment of motifs for potential collaborating factors, one of which was identified as SP1. These quantitative differences were significantly associated with AR-regulated mRNA transcription across the genome. Base-pair resolution of the DNase I cleavage profile revealed three distinct footprinting patterns associated with the AR-DNA interaction, suggesting multiple modes of AR interaction with the genome.

CONCLUSIONS

In contrast with other DNA-binding factors, AR binding to the genome does not only target regions that are accessible to DNase I cleavage prior to hormone induction. AR binding is invariably associated with an increase in chromatin accessibility and, consequently, changes in gene expression. Furthermore, we present the first in vivo evidence that a significant fraction of AR binds only to half of the full AR DNA motif. These findings indicate a dynamic quantitative relationship between chromatin structure and AR-DNA binding that impacts AR transcriptional specificity.

Chromatin immunoprecipitation: advancing analysis of nuclear hormone signaling

Recent decades have been filled with groundbreaking research in the field of endocrine hormone signaling. Pivotal events like the isolation and purification of the estrogen receptor, the cloning of glucocorticoid receptor cDNA, or dissemination of nuclear hormone receptor (NHR) DNA binding sequences are well recognized for their contributions. However, the novel genome-wide and gene-specific information obtained over the last decade describing NHR association with chromatin, cofactors, and epigenetic modifications, as well as their role in gene regulation, has been largely facilitated by the adaptation of the chromatin immunoprecipitation (ChIP) technique. Use of ChIP-based technologies has taken the field of hormone signaling from speculating about the transcription-enabling properties of acetylated chromatin and putative transcription (co-)factor genomic occupancy to demonstrating the detailed, stepwise mechanisms of factor binding and transcriptional initiation; from treating hormone-induced transcription as a steady-state event to understanding its dynamic and cyclic nature; from looking at the DNA sequences recognized by various DNA-binding domains in vitro to analyzing the cell-specific genome-wide pattern of nuclear receptor binding and interpreting its physiological implications. Not only have these events propelled hormone research, but, as some of the pioneering studies, have also contributed tremendously to the field of molecular endocrinology as a whole. In this review, we give a brief summary of some of the most important discoveries in hormone signaling using ChIP and other derivative techniques and speculate on what the future may hold.

Nucleosomes suppress spontaneous mutations base-specifically in eukaryotes

It is unknown how the composition and structure of DNA within the cell affect spontaneous mutations. Theory suggests that in eukaryotic genomes, nucleosomal DNA undergoes fewer C→T mutations because of suppressed cytosine hydrolytic deamination relative to nucleosome-depleted DNA. Comparative genomic analyses and a mutation accumulation experiment showed that nucleosome occupancy nearly eliminated cytosine deamination, resulting in an ~50% decrease of the C→T mutation rate in nucleosomal DNA. Furthermore, the rates of G→T and A→T mutations were also about twofold suppressed by nucleosomes. On the basis of these results, we conclude that nucleosome-dependent mutation spectra affect eukaryotic genome structure and evolution and may have implications for understanding the origin of mutations in cancers and in induced pluripotent stem cells.

The accessible chromatin landscape of the human genome

DNase I hypersensitive sites (DHSs) are markers of regulatory DNA and have underpinned the discovery of all classes of cis-regulatory elements including enhancers, promoters, insulators, silencers and locus control regions. Here we present the first extensive map of human DHSs identified through genome-wide profiling in 125 diverse cell and tissue types. We identify ∼2.9 million DHSs that encompass virtually all known experimentally validated cis-regulatory sequences and expose a vast trove of novel elements, most with highly cell-selective regulation. Annotating these elements using ENCODE data reveals novel relationships between chromatin accessibility, transcription, DNA methylation and regulatory factor occupancy patterns. We connect ∼580,000 distal DHSs with their target promoters, revealing systematic pairing of different classes of distal DHSs and specific promoter types. Patterning of chromatin accessibility at many regulatory regions is organized with dozens to hundreds of co-activated elements, and the transcellular DNase I sensitivity pattern at a given region can predict cell-type-specific functional behaviours. The DHS landscape shows signatures of recent functional evolutionary constraint. However, the DHS compartment in pluripotent and immortalized cells exhibits higher mutation rates than that in highly differentiated cells, exposing an unexpected link between chromatin accessibility, proliferative potential and patterns of human variation.

An expansive human regulatory lexicon encoded in transcription factor footprints

Regulatory factor binding to genomic DNA protects the underlying sequence from cleavage by DNase I, leaving nucleotide-resolution footprints. Using genomic DNase I footprinting across 41 diverse cell and tissue types, we detected 45 million transcription factor occupancy events within regulatory regions, representing differential binding to 8.4 million distinct short sequence elements. Here we show that this small genomic sequence compartment, roughly twice the size of the exome, encodes an expansive repertoire of conserved recognition sequences for DNA-binding proteins that nearly doubles the size of the human cis-regulatory lexicon. We find that genetic variants affecting allelic chromatin states are concentrated in footprints, and that these elements are preferentially sheltered from DNA methylation. High-resolution DNase I cleavage patterns mirror nucleotide-level evolutionary conservation and track the crystallographic topography of protein-DNA interfaces, indicating that transcription factor structure has been evolutionarily imprinted on the human genome sequence. We identify a stereotyped 50-base-pair footprint that precisely defines the site of transcript origination within thousands of human promoters. Finally, we describe a large collection of novel regulatory factor recognition motifs that are highly conserved in both sequence and function, and exhibit cell-selective occupancy patterns that closely parallel major regulators of development, differentiation and pluripotency.

An integrated encyclopedia of DNA elements in the human genome

The human genome encodes the blueprint of life, but the function of the vast majority of its nearly three billion bases is unknown. The Encyclopedia of DNA Elements (ENCODE) project has systematically mapped regions of transcription, transcription factor association, chromatin structure and histone modification. These data enabled us to assign biochemical functions for 80% of the genome, in particular outside of the well-studied protein-coding regions. Many discovered candidate regulatory elements are physically associated with one another and with expressed genes, providing new insights into the mechanisms of gene regulation. The newly identified elements also show a statistical correspondence to sequence variants linked to human disease, and can thereby guide interpretation of this variation. Overall, the project provides new insights into the organization and regulation of our genes and genome, and is an expansive resource of functional annotations for biomedical research.

Systematic localization of common disease-associated variation in regulatory DNA

Genome-wide association studies have identified many noncoding variants associated with common diseases and traits. We show that these variants are concentrated in regulatory DNA marked by deoxyribonuclease I (DNase I) hypersensitive sites (DHSs). Eighty-eight percent of such DHSs are active during fetal development and are enriched in variants associated with gestational exposure-related phenotypes. We identified distant gene targets for hundreds of variant-containing DHSs that may explain phenotype associations. Disease-associated variants systematically perturb transcription factor recognition sequences, frequently alter allelic chromatin states, and form regulatory networks. We also demonstrated tissue-selective enrichment of more weakly disease-associated variants within DHSs and the de novo identification of pathogenic cell types for Crohn's disease, multiple sclerosis, and an electrocardiogram trait, without prior knowledge of physiological mechanisms. Our results suggest pervasive involvement of regulatory DNA variation in common human disease and provide pathogenic insights into diverse disorders.