Anthony Shafer

Systematic Localization of Common Disease-Associated Variation in Regulatory DNA

Predictions of Genetic Disease Many genome-wide association studies (GWAS) have identified loci and variants associated with disease, but the ability to predict disease on the basis of these genetic variants remains small. Maurano et al. (p. 1190; see the Perspective by Schadt and Chang; see the cover) characterize the location of GWAS variants in the genome with respect to their proximity to regulatory DNA [marked by deoxyribonuclease I (DNase I) hypersensitive sites] by tissue type, disease, and enrichments in physiologically relevant transcription factor binding sites and networks. They found many noncoding disease associations in regulatory DNA, indicating tissue and developmental-specific regulatory roles for many common genetic variants and thus enabling links to be made between gene regulation and adult-onset disease. Genetic variants that have been associated with diseases are concentrated in regulatory regions of the genome. 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.

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.

Exonic Transcription Factor Binding Directs Codon Choice and Affects Protein Evolution

Transcription Factor Binding Sites Transcription factors (TFs) are proteins that bind to DNA to control gene transcription. Stergachis et al. (p. 1367; see the Perspective by Weatheritt and Babu) examined TF binding within the human genome in more than 80 cell types. Nearly 15% of coding regions simultaneously specify both amino acid sequence and TF recognition sites. The distribution of the TF binding sites evolutionarily constrains how codons within these regions can change, independent of encoded protein function. Thus, TF binding may represent a widespread and strong evolutionary force in coding regions. Transcription factor binding within protein-coding regions of DNA constrains how the protein can evolve. [Also see Perspective by Weatheritt and Babu] Genomes contain both a genetic code specifying amino acids and a regulatory code specifying transcription factor (TF) recognition sequences. We used genomic deoxyribonuclease I footprinting to map nucleotide resolution TF occupancy across the human exome in 81 diverse cell types. We found that ~15% of human codons are dual-use codons (“duons”) that simultaneously specify both amino acids and TF recognition sites. Duons are highly conserved and have shaped protein evolution, and TF-imposed constraint appears to be a major driver of codon usage bias. Conversely, the regulatory code has been selectively depleted of TFs that recognize stop codons. More than 17% of single-nucleotide variants within duons directly alter TF binding. Pervasive dual encoding of amino acid and regulatory information appears to be a fundamental feature of genome evolution.

High-throughput localization of functional elements by quantitative chromatin profiling

Identification of functional, noncoding elements that regulate transcription in the context of complex genomes is a major goal of modern biology. Localization of functionality to specific sequences is a requirement for genetic and computational studies. Here, we describe a generic approach, quantitative chromatin profiling, that uses quantitative analysis of in vivo chromatin structure over entire gene loci to rapidly and precisely localize cis-regulatory sequences and other functional modalities encoded by DNase I hypersensitive sites. To demonstrate the accuracy of this approach, we analyzed ∼300 kilobases of human genome sequence from diverse gene loci and cleanly delineated functional elements corresponding to a spectrum of classical cis-regulatory activities including enhancers, promoters, locus control regions and insulators as well as novel elements. Systematic, high-throughput identification of functional elements coinciding with DNase I hypersensitive sites will substantially expand our knowledge of transcriptional regulation and should simplify the search for noncoding genetic variation with phenotypic consequences.

Large-scale identification of sequence variants influencing human transcription factor occupancy in vivo

The function of human regulatory regions depends exquisitely on their local genomic environment and on cellular context, complicating experimental analysis of common disease- and trait-associated variants that localize within regulatory DNA. We use allelically resolved genomic DNase I footprinting data encompassing 166 individuals and 114 cell types to identify >60,000 common variants that directly influence transcription factor occupancy and regulatory DNA accessibility in vivo. The unprecedented scale of these data enables systematic analysis of the impact of sequence variation on transcription factor occupancy in vivo. We leverage this analysis to develop accurate models of variation affecting the recognition sites for diverse transcription factors and apply these models to discriminate nearly 500,000 common regulatory variants likely to affect transcription factor occupancy across the human genome. The approach and results provide a new foundation for the analysis and interpretation of noncoding variation in complete human genomes and for systems-level investigation of disease-associated variants.

Probing DNA shape and methylation state on a genomic scale with DNase I

DNA binding proteins find their cognate sequences within genomic DNA through recognition of specific chemical and structural features. Here we demonstrate that high-resolution DNase I cleavage profiles can provide detailed information about the shape and chemical modification status of genomic DNA. Analyzing millions of DNA backbone hydrolysis events on naked genomic DNA, we show that the intrinsic rate of cleavage by DNase I closely tracks the width of the minor groove. Integration of these DNase I cleavage data with bisulfite sequencing data for the same cell type’s genome reveals that cleavage directly adjacent to cytosine-phosphate-guanine (CpG) dinucleotides is enhanced at least eightfold by cytosine methylation. This phenomenon we show to be attributable to methylation-induced narrowing of the minor groove. Furthermore, we demonstrate that it enables simultaneous mapping of DNase I hypersensitivity and regional DNA methylation levels using dense in vivo cleavage data. Taken together, our results suggest a general mechanism by which CpG methylation can modulate protein–DNA interaction strength via the remodeling of DNA shape.

Role of DNA Methylation in Modulating Transcription Factor Occupancy

Although DNA methylation is commonly invoked as a mechanism for transcriptional repression, the extent to which it actively silences transcription factor (TF) occupancy sites in vivo is unknown. To study the role of DNA methylation in the active modulation of TF binding, we quantified the effect of DNA methylation depletion on the genomic occupancy patterns of CTCF, an abundant TF with known methylation sensitivity that is capable of autonomous binding to its target sites in chromatin. Here, we show that the vast majority (>98.5%) of the tens of thousands of unoccupied, methylated CTCF recognition sequences remain unbound upon abrogation of DNA methylation. The small fraction of sites that show methylation-dependent binding in vivo are in turn characterized by highly variable CTCF occupancy across cell types. Our results suggest that DNA methylation is not a primary groundskeeper of genomic TF landscapes, but rather a specialized mechanism for stabilizing intrinsically labile sites.

DNA methylation alone does not cause most cell-type selective transcription factor binding

Background DNA methylation at vertebrate promoters is associated with the repression of gene expression and is required for mammalian development [1]. The master genome regulator and transcription factor CTCF canonically exhibits methylation-sensitive binding in vitro ,a nd its in vivo occupancy across cell types anticorrelates with methylation at 41 % of sites genome-wide [2]. However, the genome harbors hundreds of thousands of CTCF recognition sequences, the majority of which are unbound in any cell type. Although most unbound recognition sequences harbor methylation in a given cell type, the causal role of DNA methylation in the abrogation of transcription factor binding at these potential binding sites is largely unknown [3,4]. Materials and methods Here we perform genome-wide occupancy profiling after both stable genetic and transient chemical inhibition of DNA methyltransferases in HCT116 cells lacking functional DNMT1 and DNMT3b and K562 erythroleukemia cells treated with 5-aza-2-deoxycytidine. Results We show that the vast majority of susceptible binding sites remain unoccupied upon depletion of DNA methylation. Stable loss of methylation in HCT116 cells does result in a minor increase of several thousand binding sites, largely corresponding to selective reactivation of binding sites from other cell types, especially from other malignant cell lines. Chemical inhibition in K562 cells results in a smaller but reproducible set of several hundred reactivated sites, the majority of which are corroborated by the genetic data. The narrow extent of reactivation implies that the majority of transcription factor recognition sequences are not competent for binding even in the absence of methylation, and that methylation accounts for only a modest fraction of cell-selective binding. These results offer new insight into the establishment of the global regulatory landscape, and offer an important perspective on the interpretation of disease- and trait-associated methylation differences in humans.

103 Probing DNA shape and methylation state on a genomic scale with DNase I

DNA binding proteins find their cognate sequences within genomic DNA through recognition of specific chemical and structural features. Here, we demonstrate that high-resolution DNase I cleavage profiles can provide detailed information about the shape and chemical modification status of genomic DNA. Analyzing millions of DNA-backbone hydrolysis events on naked genomic DNA, we show that the intrinsic rate of cleavage by DNase I closely tracks the width of the minor groove. Integration of these DNase I cleavage data with bisulfite sequencing data for the same cell type genome reveals that the cleavage directly adjacent to CpG dinucleotides is enhanced at least eight-fold by cytosine methylation. This phenomenon we show is attributable to methylation-induced narrowing of the minor groove. Furthermore, we demonstrate that it enables simultaneous mapping of DNase I hypersensitivity and regional DNA methylation levels using dense in vivo cleavage data. Taken together, our results suggest a general mechanism through which CpG methylation can modulate protein–DNA interaction strength via the remodeling of DNA shape.

DNA methylation modulates transcription factor occupancy chiefly at sites of high intrinsic cell-type variability

The nuclear genome of every cell harbors millions of unoccupied transcription factor (TF) recognition sequences that harbor methylated cytosines. Although DNA methylation is commonly invoked as a repressive mechanism, the extent to which it actively silences specific TF occupancy sites is unknown. To define the role of DNA methylation in modulating TF binding, we quantified the effect of DNA methyltransferase abrogation on the occupancy patterns of a ubiquitous TF capable of autonomous binding to its target sites in chromatin (CTCF). Here we show that the vast majority of unoccupied, methylated CTCF recognition sequences remain unbound upon depletion of DNA methylation. Rather, methylation-regulated binding is restricted to a small fraction of elements that exhibit high intrinsic variability in CTCF occupancy across cell types. Our results suggest that DNA methylation is not a major groundskeeper of genomic transcription factor occupancy landscapes, but rather a specialized mechanism for stabilizing epigenetically labile sites.