Search

Extrachromosomal DNA (ecDNA) is a major contributor to treatment resistance and poor outcome for patients with cancer 1,2 . Here we examine the diversity of ecDNA elements across cancer, revealing the associated tissue, genetic and mutational contexts. By analysing data from 14,778 patients with 39 tumour types from the 100,000 Genomes Project, we demonstrate that 17.1% of tumour samples contain ecDNA. We reveal a pattern highly indicative of tissue-context-based selection for ecDNAs, linking their genomic content to their tissue of origin. We show that not only is ecDNA a mechanism for amplification of driver oncogenes, but it also a mechanism that frequently amplifies immunomodulatory and inflammatory genes, such as those that modulate lymphocyte-mediated immunity and immune effector processes. Moreover, ecDNAs carrying immunomodulatory genes are associated with reduced tumour T cell infiltration. We identify ecDNAs bearing only enhancers, promoters and lncRNA elements, suggesting the combinatorial power of interactions between ecDNAs in trans. We also identify intrinsic and environmental mutational processes linked to ecDNA, including those linked to its formation, such as tobacco exposure, and progression, such as homologous recombination repair deficiency. Clinically, ecDNA detection was associated with tumour stage, more prevalent after targeted therapy and cytotoxic treatments, and associated with metastases and shorter overall survival. These results shed light on why ecDNA is a substantial clinical problem that can cooperatively drive tumour growth signals, alter transcriptional landscapes and suppress the immune system., (© 2024. The Author(s).)

Chronic inflammation and tissue fibrosis are common responses that worsen organ function, yet the molecular mechanisms governing their cross-talk are poorly understood. In diseased organs, stress-induced gene expression changes fuel maladaptive cell state transitions 1 and pathological interaction between cellular compartments. Although chronic fibroblast activation worsens dysfunction in the lungs, liver, kidneys and heart, and exacerbates many cancers 2 , the stress-sensing mechanisms initiating transcriptional activation of fibroblasts are poorly understood. Here we show that conditional deletion of the transcriptional co-activator Brd4 in infiltrating Cx3cr1 + macrophages ameliorates heart failure in mice and significantly reduces fibroblast activation. Analysis of single-cell chromatin accessibility and BRD4 occupancy in vivo in Cx3cr1 + cells identified a large enhancer proximal to interleukin-1β (IL-1β, encoded by Il1b), and a series of CRISPR-based deletions revealed the precise stress-dependent regulatory element that controls Il1b expression. Secreted IL-1β activated a fibroblast RELA-dependent (also known as p65) enhancer near the transcription factor MEOX1, resulting in a profibrotic response in human cardiac fibroblasts. In vivo, antibody-mediated IL-1β neutralization improved cardiac function and tissue fibrosis in heart failure. Systemic IL-1β inhibition or targeted Il1b deletion in Cx3cr1 + cells prevented stress-induced Meox1 expression and fibroblast activation. The elucidation of BRD4-dependent cross-talk between a specific immune cell subset and fibroblasts through IL-1β reveals how inflammation drives profibrotic cell states and supports strategies that modulate this process in heart disease and other chronic inflammatory disorders featuring tissue remodelling., Competing Interests: Competing interests D.S. is scientific co-founder, shareholder and director of Tenaya Therapeutics. S.M.H. is an executive, officer and shareholder of Amgen and is a scientific co-founder and shareholder of Tenaya Therapeutics. T.A.M. received funding from Italfarmaco for an unrelated project. K.S.P. is a shareholder of Tenaya Therapeutics., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri) 1 and African 2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Measurements of gene expression or signal transduction activity are conventionally performed using methods that require either the destruction or live imaging of a biological sample within the timeframe of interest. Here we demonstrate an alternative paradigm in which such biological activities are stably recorded to the genome. Enhancer-driven genomic recording of transcriptional activity in multiplex (ENGRAM) is based on the signal-dependent production of prime editing guide RNAs that mediate the insertion of signal-specific barcodes (symbols) into a genomically encoded recording unit. We show how this strategy can be used for multiplex recording of the cell-type-specific activities of dozens to hundreds of cis-regulatory elements with high fidelity, sensitivity and reproducibility. Leveraging signal transduction pathway-responsive cis-regulatory elements, we also demonstrate time- and concentration-dependent genomic recording of WNT, NF-κB and Tet-On activities. By coupling ENGRAM to sequential genome editing via DNA Typewriter 1 , we stably record information about the temporal dynamics of two orthogonal signalling pathways to genomic DNA. Finally we apply ENGRAM to integratively record the transient activity of nearly 100 transcription factor consensus motifs across daily windows spanning the differentiation of mouse embryonic stem cells into gastruloids, an in vitro model of early mammalian development. Although these are proof-of-concept experiments and much work remains to fully realize the possibilities, the symbolic recording of biological signals or states within cells, to the genome and over time, has broad potential to complement contemporary paradigms for how we make measurements in biological systems., (© 2024. The Author(s).)

T-lineage acute lymphoblastic leukaemia (T-ALL) is a high-risk tumour 1 that has eluded comprehensive genomic characterization, which is partly due to the high frequency of noncoding genomic alterations that result in oncogene deregulation 2,3 . Here we report an integrated analysis of genome and transcriptome sequencing of tumour and remission samples from more than 1,300 uniformly treated children with T-ALL, coupled with epigenomic and single-cell analyses of malignant and normal T cell precursors. This approach identified 15 subtypes with distinct genomic drivers, gene expression patterns, developmental states and outcomes. Analyses of chromatin topology revealed multiple mechanisms of enhancer deregulation that involve enhancers and genes in a subtype-specific manner, thereby demonstrating widespread involvement of the noncoding genome. We show that the immunophenotypically described, high-risk entity of early T cell precursor ALL is superseded by a broader category of 'early T cell precursor-like' leukaemia. This category has a variable immunophenotype and diverse genomic alterations of a core set of genes that encode regulators of hematopoietic stem cell development. Using multivariable outcome models, we show that genetic subtypes, driver and concomitant genetic alterations independently predict treatment failure and survival. These findings provide a roadmap for the classification, risk stratification and mechanistic understanding of this disease., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers 1,2 . Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations 3 . However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression., (© 2024. The Author(s).)

The current technologies to place new DNA into specific locations in plant genomes are low frequency and error-prone, and this inefficiency hampers genome-editing approaches to develop improved crops 1,2 . Often considered to be genome 'parasites', transposable elements (TEs) evolved to insert their DNA seamlessly into genomes 3-5 . Eukaryotic TEs select their site of insertion based on preferences for chromatin contexts, which differ for each TE type 6-9 . Here we developed a genome engineering tool that controls the TE insertion site and cargo delivered, taking advantage of the natural ability of the TE to precisely excise and insert into the genome. Inspired by CRISPR-associated transposases that target transposition in a programmable manner in bacteria 10-12 , we fused the rice Pong transposase protein to the Cas9 or Cas12a programmable nucleases. We demonstrated sequence-specific targeted insertion (guided by the CRISPR gRNA) of enhancer elements, an open reading frame and a gene expression cassette into the genome of the model plant Arabidopsis. We then translated this system into soybean-a major global crop in need of targeted insertion technology. We have engineered a TE 'parasite' into a usable and accessible toolkit that enables the sequence-specific targeting of custom DNA into plant genomes., (© 2024. The Author(s).)

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease 1-3 . The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb 4 . Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly 4-6 . However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies., (© 2024. The Author(s).)

Enhancers control gene expression and have crucial roles in development and homeostasis 1-3 . However, the targeted de novo design of enhancers with tissue-specific activities has remained challenging. Here we combine deep learning and transfer learning to design tissue-specific enhancers for five tissues in the Drosophila melanogaster embryo: the central nervous system, epidermis, gut, muscle and brain. We first train convolutional neural networks using genome-wide single-cell assay for transposase-accessible chromatin with sequencing (ATAC-seq) datasets and then fine-tune the convolutional neural networks with smaller-scale data from in vivo enhancer activity assays, yielding models with 13% to 76% positive predictive value according to cross-validation. We designed and experimentally assessed 40 synthetic enhancers (8 per tissue) in vivo, of which 31 (78%) were active and 27 (68%) functioned in the target tissue (100% for central nervous system and muscle). The strategy of combining genome-wide and small-scale functional datasets by transfer learning is generally applicable and should enable the design of tissue-, cell type- and cell state-specific enhancers in any system., (© 2023. The Author(s).)

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes 1 . It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models 2-6 , can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states., (© 2023. The Author(s).)

Olfactory receptor (OR) choice provides an extreme example of allelic competition for transcriptional dominance, where every olfactory neuron stably transcribes one of approximately 2,000 or more OR alleles 1,2 . OR gene choice is mediated by a multichromosomal enhancer hub that activates transcription at a single OR 3,4 , followed by OR-translation-dependent feedback that stabilizes this choice 5,6 . Here, using single-cell genomics, we show formation of many competing hubs with variable enhancer composition, only one of which retains euchromatic features and transcriptional competence. Furthermore, we provide evidence that OR transcription recruits enhancers and reinforces enhancer hub activity locally, whereas OR RNA inhibits transcription of competing ORs over distance, promoting transition to transcriptional singularity. Whereas OR transcription is sufficient to break the symmetry between equipotent enhancer hubs, OR translation stabilizes transcription at the prevailing hub, indicating that there may be sequential non-coding and coding mechanisms that are implemented by OR alleles for transcriptional prevalence. We propose that coding OR mRNAs possess non-coding functions that influence nuclear architecture, enhance their own transcription and inhibit transcription from their competitors, with generalizable implications for probabilistic cell fate decisions., (© 2023. The Author(s).)

Chromatin accessibility is essential in regulating gene expression and cellular identity, and alterations in accessibility have been implicated in driving cancer initiation, progression and metastasis 1-4 . Although the genetic contributions to oncogenic transitions have been investigated, epigenetic drivers remain less understood. Here we constructed a pan-cancer epigenetic and transcriptomic atlas using single-nucleus chromatin accessibility data (using single-nucleus assay for transposase-accessible chromatin) from 225 samples and matched single-cell or single-nucleus RNA-sequencing expression data from 206 samples. With over 1 million cells from each platform analysed through the enrichment of accessible chromatin regions, transcription factor motifs and regulons, we identified epigenetic drivers associated with cancer transitions. Some epigenetic drivers appeared in multiple cancers (for example, regulatory regions of ABCC1 and VEGFA; GATA6 and FOX-family motifs), whereas others were cancer specific (for example, regulatory regions of FGF19, ASAP2 and EN1, and the PBX3 motif). Among epigenetically altered pathways, TP53, hypoxia and TNF signalling were linked to cancer initiation, whereas oestrogen response, epithelial-mesenchymal transition and apical junction were tied to metastatic transition. Furthermore, we revealed a marked correlation between enhancer accessibility and gene expression and uncovered cooperation between epigenetic and genetic drivers. This atlas provides a foundation for further investigation of epigenetic dynamics in cancer transitions., (© 2023. The Author(s).)

Zygotic genome activation (ZGA) activates the quiescent genome to enable the maternal-to-zygotic transition 1,2 . However, the identity of transcription factors that underlie mammalian ZGA in vivo remains elusive. Here we show that OBOX, a PRD-like homeobox domain transcription factor family (OBOX1-OBOX8) 3-5 , are key regulators of mouse ZGA. Mice deficient for maternally transcribed Obox1/2/5/7 and zygotically expressed Obox3/4 had a two-cell to four-cell arrest, accompanied by impaired ZGA. The Obox knockout defects could be rescued by restoring either maternal and zygotic OBOX, which suggests that maternal and zygotic OBOX redundantly support embryonic development. Chromatin-binding analysis showed that Obox knockout preferentially affected OBOX-binding targets. Mechanistically, OBOX facilitated the 'preconfiguration' of RNA polymerase II, as the polymerase relocated from the initial one-cell binding targets to ZGA gene promoters and distal enhancers. Impaired polymerase II preconfiguration in Obox mutants was accompanied by defective ZGA and chromatin accessibility transition, as well as aberrant activation of one-cell polymerase II targets. Finally, ectopic expression of OBOX activated ZGA genes and MERVL repeats in mouse embryonic stem cells. These data thus demonstrate that OBOX regulates mouse ZGA and early embryogenesis., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters 1-4 . However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs 5,6 . Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

For cells to initiate and sustain a differentiated state, it is necessary that a 'memory' of this state is transmitted through mitosis to the daughter cells 1-3 . Mammalian switch/sucrose non-fermentable (SWI/SNF) complexes (also known as Brg1/Brg-associated factors, or BAF) control cell identity by modulating chromatin architecture to regulate gene expression 4-7 , but whether they participate in cell fate memory is unclear. Here we provide evidence that subunits of SWI/SNF act as mitotic bookmarks to safeguard cell identity during cell division. The SWI/SNF core subunits SMARCE1 and SMARCB1 are displaced from enhancers but are bound to promoters during mitosis, and we show that this binding is required for appropriate reactivation of bound genes after mitotic exit. Ablation of SMARCE1 during a single mitosis in mouse embryonic stem cells is sufficient to disrupt gene expression, impair the occupancy of several established bookmarks at a subset of their targets and cause aberrant neural differentiation. Thus, SWI/SNF subunit SMARCE1 has a mitotic bookmarking role and is essential for heritable epigenetic fidelity during transcriptional reprogramming., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Transcriptional enhancer elements are responsible for orchestrating the temporal and spatial control over gene expression that is crucial for programming cell identity during development 1-3 . Here we describe a novel enhancer element that is important for regulating the expression of Prox1 in lymphatic endothelial cells. This evolutionarily conserved enhancer is bound by key lymphatic transcriptional regulators including GATA2, FOXC2, NFATC1 and PROX1. Genome editing of the enhancer to remove five nucleotides encompassing the GATA2-binding site resulted in perinatal death of homozygous mutant mice due to profound lymphatic vascular defects. Lymphatic endothelial cells in enhancer mutant mice exhibited reduced expression of genes characteristic of lymphatic endothelial cell identity and increased expression of genes characteristic of haemogenic endothelium, and acquired the capacity to generate haematopoietic cells. These data not only reveal a transcriptional enhancer element important for regulating Prox1 expression and lymphatic endothelial cell identity but also demonstrate that the lymphatic endothelium has haemogenic capacity, ordinarily repressed by Prox1., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences 1,2 . Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer 3 . However, only a small minority of SVs are associated with altered gene expression 4,5 , and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Diffuse large B cell lymphoma (DLBCL) is the most common B cell non-Hodgkin lymphoma and remains incurable in around 40% of patients. Efforts to sequence the coding genome identified several genes and pathways that are altered in this disease, including potential therapeutic targets 1-5 . However, the non-coding genome of DLBCL remains largely unexplored. Here we show that active super-enhancers are highly and specifically hypermutated in 92% of samples from individuals with DLBCL, display signatures of activation-induced cytidine deaminase activity, and are linked to genes that encode B cell developmental regulators and oncogenes. As evidence of oncogenic relevance, we show that the hypermutated super-enhancers linked to the BCL6, BCL2 and CXCR4 proto-oncogenes prevent the binding and transcriptional downregulation of the corresponding target gene by transcriptional repressors, including BLIMP1 (targeting BCL6) and the steroid receptor NR3C1 (targeting BCL2 and CXCR4). Genetic correction of selected mutations restored repressor DNA binding, downregulated target gene expression and led to the counter-selection of cells containing corrected alleles, indicating an oncogenic dependency on the super-enhancer mutations. This pervasive super-enhancer mutational mechanism reveals a major set of genetic lesions deregulating gene expression, which expands the involvement of known oncogenes in DLBCL pathogenesis and identifies new deregulated gene targets of therapeutic relevance., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Gene regulation in the human genome is controlled by distal enhancers that activate specific nearby promoters 1 . A proposed model for this specificity is that promoters have sequence-encoded preferences for certain enhancers, for example, mediated by interacting sets of transcription factors or cofactors 2 . This 'biochemical compatibility' model has been supported by observations at individual human promoters and by genome-wide measurements in Drosophila 3-9 . However, the degree to which human enhancers and promoters are intrinsically compatible has not yet been systematically measured, and how their activities combine to control RNA expression remains unclear. Here we design a high-throughput reporter assay called enhancer × promoter self-transcribing active regulatory region sequencing (ExP STARR-seq) and applied it to examine the combinatorial compatibilities of 1,000 enhancer and 1,000 promoter sequences in human K562 cells. We identify simple rules for enhancer-promoter compatibility, whereby most enhancers activate all promoters by similar amounts, and intrinsic enhancer and promoter activities multiplicatively combine to determine RNA output (R 2  = 0.82). In addition, two classes of enhancers and promoters show subtle preferential effects. Promoters of housekeeping genes contain built-in activating motifs for factors such as GABPA and YY1, which decrease the responsiveness of promoters to distal enhancers. Promoters of variably expressed genes lack these motifs and show stronger responsiveness to enhancers. Together, this systematic assessment of enhancer-promoter compatibility suggests a multiplicative model tuned by enhancer and promoter class to control gene transcription in the human genome., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The divergence of the common dendritic cell progenitor 1-3 (CDP) into the conventional type 1 and type 2 dendritic cell (cDC1 and cDC2, respectively) lineages 4,5 is poorly understood. Some transcription factors act in the commitment of already specified progenitors-such as BATF3, which stabilizes Irf8 autoactivation at the +32 kb Irf8 enhancer 4,6 -but the mechanisms controlling the initial divergence of CDPs remain unknown. Here we report the transcriptional basis of CDP divergence and describe the first requirements for pre-cDC2 specification. Genetic epistasis analysis 7 suggested that Nfil3 acts upstream of Id2, Batf3 and Zeb2 in cDC1 development but did not reveal its mechanism or targets. Analysis of newly generated NFIL3 reporter mice showed extremely transient NFIL3 expression during cDC1 specification. CUT&RUN and chromatin immunoprecipitation followed by sequencing identified endogenous NFIL3 binding in the -165 kb Zeb2 enhancer 8 at three sites that also bind the CCAAT-enhancer-binding proteins C/EBPα and C/EBPβ. In vivo mutational analysis using CRISPR-Cas9 targeting showed that these NFIL3-C/EBP sites are functionally redundant, with C/EBPs supporting and NFIL3 repressing Zeb2 expression at these sites. A triple mutation of all three NFIL3-C/EBP sites ablated Zeb2 expression in myeloid, but not lymphoid progenitors, causing the complete loss of pre-cDC2 specification and mature cDC2 development in vivo. These mice did not generate T helper 2 (T H 2) cell responses against Heligmosomoides polygyrus infection, consistent with cDC2 supporting T H 2 responses to helminths 9-11 . Thus, CDP divergence into cDC1 or cDC2 is controlled by competition between NFIL3 and C/EBPs at the -165 kb Zeb2 enhancer., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

All multicellular organisms rely on differential gene transcription regulated by genomic enhancers, which function through cofactors that are recruited by transcription factors 1,2 . Emerging evidence suggests that not all cofactors are required at all enhancers 3-5 , yet whether these observations reflect more general principles or distinct types of enhancers remained unknown. Here we categorized human enhancers by their cofactor dependencies and show that these categories provide a framework to understand the sequence and chromatin diversity of enhancers and their roles in different gene-regulatory programmes. We quantified enhancer activities along the entire human genome using STARR-seq 6 in HCT116 cells, following the rapid degradation of eight cofactors. This analysis identified different types of enhancers with distinct cofactor requirements, sequences and chromatin properties. Some enhancers were insensitive to the depletion of the core Mediator subunit MED14 or the bromodomain protein BRD4 and regulated distinct transcriptional programmes. In particular, canonical Mediator 7 seemed dispensable for P53-responsive enhancers, and MED14-depleted cells induced endogenous P53 target genes. Similarly, BRD4 was not required for the transcription of genes that bear CCAAT boxes and a TATA box (including histone genes and LTR12 retrotransposons) or for the induction of heat-shock genes. This categorization of enhancers through cofactor dependencies reveals distinct enhancer types that can bypass broadly utilized cofactors, which illustrates how alternative ways to activate transcription separate gene expression programmes and provide a conceptual framework to understand enhancer function and regulatory specificity., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The prevailing view of metazoan gene regulation is that individual genes are independently regulated by their own dedicated sets of transcriptional enhancers. Past studies have reported long-range gene-gene associations 1-3 , but their functional importance in regulating transcription remains unclear. Here we used quantitative single-cell live imaging methods to provide a demonstration of co-dependent transcriptional dynamics of genes separated by large genomic distances in living Drosophila embryos. We find extensive physical and functional associations of distant paralogous genes, including co-regulation by shared enhancers and co-transcriptional initiation over distances of nearly 250 kilobases. Regulatory interconnectivity depends on promoter-proximal tethering elements, and perturbations in these elements uncouple transcription and alter the bursting dynamics of distant genes, suggesting a role of genome topology in the formation and stability of co-transcriptional hubs. Transcriptional coupling is detected throughout the fly genome and encompasses a broad spectrum of conserved developmental processes, suggesting a general strategy for long-range integration of gene activity., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Chromosome structure in mammals is thought to regulate transcription by modulating three-dimensional interactions between enhancers and promoters, notably through CTCF-mediated loops and topologically associating domains (TADs) 1-4 . However, how chromosome interactions are actually translated into transcriptional outputs remains unclear. Here, to address this question, we use an assay to position an enhancer at large numbers of densely spaced chromosomal locations relative to a fixed promoter, and measure promoter output and interactions within a genomic region with minimal regulatory and structural complexity. A quantitative analysis of hundreds of cell lines reveals that the transcriptional effect of an enhancer depends on its contact probabilities with the promoter through a nonlinear relationship. Mathematical modelling suggests that nonlinearity might arise from transient enhancer-promoter interactions being translated into slower promoter bursting dynamics in individual cells, therefore uncoupling the temporal dynamics of interactions from those of transcription. This uncovers a potential mechanism of how distal enhancers act from large genomic distances, and of how topologically associating domain boundaries block distal enhancers. Finally, we show that enhancer strength also determines absolute transcription levels as well as the sensitivity of a promoter to CTCF-mediated transcriptional insulation. Our measurements establish general principles for the context-dependent role of chromosome structure in long-range transcriptional regulation., (© 2022. The Author(s).)

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing 1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data 3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments 4 . By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum., (© 2021. The Author(s).)

The similarities and differences between nervous systems of various species result from developmental constraints and specific adaptations 1-4 . Comparative analyses of the prefrontal cortex (PFC), a cerebral cortex region involved in higher-order cognition and complex social behaviours, have identified true and potential human-specific structural and molecular specializations 4-8 , such as an exaggerated PFC-enriched anterior-posterior dendritic spine density gradient 5 . These changes are probably mediated by divergence in spatiotemporal gene regulation 9-17 , which is particularly prominent in the midfetal human cortex 15,18-20 . Here we analysed human and macaque transcriptomic data 15,20 and identified a transient PFC-enriched and laminar-specific upregulation of cerebellin 2 (CBLN2), a neurexin (NRXN) and glutamate receptor-δ GRID/GluD-associated synaptic organizer 21-27 , during midfetal development that coincided with the initiation of synaptogenesis. Moreover, we found that species differences in level of expression and laminar distribution of CBLN2 are, at least in part, due to Hominini-specific deletions containing SOX5-binding sites within a retinoic acid-responsive CBLN2 enhancer. In situ genetic humanization of the mouse Cbln2 enhancer drives increased and ectopic laminar Cbln2 expression and promotes PFC dendritic spine formation. These findings suggest a genetic and molecular basis for the anterior-posterior cortical gradient and disproportionate increase in the Hominini PFC of dendritic spines and a developmental mechanism that may link dysfunction of the NRXN-GRID-CBLN2 complex to the pathogenesis of neuropsychiatric disorders., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

In higher eukaryotes, many genes are regulated by enhancers that are 10 4 -10 6  base pairs (bp) away from the promoter. Enhancers contain transcription-factor-binding sites (which are typically around 7-22 bp), and physical contact between the promoters and enhancers is thought to be required to modulate gene expression. Although chromatin architecture has been mapped extensively at resolutions of 1 kilobase and above; it has not been possible to define physical contacts at the scale of the proteins that determine gene expression. Here we define these interactions in detail using a chromosome conformation capture method (Micro-Capture-C) that enables the physical contacts between different classes of regulatory elements to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show that transcription factors have an important role in the maintenance of the contacts between enhancers and promoters. Our data show that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. This supports a model in which chromatin loop extrusion 1 is dependent on cohesin loading at active promoters and enhancers, which explains the formation of tissue-specific chromatin domains without changes in CTCF binding.

Defects in DNA repair frequently lead to neurodevelopmental and neurodegenerative diseases, underscoring the particular importance of DNA repair in long-lived post-mitotic neurons 1,2 . The cellular genome is subjected to a constant barrage of endogenous DNA damage, but surprisingly little is known about the identity of the lesion(s) that accumulate in neurons and whether they accrue throughout the genome or at specific loci. Here we show that post-mitotic neurons accumulate unexpectedly high levels of DNA single-strand breaks (SSBs) at specific sites within the genome. Genome-wide mapping reveals that SSBs are located within enhancers at or near CpG dinucleotides and sites of DNA demethylation. These SSBs are repaired by PARP1 and XRCC1-dependent mechanisms. Notably, deficiencies in XRCC1-dependent short-patch repair increase DNA repair synthesis at neuronal enhancers, whereas defects in long-patch repair reduce synthesis. The high levels of SSB repair in neuronal enhancers are therefore likely to be sustained by both short-patch and long-patch processes. These data provide the first evidence of site- and cell-type-specific SSB repair, revealing unexpected levels of localized and continuous DNA breakage in neurons. In addition, they suggest an explanation for the neurodegenerative phenotypes that occur in patients with defective SSB repair.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease 1 . Many of the underlying causal variants may affect enhancers 2,3 , but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types 4 . Here we apply this ABC model to create enhancer-gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions.

Annotating the molecular basis of human disease remains an unsolved challenge, as 93% of disease loci are non-coding and gene-regulatory annotations are highly incomplete 1-3 . Here we present EpiMap, a compendium comprising 10,000 epigenomic maps across 800 samples, which we used to define chromatin states, high-resolution enhancers, enhancer modules, upstream regulators and downstream target genes. We used this resource to annotate 30,000 genetic loci that were associated with 540 traits 4 , predicting trait-relevant tissues, putative causal nucleotide variants in enriched tissue enhancers and candidate tissue-specific target genes for each. We partitioned multifactorial traits into tissue-specific contributing factors with distinct functional enrichments and disease comorbidity patterns, and revealed both single-factor monotropic and multifactor pleiotropic loci. Top-scoring loci frequently had multiple predicted driver variants, converging through multiple enhancers with a common target gene, multiple genes in common tissues, or multiple genes and multiple tissues, indicating extensive pleiotropy. Our results demonstrate the importance of dense, rich, high-resolution epigenomic annotations for the investigation of complex traits.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species 1 . However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Changes in gene regulation underlie much of phenotypic evolution 1 . However, our understanding of the potential for regulatory evolution is biased, because most evidence comes from either natural variation or limited experimental perturbations 2 . Using an automated robotics pipeline, we surveyed an unbiased mutation library for a developmental enhancer in Drosophila melanogaster. We found that almost all mutations altered gene expression and that parameters of gene expression-levels, location, and state-were convolved. The widespread pleiotropic effects of most mutations may constrain the evolvability of developmental enhancers. Consistent with these observations, comparisons of diverse Drosophila larvae revealed apparent biases in the phenotypes influenced by the enhancer. Developmental enhancers may encode a higher density of regulatory information than has been appreciated previously, imposing constraints on regulatory evolution.

Transcription factors are DNA-binding proteins that have key roles in gene regulation 1,2 . Genome-wide occupancy maps of transcriptional regulators are important for understanding gene regulation and its effects on diverse biological processes 3-6 . However, only a minority of the more than 1,600 transcription factors encoded in the human genome has been assayed. Here we present, as part of the ENCODE (Encyclopedia of DNA Elements) project, data and analyses from chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) experiments using the human HepG2 cell line for 208 chromatin-associated proteins (CAPs). These comprise 171 transcription factors and 37 transcriptional cofactors and chromatin regulator proteins, and represent nearly one-quarter of CAPs expressed in HepG2 cells. The binding profiles of these CAPs form major groups associated predominantly with promoters or enhancers, or with both. We confirm and expand the current catalogue of DNA sequence motifs for transcription factors, and describe motifs that correspond to other transcription factors that are co-enriched with the primary ChIP target. For example, FOX family motifs are enriched in ChIP-seq peaks of 37 other CAPs. We show that motif content and occupancy patterns can distinguish between promoters and enhancers. This catalogue reveals high-occupancy target regions at which many CAPs associate, although each contains motifs for only a minority of the numerous associated transcription factors. These analyses provide a more complete overview of the gene regulatory networks that define this cell type, and demonstrate the usefulness of the large-scale production efforts of the ENCODE Consortium.

Genetic variations underlying susceptibility to complex autoimmune and allergic diseases are concentrated within noncoding regulatory elements termed enhancers 1 . The functions of a large majority of disease-associated enhancers are unknown, in part owing to their distance from the genes they regulate, a lack of understanding of the cell types in which they operate, and our inability to recapitulate the biology of immune diseases in vitro. Here, using shared synteny to guide loss-of-function analysis of homologues of human enhancers in mice, we show that the prominent autoimmune and allergic disease risk locus at chromosome 11q13.5 2-7 contains a distal enhancer that is functional in CD4 + regulatory T (T reg ) cells and required for T reg -mediated suppression of colitis. The enhancer recruits the transcription factors STAT5 and NF-κB to mediate signal-driven expression of Lrrc32, which encodes the protein glycoprotein A repetitions predominant (GARP). Whereas disruption of the Lrrc32 gene results in early lethality, mice lacking the enhancer are viable but lack GARP expression in Foxp3 + T reg cells, which are unable to control colitis in a cell-transfer model of the disease. In human T reg cells, the enhancer forms conformational interactions with the promoter of LRRC32 and enhancer risk variants are associated with reduced histone acetylation and GARP expression. Finally, functional fine-mapping of 11q13.5 using CRISPR-activation (CRISPRa) identifies a CRISPRa-responsive element in the vicinity of risk variant rs11236797 capable of driving GARP expression. These findings provide a mechanistic basis for association of the 11q13.5 risk locus with immune-mediated diseases and identify GARP as a potential target in their therapy.

Physical interactions between distal regulatory elements have a key role in regulating gene expression, but the extent to which these interactions vary between cell types and contribute to cell-type-specific gene expression remains unclear. Here, to address these questions as part of phase III of the Encyclopedia of DNA Elements (ENCODE), we mapped cohesin-mediated chromatin loops, using chromatin interaction analysis by paired-end tag sequencing (ChIA-PET), and analysed gene expression in 24 diverse human cell types, including core ENCODE cell lines. Twenty-eight per cent of all chromatin loops vary across cell types; these variations modestly correlate with changes in gene expression and are effective at grouping cell types according to their tissue of origin. The connectivity of genes corresponds to different functional classes, with housekeeping genes having few contacts, and dosage-sensitive genes being more connected to enhancer elements. This atlas of chromatin loops complements the diverse maps of regulatory architecture that comprise the ENCODE Encyclopedia, and will help to support emerging analyses of genome structure and function.

Cytosine DNA methylation is essential for mammalian development but understanding of its spatiotemporal distribution in the developing embryo remains limited 1,2 . Here, as part of the mouse Encyclopedia of DNA Elements (ENCODE) project, we profiled 168 methylomes from 12 mouse tissues or organs at 9 developmental stages from embryogenesis to adulthood. We identified 1,808,810 genomic regions that showed variations in CG methylation by comparing the methylomes of different tissues or organs from different developmental stages. These DNA elements predominantly lose CG methylation during fetal development, whereas the trend is reversed after birth. During late stages of fetal development, non-CG methylation accumulated within the bodies of key developmental transcription factor genes, coinciding with their transcriptional repression. Integration of genome-wide DNA methylation, histone modification and chromatin accessibility data enabled us to predict 461,141 putative developmental tissue-specific enhancers, the human orthologues of which were enriched for disease-associated genetic variants. These spatiotemporal epigenome maps provide a resource for studies of gene regulation during tissue or organ progression, and a starting point for investigating regulatory elements that are involved in human developmental disorders.

The Encyclopedia of DNA Elements (ENCODE) project has established a genomic resource for mammalian development, profiling a diverse panel of mouse tissues at 8 developmental stages from 10.5 days after conception until birth, including transcriptomes, methylomes and chromatin states. Here we systematically examined the state and accessibility of chromatin in the developing mouse fetus. In total we performed 1,128 chromatin immunoprecipitation with sequencing (ChIP-seq) assays for histone modifications and 132 assay for transposase-accessible chromatin using sequencing (ATAC-seq) assays for chromatin accessibility across 72 distinct tissue-stages. We used integrative analysis to develop a unified set of chromatin state annotations, infer the identities of dynamic enhancers and key transcriptional regulators, and characterize the relationship between chromatin state and accessibility during developmental gene regulation. We also leveraged these data to link enhancers to putative target genes and demonstrate tissue-specific enrichments of sequence variants associated with disease in humans. The mouse ENCODE data sets provide a compendium of resources for biomedical researchers and achieve, to our knowledge, the most comprehensive view of chromatin dynamics during mammalian fetal development to date.

Highly structured RNA molecules usually interact with each other, and associate with various RNA-binding proteins, to regulate critical biological processes. However, RNA structures and interactions in intact cells remain largely unknown. Here, by coupling proximity ligation mediated by RNA-binding proteins with deep sequencing, we report an RNA in situ conformation sequencing (RIC-seq) technology for the global profiling of intra- and intermolecular RNA-RNA interactions. This technique not only recapitulates known RNA secondary structures and tertiary interactions, but also facilitates the generation of three-dimensional (3D) interaction maps of RNA in human cells. Using these maps, we identify noncoding RNA targets globally, and discern RNA topological domains and trans-interacting hubs. We reveal that the functional connectivity of enhancers and promoters can be assigned using their pairwise-interacting RNAs. Furthermore, we show that CCAT1-5L-a super-enhancer hub RNA-interacts with the RNA-binding protein hnRNPK, as well as RNA derived from the MYC promoter and enhancer, to boost MYC transcription by modulating chromatin looping. Our study demonstrates the power and applicability of RIC-seq in discovering the 3D structures, interactions and regulatory roles of RNA.

Proteins of the bromodomain and extra-terminal (BET) domain family are epigenetic readers that bind acetylated histones through their bromodomains to regulate gene transcription. Dual-bromodomain BET inhibitors (DbBi) that bind with similar affinities to the first (BD1) and second (BD2) bromodomains of BRD2, BRD3, BRD4 and BRDt have displayed modest clinical activity in monotherapy cancer trials. A reduced number of thrombocytes in the blood (thrombocytopenia) as well as symptoms of gastrointestinal toxicity are dose-limiting adverse events for some types of DbBi 1-5 . Given that similar haematological and gastrointestinal defects were observed after genetic silencing of Brd4 in mice 6 , the platelet and gastrointestinal toxicities may represent on-target activities associated with BET inhibition. The two individual bromodomains in BET family proteins may have distinct functions 7-9 and different cellular phenotypes after pharmacological inhibition of one or both bromodomains have been reported 10,11 , suggesting that selectively targeting one of the bromodomains may result in a different efficacy and tolerability profile compared with DbBi. Available compounds that are selective to individual domains lack sufficient potency and the pharmacokinetics properties that are required for in vivo efficacy and tolerability assessment 10-13 . Here we carried out a medicinal chemistry campaign that led to the discovery of ABBV-744, a highly potent and selective inhibitor of the BD2 domain of BET family proteins with drug-like properties. In contrast to the broad range of cell growth inhibition induced by DbBi, the antiproliferative activity of ABBV-744 was largely, but not exclusively, restricted to cell lines of acute myeloid leukaemia and prostate cancer that expressed the full-length androgen receptor (AR). ABBV-744 retained robust activity in prostate cancer xenografts, and showed fewer platelet and gastrointestinal toxicities than the DbBi ABBV-075 14 . Analyses of RNA expression and chromatin immunoprecipitation followed by sequencing revealed that ABBV-744 displaced BRD4 from AR-containing super-enhancers and inhibited AR-dependent transcription, with less impact on global transcription compared with ABBV-075. These results underscore the potential value of selectively targeting the BD2 domain of BET family proteins for cancer therapy.

Xist represents a paradigm for the function of long non-coding RNA in epigenetic regulation, although how it mediates X-chromosome inactivation (XCI) remains largely unexplained. Several proteins that bind to Xist RNA have recently been identified, including the transcriptional repressor SPEN 1-3 , the loss of which has been associated with deficient XCI at multiple loci 2-6 . Here we show in mice that SPEN is a key orchestrator of XCI in vivo and we elucidate its mechanism of action. We show that SPEN is essential for initiating gene silencing on the X chromosome in preimplantation mouse embryos and in embryonic stem cells. SPEN is dispensable for maintenance of XCI in neural progenitors, although it significantly decreases the expression of genes that escape XCI. We show that SPEN is immediately recruited to the X chromosome upon the upregulation of Xist, and is targeted to enhancers and promoters of active genes. SPEN rapidly disengages from chromatin upon gene silencing, suggesting that active transcription is required to tether SPEN to chromatin. We define the SPOC domain as a major effector of the gene-silencing function of SPEN, and show that tethering SPOC to Xist RNA is sufficient to mediate gene silencing. We identify the protein partners of SPOC, including NCoR/SMRT, the m 6 A RNA methylation machinery, the NuRD complex, RNA polymerase II and factors involved in the regulation of transcription initiation and elongation. We propose that SPEN acts as a molecular integrator for the initiation of XCI, bridging Xist RNA with the transcription machinery-as well as with nucleosome remodellers and histone deacetylases-at active enhancers and promoters.

Formation of the three primary germ layers during gastrulation is an essential step in the establishment of the vertebrate body plan and is associated with major transcriptional changes 1-5 . Global epigenetic reprogramming accompanies these changes 6-8 , but the role of the epigenome in regulating early cell-fate choice remains unresolved, and the coordination between different molecular layers is unclear. Here we describe a single-cell multi-omics map of chromatin accessibility, DNA methylation and RNA expression during the onset of gastrulation in mouse embryos. The initial exit from pluripotency coincides with the establishment of a global repressive epigenetic landscape, followed by the emergence of lineage-specific epigenetic patterns during gastrulation. Notably, cells committed to mesoderm and endoderm undergo widespread coordinated epigenetic rearrangements at enhancer marks, driven by ten-eleven translocation (TET)-mediated demethylation and a concomitant increase of accessibility. By contrast, the methylation and accessibility landscape of ectodermal cells is already established in the early epiblast. Hence, regulatory elements associated with each germ layer are either epigenetically primed or remodelled before cell-fate decisions, providing the molecular framework for a hierarchical emergence of the primary germ layers.

Epigenetic aberrations are widespread in cancer, yet the underlying mechanisms and causality remain poorly understood 1-3 . A subset of gastrointestinal stromal tumours (GISTs) lack canonical kinase mutations but instead have succinate dehydrogenase (SDH) deficiency and global DNA hyper-methylation 4,5 . Here, we associate this hyper-methylation with changes in genome topology that activate oncogenic programs. To investigate epigenetic alterations systematically, we mapped DNA methylation, CTCF insulators, enhancers, and chromosome topology in KIT-mutant, PDGFRA-mutant and SDH-deficient GISTs. Although these respective subtypes shared similar enhancer landscapes, we identified hundreds of putative insulators where DNA methylation replaced CTCF binding in SDH-deficient GISTs. We focused on a disrupted insulator that normally partitions a core GIST super-enhancer from the FGF4 oncogene. Recurrent loss of this insulator alters locus topology in SDH-deficient GISTs, allowing aberrant physical interaction between enhancer and oncogene. CRISPR-mediated excision of the corresponding CTCF motifs in an SDH-intact GIST model disrupted the boundary between enhancer and oncogene, and strongly upregulated FGF4 expression. We also identified a second recurrent insulator loss event near the KIT oncogene, which is also highly expressed across SDH-deficient GISTs. Finally, we established a patient-derived xenograft (PDX) from an SDH-deficient GIST that faithfully maintains the epigenetics of the parental tumour, including hypermethylation and insulator defects. This PDX model is highly sensitive to FGF receptor (FGFR) inhibition, and more so to combined FGFR and KIT inhibition, validating the functional significance of the underlying epigenetic lesions. Our study reveals how epigenetic alterations can drive oncogenic programs in the absence of canonical kinase mutations, with implications for mechanistic targeting of aberrant pathways in cancers.

The regulated transcription of genes determines cell identity and function. Recent structural studies have elucidated mechanisms that govern the regulation of transcription by RNA polymerases during the initiation and elongation phases. Microscopy studies have revealed that transcription involves the condensation of factors in the cell nucleus. A model is emerging for the transcription of protein-coding genes in which distinct transient condensates form at gene promoters and in gene bodies to concentrate the factors required for transcription initiation and elongation, respectively. The transcribing enzyme RNA polymerase II may shuttle between these condensates in a phosphorylation-dependent manner. Molecular principles are being defined that rationalize transcriptional organization and regulation, and that will guide future investigations.

The synthesis of pre-mRNA by RNA polymerase II (Pol II) involves the formation of a transcription initiation complex, and a transition to an elongation complex 1-4 . The large subunit of Pol II contains an intrinsically disordered C-terminal domain that is phosphorylated by cyclin-dependent kinases during the transition from initiation to elongation, thus influencing the interaction of the C-terminal domain with different components of the initiation or the RNA-splicing apparatus 5,6 . Recent observations suggest that this model provides only a partial picture of the effects of phosphorylation of the C-terminal domain 7-12 . Both the transcription-initiation machinery and the splicing machinery can form phase-separated condensates that contain large numbers of component molecules: hundreds of molecules of Pol II and mediator are concentrated in condensates at super-enhancers 7,8 , and large numbers of splicing factors are concentrated in nuclear speckles, some of which occur at highly active transcription sites 9-12 . Here we investigate whether the phosphorylation of the Pol II C-terminal domain regulates the incorporation of Pol II into phase-separated condensates that are associated with transcription initiation and splicing. We find that the hypophosphorylated C-terminal domain of Pol II is incorporated into mediator condensates and that phosphorylation by regulatory cyclin-dependent kinases reduces this incorporation. We also find that the hyperphosphorylated C-terminal domain is preferentially incorporated into condensates that are formed by splicing factors. These results suggest that phosphorylation of the Pol II C-terminal domain drives an exchange from condensates that are involved in transcription initiation to those that are involved in RNA processing, and implicates phosphorylation as a mechanism that regulates condensate preference.

Large-scale genome sequencing is poised to provide a substantial increase in the rate of discovery of disease-associated mutations, but the functional interpretation of such mutations remains challenging. Here we show that deletions of a sequence on human chromosome 16 that we term the intestine-critical region (ICR) cause intractable congenital diarrhoea in infants 1,2 . Reporter assays in transgenic mice show that the ICR contains a regulatory sequence that activates transcription during the development of the gastrointestinal system. Targeted deletion of the ICR in mice caused symptoms that recapitulated the human condition. Transcriptome analysis revealed that an unannotated open reading frame (Percc1) flanks the regulatory sequence, and the expression of this gene was lost in the developing gut of mice that lacked the ICR. Percc1-knockout mice displayed phenotypes similar to those observed upon ICR deletion in mice and patients, whereas an ICR-driven Percc1 transgene was sufficient to rescue the phenotypes found in mice that lacked the ICR. Together, our results identify a gene that is critical for intestinal function and underscore the need for targeted in vivo studies to interpret the growing number of clinical genetic findings that do not affect known protein-coding genes.

Transcriptional cofactors (COFs) communicate regulatory cues from enhancers to promoters and are central effectors of transcription activation and gene expression 1 . Although some COFs have been shown to prefer certain promoter types 2-5 over others (for example, see refs 6,7 ), the extent to which different COFs display intrinsic specificities for distinct promoters is unclear. Here we use a high-throughput promoter-activity assay in Drosophila melanogaster S2 cells to screen 23 COFs for their ability to activate 72,000 candidate core promoters (CPs). We observe differential activation of CPs, indicating distinct regulatory preferences or 'compatibilities' 8,9 between COFs and specific types of CPs. These functionally distinct CP types are differentially enriched for known sequence elements 2,4 , such as the TATA box, downstream promoter element (DPE) or TCT motif, and display distinct chromatin properties at endogenous loci. Notably, the CP types differ in their relative abundance of H3K4me3 and H3K4me1 marks (see also refs 10-12 ), suggesting that these histone modifications might distinguish trans-regulatory factors rather than promoter- versus enhancer-type cis-regulatory elements. We confirm the existence of distinct COF-CP compatibilities in two additional Drosophila cell lines and in human cells, for which we find COFs that prefer TATA-box or CpG-island promoters, respectively. Distinct compatibilities between COFs and promoters can explain how different enhancers specifically activate distinct sets of genes 9 , alternative promoters within the same genes, and distinct transcription start sites within the same promoter 13 . Thus, COF-promoter compatibilities may underlie distinct transcriptional programs in species as divergent as flies and humans.

Neuronal-activity-dependent transcription couples sensory experience to adaptive responses of the brain including learning and memory. Mechanisms of activity-dependent gene expression including alterations of the epigenome have been characterized 1-8 . However, the fundamental question of whether sensory experience remodels chromatin architecture in the adult brain in vivo to induce neural code transformations and learning and memory remains to be addressed. Here we use in vivo calcium imaging, optogenetics and pharmacological approaches to show that granule neuron activation in the anterior dorsal cerebellar vermis has a crucial role in a delay tactile startle learning paradigm in mice. Of note, using large-scale transcriptome and chromatin profiling, we show that activation of the motor-learning-linked granule neuron circuit reorganizes neuronal chromatin including through long-distance enhancer-promoter and transcriptionally active compartment interactions to orchestrate distinct granule neuron gene expression modules. Conditional CRISPR knockout of the chromatin architecture regulator cohesin in anterior dorsal cerebellar vermis granule neurons in adult mice disrupts enhancer-promoter interactions, activity-dependent transcription and motor learning. These findings define how sensory experience patterns chromatin architecture and neural circuit coding in the brain to drive motor learning.

Precision oncology hinges on linking tumour genotype with molecularly targeted drugs 1 ; however, targeting the frequently dysregulated metabolic landscape of cancer has proven to be a major challenge 2 . Here we show that tissue context is the major determinant of dependence on the nicotinamide adenine dinucleotide (NAD) metabolic pathway in cancer. By analysing more than 7,000 tumours and 2,600 matched normal samples of 19 tissue types, coupled with mathematical modelling and extensive in vitro and in vivo analyses, we identify a simple and actionable set of 'rules'. If the rate-limiting enzyme of de novo NAD synthesis, NAPRT, is highly expressed in a normal tissue type, cancers that arise from that tissue will have a high frequency of NAPRT amplification and be completely and irreversibly dependent on NAPRT for survival. By contrast, tumours that arise from normal tissues that do not express NAPRT highly are entirely dependent on the NAD salvage pathway for survival. We identify the previously unknown enhancer that underlies this dependence. Amplification of NAPRT is shown to generate a pharmacologically actionable tumour cell dependence for survival. Dependence on another rate-limiting enzyme of the NAD synthesis pathway, NAMPT, as a result of enhancer remodelling is subject to resistance by NMRK1-dependent synthesis of NAD. These results identify a central role for tissue context in determining the choice of NAD biosynthetic pathway, explain the failure of NAMPT inhibitors, and pave the way for more effective treatments.