DNaseI 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 choreographed with dozens to hundreds of co-activated elements, and the trans-cellular DNaseI sensitivity pattern at a given region can predict cell type-specific functional behaviors. 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.

Epigenome editing with the CRISPR/Cas9 platform is a promising technology to modulate gene expression to direct cell phenotype and to dissect the causal epigenetic mechanisms of gene regulation. Fusions of the nuclease-inactive dCas9 to the KRAB repressor (dCas9-KRAB) can silence target gene expression, but the genome-wide specificity and the extent of heterochromatin formation catalyzed by dCas9-KRAB is not known. We targeted dCas9-KRAB to the HS2 enhancer, a distal regulatory element that orchestrates expression of multiple globin genes. Genome-wide analyses demonstrated that localization of dCas9-KRAB to HS2 specifically induced H3K9 tri-methylation (H3K9me3) at the enhancer and reduced the chromatin accessibility of both the enhancer and its promoter targets. Targeted epigenetic modification of HS2 silenced the expression of multiple globin genes, with minimal off-target changes in gene expression. These results demonstrate that repression mediated by dCas9-KRAB is sufficiently specific to disrupt the activity of individual enhancers via local modification of the epigenome.

Large genome-mapping consortia and thousands of genome-wide association studies have identified non-protein coding elements in the genome as a having a central role in various biological processes. However, decoding the function of the millions of putative regulatory elements discovered in these studies remains challenging. CRISPR–Cas9-based epigenome editing technologies have enabled precise perturbation of the activity of specific regulatory elements. Here we describe CRISPR–Cas9-based epigenomic regulatory element screening (CERES) for improved high-throughput screening of regulatory element activity within the native genomic context. Using dCas9KRAB repressor and dCas9p300 activator constructs and lentiviral sgRNA libraries targeting DNase I hypersensitive sites surrounding a gene of interest, we perform both loss- and gain-of-function screens to identify regulatory elements for the β-globin and the HER2 loci in human cells. CERES readily identified known and novel regulatory elements, some of which were dependent on cell type or direction of perturbation. This technology allows the high-throughput functional annotation of putative regulatory elements in their native chromosomal context.

Genome-wide association studies (GWAS) have identified over 100 risk loci for schizophrenia, but the causal mechanisms remain largely unknown. We performed a transcriptome-wide association study (TWAS) integrating a schizophrenia GWAS of 79,845 individuals from the Psychiatric Genomics Consortium with expression data from brain, blood, and adipose tissues across 3,693 primarily control individuals. We identified 157 TWAS significant genes, of which 35 did not overlap a known GWAS locus. 42/157 genes were associated to specific chromatin features measured in independent samples, highlighting potential regulatory targets for follow-up. Suppression of one identified susceptibility gene, MAPK3, in zebrafish showed a significant effect on neurodevelopmental phenotypes. Expression and splicing from brain captured the majority of the TWAS effect across all genes. This large-scale connection of associations to target genes, tissues, and regulatory features is an essential step in moving towards a mechanistic understanding of GWAS.

Summary Most human transcription factors bind a small subset of potential genomic sites and often use different subsets in different cell types. To identify mechanisms that govern cell type-specific transcription factor binding, we used an integrative approach to study estrogen receptor α (ER). We found that ER exhibits two distinct modes of binding. Shared sites, bound in multiple cell types, are characterized by high affinity estrogen response elements (EREs), inaccessible chromatin and a lack of DNA methylation, while cell-specific sites are characterized by a lack of EREs, co-occurrence with other transcription factors and cell type-specific chromatin accessibility and DNA methylation. These observations enabled accurate quantitative models of ER binding that suggest tethering of ER to one-third of cell-specific sites. The distinct properties of cell-specific binding were also observed with glucocorticoid receptor and for ER in primary mouse tissues, representing an elegant genomic encoding scheme for generating cell type-specific gene regulation.