Polycomb-mediated genome architecture enables long-range spreading of H3K27 methylation.

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding­deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

ecDNA hubs drive cooperative intermolecular oncogene expression.

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs-clusters of around 10-100 ecDNAs within the nucleus-enable intermolecular enhancer-gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer-gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy.

LKB1 inactivation modulates chromatin accessibility to drive metastatic progression.

Metastasis is the leading cause of cancer-related deaths and enables cancer cells to compromise organ function by expanding in secondary sites. Since primary tumours and metastases often share the same constellation of driver mutations, the mechanisms that drive their distinct phenotypes are unclear. Here we show that inactivation of the frequently mutated tumour suppressor gene LKB1 (encoding liver kinase B1) has evolving effects throughout the progression of lung cancer, which leads to the differential epigenetic re-programming of early-stage primary tumours compared with late-stage metastases. By integrating genome-scale CRISPR-Cas9 screening with bulk and single-cell multi-omic analyses, we unexpectedly identify LKB1 as a master regulator of chromatin accessibility in lung adenocarcinoma primary tumours. Using an in vivo model of metastatic progression, we further show that loss of LKB1 activates the early endoderm transcription factor SOX17 in metastases and a metastatic-like sub-population of cancer cells within primary tumours. The expression of SOX17 is necessary and sufficient to drive a second wave of epigenetic changes in LKB1-deficient cells that enhances metastatic ability. Overall, our study demonstrates how the downstream effects of an individual driver mutation can change throughout cancer development, with implications for stage-specific therapeutic resistance mechanisms and the gene regulatory underpinnings of metastatic evolution.

High-throughput single-cell chromatin accessibility CRISPR screens enable unbiased identification of regulatory networks in cancer.

Chromatin accessibility profiling can identify putative regulatory regions genome wide; however, pooled single-cell methods for assessing the effects of regulatory perturbations on accessibility are limited. Here, we report a modified droplet-based single-cell ATAC-seq protocol for perturbing and evaluating dynamic single-cell epigenetic states. This method (Spear-ATAC) enables simultaneous read-out of chromatin accessibility profiles and integrated sgRNA spacer sequences from thousands of individual cells at once. Spear-ATAC profiling of 104,592 cells representing 414 sgRNA knock-down populations reveals the temporal dynamics of epigenetic responses to regulatory perturbations in cancer cells and the associations between transcription factor binding profiles.

ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis.

The advent of single-cell chromatin accessibility profiling has accelerated the ability to map gene regulatory landscapes but has outpaced the development of scalable software to rapidly extract biological meaning from these data. Here we present a software suite for single-cell analysis of regulatory chromatin in R (ArchR; https://www.archrproject.com/ ) that enables fast and comprehensive analysis of single-cell chromatin accessibility data. ArchR provides an intuitive, user-focused interface for complex single-cell analyses, including doublet removal, single-cell clustering and cell type identification, unified peak set generation, cellular trajectory identification, DNA element-to-gene linkage, transcription factor footprinting, mRNA expression level prediction from chromatin accessibility and multi-omic integration with single-cell RNA sequencing (scRNA-seq). Enabling the analysis of over 1.2 million single cells within 8 h on a standard Unix laptop, ArchR is a comprehensive software suite for end-to-end analysis of single-cell chromatin accessibility that will accelerate the understanding of gene regulation at the resolution of individual cells.

Chromatin accessibility landscapes of skin cells in systemic sclerosis nominate dendritic cells in disease pathogenesis.

Systemic sclerosis (SSc) is a disease at the intersection of autoimmunity and fibrosis. However, the epigenetic regulation and the contributions of diverse cell types to SSc remain unclear. Here we survey, using ATAC-seq, the active DNA regulatory elements of eight types of primary cells in normal skin from healthy controls, as well as clinically affected and unaffected skin from SSc patients. We find that accessible DNA elements in skin-resident dendritic cells (DCs) exhibit the highest enrichment of SSc-associated single-nucleotide polymorphisms (SNPs) and predict the degrees of skin fibrosis in patients. DCs also have the greatest disease-associated changes in chromatin accessibility and the strongest alteration of cell-cell interactions in SSc lesions. Lastly, data from an independent cohort of patients with SSc confirm a significant increase of DCs in lesioned skin. Thus, the DCs epigenome links inherited susceptibility and clinically apparent fibrosis in SSc skin, and can be an important driver of SSc pathogenesis.

Single-cell epigenomic analyses implicate candidate causal variants at inherited risk loci for Alzheimer's and Parkinson's diseases.

Genome-wide association studies of neurological diseases have identified thousands of variants associated with disease phenotypes. However, most of these variants do not alter coding sequences, making it difficult to assign their function. Here, we present a multi-omic epigenetic atlas of the adult human brain through profiling of single-cell chromatin accessibility landscapes and three-dimensional chromatin interactions of diverse adult brain regions across a cohort of cognitively healthy individuals. We developed a machine-learning classifier to integrate this multi-omic framework and predict dozens of functional SNPs for Alzheimer's and Parkinson's diseases, nominating target genes and cell types for previously orphaned loci from genome-wide association studies. Moreover, we dissected the complex inverted haplotype of the MAPT (encoding tau) Parkinson's disease risk locus, identifying putative ectopic regulatory interactions in neurons that may mediate this disease association. This work expands understanding of inherited variation and provides a roadmap for the epigenomic dissection of causal regulatory variation in disease.

Chromatin Landscape Underpinning Human Dendritic Cell Heterogeneity.

Human dendritic cells (DCs) comprise subsets with distinct phenotypic and functional characteristics, but the transcriptional programs that dictate their identity remain elusive. Here, we analyze global chromatin accessibility profiles across resting and stimulated human DC subsets by means of the assay for transposase-accessible chromatin using sequencing (ATAC-seq). We uncover specific regions of chromatin accessibility for each subset and transcriptional regulators of DC function. By comparing plasmacytoid DC responses to IFN-I-producing and non-IFN-I-producing conditions, we identify genetic programs related to their function. Finally, by intersecting chromatin accessibility with genome-wide association studies, we recognize DC subset-specific enrichment of heritability in autoimmune diseases. Our results unravel the basis of human DC subset heterogeneity and provide a framework for their analysis in disease pathogenesis.

c-Jun overexpression in CAR T cells induces exhaustion resistance.

Chimeric antigen receptor (CAR) T cells mediate anti-tumour effects in a small subset of patients with cancer1-3, but dysfunction due to T cell exhaustion is an important barrier to progress4-6. To investigate the biology of exhaustion in human T cells expressing CAR receptors, we used a model system with a tonically signaling CAR, which induces hallmark features of exhaustion6. Exhaustion was associated with a profound defect in the production of IL-2, along with increased chromatin accessibility of AP-1 transcription factor motifs and overexpression of the bZIP and IRF transcription factors that have been implicated in mediating dysfunction in exhausted T cells7-10. Here we show that CAR T cells engineered to overexpress the canonical AP-1 factor c-Jun have enhanced expansion potential, increased functional capacity, diminished terminal differentiation and improved anti-tumour potency in five different mouse tumour models in vivo. We conclude that a functional deficiency in c-Jun mediates dysfunction in exhausted human T cells, and that engineering CAR T cells to overexpress c-Jun renders them resistant to exhaustion, thereby addressing a major barrier to progress for this emerging class of therapeutic agents.

Single-cell multiomic analysis identifies regulatory programs in mixed-phenotype acute leukemia.

Identifying the causes of human diseases requires deconvolution of abnormal molecular phenotypes spanning DNA accessibility, gene expression and protein abundance1-3. We present a single-cell framework that integrates highly multiplexed protein quantification, transcriptome profiling and analysis of chromatin accessibility. Using this approach, we establish a normal epigenetic baseline for healthy blood development, which we then use to deconvolve aberrant molecular features within blood from patients with mixed-phenotype acute leukemia4,5. Despite widespread epigenetic heterogeneity within the patient cohort, we observe common malignant signatures across patients as well as patient-specific regulatory features that are shared across phenotypic compartments of individual patients. Integrative analysis of transcriptomic and chromatin-accessibility maps identified 91,601 putative peak-to-gene linkages and transcription factors that regulate leukemia-specific genes, such as RUNX1-linked regulatory elements proximal to the marker gene CD69. These results demonstrate how integrative, multiomic analysis of single cells within the framework of normal development can reveal both distinct and shared molecular mechanisms of disease from patient samples.

GWAS for systemic sclerosis identifies multiple risk loci and highlights fibrotic and vasculopathy pathways.

Systemic sclerosis (SSc) is an autoimmune disease that shows one of the highest mortality rates among rheumatic diseases. We perform a large genome-wide association study (GWAS), and meta-analysis with previous GWASs, in 26,679 individuals and identify 27 independent genome-wide associated signals, including 13 new risk loci. The novel associations nearly double the number of genome-wide hits reported for SSc thus far. We define 95% credible sets of less than 5 likely causal variants in 12 loci. Additionally, we identify specific SSc subtype-associated signals. Functional analysis of high-priority variants shows the potential function of SSc signals, with the identification of 43 robust target genes through HiChIP. Our results point towards molecular pathways potentially involved in vasculopathy and fibrosis, two main hallmarks in SSc, and highlight the spectrum of critical cell types for the disease. This work supports a better understanding of the genetic basis of SSc and provides directions for future functional experiments.

Circular ecDNA promotes accessible chromatin and high oncogene expression.

Oncogenes are commonly amplified on particles of extrachromosomal DNA (ecDNA) in cancer1,2, but our understanding of the structure of ecDNA and its effect on gene regulation is limited. Here, by integrating ultrastructural imaging, long-range optical mapping and computational analysis of whole-genome sequencing, we demonstrate the structure of circular ecDNA. Pan-cancer analyses reveal that oncogenes encoded on ecDNA are among the most highly expressed genes in the transcriptome of the tumours, linking increased copy number with high transcription levels. Quantitative assessment of the chromatin state reveals that although ecDNA is packaged into chromatin with intact domain structure, it lacks higher-order compaction that is typical of chromosomes and displays significantly enhanced chromatin accessibility. Furthermore, ecDNA is shown to have a significantly greater number of ultra-long-range interactions with active chromatin, which provides insight into how the structure of circular ecDNA affects oncogene function, and connects ecDNA biology with modern cancer genomics and epigenetics.

Cryptic activation of an Irf8 enhancer governs cDC1 fate specification.

Induction of the transcription factor Irf8 in the common dendritic cell progenitor (CDP) is required for classical type 1 dendritic cell (cDC1) fate specification, but the mechanisms controlling this induction are unclear. In the present study Irf8 enhancers were identified via chromatin profiling of dendritic cells and CRISPR/Cas9 genome editing was used to assess their roles in Irf8 regulation. An enhancer 32 kilobases (kb) downstream of the Irf8 transcriptional start site (+32-kb Irf8) that was active in mature cDC1s was required for the development of this lineage, but not for its specification. Instead, a +41-kb Irf8 enhancer, previously thought to be active only in plasmacytoid dendritic cells, was found to also be transiently accessible in cDC1 progenitors, and deleting this enhancer prevented the induction of Irf8 in CDPs and abolished cDC1 specification. Thus, cryptic activation of the +41-kb Irf8 enhancer in dendritic cell progenitors is responsible for cDC1 fate specification.

An Nfil3-Zeb2-Id2 pathway imposes Irf8 enhancer switching during cDC1 development.

Classical type 1 dendritic cells (cDC1s) are required for antiviral and antitumor immunity, which necessitates an understanding of their development. Development of the cDC1 progenitor requires an E-protein-dependent enhancer located 41 kilobases downstream of the transcription start site of the transcription factor Irf8 (+41-kb Irf8 enhancer), but its maturation instead requires the Batf3-dependent +32-kb Irf8 enhancer. To understand this switch, we performed single-cell RNA sequencing of the common dendritic cell progenitor (CDP) and identified a cluster of cells that expressed transcription factors that influence cDC1 development, such as Nfil3, Id2 and Zeb2. Genetic epistasis among these factors revealed that Nfil3 expression is required for the transition from Zeb2hi and Id2lo CDPs to Zeb2lo and Id2hi CDPs, which represent the earliest committed cDC1 progenitors. This genetic circuit blocks E-protein activity to exclude plasmacytoid dendritic cell potential and explains the switch in Irf8 enhancer usage during cDC1 development.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Clonal replacement of tumor-specific T cells following PD-1 blockade.

Immunotherapies that block inhibitory checkpoint receptors on T cells have transformed the clinical care of patients with cancer1. However, whether the T cell response to checkpoint blockade relies on reinvigoration of pre-existing tumor-infiltrating lymphocytes or on recruitment of novel T cells remains unclear2-4. Here we performed paired single-cell RNA and T cell receptor sequencing on 79,046 cells from site-matched tumors from patients with basal or squamous cell carcinoma before and after anti-PD-1 therapy. Tracking T cell receptor clones and transcriptional phenotypes revealed coupling of tumor recognition, clonal expansion and T cell dysfunction marked by clonal expansion of CD8+CD39+ T cells, which co-expressed markers of chronic T cell activation and exhaustion. However, the expansion of T cell clones did not derive from pre-existing tumor-infiltrating T lymphocytes; instead, the expanded clones consisted of novel clonotypes that had not previously been observed in the same tumor. Clonal replacement of T cells was preferentially observed in exhausted CD8+ T cells and evident in patients with basal or squamous cell carcinoma. These results demonstrate that pre-existing tumor-specific T cells may have limited reinvigoration capacity, and that the T cell response to checkpoint blockade derives from a distinct repertoire of T cell clones that may have just recently entered the tumor.

HiChIRP reveals RNA-associated chromosome conformation.

Modular domains of long non-coding RNAs can serve as scaffolds to bring distant regions of the linear genome into spatial proximity. Here, we present HiChIRP, a method leveraging bio-orthogonal chemistry and optimized chromosome conformation capture conditions, which enables interrogation of chromatin architecture focused around a specific RNA of interest down to approximately ten copies per cell. HiChIRP of three nuclear RNAs reveals insights into promoter interactions (7SK), telomere biology (telomerase RNA component) and inflammatory gene regulation (lincRNA-EPS).

Enhancer Connectome Nominates Target Genes of Inherited Risk Variants from Inflammatory Skin Disorders.

The vast majority of polymorphisms for human dermatologic diseases fall in noncoding DNA regions, leading to difficulty interpreting their functional significance. Recent work using chromosome conformation capture technology in combination with chromatin immunoprecipitation (ChIP) has provided a systematic means of linking noncoding variants in active enhancer loci to putative gene targets. Here, we apply H3K27ac HiChIP high-resolution contact maps, generated from primary human T-cell subsets (CD4+ naïve, T helper type 17, and regulatory T cells), to 21 dermatologic conditions associated with single nucleotide polymorphisms from 106 genome-wide association studies. This "enhancer connectome" identified 1,492 HiChIP gene targets from 542 noncoding SNPs (P ≤ 5.0 × 10-8). SNP-containing enhancers from inflammatory skin conditions were significantly enriched at the HLA locus and also targeted several key factors from the JAK-STAT signaling pathway, but nonimmune conditions were not. A focused profiling of systemic lupus erythematosus HiChIP genes identified enhancer interactions with factors important for effector CD4+ T-cell differentiation and function, including IRF8 and members of the Ikaros family of zinc-finger proteins. Our results show the ability of the enhancer connectome to nominate functionally relevant candidates from genome-wide association study-identified variants, representing a powerful tool to guide future studies into the genomic regulatory mechanisms underlying dermatologic diseases.