R. Scott Hansen

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.

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.

Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

Faithful transmission of genetic material to daughter cells involves a characteristic temporal order of DNA replication, which may play a significant role in the inheritance of epigenetic states. We developed a genome-scale approach—Repli Seq—to map temporally ordered replicating DNA using massively parallel sequencing and applied it to study regional variation in human DNA replication time across multiple human cell types. The method requires as few as 8,000 cytometry-fractionated cells for a single analysis, and provides high-resolution DNA replication patterns with respect to both cell-cycle time and genomic position. We find that different cell types exhibit characteristic replication signatures that reveal striking plasticity in regional replication time patterns covering at least 50% of the human genome. We also identified autosomal regions with marked biphasic replication timing that include known regions of monoallelic expression as well as many previously uncharacterized domains. Comparison with high-resolution genome-wide profiles of DNaseI sensitivity revealed that DNA replication typically initiates within foci of accessible chromatin comprising clustered DNaseI hypersensitive sites, and that replication time is better correlated with chromatin accessibility than with gene expression. The data collectively provide a unique, genome-wide picture of the epigenetic compartmentalization of the human genome and suggest that cell-lineage specification involves extensive reprogramming of replication timing patterns.

Topologically associating domains are stable units of replication-timing regulation

Eukaryotic chromosomes replicate in a temporal order known as the replication-timing program. In mammals, replication timing is cell-type-specific with at least half the genome switching replication timing during development, primarily in units of 400–800 kilobases (‘replication domains’), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements. Early and late replication correlate, respectively, with open and closed three-dimensional chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, late replication correlates with lamina-associated domains (LADs). Recent Hi-C mapping has unveiled substructure within chromatin compartments called topologically associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to replication domains. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure. Here we localize boundaries of replication domains to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, replication domain boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure replication domain boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type-specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell-type-specific sub-nuclear compartmentalization and replication timing with developmentally stable structural domains and offer a unified model for large-scale chromosome structure and function.

Cell-type-specific replication initiation programs set fragility of the FRA3B fragile site

Common fragile sites have long been identified by cytogeneticists as chromosomal regions prone to breakage upon replication stress. They are increasingly recognized to be preferential targets for oncogene-induced DNA damage in pre-neoplastic lesions and hotspots for chromosomal rearrangements in various cancers. Common fragile site instability was attributed to the fact that they contain sequences prone to form secondary structures that may impair replication fork movement, possibly leading to fork collapse resulting in DNA breaks. Here we show, in contrast to this view, that the fragility of FRA3B—the most active common fragile site in human lymphocytes—does not rely on fork slowing or stalling but on a paucity of initiation events. Indeed, in lymphoblastoid cells, but not in fibroblasts, initiation events are excluded from a FRA3B core extending approximately 700 kilobases, which forces forks coming from flanking regions to cover long distances in order to complete replication. We also show that origins of the flanking regions fire in mid-S phase, leaving the site incompletely replicated upon fork slowing. Notably, FRA3B instability is specific to cells showing this particular initiation pattern. The fact that both origin setting and replication timing are highly plastic in mammalian cells explains the tissue specificity of common fragile site instability we observed. Thus, we propose that common fragile sites correspond to the latest initiation-poor regions to complete replication in a given cell type. For historical reasons, common fragile sites have been essentially mapped in lymphocytes. Therefore, common fragile site contribution to chromosomal rearrangements in tumours should be reassessed after mapping fragile sites in the cell type from which each tumour originates.

Conservation of trans-acting circuitry during mammalian regulatory evolution

The basic body plan and major physiological axes have been highly conserved during mammalian evolution, yet only a small fraction of the human genome sequence appears to be subject to evolutionary constraint. To quantify cis- versus trans-acting contributions to mammalian regulatory evolution, we performed genomic DNase I footprinting of the mouse genome across 25 cell and tissue types, collectively defining ∼8.6 million transcription factor (TF) occupancy sites at nucleotide resolution. Here we show that mouse TF footprints conjointly encode a regulatory lexicon that is ∼95% similar with that derived from human TF footprints. However, only ∼20% of mouse TF footprints have human orthologues. Despite substantial turnover of the cis-regulatory landscape, nearly half of all pairwise regulatory interactions connecting mouse TF genes have been maintained in orthologous human cell types through evolutionary innovation of TF recognition sequences. Furthermore, the higher-level organization of mouse TF-to-TF connections into cellular network architectures is nearly identical with human. Our results indicate that evolutionary selection on mammalian gene regulation is targeted chiefly at the level of trans-regulatory circuitry, enabling and potentiating cis-regulatory plasticity.

Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution

To study the evolutionary dynamics of regulatory DNA, we mapped >1.3 million deoxyribonuclease I–hypersensitive sites (DHSs) in 45 mouse cell and tissue types, and systematically compared these with human DHS maps from orthologous compartments. We found that the mouse and human genomes have undergone extensive cis-regulatory rewiring that combines branch-specific evolutionary innovation and loss with widespread repurposing of conserved DHSs to alternative cell fates, and that this process is mediated by turnover of transcription factor (TF) recognition elements. Despite pervasive evolutionary remodeling of the location and content of individual cis-regulatory regions, within orthologous mouse and human cell types the global fraction of regulatory DNA bases encoding recognition sites for each TF has been strictly conserved. Our findings provide new insights into the evolutionary forces shaping mammalian regulatory DNA landscapes. Mouse-to-human genomic comparisons illuminate conserved transcriptional programs despite regulatory rewiring. Rewiring the gene regulatory landscape DNAse I hypersensitive sites (DHSs) correlate with genomic locations that control where messenger RNA is to be produced. DHSs differ, depending on the cell type, developmental stage, and species. Viestra et al. compared mouse and human genome-wide DHS maps. Approximately one-third of the DHSs are conserved between the species, which separated approximately 550 million years ago. Most DHSs fell into tissue-specific cohorts; however, these were generally not conserved between the human and mouse. It seems that the majority of DHSs evolve because of changes in the sequence that gradually change how the region is regulated. Science, this issue p. 1007

Genetic variation in ICF syndrome: evidence for genetic heterogeneity.

ICF syndrome is a rare autosomal recessive immunoglobulin deficiency, sometimes combined with defective cellular immunity. Other features that are frequently observed in ICF syndrome patients include facial dysmorphism, developmental delay, and recurrent infections. The most diagnostic feature of ICF syndrome is the branching of chromosomes 1, 9, and 16 due to pericentromeric instability. Positional candidate cloning recently discovered the de novo DNA methyltransferase 3B (DNMT3B) as the responsible gene by identifying seven different mutations in nine ICF patients. DNMT3B specifically methylates repeat sequences adjacent to the centromeres of chromosome 1, 9, and 16. Our panel of 14 ICF patients was subjected to mutation analysis in the DNMT3B gene. Mutations in DNMT3B were discovered in only nine of our 14 ICF patients. Moreover, two ICF patients from consanguineous families who did not show autozygosity (i.e. homozygosity by descent) for the DNMT3B locus did not reveal DNMT3B mutations, suggesting genetic heterogeneity for this disease. Mutation analysis revealed 11 different mutations, including seven novel ones: eight different missense mutations, two different nonsense mutations, and a splice‐site mutation leading to the insertion of three aa’s. The missense mutations occurred in or near the catalytic domain of DNMT3B protein, indicating a possible interference with the normal functioning of the enzyme. However, none of the ICF patients was homozygous for a nonsense allele, suggesting that absence of this enzyme is not compatible with life. Compound heterozygosity for a missense and a nonsense mutation did not seem to correlate with a more severe phenotype. Hum Mutat 16:509–517, 2000.

Encoding PCR Products with Batch-stamps and Barcodes

Polymerase chain reaction (PCR) has become the mainstay of DNA sequence analysis. Yet there is always uncertainty concerning the source of the template DNA that gave rise to a particular PCR product. The risks of contamination, biased amplification, and product redundancy are especially high when limited amounts of template DNA are used. We have developed and applied molecular encoding principles to solve this source-uncertainty problem for DNA sequences generated by standard PCR. Batch-stamps specify the date and sample identity, and barcodes detect template redundancy. Our approach thus enables classification of each PCR-derived sequence as valid, contaminant, or redundant, and provides a measure of sequence diversity. We recommend that batch-stamps and barcodes be used when amplifying irreplaceable DNAs and cDNAs available for forensic, clinical, single cell, and ancient DNA analyses.