Brentuximab Vedotin-Driven Microtubule Disruption Results in Endoplasmic Reticulum Stress Leading to Immunogenic Cell Death and Antitumor Immunity.

Brentuximab vedotin, a CD30-directed antibody-drug conjugate (ADC), is approved for clinical use in multiple CD30-expressing lymphomas. The cytotoxic payload component of brentuximab vedotin is monomethyl auristatin E (MMAE), a highly potent microtubule-disrupting agent. Preclinical results provided here demonstrate that treatment of cancer cells with brentuximab vedotin or free MMAE leads to a catastrophic disruption of the microtubule network eliciting a robust endoplasmic reticulum (ER) stress response that culminates in the induction of the classic hallmarks of immunogenic cell death (ICD). In accordance with the induction of ICD, brentuximab vedotin-killed lymphoma cells drove innate immune cell activation in vitro and in vivo. In the "gold-standard" test of ICD, vaccination of mice with brentuximab vedotin or free MMAE-killed tumor cells protected animals from tumor rechallenge; in addition, T cells transferred from previously vaccinated animals slowed tumor growth in immunodeficient mice. Immunity acquired from killed tumor cell vaccination was further amplified by the addition of PD-1 blockade. In a humanized model of CD30+ B-cell tumors, treatment with brentuximab vedotin drove the expansion and recruitment of autologous Epstein-Barr virus-reactive CD8+ T cells potentiating the activity of anti-PD-1 therapy. Together, these data support the ability of brentuximab vedotin and MMAE to drive ICD in tumor cells resulting in the activation of antigen-presenting cells and augmented T-cell immunity. These data provide a strong rationale for the clinical combination of brentuximab vedotin and other MMAE-based ADCs with checkpoint inhibitors.

Nonfucosylation of an anti-TIGIT antibody enhances FcγR engagement, driving innate immune activation and antitumor activity.

TIGIT is an immune checkpoint receptor expressed on activated and memory T cells, immunosuppressive T regulatory cells, and natural killer (NK) cells. TIGIT has emerged as an attractive target for antitumor therapies, due to its proposed immunosuppressive effects on lymphocyte function and T cell activation. We generated an anti-TIGIT monoclonal antibody (mAb) that binds with high affinity to human, non-human primate, and murine TIGIT and through multiple experimental methodologies demonstrated that checkpoint blockade alone is insufficient for antitumor activity. Generating anti-TIGIT mAbs with various Fc backbones we show that muting the Fc-Fcγ receptor (FcγR) interaction failed to drive antitumor activity, while mAbs with Fc functional backbones demonstrate substantial antitumor activity, mediated through activation of antigen-presenting cells (APCs), T cell priming, and NK-mediated depletion of suppressive Tregs and exhausted T cells. Further, nonfucosylation of the Fc backbone resulted in enhanced immune responses and antitumor activity relative to the intact IgG1 backbone. The improved activity correlated with the biased FcγR interaction profile of the nonfucosylated anti-TIGIT mAb, which supports that FcγRIIIa binding with decreased FcγRIIb binding favorably activates APCs and enhances tumor-specific CD8+ T cell responses. The anti-TIGIT mAbs with intact FcγR interacting backbones also demonstrated synergistic enhancement of other standard antitumor treatments, including anti-PD-1 treatment and a model monomethyl auristatin E antibody-drug conjugate. These findings highlight the importance of the anti-TIGIT mAb's Fc backbone to its antitumor activity and the extent to which this activity can be enhanced through nonfucosylation of the backbone.

SGN-B7H4V, an investigational vedotin ADC directed to the immune checkpoint ligand B7-H4, shows promising activity in preclinical models.

BACKGROUND: SGN-B7H4V is a novel investigational vedotin antibody-drug conjugate (ADC) comprising a B7-H4-directed human monoclonal antibody conjugated to the cytotoxic payload monomethyl auristatin E (MMAE) via a protease-cleavable maleimidocaproyl valine citrulline (mc-vc) linker. This vedotin linker-payload system has been clinically validated in multiple Food and Drug Administration approved agents including brentuximab vedotin, enfortumab vedotin, and tisotumab vedotin. B7-H4 is an immune checkpoint ligand with elevated expression on a variety of solid tumors, including breast, ovarian, and endometrial tumors, and limited normal tissue expression. SGN-B7H4V is designed to induce direct cytotoxicity against target cells by binding to B7-H4 on the surface of target cells and releasing the cytotoxic payload MMAE upon internalization of the B7-H4/ADC complex. METHODS: B7-H4 expression was characterized by immunohistochemistry across multiple solid tumor types. The ability of SGN-B7H4V to kill B7-H4-expressing tumor cells in vitro and in vivo in a variety of xenograft tumor models was also evaluated. Finally, the antitumor activity of SGN-B7H4V as monotherapy and in combination with an anti-programmed cell death-1 (PD-1) agent was evaluated using an immunocompetent murine B7-H4-expressing Renca tumor model. RESULTS: Immunohistochemistry confirmed B7-H4 expression across multiple solid tumors, with the highest prevalence in breast, endometrial, and ovarian tumors. In vitro, SGN-B7H4V killed B7-H4-expressing tumor cells by MMAE-mediated direct cytotoxicity and antibody-mediated effector functions including antibody-dependent cellular cytotoxicity and antibody-dependent cellular phagocytosis. In vivo, SGN-B7H4V demonstrated strong antitumor activity in multiple xenograft models of breast and ovarian cancer, including xenograft tumors with heterogeneous B7-H4 expression, consistent with the ability of vedotin ADCs to elicit a bystander effect. In an immunocompetent murine B7-H4-expressing tumor model, SGN-B7H4V drove robust antitumor activity as a monotherapy that was enhanced when combined with an anti-PD-1 agent. CONCLUSION: The immune checkpoint ligand B7-H4 is a promising molecular target expressed by multiple solid tumors. SGN-B7H4V demonstrates robust antitumor activity in preclinical models through multiple potential mechanisms. Altogether, these preclinical data support the evaluation of SGN-B7H4V as a monotherapy in the ongoing phase 1 study of SGN-B7H4V in advanced solid tumors (NCT05194072) and potential future clinical combinations with immunotherapies.

HER2-Selective and Reversible Tyrosine Kinase Inhibitor Tucatinib Potentiates the Activity of T-DM1 in Preclinical Models of HER2-positive Breast Cancer.

The oncogenic receptor HER2 is overexpressed in many cancers, including up to 20% of breast cancers. Despite the availability of HER2-targeted treatments, patients' disease often progresses during therapy, underscoring the need for novel treatment strategies. The addition of tucatinib, a reversible, highly selective HER2 tyrosine kinase inhibitor (TKI), to treatment with trastuzumab and capecitabine significantly improved survival outcomes of patients with HER2-positive metastatic breast cancer, including those with active brain metastases. We rationalized that combining tucatinib with other HER2-targeting agents with complementary mechanisms of action would further increase efficacy against tumors. We characterized the activity of tucatinib with the antibody­drug conjugate T-DM1 in preclinical models of breast cancer, including HER2-positive breast cancer cells and patient-derived xenograft (PDX) models. Mechanistic details on tucatinib activity were obtained in internalization and catabolism studies. In combination, tucatinib and T-DM1 showed an enhanced, often synergistic, cytotoxic response and demonstrated improved antitumor activity in vivo, including in PDX models refractory to T-DM1 single-agent activity. Mechanistically, tucatinib mediated an increase in inactive HER2 molecules at the cell surface through inhibition of HER2 ubiquitination, resulting in increased internalization and catabolism of T-DM1. The combination was correlated with enhanced HER2 pathway inhibition, decreased proliferation, and increased apoptosis. In a xenograft model of brain metastasis, tucatinib penetrated intracranial tumor tissues, inhibiting tumor growth and improving survival. These results suggest that tucatinib may be the optimal TKI partner for HER2-targeted therapies and support clinical studies of its combination with T-DM1, including in patients with brain metastases. SIGNIFICANCE: The preclinical findings in breast cancer models presented here demonstrate that combining tucatinib with T-DM1 enhances the antitumor activity of either agent alone, supporting clinical studies of the combination in HER2-positive breast cancer, including in patients with brain metastases, which remains an important unmet medical need.

SGN-CD228A Is an Investigational CD228-Directed Antibody-Drug Conjugate with Potent Antitumor Activity across a Wide Spectrum of Preclinical Solid Tumor Models.

SGN-CD228A is an investigational antibody-drug conjugate (ADC) directed to melanotransferrin (CD228, MELTF, MFI2, p97), a cell-surface protein first identified in melanoma. SGN-CD228A consists of a humanized antibody, hL49, with high specificity and affinity for CD228 that is stably conjugated to 8 molecules of the clinically validated microtubule-disrupting agent monomethyl auristatin E (MMAE) via a novel glucuronide linker. We performed comprehensive IHC studies, which corroborated published RNA sequencing data and confirmed low CD228 expression in normal tissues and high expression in several cancers, including melanoma, squamous non-small cell lung cancer (NSCLC), triple-negative breast cancer (TNBC), colorectal cancer, and pancreatic cancer. SGN-CD228A was efficiently internalized in various tumor cell types, and its cytotoxic activity was dependent on CD228 expression and internalization and intrinsic sensitivity to the MMAE payload. Compared with the valine-citrulline dipeptide linker, the novel glucuronide linker increased the cellular retention of MMAE in vitro and conferred improved antitumor activity against melanoma cell lines in vitro and in vivo. In addition, SGN-CD228A was active across melanoma, TNBC, and NSCLC cell line- and patient-derived xenograft models with heterogeneous antigen expression. In vivo, CD228 expression was important for response to SGN-CD228A but was not well correlated across all tumor types, suggesting that other factors associated with ADC activity are important. Overall, SGN-CD228A is a CD228-directed, investigational ADC that employs innovative technology and has compelling preclinical antitumor activity. SGN-CD228A is investigated in a Phase I clinical trial (NCT04042480) in patients with advanced solid tumors.

Mapping and Dynamics of Regulatory DNA in Maturing Arabidopsis thaliana Siliques.

The genome is reprogrammed during development to produce diverse cell types, largely through altered expression and activity of key transcription factors. The accessibility and critical functions of epidermal cells have made them a model for connecting transcriptional events to development in a range of model systems. In Arabidopsis thaliana and many other plants, fertilization triggers differentiation of specialized epidermal seed coat cells that have a unique morphology caused by large extracellular deposits of polysaccharides. Here, we used DNase I-seq to generate regulatory landscapes of A. thaliana seeds at two critical time points in seed coat maturation (4 and 7 DPA), enriching for seed coat cells with the INTACT method. We found over 3,000 developmentally dynamic regulatory DNA elements and explored their relationship with nearby gene expression. The dynamic regulatory elements were enriched for motifs for several transcription factors families; most notably the TCP family at the earlier time point and the MYB family at the later one. To assess the extent to which the observed regulatory sites in seeds added to previously known regulatory sites in A. thaliana, we compared our data to 11 other data sets generated with 7-day-old seedlings for diverse tissues and conditions. Surprisingly, over a quarter of the regulatory, i.e. accessible, bases observed in seeds were novel. Notably, plant regulatory landscapes from different tissues, cell types, or developmental stages were more dynamic than those generated from bulk tissue in response to environmental perturbations, highlighting the importance of extending studies of regulatory DNA to single tissues and cell types during development.

DNase I hypersensitivity analysis of the mouse brain and retina identifies region-specific regulatory elements.

BACKGROUND: The brain, spinal cord, and neural retina comprise the central nervous system (CNS) of vertebrates. Understanding the regulatory mechanisms that underlie the enormous cell-type diversity of the CNS is a significant challenge. Whole-genome mapping of DNase I-hypersensitive sites (DHSs) has been used to identify cis-regulatory elements in many tissues. We have applied this approach to the mouse CNS, including developing and mature neural retina, whole brain, and two well-characterized brain regions, the cerebellum and the cerebral cortex. RESULTS: For the various regions and developmental stages of the CNS that we analyzed, there were approximately the same number of DHSs; however, there were many DHSs unique to each CNS region and developmental stage. Many of the DHSs are likely to mark enhancers that are specific to the specific CNS region and developmental stage. We validated the DNase I mapping approach for identification of CNS enhancers using the existing VISTA Browser database and with in vivo and in vitro electroporation of the retina. Analysis of transcription factor consensus sites within the DHSs shows distinct region-specific profiles of transcriptional regulators particular to each region. Clustering developmentally dynamic DHSs in the retina revealed enrichment of developmental stage-specific transcriptional regulators. Additionally, we found reporter gene activity in the retina driven from several previously uncharacterized regulatory elements surrounding the neurodevelopmental gene Otx2. Identification of DHSs shared between mouse and human showed region-specific differences in the evolution of cis-regulatory elements. CONCLUSIONS: Overall, our results demonstrate the potential of genome-wide DNase I mapping to cis-regulatory questions regarding the regional diversity within the CNS. These data represent an extensive catalogue of potential cis-regulatory elements within the CNS that display region and temporal specificity, as well as a set of DHSs common to CNS tissues. Further examination of evolutionary conservation of DHSs between CNS regions and different species may reveal important cis-regulatory elements in the evolution of the mammalian CNS.

Cell-of-origin chromatin organization shapes the mutational landscape of cancer.

Cancer is a disease potentiated by mutations in somatic cells. Cancer mutations are not distributed uniformly along the human genome. Instead, different human genomic regions vary by up to fivefold in the local density of cancer somatic mutations, posing a fundamental problem for statistical methods used in cancer genomics. Epigenomic organization has been proposed as a major determinant of the cancer mutational landscape. However, both somatic mutagenesis and epigenomic features are highly cell-type-specific. We investigated the distribution of mutations in multiple independent samples of diverse cancer types and compared them to cell-type-specific epigenomic features. Here we show that chromatin accessibility and modification, together with replication timing, explain up to 86% of the variance in mutation rates along cancer genomes. The best predictors of local somatic mutation density are epigenomic features derived from the most likely cell type of origin of the corresponding malignancy. Moreover, we find that cell-of-origin chromatin features are much stronger determinants of cancer mutation profiles than chromatin features of matched cancer cell lines. Furthermore, we show that the cell type of origin of a cancer can be accurately determined based on the distribution of mutations along its genome. Thus, the DNA sequence of a cancer genome encompasses a wealth of information about the identity and epigenomic features of its cell of origin.

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.

A comparative encyclopedia of DNA elements in the mouse genome.

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

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.

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.

Mapping and dynamics of regulatory DNA and transcription factor networks in A. thaliana.

Our understanding of gene regulation in plants is constrained by our limited knowledge of plant cis-regulatory DNA and its dynamics. We mapped DNase I hypersensitive sites (DHSs) in A. thaliana seedlings and used genomic footprinting to delineate ∼ 700,000 sites of in vivo transcription factor (TF) occupancy at nucleotide resolution. We show that variation associated with 72 diverse quantitative phenotypes localizes within DHSs. TF footprints encode an extensive cis-regulatory lexicon subject to recent evolutionary pressures, and widespread TF binding within exons may have shaped codon usage patterns. The architecture of A. thaliana TF regulatory networks is strikingly similar to that of animals in spite of diverged regulatory repertoires. We analyzed regulatory landscape dynamics during heat shock and photomorphogenesis, disclosing thousands of environmentally sensitive elements and enabling mapping of key TF regulatory circuits underlying these fundamental responses. Our results provide an extensive resource for the study of A. thaliana gene regulation and functional biology.

Domains of genome-wide gene expression dysregulation in Down's syndrome.

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

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 components of the nucleotide excision repair system. Together, our results connect chromatin structure, gene regulation and cancer-associated somatic mutation.

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.

Genome-scale mapping of DNase I hypersensitivity.

DNase I-seq is a global and high-resolution method that uses the nonspecific endonuclease DNase I to map chromatin accessibility. These accessible regions, designated as DNase I hypersensitive sites (DHSs), define the regulatory features, (e.g., promoters, enhancers, insulators, and locus control regions) of complex genomes. In this unit, methods are described for nuclei isolation, digestion of nuclei with limiting concentrations of DNase I, and the biochemical fractionation of DNase I hypersensitive sites in preparation for high-throughput sequencing. DNase I-seq is an unbiased and robust method that is not predicated on an a priori understanding of regulatory patterns or chromatin features.

Patterns of regulatory activity across diverse human cell types predict tissue identity, transcription factor binding, and long-range interactions.

Regulatory elements recruit transcription factors that modulate gene expression distinctly across cell types, but the relationships among these remains elusive. To address this, we analyzed matched DNase-seq and gene expression data for 112 human samples representing 72 cell types. We first defined more than 1800 clusters of DNase I hypersensitive sites (DHSs) with similar tissue specificity of DNase-seq signal patterns. We then used these to uncover distinct associations between DHSs and promoters, CpG islands, conserved elements, and transcription factor motif enrichment. Motif analysis within clusters identified known and novel motifs in cell-type-specific and ubiquitous regulatory elements and supports a role for AP-1 regulating open chromatin. We developed a classifier that accurately predicts cell-type lineage based on only 43 DHSs and evaluated the tissue of origin for cancer cell types. A similar classifier identified three sex-specific loci on the X chromosome, including the XIST lincRNA locus. By correlating DNase I signal and gene expression, we predicted regulated genes for more than 500K DHSs. Finally, we introduce a web resource to enable researchers to use these results to explore these regulatory patterns and better understand how expression is modulated within and across human cell types.

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.