LHX2- and LDB1-mediated trans interactions regulate olfactory receptor choice

Planarian LDB and SSDP proteins scaffold transcriptional complexes for regeneration and patterning

Sequence-specific transcription factors often function as components of large regulatory complexes. LIM-domain binding protein (LDB) and single-stranded DNA-binding protein (SSDP) function as core scaffolds of transcriptional complexes in animals and plants. Little is known about potential partners and functions for LDB/SSDP complexes in the context of tissue regeneration. In this work, we find that planarian LDB1 and SSDP2 promote tissue regeneration, with a particular function in anterior regeneration and mediolateral polarity reestablishment. We find that LDB1 and SSDP2 interact with one another and with characterized planarian LIM-HD proteins Arrowhead, Islet1, and Lhx1/5-1. We also show that SSDP2 and LDB1 function with islet1 in polarity reestablishment and with lhx1/5-1 in serotonergic neuron maturation. Finally, we find new roles for LDB1 and SSDP2 in regulating gene expression in the planarian intestine and parenchyma; these functions are likely LIM-HD-independent. Together, our work provides insight into LDB/SSDP complexes in a highly regenerative organism. Further, our work provides a strong starting point for identifying and characterizing potential binding partners of LDB1 and SSDP2 and for exploring roles for these proteins in diverse aspects of planarian physiology.

Epigenetic and Genetic Profiling of Comorbidity Patterns among Substance Dependence Diagnoses

Objective

This study investigated the genetic and epigenetic mechanisms underlying the comorbidity patterns of five substance dependence diagnoses (SDs; alcohol, AD; cannabis, CaD; cocaine, CoD; opioid, OD; tobacco, TD).

Methods

A latent class analysis (LCA) was performed on 31,197 individuals (average age 42±11 years; 49% females) from six cohorts to identify comorbid DSM-IV SD patterns. In subsets of this sample, we tested SD-latent classes with respect to polygenic burden of psychiatric and behavioral traits and epigenome-wide changes in three population groups.

Results

An LCA identified four latent classes related to SD comorbidities: AD+TD, CoD+TD, AD+CoD+OD+TD (i.e., polysubstance use, PSU), and TD. In the epigenome-wide association analysis, SPATA4 cg02833127 was associated with CoD+TD, AD+TD, and PSU latent classes. AD+TD latent class was also associated with CpG sites located on ARID1B , NOTCH1 , SERTAD4, and SIN3B , while additional epigenome-wide significant associations with CoD+TD latent class were observed in ANO6 and MOV10 genes. PSU-latent class was also associated with a differentially methylated region in LDB1 . We also observed shared polygenic score (PGS) associations for PSU, AD+TD, and CoD+TD latent classes (i.e., attention-deficit hyperactivity disorder, anxiety, educational attainment, and schizophrenia PGS). In contrast, TD-latent class was exclusively associated with posttraumatic stress disorder-PGS. Other specific associations were observed for PSU-latent class (subjective wellbeing-PGS and neuroticism-PGS) and AD+TD-latent class (bipolar disorder-PGS).

Conclusions

We identified shared and unique genetic and epigenetic mechanisms underlying SD comorbidity patterns. These findings highlight the importance of modeling the co-occurrence of SD diagnoses when investigating the molecular basis of addiction-related traits.

The power of proximity: mechanisms and biological roles of transvection

The phenomenon of transvection, defined as a proximity-dependent interallelic interaction, has been observed in the context of complementation between mutant alleles for numerous Drosophila genes. Cases of transvection-like phenomena have also been observed in other species, including mammals. However, the potential contribution of transvection to wild-type gene regulation and the underlying mechanisms remain uncertain. Here, I review recent evidence demonstrating the relevance of transvection in physiological contexts. These findings suggest that transvection represents an additional layer of gene regulation that allows cells to fine-tune gene expression based on the proximity of homologous alleles. In addition, recent studies have measured the physical distance between interacting alleles, revealing unexpectedly large and variable distances. I will discuss how these distances are compatible with the 'hub' model of transcriptional regulation.

Spandrels of the cell nucleus

S.J. Gould and R. Lewontin in their famous "Spandrels paper" (1979) argued that many anatomical elements arise in evolution not due to their "current utility" but rather due to other "reasons for origin", such as other developmental processes, physical constraints and mechanical forces. Here, in the same spirit, we argue that a variety of molecular processes, physical constraints, and mechanical forces, alone or together, generate structures that are detectable in the cell nucleus, yet these structures themselves may not carry any specific function, being a mere reflection of processes that produced them.

Emerging methods and applications in 3D genomics

Since the advent of Hi-C in 2009, a plethora of high-throughput sequencing methods have emerged to profile the three-dimensional (3D) organization of eukaryotic genomes, igniting the era of 3D genomics. In recent years, the genomic resolution achievable by these approaches has dramatically increased and several single-cell versions of Hi-C have been developed. Moreover, a new repertoire of tools not based on proximity ligation of digested chromatin has emerged, enabling the investigation of the higher-order organization of chromatin in the nucleus. In this review, we summarize the expanding portfolio of 3D genomic technologies, highlighting recent developments and applications from the past three years. Lastly, we present an outlook of where this technology-driven field might be headed.

Regulatory mechanisms orchestrating cellular diversity of Cd36+ olfactory sensory neurons revealed by scRNA-seq and scATAC-seq analysis

Recent discoveries have revealed remarkable complexity within olfactory sensory neurons (OSNs), including the existence of two OSN populations based on the expression of Cd36. However, the regulatory mechanisms governing this cellular diversity in the same cell type remain elusive. Here, we show the preferential expression of 79 olfactory receptors in Cd36+ OSNs and the anterior projection characteristics of Cd36+ OSNs, indicating the non-randomness of Cd36 expression. The integrated analysis of single-cell RNA sequencing (scRNA-seq) and scATAC-seq reveals that the differences in Cd36+/- OSNs occur at the immature OSN stage, with Mef2a and Hdac9 being important regulators of developmental divergence. We hypothesize that the absence of Hdac9 may affect the activation of Mef2a, leading to the up-regulation of Mef2a target genes, including teashirt zinc finger family member 1 (Tshz1), in the Cd36+ OSN lineage. We validate that Tshz1 directly promotes Cd36 expression through enhancer bindings. Our study unravels the intricate regulatory landscape and principles governing cellular diversity in the olfactory system.

Prediction and functional interpretation of inter-chromosomal genome architecture from DNA sequence with TwinC

Three-dimensional nuclear DNA architecture comprises well-studied intra-chromosomal (cis) folding and less characterized inter-chromosomal (trans) interfaces. Current predictive models of 3D genome folding can effectively infer pairwise cis-chromatin interactions from the primary DNA sequence but generally ignore trans contacts. There is an unmet need for robust models of trans-genome organization that provide insights into their underlying principles and functional relevance. We present TwinC, an interpretable convolutional neural network model that reliably predicts trans contacts measurable through genome-wide chromatin conformation capture (Hi-C). TwinC uses a paired sequence design from replicate Hi-C experiments to learn single base pair relevance in trans interactions across two stretches of DNA. The method achieves high predictive accuracy (AUROC=0.80) on a cross-chromosomal test set from Hi-C experiments in heart tissue. Mechanistically, the neural network learns the importance of compartments, chromatin accessibility, clustered transcription factor binding and G-quadruplexes in forming trans contacts. In summary, TwinC models and interprets trans genome architecture, shedding light on this poorly understood aspect of gene regulation.

Extrachromosomal DNA Associates with Nuclear Condensates and Reorganizes Chromatin Structures to Enhance Oncogenic Transcription

Extrachromosomal, circular DNA (ecDNA) is a prevalent oncogenic alteration in cancer genomes, often associated with aggressive tumor behavior and poor patient outcome. While previous studies proposed a chromatin-based mobile enhancer model for ecDNA-driven oncogenesis, its precise mechanism and impact remains unclear across diverse cancer types. Our study, utilizing advanced multi-omics profiling, epigenetic editing, and imaging approaches in three cancer models, reveals that ecDNA hubs are an integrated part of nuclear condensates and exhibit cancer-type specific chromatin connectivity. Epigenetic silencing of the ecDNA-specific regulatory modules or chemically disrupting liquid-liquid phase separation breaks down ecDNA hubs, displaces MED1 co-activator binding, inhibits oncogenic transcription, and promotes cell death. These findings substantiate the trans -activator function of ecDNA and underscore a structural mechanism driving oncogenesis. This refined understanding expands our views of oncogene regulation and opens potential avenues for novel therapeutic strategies in cancer treatment.

Molecular, Cellular, and Developmental Organization of the Mouse Vomeronasal organ at Single Cell Resolution

We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors and a population of canonical olfactory sensory neurons in the VNO. High resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and vomeronasal receptors, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.

Chromosome compartmentalization: causes, changes, consequences, and conundrums

The spatial segregation of the genome into compartments is a major feature of 3D genome organization. New data on mammalian chromosome organization across different conditions reveal important information about how and why these compartments form and change. A combination of epigenetic state, nuclear body tethering, physical forces, gene expression, and replication timing (RT) can all influence the establishment and alteration of chromosome compartments. We review the causes and implications of genomic regions undergoing a 'compartment switch' that changes their physical associations and spatial location in the nucleus. About 20-30% of genomic regions change compartment during cell differentiation or cancer progression, whereas alterations in response to a stimulus within a cell type are usually much more limited. However, even a change in 1-2% of genomic bins may have biologically relevant implications. Finally, we review the effects of compartment changes on gene regulation, DNA damage repair, replication, and the physical state of the cell.

YY1-controlled regulatory connectivity and transcription are influenced by the cell cycle

Few transcription factors have been examined for their direct roles in physically connecting enhancers and promoters. Here acute degradation of Yin Yang 1 (YY1) in erythroid cells revealed its requirement for the maintenance of numerous enhancer-promoter loops, but not compartments or domains. Despite its reported ability to interact with cohesin, the formation of YY1-dependent enhancer-promoter loops does not involve stalling of cohesin-mediated loop extrusion. Integrating mitosis-to-G1-phase dynamics, we observed partial retention of YY1 on mitotic chromatin, predominantly at gene promoters, followed by rapid rebinding during mitotic exit, coinciding with enhancer-promoter loop establishment. YY1 degradation during the mitosis-to-G1-phase interval revealed a set of enhancer-promoter loops that require YY1 for establishment during G1-phase entry but not for maintenance in interphase, suggesting that cell cycle stage influences YY1's architectural function. Thus, as revealed here for YY1, chromatin architectural functions of transcription factors can vary in their interplay with CTCF and cohesin as well as by cell cycle stage.

LDB1 establishes multi-enhancer networks to regulate gene expression

How specific enhancer-promoter pairing is established is still mostly unclear. Besides the CTCF/cohesin machinery, only a few nuclear factors have been studied for a direct role in physically connecting regulatory elements. Here, we show via acute degradation experiments that LDB1 directly and broadly promotes enhancer-promoter loops. Most LDB1-mediated contacts, even those spanning hundreds of kb, can form in the absence of CTCF, cohesin, or YY1 as determined via the use of multiple degron systems. Moreover, an engineered LDB1-driven chromatin loop is cohesin independent. Cohesin-driven loop extrusion does not stall at LDB1 occupied sites but may aid the formation of a subset of LDB1 anchored loops. Leveraging the dynamic reorganization of nuclear architecture during the transition from mitosis to G1-phase, we establish a relationship between LDB1-dependent interactions in the context of TAD organization and gene activation. Lastly, Tri-C and Region Capture Micro-C reveal that LDB1 organizes multi-enhancer networks to activate transcription. This establishes LDB1 as a direct driver of regulatory network inter-connectivity.

Barcoding distinct neurons

The genomic landscape of a cell surface protein reveals how neuron identity is displayed.

Cis-regulatory evolution of the recently expanded Ly49 gene family

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

The influence of olfactory experience on the birthrates of olfactory sensory neurons with specific odorant receptor identities

Olfactory sensory neurons (OSNs) are one of a few neuron types that are generated continuously throughout life in mammals. The persistence of olfactory sensory neurogenesis beyond early development has long been thought to function simply to replace neurons that are lost or damaged through exposure to environmental insults. The possibility that olfactory sensory neurogenesis may also serve an adaptive function has received relatively little consideration, largely due to the assumption that the generation of new OSNs is stochastic with respect to OSN subtype, as defined by the single odorant receptor gene that each neural precursor stochastically chooses for expression out of hundreds of possibilities. Accordingly, the relative birthrates of different OSN subtypes are predicted to be constant and impervious to olfactory experience. This assumption has been called into question, however, by evidence that the birthrates of specific OSN subtypes can be selectively altered by manipulating olfactory experience through olfactory deprivation, enrichment, and conditioning paradigms. Moreover, studies of recovery of the OSN population following injury provide further evidence that olfactory sensory neurogenesis may not be strictly stochastic with respect to subtype. Here we review this evidence and consider mechanistic and functional implications of the prospect that specific olfactory experiences can regulate olfactory sensory neurogenesis rates in a subtype-selective manner.

Simultaneous single-cell three-dimensional genome and gene expression profiling uncovers dynamic enhancer connectivity underlying olfactory receptor choice

The simultaneous measurement of three-dimensional (3D) genome structure and gene expression of individual cells is critical for understanding a genome's structure-function relationship, yet this is challenging for existing methods. Here we present 'Linking mRNA to Chromatin Architecture (LiMCA)', which jointly profiles the 3D genome and transcriptome with exceptional sensitivity and from low-input materials. Combining LiMCA and our high-resolution scATAC-seq assay, METATAC, we successfully characterized chromatin accessibility, as well as paired 3D genome structures and gene expression information, of individual developing olfactory sensory neurons. We expanded the repertoire of known olfactory receptor (OR) enhancers and discovered unexpected rules of their dynamics: OR genes and their enhancers are most accessible during early differentiation. Furthermore, we revealed the dynamic spatial relationship between ORs and enhancers behind stepwise OR expression. These findings offer valuable insights into how 3D connectivity of ORs and enhancers dynamically orchestrate the 'one neuron-one receptor' selection process.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Establishment and maintenance of random monoallelic expression

This Review elucidates the regulatory principles of random monoallelic expression by focusing on two well-studied examples: the X-chromosome inactivation regulator Xist and the olfactory receptor gene family. Although the choice of a single X chromosome or olfactory receptor occurs in different developmental contexts, common gene regulatory principles guide monoallelic expression in both systems. In both cases, an event breaks the symmetry between genetically and epigenetically identical copies of the gene, leading to the expression of one single random allele, stabilized through negative feedback control. Although many regulatory steps that govern the establishment and maintenance of monoallelic expression have been identified, key pieces of the puzzle are still missing. We provide an overview of the current knowledge and models for the monoallelic expression of Xist and olfactory receptors. We discuss their similarities and differences, and highlight open questions and approaches that could guide the study of other monoallelically expressed genes.

Potential roles of inter-chromosomal interactions in cell fate determination

Mammalian genomic DNA is packed in a small nucleus, and its folding and organization in the nucleus are critical for gene regulation and cell fate determination. In interphase, chromosomes are compartmentalized into certain nuclear spaces and territories that are considered incompatible with each other. The regulation of gene expression is influenced by the epigenetic characteristics of topologically associated domains and A/B compartments within chromosomes (intrachromosomal). Previously, interactions among chromosomes detected via chromosome conformation capture-based methods were considered noise or artificial errors. However, recent studies based on newly developed ligation-independent methods have shown that inter-chromosomal interactions play important roles in gene regulation. This review summarizes the recent understanding of spatial genomic organization in mammalian interphase nuclei and discusses the potential mechanisms that determine cell identity. In addition, this review highlights the potential role of inter-chromosomal interactions in early mouse development.

Transcriptional priming and chromatin regulation during stochastic cell fate specification

Stochastic cell fate specification, in which a cell chooses between two or more fates with a set probability, diversifies cell subtypes in development. Although this is a vital process across species, a common mechanism for these cell fate decisions remains elusive. This review examines two well-characterized stochastic cell fate decisions to identify commonalities between their developmental programmes. In the fly eye, two subtypes of R7 photoreceptors are specified by the stochastic ON/OFF expression of a transcription factor, spineless. In the mouse olfactory system, olfactory sensory neurons (OSNs) randomly select to express one copy of an olfactory receptor (OR) gene out of a pool of 2800 alleles. Despite the differences in these sensory systems, both stochastic fate choices rely on the dynamic interplay between transcriptional priming, chromatin regulation and terminal gene expression. The coupling of transcription and chromatin modifications primes gene loci in undifferentiated neurons, enabling later expression during terminal differentiation. Here, we compare these mechanisms, examine broader implications for gene regulation during development and posit key challenges moving forward. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Epigenetic programming of stochastic olfactory receptor choice

The mammalian sense of smell relies upon a vast array of receptor proteins to detect odorant compounds present in the environment. The proper deployment of these receptor proteins in olfactory sensory neurons is orchestrated by a suite of epigenetic processes that remodel the olfactory genes in differentiating neuronal progenitors. The goal of this review is to elucidate the central role of gene regulatory processes acting in neuronal progenitors of olfactory sensory neurons that lead to a singular expression of an odorant receptor in mature olfactory sensory neurons. We begin by describing the principal features of odorant receptor gene expression in mature olfactory sensory neurons. Next, we delineate our current understanding of how these features emerge from multiple gene regulatory mechanisms acting in neuronal progenitors. Finally, we close by discussing the key gaps in our understanding of how these regulatory mechanisms work and how they interact with each other over the course of differentiation.

Insertion of a neomycin selection cassette in the Amigo1 locus alters gene expression in the olfactory epithelium leading to region-specific defects in olfactory receptor neuron development

During development of the nervous system, neurons connect to one another in a precisely organized manner. Sensory systems provide a good example of this organization, whereby the composition of the outside world is represented in the brain by neuronal maps. Establishing correct patterns of neural circuitry is crucial, as inaccurate map formation can lead to severe disruptions in sensory processing. In rodents, olfactory stimuli modulate a wide variety of behaviors essential for survival. The formation of the olfactory glomerular map is dependent on molecular cues that guide olfactory receptor neuron axons to broad regions of the olfactory bulb and on cell adhesion molecules that promote axonal sorting into specific synaptic units in this structure. Here, we demonstrate that the cell adhesion molecule Amigo1 is expressed in a subpopulation of olfactory receptor neurons, and we investigate its role in the precise targeting of olfactory receptor neuron axons to the olfactory bulb using a genetic loss-of-function approach in mice. While ablation of Amigo1 did not lead to alterations in olfactory sensory neuron axonal targeting, our experiments revealed that the presence of a neomycin resistance selection cassette in the Amigo1 locus can lead to off-target effects that are not due to loss of Amigo1 expression, including unexpected altered gene expression in olfactory receptor neurons and reduced glomerular size in the ventral region of the olfactory bulb. Our results demonstrate that insertion of a neomycin selection cassette into the mouse genome can have specific deleterious effects on the development of the olfactory system and highlight the importance of removing antibiotic resistance cassettes from genetic loss-of-function mouse models when studying olfactory system development.

Molecular mechanisms of differentiation and class choice of olfactory sensory neurons

The sense of smell is intricately linked to essential animal behaviors necessary for individual survival and species preservation. During vertebrate evolution, odorant receptors (ORs), responsible for detecting odor molecules, have evolved to adapt to changing environments, transitioning from aquatic to terrestrial habitats and accommodating increasing complex chemical environments. These evolutionary pressures have given rise to the largest gene family in vertebrate genomes. Vertebrate ORs are phylogenetically divided into two major classes; class I and class II. Class I OR genes, initially identified in fish and frog, have persisted across vertebrate species. On the other hand, class II OR genes are unique to terrestrial animals, accounting for ~90% of mammalian OR genes. In mice, each olfactory sensory neuron (OSN) expresses a single functional allele of a single OR gene from either the class I or class II OR repertoire. This one neuron-one receptor rule is established through two sequential steps: specification of OR class and subsequent exclusive OR expression from the corresponding OR class. Consequently, OSNs acquire diverse neuronal identities during the process of OSN differentiation, enabling animals to detect a wide array of odor molecules. This review provides an overview of the OSN differentiation process through which OSN diversity is achieved, primarily using the mouse as a model animal.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

Nescient helix-loop-helix 1 (Nhlh1) is a novel activating transcription factor 5 (ATF5) target gene in olfactory and vomeronasal sensory neurons in mice

Activating transcription factor 5 (ATF5) is a transcription factor that belongs to the cAMP-response element-binding protein/ATF family and is essential for the differentiation and survival of sensory neurons in mouse olfactory organs. However, transcriptional target genes for ATF5 have yet to be identified. In the present study, chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) experiments were performed to verify ATF5 target genes in the main olfactory epithelium and vomeronasal organ in the postnatal pups. ChIP-qPCR was conducted using hemagglutinin (HA)-tagged ATF5 knock-in olfactory organs. The results obtained demonstrated that ATF5-HA fusion proteins bound to the CCAAT/enhancer-binding protein-ATF response element (CARE) site in the enhancer region of nescient helix-loop-helix 1 (Nhlh1), a transcription factor expressed in differentiating olfactory and vomeronasal sensory neurons. Nhlh1 mRNA expression was downregulated in ATF5-deficient (ATF5-/-) olfactory organs. The LIM/homeobox protein transcription factor Lhx2 co-localized with ATF5 in the nuclei of olfactory and vomeronasal sensory neurons and bound to the homeodomain site proximal to the CARE site in the Nhlh1 gene. The CARE region of the Nhlh1 gene was enriched by the active enhancer marker, acetyl-histone H3 (Lys27). The present study identified Nhlh1 as a novel target gene for ATF5 in murine olfactory organs. ATF5 may upregulate Nhlh1 expression in concert with Lhx2, thereby promoting the differentiation of olfactory and vomeronasal sensory neurons.

Identification and Validation of eRNA as a Prognostic Indicator for Cervical Cancer

The survival of CESC patients is closely related to the expression of enhancer RNA (eRNA). In this work, we downloaded eRNA expression, clinical, and gene expression data from the TCeA and TCGA portals. A total of 7936 differentially expressed eRNAs were discovered by limma analysis, and the relationship between these eRNAs and survival was analyzed by univariate Cox hazard analysis, LASSO regression, and multivariate Cox hazard analysis to obtain an 8-eRNA model. Risk score heat maps, KM curves, ROC analysis, robustness analysis, and nomograms further indicate that this 8-eRNA model is a novel indicator with high prognostic performance independent of clinicopathological classification. The model divided patients into high-risk and low-risk groups, compared pathway diversity between the two groups through GSEA analysis, and provided potential therapeutic agents for high-risk patients.

Chromatin gene-gene loops support the cross-regulation of genes with related function

Chromatin loops between gene pairs have been observed in diverse contexts in both flies and vertebrates. Combining high-resolution Capture-C, DNA fluorescence in situ hybridization, and genetic perturbations, we dissect the functional role of three loops between genes with related function during Drosophila embryogenesis. By mutating the loop anchor (but not the gene) or the gene (but not loop anchor), we disentangle loop formation and gene expression and show that the 3D proximity of paralogous gene loci supports their co-regulation. Breaking the loop leads to either an attenuation or enhancement of expression and perturbs their relative levels of expression and cross-regulation. Although many loops appear constitutive across embryogenesis, their function can change in different developmental contexts. Taken together, our results indicate that chromatin gene-gene loops act as architectural scaffolds that can be used in different ways in different contexts to fine-tune the coordinated expression of genes with related functions and sustain their cross-regulation.

A noncoding role of coding mRNA in monogenic olfactory receptor choice

In a recent study, Pourmorady and colleagues uncovered a noncoding role for olfactory receptor (OR)-coding mRNA in mediating nuclear architecture and singular OR choice. The OR mRNAs reinforce the prevailing enhancer hub and inhibit other competitors, facilitating transition from polygenic to singular OR expression.

SpecLoop predicts cell type-specific chromatin loop via transcription factor cooperation

Cell-type-Specific Chromatin Loops (CSCLs) are crucial for gene regulation and cell fate determination. However, the mechanisms governing their establishment remain elusive. Here, we present SpecLoop, a network regularization-based machine learning framework, to investigate the role of transcription factors (TFs) cooperation in CSCL formation. SpecLoop integrates multi-omics data, including gene expression, chromatin accessibility, sequence, protein-protein interaction, and TF binding motif data, to predict CSCLs and identify TF cooperations. Using high resolution Hi-C data as the gold standard, SpecLoop accurately predicts CSCL in GM12878, IMR90, HeLa-S3, K562, HUVEC, HMEC, and NHEK seven cell types, with the AUROC values ranging from 0.8645 to 0.9852 and AUPR values ranging from 0.8654 to 0.9734. Notably SpecLoop demonstrates improved accuracy in predicting long-distance CSCLs and identifies TF complexes with strong predictive ability. Our study systematically explores the TFs and TF pairs associated with CSCL through effective integration of diverse omics data. SpecLoop is freely available at https://github.com/AMSSwanglab/SpecLoop.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Maintenance of neuronal identity in C. elegans and beyond: Lessons from transcription and chromatin factors

Neurons are remarkably long-lived, non-dividing cells that must maintain their functional features (e.g., electrical properties, chemical signaling) for extended periods of time - decades in humans. How neurons accomplish this incredible feat is poorly understood. Here, we review recent advances, primarily in the nematode C. elegans, that have enhanced our understanding of the molecular mechanisms that enable post-mitotic neurons to maintain their functionality across different life stages. We begin with "terminal selectors" - transcription factors necessary for the establishment and maintenance of neuronal identity. We highlight new findings on five terminal selectors (CHE-1 [Glass], UNC-3 [Collier/Ebf1-4], LIN-39 [Scr/Dfd/Hox4-5], UNC-86 [Acj6/Brn3a-c], AST-1 [Etv1/ER81]) from different transcription factor families (ZNF, COE, HOX, POU, ETS). We compare the functions of these factors in specific neuron types of C. elegans with the actions of their orthologs in other invertebrate (D. melanogaster) and vertebrate (M. musculus) systems, highlighting remarkable functional conservation. Finally, we reflect on recent findings implicating chromatin-modifying proteins, such as histone methyltransferases and Polycomb proteins, in the control of neuronal terminal identity. Altogether, these new studies on transcription factors and chromatin modifiers not only shed light on the fundamental problem of neuronal identity maintenance, but also outline mechanistic principles of gene regulation that may operate in other long-lived, post-mitotic cell types.

RNA-mediated symmetry breaking enables singular olfactory receptor choice

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 alleles1,2. OR gene choice is mediated by a multichromosomal enhancer hub that activates transcription at a single OR3,4, followed by OR-translation-dependent feedback that stabilizes this choice5,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.

Opposing, spatially-determined epigenetic forces impose restrictions on stochastic olfactory receptor choice

Olfactory receptor (OR) choice represents an example of genetically hardwired stochasticity, where every olfactory neuron expresses one out of ~2000 OR alleles in the mouse genome in a probabilistic, yet stereotypic fashion. Here, we propose that topographic restrictions in OR expression are established in neuronal progenitors by two opposing forces: polygenic transcription and genomic silencing, both of which are influenced by dorsoventral gradients of transcription factors NFIA, B, and X. Polygenic transcription of OR genes may define spatially constrained OR repertoires, among which one OR allele is selected for singular expression later in development. Heterochromatin assembly and genomic compartmentalization of OR alleles also vary across the axes of the olfactory epithelium and may preferentially eliminate ectopically expressed ORs with more dorsal expression destinations from this 'privileged' repertoire. Our experiments identify early transcription as a potential 'epigenetic' contributor to future developmental patterning and reveal how two spatially responsive probabilistic processes may act in concert to establish deterministic, precise, and reproducible territories of stochastic gene expression.

An allele-selective inter-chromosomal protein bridge supports monogenic antigen expression in the African trypanosome

UPF1-like helicases play roles in telomeric heterochromatin formation and X-chromosome inactivation, and also in monogenic variant surface glycoprotein (VSG) expression via VSG exclusion-factor-2 (VEX2), a UPF1-related protein in the African trypanosome. We show that VEX2 associates with chromatin specifically at the single active VSG expression site on chromosome 6, forming an allele-selective connection, via VEX1, to the trans-splicing locus on chromosome 9, physically bridging two chromosomes and the VSG transcription and splicing compartments. We further show that the VEX-complex is multimeric and self-regulates turnover to tightly control its abundance. Using single cell transcriptomics following VEX2-depletion, we observed simultaneous derepression of many other telomeric VSGs and multi-allelic VSG expression in individual cells. Thus, an allele-selective, inter-chromosomal, and self-limiting VEX1-2 bridge supports monogenic VSG expression and multi-allelic VSG exclusion.

Joint epigenome profiling reveals cell-type-specific gene regulatory programmes in human cortical organoids

Gene expression is regulated by multiple epigenetic mechanisms, which are coordinated in development and disease. However, current multiomics methods are frequently limited to one or two modalities at a time, making it challenging to obtain a comprehensive gene regulatory signature. Here, we describe a method-3D genome, RNA, accessibility and methylation sequencing (3DRAM-seq)-that simultaneously interrogates spatial genome organization, chromatin accessibility and DNA methylation genome-wide and at high resolution. We combine 3DRAM-seq with immunoFACS and RNA sequencing in cortical organoids to map the cell-type-specific regulatory landscape of human neural development across multiple epigenetic layers. Finally, we apply a massively parallel reporter assay to profile cell-type-specific enhancer activity in organoids and to functionally assess the role of key transcription factors for human enhancer activation and function. More broadly, 3DRAM-seq can be used to profile the multimodal epigenetic landscape in rare cell types and different tissues.

Lineage-specific 3D genome organization is assembled at multiple scales by IKAROS

A generic level of chromatin organization generated by the interplay between cohesin and CTCF suffices to limit promiscuous interactions between regulatory elements, but a lineage-specific chromatin assembly that supersedes these constraints is required to configure the genome to guide gene expression changes that drive faithful lineage progression. Loss-of-function approaches in B cell precursors show that IKAROS assembles interactions across megabase distances in preparation for lymphoid development. Interactions emanating from IKAROS-bound enhancers override CTCF-imposed boundaries to assemble lineage-specific regulatory units built on a backbone of smaller invariant topological domains. Gain of function in epithelial cells confirms IKAROS' ability to reconfigure chromatin architecture at multiple scales. Although the compaction of the Igκ locus required for genome editing represents a function of IKAROS unique to lymphocytes, the more general function to preconfigure the genome to support lineage-specific gene expression and suppress activation of extra-lineage genes provides a paradigm for lineage restriction.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Systematic identification of inter-chromosomal interaction networks supports the existence of RNA factories

Most studies of genome organization have focused on intra-chromosomal (cis) contacts because they harbor key features such as DNA loops and topologically associating domains. Inter-chromosomal (trans) contacts have received much less attention, and tools for interrogating potential biologically relevant trans structures are lacking. Here, we develop a computational framework to identify sets of loci that jointly interact in trans from Hi-C data. This method, trans-C, initiates probabilistic random walks with restarts from a set of seed loci to traverse an input Hi-C contact network, thereby identifying sets of trans-contacting loci. We validate trans-C in three increasingly complex models of established trans contacts: the Plasmodium falciparum var genes, the mouse olfactory receptor "Greek islands", and the human RBM20 cardiac splicing factory. We then apply trans-C to systematically test the hypothesis that genes co-regulated by the same trans-acting element (i.e., a transcription or splicing factor) co-localize in three dimensions to form "RNA factories" that maximize the efficiency and accuracy of RNA biogenesis. We find that many loci with multiple binding sites of the same transcription factor interact with one another in trans, especially those bound by transcription factors with intrinsically disordered domains. Similarly, clustered binding of a subset of RNA binding proteins correlates with trans interaction of the encoding loci. These findings support the existence of trans interacting chromatin domains (TIDs) driven by RNA biogenesis. Trans-C provides an efficient computational framework for studying these and other types of trans interactions, empowering studies of a poorly understood aspect of genome architecture.

Opposing, spatially-determined epigenetic forces impose restrictions on stochastic olfactory receptor choice

Olfactory receptor (OR) choice represents an example of genetically hardwired stochasticity, where every olfactory neuron expresses one out of ~2000 OR alleles in a probabilistic, yet stereotypic fashion. Here, we propose that topographic restrictions in OR expression are established in neuronal progenitors by two opposing forces: polygenic transcription and genomic silencing, both of which are influenced by dorsoventral gradients of transcription factors NFIA, B, and X. Polygenic transcription of OR genes may define spatially constrained OR repertoires, among which one OR allele is selected for singular expression later in development. Heterochromatin assembly and genomic compartmentalization of OR alleles also vary across the axes of the olfactory epithelium and may preferentially eliminate ectopically expressed ORs with more dorsal expression destinations from this "privileged" repertoire. Our experiments identify early transcription as a potential "epigenetic" contributor to future developmental patterning and reveal how two spatially responsive probabilistic processes may act in concert to establish deterministic, precise, and reproducible territories of stochastic gene expression.

Lifelong restructuring of 3D genome architecture in cerebellar granule cells

The cerebellum contains most of the neurons in the human brain and exhibits distinctive modes of development and aging. In this work, by developing our single-cell three-dimensional (3D) genome assay-diploid chromosome conformation capture, or Dip-C-into population-scale (Pop-C) and virus-enriched (vDip-C) modes, we resolved the first 3D genome structures of single cerebellar cells, created life-spanning 3D genome atlases for both humans and mice, and jointly measured transcriptome and chromatin accessibility during development. We found that although the transcriptome and chromatin accessibility of cerebellar granule neurons mature in early postnatal life, 3D genome architecture gradually remodels throughout life, establishing ultra-long-range intrachromosomal contacts and specific interchromosomal contacts that are rarely seen in neurons. These results reveal unexpected evolutionarily conserved molecular processes that underlie distinctive features of neural development and aging across the mammalian life span.

Total-Activity Conservation Analysis and Design of Boolean Networks

The law of conservation of mass, represented in Boolean networks (BNs) as total-activity conservation, is one of the typical properties of biological networks. This article analyzes the total-activity conservation of BNs based on the algebraic state-space representation (ASSR) approach. First, the total-activity-conservative matrix is defined and a matrix-based criterion is proposed to verify the total-activity conservation of BNs. Meanwhile, when function perturbation is considered, robust total-activity conservation is investigated. Second, by means of the pseudo-Boolean function generated by total-activity conservation, a constructive design procedure of the Boolean dynamics is given to achieve the total-activity conservation. Third, the total-activity conservation of switched Boolean networks (SBNs) under arbitrary switching signal is studied, which together with network aggregation achieves the total-activity conservation of large-scale BNs.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Chromatin loop dynamics during cellular differentiation are associated with changes to both anchor and internal regulatory features

Three-dimensional (3D) chromatin structure has been shown to play a role in regulating gene transcription during biological transitions. Although our understanding of loop formation and maintenance is rapidly improving, much less is known about the mechanisms driving changes in looping and the impact of differential looping on gene transcription. One limitation has been a lack of well-powered differential looping data sets. To address this, we conducted a deeply sequenced Hi-C time course of megakaryocyte development comprising four biological replicates and 6 billion reads per time point. Statistical analysis revealed 1503 differential loops. Gained loop anchors were enriched for AP-1 occupancy and were characterized by large increases in histone H3K27ac (over 11-fold) but relatively small increases in CTCF and RAD21 binding (1.26- and 1.23-fold, respectively). Linear modeling revealed that changes in histone H3K27ac, chromatin accessibility, and JUN binding were better correlated with changes in looping than RAD21 and almost as well correlated as CTCF. Changes to epigenetic features between-rather than at-boundaries were highly predictive of changes in looping. Together these data suggest that although CTCF and RAD21 may be the core machinery dictating where loops form, other features (both at the anchors and within the loop boundaries) may play a larger role than previously anticipated in determining the relative loop strength across cell types and conditions.

Mega-Enhancer Bodies Organize Neuronal Long Genes in the Cerebellum

Dynamic regulation of gene expression plays a key role in establishing the diverse neuronal cell types in the brain. Recent findings in genome biology suggest that three-dimensional (3D) genome organization has important, but mechanistically poorly understood functions in gene transcription. Beyond local genomic interactions between promoters and enhancers, we find that cerebellar granule neurons undergoing differentiation in vivo exhibit striking increases in long-distance genomic interactions between transcriptionally active genomic loci, which are separated by tens of megabases within a chromosome or located on different chromosomes. Among these interactions, we identify a nuclear subcompartment enriched for near-megabase long enhancers and their associated neuronal long genes encoding synaptic or signaling proteins. Neuronal long genes are differentially recruited to this enhancer-dense subcompartment to help shape the transcriptional identities of granule neuron subtypes in the cerebellum. SPRITE analyses of higher-order genomic interactions, together with IGM-based 3D genome modeling and imaging approaches, reveal that the enhancer-dense subcompartment forms prominent nuclear structures, which we term mega-enhancer bodies. These novel nuclear bodies reside in the nuclear periphery, away from other transcriptionally active structures, including nuclear speckles located in the nuclear interior. Together, our findings define additional layers of higher-order 3D genome organization closely linked to neuronal maturation and identity in the brain.

LHX2 haploinsufficiency causes a variable neurodevelopmental disorder

PURPOSE

LHX2 encodes the LIM homeobox 2 transcription factor (LHX2), which is highly expressed in brain and well conserved across species, but it has not been clearly linked to neurodevelopmental disorders (NDDs) to date.

METHODS

Through international collaboration, we identified 19 individuals from 18 families with variable neurodevelopmental phenotypes, carrying a small chromosomal deletion, likely gene-disrupting or missense variants in LHX2. Functional consequences of missense variants were investigated in cellular systems.

RESULTS

Affected individuals presented with developmental and/or behavioral abnormalities, autism spectrum disorder, variable intellectual disability, and microcephaly. We observed nucleolar accumulation for 2 missense variants located within the DNA-binding HOX domain, impaired interaction with co-factor LDB1 for another variant located in the protein-protein interaction-mediating LIM domain, and impaired transcriptional activation by luciferase assay for 4 missense variants.

CONCLUSION

We implicate LHX2 haploinsufficiency by deletion and likely gene-disrupting variants as causative for a variable NDD. Our findings suggest a loss-of-function mechanism also for likely pathogenic LHX2 missense variants. Together, our observations underscore the importance of LHX2 in the nervous system and for variable neurodevelopmental phenotypes.

Single cell Hi-C identifies plastic chromosome conformations underlying the gastrulation enhancer landscape

Embryonic development involves massive proliferation and differentiation of cell lineages. This must be supported by chromosome replication and epigenetic reprogramming, but how proliferation and cell fate acquisition are balanced in this process is not well understood. Here we use single cell Hi-C to map chromosomal conformations in post-gastrulation mouse embryo cells and study their distributions and correlations with matching embryonic transcriptional atlases. We find that embryonic chromosomes show a remarkably strong cell cycle signature. Despite that, replication timing, chromosome compartment structure, topological associated domains (TADs) and promoter-enhancer contacts are shown to be variable between distinct epigenetic states. About 10% of the nuclei are identified as primitive erythrocytes, showing exceptionally compact and organized compartment structure. The remaining cells are broadly associated with ectoderm and mesoderm identities, showing only mild differentiation of TADs and compartment structures, but more specific localized contacts in hundreds of ectoderm and mesoderm promoter-enhancer pairs. The data suggest that while fully committed embryonic lineages can rapidly acquire specific chromosomal conformations, most embryonic cells are showing plastic signatures driven by complex and intermixed enhancer landscapes.