Jason D. Lieb

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.

The DNA-encoded nucleosome organization of a eukaryotic genome

Nucleosome organization is critical for gene regulation. In living cells this organization is determined by multiple factors, including the action of chromatin remodellers, competition with site-specific DNA-binding proteins, and the DNA sequence preferences of the nucleosomes themselves. However, it has been difficult to estimate the relative importance of each of these mechanisms in vivo, because in vivo nucleosome maps reflect the combined action of all influencing factors. Here we determine the importance of nucleosome DNA sequence preferences experimentally by measuring the genome-wide occupancy of nucleosomes assembled on purified yeast genomic DNA. The resulting map, in which nucleosome occupancy is governed only by the intrinsic sequence preferences of nucleosomes, is similar to in vivo nucleosome maps generated in three different growth conditions. In vitro, nucleosome depletion is evident at many transcription factor binding sites and around gene start and end sites, indicating that nucleosome depletion at these sites in vivo is partly encoded in the genome. We confirm these results with a micrococcal nuclease-independent experiment that measures the relative affinity of nucleosomes for ∼40,000 double-stranded 150-base-pair oligonucleotides. Using our in vitro data, we devise a computational model of nucleosome sequence preferences that is significantly correlated with in vivo nucleosome occupancy in Caenorhabditis elegans. Our results indicate that the intrinsic DNA sequence preferences of nucleosomes have a central role in determining the organization of nucleosomes in vivo.

ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia

Chromatin immunoprecipitation (ChIP) followed by high-throughput DNA sequencing (ChIP-seq) has become a valuable and widely used approach for mapping the genomic location of transcription-factor binding and histone modifications in living cells. Despite its widespread use, there are considerable differences in how these experiments are conducted, how the results are scored and evaluated for quality, and how the data and metadata are archived for public use. These practices affect the quality and utility of any global ChIP experiment. Through our experience in performing ChIP-seq experiments, the ENCODE and modENCODE consortia have developed a set of working standards and guidelines for ChIP experiments that are updated routinely. The current guidelines address antibody validation, experimental replication, sequencing depth, data and metadata reporting, and data quality assessment. We discuss how ChIP quality, assessed in these ways, affects different uses of ChIP-seq data. All data sets used in the analysis have been deposited for public viewing and downloading at the ENCODE (http://encodeproject.org/ENCODE/) and modENCODE (http://www.modencode.org/) portals.

Evidence for nucleosome depletion at active regulatory regions genome-wide

The identification of nuclease-hypersensitive sites in an active globin gene and in the 5′ regions of fruit fly heat shock genes first suggested that chromatin changes accompany gene regulation in vivo. Here we present evidence that the basic repeating units of eukaryotic chromatin, nucleosomes, are depleted from active regulatory elements throughout the Saccharomyces cerevisiae genome in vivo. We found that during rapid mitotic growth, the level of nucleosome occupancy is inversely proportional to the transcriptional initiation rate at the promoter. We also observed a partial loss of histone H3 and H4 tetramers from the coding regions of the most heavily transcribed genes. Alterations in the global transcriptional program caused by heat shock or a change in carbon source resulted in an increased nucleosome occupancy at repressed promoters, and a decreased nucleosome occupancy at promoters that became active. Nuclease-hypersensitive sites occur in species from yeast to humans and result from chromatin perturbation. Given the conservation of sequence and function among components of both chromatin and the transcriptional machinery, nucleosome depletion at promoters may be a fundamental feature of eukaryotic transcriptional regulation.

Promoter-specific binding of Rap1 revealed by genome-wide maps of protein-DNA association

We determined the distribution of repressor-activator protein 1 (Rap1) and the accessory silencing proteins Sir2, Sir3 and Sir4 in vivo on the entire yeast genome, at a resolution of 2 kb. Rap1 is central to the cellular economy during rapid growth, targeting 294 loci, about 5% of yeast genes, and participating in the activation of 37% of all RNA polymerase II initiation events in exponentially growing cells. Although the DNA sequence recognized by Rap1 is found in both coding and intergenic sequences, the binding of Rap1 to the genome was highly specific to intergenic regions with the potential to act as promoters. This global phenomenon, which may be a general characteristic of sequence-specific transcriptional factors, indicates the existence of a genome-wide molecular mechanism for marking promoter regions.

Unlocking the secrets of the genome

Despite the successes of genomics, little is known about how genetic information produces complex organisms. A look at the crucial functional elements of fly and worm genomes could change that. The National Human Genome Research Institutes modENCODE project (the model organism ENCyclopedia Of DNA Elements) was set up in 2007 with the goal of identifying all the sequence-based functional elements in the genomes of two important experimental organisms, Caenorhabditis elegans and Drosophila melanogaster. Armed with modENCODE data, geneticists will be able to undertake the comprehensive molecular studies of regulatory networks that hold the key to how complex multicellular organisms arise from the list of instructions coded in the genome. In this issue, modENCODE team members outline their plan of campaign. Data from the project are to be made available on http://www.modencode.org and elsewhere as the work progresses.

A map of open chromatin in human pancreatic islets

Tissue-specific transcriptional regulation is central to human disease. To identify regulatory DNA active in human pancreatic islets, we profiled chromatin by formaldehyde-assisted isolation of regulatory elements coupled with high-throughput sequencing (FAIRE-seq). We identified ∼80,000 open chromatin sites. Comparison of FAIRE-seq data from islets to that from five non-islet cell lines revealed ∼3,300 physically linked clusters of islet-selective open chromatin sites, which typically encompassed single genes that have islet-specific expression. We mapped sequence variants to open chromatin sites and found that rs7903146, a TCF7L2 intronic variant strongly associated with type 2 diabetes, is located in islet-selective open chromatin. We found that human islet samples heterozygous for rs7903146 showed allelic imbalance in islet FAIRE signals and that the variant alters enhancer activity, indicating that genetic variation at this locus acts in cis with local chromatin and regulatory changes. These findings illuminate the tissue-specific organization of cis-regulatory elements and show that FAIRE-seq can guide the identification of regulatory variants underlying disease susceptibility.

Nucleosome Dynamics Define Transcriptional Enhancers

Chromatin plays a central role in eukaryotic gene regulation. We performed genome-wide mapping of epigenetically marked nucleosomes to determine their position both near transcription start sites and at distal regulatory elements, including enhancers. In prostate cancer cells, where androgen receptor binds primarily to enhancers, we found that androgen treatment dismisses a central nucleosome present at androgen receptor binding sites that is flanked by a pair of marked nucleosomes. A new quantitative model built on the behavior of such nucleosome pairs correctly identified regions bound by the regulators of the immediate androgen response, including androgen receptor and FOXA1. More importantly, this model also correctly predicted previously unidentified binding sites for other transcription factors present after prolonged androgen stimulation, including OCT1 and NKX3-1. Therefore, quantitative modeling of enhancer structure provides a powerful predictive method to infer the identity of transcription factors involved in cellular responses to specific stimuli.

A Library of Yeast Transcription Factor Motifs Reveals a Widespread Function for Rsc3 in Targeting Nucleosome Exclusion at Promoters

The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially approximately 100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.

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.