Enhancer redundancy in development and disease

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

Predicting genome-wide tissue-specific enhancers via combinatorial transcription factor genomic occupancy analysis

Enhancers are non-coding cis-regulatory elements crucial for transcriptional regulation. Mutations in enhancers can disrupt gene regulation, leading to disease phenotypes. Identifying enhancers and their tissue-specific activity is challenging due to their lack of stereotyped sequences. This study presents a sequence-based computational model that uses combinatorial transcription factor (TF) genomic occupancy to predict tissue-specific enhancers. Trained on diverse datasets, including ENCODE and Vista enhancer browser data, the model predicted 25 000 forebrain-specific cis-regulatory modules (CRMs) in the human genome. Validation using biochemical features, disease-associated SNPs, and in vivo zebrafish analysis confirmed its effectiveness. This model aids in predicting enhancers lacking well-characterized chromatin features, complementing experimental approaches in tissue-specific enhancer discovery.

EnhancerNet: a predictive model of cell identity dynamics through enhancer selection

Understanding how cell identity is encoded by the genome and acquired during differentiation is a central challenge in cell biology. I have developed a theoretical framework called EnhancerNet, which models the regulation of cell identity through the lens of transcription factor-enhancer interactions. I demonstrate that autoregulation in these interactions imposes a constraint on the model, resulting in simplified dynamics that can be parameterized from observed cell identities. Despite its simplicity, EnhancerNet recapitulates a broad range of experimental observations on cell identity dynamics, including enhancer selection, cell fate induction, hierarchical differentiation through multipotent progenitor states and direct reprogramming by transcription factor overexpression. The model makes specific quantitative predictions, reproducing known reprogramming recipes and the complex haematopoietic differentiation hierarchy without fitting unobserved parameters. EnhancerNet provides insights into how new cell types could evolve and highlights the functional importance of distal regulatory elements with dynamic chromatin in multicellular evolution.

A Ctnnb1 enhancer transcriptionally regulates Wnt signaling dosage to balance homeostasis and tumorigenesis of intestinal epithelia

The β-catenin-dependent canonical Wnt signaling is pivotal in organ development, tissue homeostasis, and cancer. Here, we identified an upstream enhancer of Ctnnb1 - the coding gene for β-catenin, named ieCtnnb1 (intestinal enhancer of Ctnnb1), which is crucial for intestinal homeostasis. ieCtnnb1 is predominantly active in the base of small intestinal crypts and throughout the epithelia of large intestine. Knockout of ieCtnnb1 led to a reduction in Ctnnb1 transcription, compromising the canonical Wnt signaling in intestinal crypts. Single-cell sequencing revealed that ieCtnnb1 knockout altered epithelial compositions and potentially compromised functions of small intestinal crypts. While deletion of ieCtnnb1 hampered epithelial turnovers in physiologic conditions, it prevented occurrence and progression of Wnt/β-catenin-driven colorectal cancers. Human ieCTNNB1 drove reporter gene expression in a pattern highly similar to mouse ieCtnnb1. ieCTNNB1 contains a single-nucleotide polymorphism associated with CTNNB1 expression levels in human gastrointestinal epithelia. The enhancer activity of ieCTNNB1 in colorectal cancer tissues was stronger than that in adjacent normal tissues. HNF4α and phosphorylated CREB1 were identified as key trans-factors binding to ieCTNNB1 and regulating CTNNB1 transcription. Together, these findings unveil an enhancer-dependent mechanism controlling the dosage of Wnt signaling and homeostasis in intestinal epithelia.

Interpreting cis-regulatory interactions from large-scale deep neural networks

The rise of large-scale, sequence-based deep neural networks (DNNs) for predicting gene expression has introduced challenges in their evaluation and interpretation. Current evaluations align DNN predictions with orthogonal experimental data, providing insights into generalization but offering limited insights into their decision-making process. Existing model explainability tools focus mainly on motif analysis, which becomes complex when interpreting longer sequences. Here we present cis-regulatory element model explanations (CREME), an in silico perturbation toolkit that interprets the rules of gene regulation learned by a genomic DNN. Applying CREME to Enformer, a state-of-the-art DNN, we identify cis-regulatory elements that enhance or silence gene expression and characterize their complex interactions. CREME can provide interpretations across multiple scales of genomic organization, from cis-regulatory elements to fine-mapped functional sequence elements within them, offering high-resolution insights into the regulatory architecture of the genome. CREME provides a powerful toolkit for translating the predictions of genomic DNNs into mechanistic insights of gene regulation.

The Oct4-related PouV gene, pou5f3, mediates isthmus development in zebrafish by directly and dynamically regulating pax2a

Using a transgenic zebrafish line harboring a heat-inducible dominant-interference pou5f3 gene (en-pou5f3), we reported that this PouV gene is involved in isthmus development at the midbrain-hindbrain boundary (MHB), which patterns the midbrain and cerebellum. Importantly, the functions of pou5f3 reportedly differ before and after the end of gastrulation. In the present study, we examined in detail the effects of en-pou5f3 induction on isthmus development during embryogenesis. When en-pou5f3 was induced around the end of gastrulation (bud stage), the isthmus was abrogated or deformed by the end of somitogenesis (24 hours post-fertilization). At this stage, the expression of MHB markers -- such as pax2a, fgf8a, wnt1, and gbx2 -- was absent in embryos lacking the isthmus structure, whereas it was present, although severely distorted, in embryos with a deformed isthmus. We further found that, after en-pou5f3 induction at late gastrulation, pax2a, fgf8a, and wnt1 were immediately and irreversibly downregulated, whereas the expression of en2a and gbx2 was reduced only weakly and slowly. Induction of en-pou5f3 at early somite stages also immediately downregulated MHB genes, particularly pax2a, but their expression was restored later. Overall, the data suggested that pou5f3 directly upregulates at least pax2a and possibly fgf8a and wnt1, which function in parallel in establishing the MHB, and that the role of pou5f3 dynamically changes around the end of gastrulation. We next examined the transcriptional regulation of pax2a using both in vitro and in vivo reporter analyses; the results showed that two upstream 1.0-kb regions with sequences conserved among vertebrates specifically drove transcription at the MHB. These reporter analyses confirmed that development of the isthmic organizer is regulated by PouV through direct regulation of pax2/pax2a in vertebrate embryos.

Integration of chromosome locations and functional aspects of enhancers and topologically associating domains in knowledge graphs enables versatile queries about gene regulation

Knowledge about transcription factor binding and regulation, target genes, cis-regulatory modules and topologically associating domains is not only defined by functional associations like biological processes or diseases but also has a determinative genome location aspect. Here, we exploit these location and functional aspects together to develop new strategies to enable advanced data querying. Many databases have been developed to provide information about enhancers, but a schema that allows the standardized representation of data, securing interoperability between resources, has been lacking. In this work, we use knowledge graphs for the standardized representation of enhancers and topologically associating domains, together with data about their target genes, transcription factors, location on the human genome, and functional data about diseases and gene ontology annotations. We used this schema to integrate twenty-five enhancer datasets and two domain datasets, creating the most powerful integrative resource in this field to date. The knowledge graphs have been implemented using the Resource Description Framework and integrated within the open-access BioGateway knowledge network, generating a resource that contains an interoperable set of knowledge graphs (enhancers, TADs, genes, proteins, diseases, GO terms, and interactions between domains). We show how advanced queries, which combine functional and location restrictions, can be used to develop new hypotheses about functional aspects of gene expression regulation.

H3K4me2 distinguishes a distinct class of enhancers during the maternal-to-zygotic transition

After egg fertilization, an initially silent embryonic genome is transcriptionally activated during the maternal-to-zygotic transition. In zebrafish, maternal vertebrate pluripotency factors Nanog, Pou5f3 (OCT4 homolog), and Sox19b (SOX2 homolog) (NPS) play essential roles in orchestrating embryonic genome activation, acting as "pioneers" that open condensed chromatin and mediate acquisition of activating histone modifications. However, some embryonic gene transcription still occurs in the absence of these factors, suggesting the existence of other mechanisms regulating genome activation. To identify chromatin signatures of these unknown pathways, we profiled the histone modification landscape of zebrafish embryos using CUT&RUN. Our regulatory map revealed two subclasses of enhancers distinguished by presence or absence of H3K4me2. Enhancers lacking H3K4me2 tend to require NPS factors for de novo activation, while enhancers bearing H3K4me2 are epigenetically bookmarked by DNA hypomethylation to recapitulate gamete activity in the embryo, independent of NPS pioneering. Thus, parallel enhancer activation pathways combine to induce transcriptional reprogramming to pluripotency in the early embryo.

Notch-Dependent Expression of the Drosophila Hey Gene Is Supported by a Pair of Enhancers with Overlapping Activities

Drosophila Hey is a basic helix-loop-helix-orange (bHLH-O) protein with an important role in the establishment of distinct identities of postmitotic cells. We have previously identified Hey as a transcriptional target and effector of Notch signalling during the asymmetric division of neuronal progenitors, generating neurons of two types, and we have shown that Notch-dependent expression of Hey also marks a subpopulation of the newborn enteroendocrine (EE) cells in the midgut primordium of the embryo. Here, we investigate the transcriptional regulation of Hey in neuronal and intestinal tissues. We isolated two genomic regions upstream of the promoter (HeyUP) and in the second intron (HeyIN2) of the Hey gene, based on the presence of binding motifs for Su(H), the transcription factor that mediates Notch activity. We found that both regions can direct the overlapping expression patterns of reporter transgenes recapitulating endogenous Hey expression. Moreover, we showed that while HeyIN2 represents a Notch-dependent enhancer, HeyUP confers both Notch-dependent and independent transcriptional regulation. We induced mutations that removed the Su(H) binding motifs in either region and then studied the enhancer functionality in the respective Hey mutant lines. Our results provide direct evidence that although both enhancers support Notch-dependent regulation of the Hey gene, their role is redundant, as a Hey loss-of-function lethal phenotype is observed only after deletion of all their Su(H) binding motifs by CRISPR/Cas9.

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

BACKGROUND

Increasing evidence suggests that a substantial proportion of disease-associated mutations occur in enhancers, regions of non-coding DNA essential to gene regulation. Understanding the structures and mechanisms of the regulatory programs this variation affects can shed light on the apparatuses of human diseases.

RESULTS

We collect epigenetic and gene expression datasets from seven early time points during neural differentiation. Focusing on this model system, we construct networks of enhancer-promoter interactions, each at an individual stage of neural induction. These networks serve as the base for a rich series of analyses, through which we demonstrate their temporal dynamics and enrichment for various disease-associated variants. We apply the Girvan-Newman clustering algorithm to these networks to reveal biologically relevant substructures of regulation. Additionally, we demonstrate methods to validate predicted enhancer-promoter interactions using transcription factor overexpression and massively parallel reporter assays.

CONCLUSIONS

Our findings suggest a generalizable framework for exploring gene regulatory programs and their dynamics across developmental processes; this includes a comprehensive approach to studying the effects of disease-associated variation on transcriptional networks. The techniques applied to our networks have been published alongside our findings as a computational tool, E-P-INAnalyzer. Our procedure can be utilized across different cellular contexts and disorders.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6 weeks; n = 10) and mouse (post-natal day [P]0; n = 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

Gapped-kmer sequence modeling robustly identifies regulatory vocabularies and distal enhancers conserved between evolutionarily distant mammals

Gene regulatory elements drive complex biological phenomena and their mutations are associated with common human diseases. The impacts of human regulatory variants are often tested using model organisms such as mice. However, mapping human enhancers to conserved elements in mice remains a challenge, due to both rapid enhancer evolution and limitations of current computational methods. We analyze distal enhancers across 45 matched human/mouse cell/tissue pairs from a comprehensive dataset of DNase-seq experiments, and show that while cell-specific regulatory vocabulary is conserved, enhancers evolve more rapidly than promoters and CTCF binding sites. Enhancer conservation rates vary across cell types, in part explainable by tissue specific transposable element activity. We present an improved genome alignment algorithm using gapped-kmer features, called gkm-align, and make genome wide predictions for 1,401,803 orthologous regulatory elements. We show that gkm-align discovers 23,660 novel human/mouse conserved enhancers missed by previous algorithms, with strong evidence of conserved functional activity.

HyperCas12a enables highly-multiplexed epigenome editing screens

Interactions between multiple genes or cis-regulatory elements (CREs) underlie a wide range of biological processes in both health and disease. High-throughput screens using dCas9 fused to epigenome editing domains have allowed researchers to assess the impact of activation or repression of both coding and non-coding genomic regions on a phenotype of interest, but assessment of genetic interactions between those elements has been limited to pairs. Here, we combine a hyper-efficient version of Lachnospiraceae bacterium dCas12a (dHyperLbCas12a) with RNA Polymerase II expression of long CRISPR RNA (crRNA) arrays to enable efficient highly-multiplexed epigenome editing. We demonstrate that this system is compatible with several activation and repression domains, including the P300 histone acetyltransferase domain and SIN3A interacting domain (SID). We also show that the dCas12a platform can perform simultaneous activation and repression using a single crRNA array via co-expression of multiple dCas12a orthologues. Lastly, demonstrate that the dCas12a system is highly effective for high-throughput screens. We use dHyperLbCas12a-KRAB and a ~19,000-member barcoded library of crRNA arrays containing six crRNAs each to dissect the independent and combinatorial contributions of CREs to the dose-dependent control of gene expression at a glucocorticoid-responsive locus. The tools and methods introduced here create new possibilities for highly multiplexed control of gene expression in a wide variety of biological systems.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Chromatin profiling identifies chondrocyte-specific Sox9 enhancers important for skeletal development

The transcription factor SRY-related HMG box 9 (Sox9) is essential for chondrogenesis. Mutations in and around SOX9 cause campomelic dysplasia (CD) characterized by skeletal malformations. Although the function of Sox9 in this context is well studied, the mechanisms that regulate Sox9 expression in chondrocytes remain to be elucidated. Here, we have used genome-wide profiling to identify 2 Sox9 enhancers located in a proximal breakpoint cluster responsible for CD. Enhancer activity of E308 (located 308 kb 5' upstream) and E160 (located 160 kb 5' upstream) correlated with Sox9 expression levels, and both enhancers showed a synergistic effect in vitro. While single deletions in mice had no apparent effect, simultaneous deletion of both E308 and E160 caused a dwarf phenotype, concomitant with a reduction of Sox9 expression in chondrocytes. Moreover, bone morphogenetic protein 2-dependent chondrocyte differentiation of limb bud mesenchymal cells was severely attenuated in E308/E160 deletion mice. Finally, we found that an open chromatin region upstream of the Sox9 gene was reorganized in the E308/E160 deletion mice to partially compensate for the loss of E308 and E160. In conclusion, our findings reveal a mechanism of Sox9 gene regulation in chondrocytes that might aid in our understanding of the pathophysiology of skeletal disorders.

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

A promoter-dependent upstream activator augments CFTR expression in diverse epithelial cell types

The cystic fibrosis transmembrane conductance regulator (CFTR) gene encodes an anion-selective channel found in epithelial cell membranes. Mutations in CFTR cause cystic fibrosis (CF), an inherited disorder that impairs epithelial function in multiple organs. Most men with CF are infertile due to loss of intact genital ducts. Here we investigated a novel epididymis-selective cis-regulatory element (CRE), located within a peak of open chromatin at -9.5 kb 5' to the CFTR gene promoter. Activation of the -9.5 kb CRE alone by CRISPRa had no impact on CFTR gene expression. However, CRISPRa co-activation of the -9.5 kb CRE and the CFTR gene promoter in epididymis cells significantly augmented CFTR mRNA and protein expression when compared to promoter activation alone. This increase was accompanied by enhanced chromatin accessibility at both sites. Furthermore, the combined CRISPRa strategy activated CFTR expression in other epithelial cells that lack open chromatin at the -9.5 kb site and in which the locus is normally inactive. However, the -9.5 kb CRE does not function as a classical enhancer of the CFTR promoter in transient reporter gene assays. These data provide a novel mechanism for activating/augmenting CFTR expression, which may have therapeutic utility for mutations that perturb CFTR transcription.

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

Background

Increasing evidence suggests that a substantial proportion of disease-associated mutations occur in enhancers, regions of non-coding DNA essential to gene regulation. Understanding the structures and mechanisms of regulatory programs this variation affects can shed light on the apparatuses of human diseases.

Results

We collected epigenetic and gene expression datasets from seven early time points during neural differentiation. Focusing on this model system, we constructed networks of enhancer-promoter interactions, each at an individual stage of neural induction. These networks served as the base for a rich series of analyses, through which we demonstrated their temporal dynamics and enrichment for various disease-associated variants. We applied the Girvan-Newman clustering algorithm to these networks to reveal biologically relevant substructures of regulation. Additionally, we demonstrated methods to validate predicted enhancer-promoter interactions using transcription factor overexpression and massively parallel reporter assays.

Conclusions

Our findings suggest a generalizable framework for exploring gene regulatory programs and their dynamics across developmental processes. This includes a comprehensive approach to studying the effects of disease-associated variation on transcriptional networks. The techniques applied to our networks have been published alongside our findings as a computational tool, E-P-INAnalyzer. Our procedure can be utilized across different cellular contexts and disorders.

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

Transcription regulation by long non-coding RNAs: mechanisms and disease relevance

Long non-coding RNAs (lncRNAs) outnumber protein-coding transcripts, but their functions remain largely unknown. In this Review, we discuss the emerging roles of lncRNAs in the control of gene transcription. Some of the best characterized lncRNAs have essential transcription cis-regulatory functions that cannot be easily accomplished by DNA-interacting transcription factors, such as XIST, which controls X-chromosome inactivation, or imprinted lncRNAs that direct allele-specific repression. A growing number of lncRNA transcription units, including CHASERR, PVT1 and HASTER (also known as HNF1A-AS1) act as transcription-stabilizing elements that fine-tune the activity of dosage-sensitive genes that encode transcription factors. Genetic experiments have shown that defects in such transcription stabilizers often cause severe phenotypes. Other lncRNAs, such as lincRNA-p21 (also known as Trp53cor1) and Maenli (Gm29348) contribute to local activation of gene transcription, whereas distinct lncRNAs influence gene transcription in trans. We discuss findings of lncRNAs that elicit a function through either activation of their transcription, transcript elongation and processing or the lncRNA molecule itself. We also discuss emerging evidence of lncRNA involvement in human diseases, and their potential as therapeutic targets.

Multiple cis-regulatory elements control prox1a expression in distinct lymphatic vascular beds

During embryonic development, lymphatic endothelial cell (LEC) precursors are distinguished from blood endothelial cells by the expression of Prospero-related homeobox 1 (Prox1), which is essential for lymphatic vasculature formation in mouse and zebrafish. Prox1 expression initiation precedes LEC sprouting and migration, serving as the marker of specified LECs. Despite its crucial role in lymphatic development, Prox1 upstream regulation in LECs remains to be uncovered. SOX18 and COUP-TFII are thought to regulate Prox1 in mice by binding its promoter region. However, the specific regulation of Prox1 expression in LECs remains to be studied in detail. Here, we used evolutionary conservation and chromatin accessibility to identify enhancers located in the proximity of zebrafish prox1a active in developing LECs. We confirmed the functional role of the identified sequences through CRISPR/Cas9 mutagenesis of a lymphatic valve enhancer. The deletion of this region results in impaired valve morphology and function. Overall, our results reveal an intricate control of prox1a expression through a collection of enhancers. Ray-finned fish-specific distal enhancers drive pan-lymphatic expression, whereas vertebrate-conserved proximal enhancers refine expression in functionally distinct subsets of lymphatic endothelium.

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

Conserved and novel enhancers in the Aedes aegypti single-minded locus recapitulate embryonic ventral midline gene expression

Transcriptional cis-regulatory modules, e.g., enhancers, control the time and location of metazoan gene expression. While changes in enhancers can provide a powerful force for evolution, there is also significant deep conservation of enhancers for developmentally important genes, with function and sequence characteristics maintained over hundreds of millions of years of divergence. Not well understood, however, is how the overall regulatory composition of a locus evolves, with important outstanding questions such as how many enhancers are conserved vs. novel, and to what extent are the locations of conserved enhancers within a locus maintained? We begin here to address these questions with a comparison of the respective single-minded (sim) loci in the two dipteran species Drosophila melanogaster (fruit fly) and Aedes aegypti (mosquito). sim encodes a highly conserved transcription factor that mediates development of the arthropod embryonic ventral midline. We identify two enhancers in the A. aegypti sim locus and demonstrate that they function equivalently in both transgenic flies and transgenic mosquitoes. One A. aegypti enhancer is highly similar to known Drosophila counterparts in its activity, location, and autoregulatory capability. The other differs from any known Drosophila sim enhancers with a novel location, failure to autoregulate, and regulation of expression in a unique subset of midline cells. Our results suggest that the conserved pattern of sim expression in the two species is the result of both conserved and novel regulatory sequences. Further examination of this locus will help to illuminate how the overall regulatory landscape of a conserved developmental gene evolves.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Interpreting Cis-Regulatory Interactions from Large-Scale Deep Neural Networks for Genomics

The rise of large-scale, sequence-based deep neural networks (DNNs) for predicting gene expression has introduced challenges in their evaluation and interpretation. Current evaluations align DNN predictions with experimental perturbation assays, which provides insights into the generalization capabilities within the studied loci but offers a limited perspective of what drives their predictions. Moreover, existing model explainability tools focus mainly on motif analysis, which becomes complex when interpreting longer sequences. Here we introduce CREME, an in silico perturbation toolkit that interrogates large-scale DNNs to uncover rules of gene regulation that it learns. Using CREME, we investigate Enformer, a prominent DNN in gene expression prediction, revealing cis-regulatory elements (CREs) that directly enhance or silence target genes. We explore the intricate complexity of higher-order CRE interactions, the relationship between CRE distance from transcription start sites on gene expression, as well as the biochemical features of enhancers and silencers learned by Enformer. Moreover, we demonstrate the flexibility of CREME to efficiently uncover a higher-resolution view of functional sequence elements within CREs. This work demonstrates how CREME can be employed to translate the powerful predictions of large-scale DNNs to study open questions in gene regulation.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

The pleiotropic enhancer enh9 promotes cell proliferation and migration in non-small cell lung cancer via ERMP1 and PD-L1

Enhancers, cis-acting DNA elements for transcriptional regulation, are important regulators of cell identity and disease. However, of the hundreds of thousands of enhancers annotated in the human genome, only a few have been studied for their regulatory mechanisms and functions in cancer progression and therapeutic resistance. Here, we report the pleiotropy of one enhancer (named enh9) in both cell proliferation and migration in non-small cell lung cancer (NSCLC) cells. By integrating multi-genomic data, ERMP1 and PD-L1 were screened out as potential targets of enh9. CUT&Tag sequencing demonstrated that enh9 was involved in the genomic interactions between the transcription factor RELA and the promoters of ERMP1 and PD-L1. In addition, ERMP1 and PD-L1 were validated to be involved in cell proliferation and migration, respectively. Our study fully elucidated the function and transcriptional regulation mechanisms of enh9 in NSCLC. The exploration on enhancers is promising to provide new insights for cancer diagnosis and therapy.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Higher-order combinatorial chromatin perturbations by engineered CRISPR-Cas12a for functional genomics

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting 1-3 genomic sites per cell. To develop a tool for higher-order ( > 3) combinatorial targeting of genomic sites with CRISPRi in functional genomics screens, we engineered an Acidaminococcus Cas12a variant -- referred to as mul tiplexed transcriptional interference AsCas12a (multiAsCas12a). multiAsCas12a incorporates a key mutation, R1226A, motivated by the hypothesis of nicking-induced stabilization of the ribonucleoprotein:DNA complex for improving CRISPRi activity. multiAsCas12a significantly outperforms prior state-of-the-art Cas12a variants in combinatorial CRISPRi targeting using high-order multiplexed arrays of lentivirally transduced CRISPR RNAs (crRNA), including in high-throughput pooled screens using 6-plex crRNA array libraries. Using multiAsCas12a CRISPRi, we discover new enhancer elements and dissect the combinatorial function of cis-regulatory elements. These results instantiate a group testing framework for efficiently surveying potentially numerous combinations of chromatin perturbations for biological discovery and engineering.

Decoding the Epigenetics and Chromatin Loop Dynamics of Androgen Receptor-Mediated Transcription

Androgen receptor (AR)-mediated transcription plays a critical role in normal prostate development and prostate cancer growth. AR drives gene expression by binding to thousands of cis-regulatory elements (CRE) that loop to hundreds of target promoters. With multiple CREs interacting with a single promoter, it remains unclear how individual AR bound CREs contribute to gene expression. To characterize the involvement of these CREs, we investigated the AR-driven epigenetic and chromosomal chromatin looping changes. We collected a kinetic multi-omic dataset comprised of steady-state mRNA, chromatin accessibility, transcription factor binding, histone modifications, chromatin looping, and nascent RNA. Using an integrated regulatory network, we found that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we showed that binding of AR does not result in a substantial rewiring of chromatin loops, but instead increases the contact frequency of pre-existing loops to target promoters. Our results show that gene expression strongly correlates to the changes in contact frequency. We then proposed and experimentally validated an unbalanced multi-enhancer model where the impact on gene expression of AR-bound enhancers is heterogeneous, and is proportional to their contact frequency with target gene promoters. Overall, these findings provide new insight into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Transposable element-initiated enhancer-like elements generate the subgenome-biased spike specificity of polyploid wheat

Transposable elements (TEs) comprise ~85% of the common wheat genome, which are highly diverse among subgenomes, possibly contribute to polyploid plasticity, but the causality is only assumed. Here, by integrating data from gene expression cap analysis and epigenome profiling via hidden Markov model in common wheat, we detect a large proportion of enhancer-like elements (ELEs) derived from TEs producing nascent noncoding transcripts, namely ELE-RNAs, which are well indicative of the regulatory activity of ELEs. Quantifying ELE-RNA transcriptome across typical developmental stages reveals that TE-initiated ELE-RNAs are mainly from RLG_famc7.3 specifically expanded in subgenome A. Acquisition of spike-specific transcription factor binding likely confers spike-specific expression of RLG_famc7.3-initiated ELE-RNAs. Knockdown of RLG_famc7.3-initiated ELE-RNAs resulted in global downregulation of spike-specific genes and abnormal spike development. These findings link TE expansion to regulatory specificity and polyploid developmental plasticity, highlighting the functional impact of TE-driven regulatory innovation on polyploid evolution.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal discs; in contrast, small IDE-reporter transgenes are expressed mainly in the anterior compartment. We found that En binds to the IDEs and suggest that it may directly repress IDE function and modulate En expression levels. We identified two en IDEs, O and S. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from an endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data suggest that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic memory set by the embryonic enhancers. This epigenetic regulation is similar to that of the Ultrabithorax IDEs and we suggest that the activity of late-acting enhancers in other genes may be similarly regulated.

Unique activities of two overlapping PAX6 retinal enhancers

Enhancers play a critical role in development by precisely modulating spatial, temporal, and cell type-specific gene expression. Sequence variants in enhancers have been implicated in diseases; however, establishing the functional consequences of these variants is challenging because of a lack of understanding of precise cell types and developmental stages where the enhancers are normally active. PAX6 is the master regulator of eye development, with a regulatory landscape containing multiple enhancers driving the expression in the eye. Whether these enhancers perform additive, redundant or distinct functions is unknown. Here, we describe the precise cell types and regulatory activity of two PAX6 retinal enhancers, HS5 and NRE. Using a unique combination of live imaging and single-cell RNA sequencing in dual enhancer-reporter zebrafish embryos, we uncover differences in the spatiotemporal activity of these enhancers. Our results show that although overlapping, these enhancers have distinct activities in different cell types and therefore likely nonredundant functions. This work demonstrates that unique cell type-specific activities can be uncovered for apparently similar enhancers when investigated at high resolution in vivo.

Conserved enhancers control notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Iron Oxide Nanoparticles with and without Cobalt Functionalization Provoke Changes in the Transcription Profile via Epigenetic Modulation of Enhancer Activity

Despite the progress in the field of nanotoxicology, much about the cellular mechanisms that mediate the adverse effects of nanoparticles (NPs) and, in particular, the possible role of epigenetics in nanotoxicity, remains to be clarified. Therefore, we studied the changes occurring in the genome-wide distribution of H3K27ac, H3K4me1, H3K9me2, and H3K27me3 histone modifications and compared them with the transcriptome after exposing NIH3T3 cells to iron-based magnetic NPs (i.e., Fe2O3 and Fe2O3@Co NPs). We found that the transcription response is mainly due to changes in the genomic distribution of H3K27ac that can modulate the activity of enhancers. We propose that alteration of the epigenetic landscape is a key mechanism in defining the gene expression program changes resulting in nanotoxicity. With this approach, it is possible to construct a data set of genomic regions that could be useful for defining toxicity in a manner that is more comprehensive than what is possible with the present toxicology assays.

Transposable elements in mammalian chromatin organization

Transposable elements (TEs) are mobile DNA elements that comprise almost 50% of mammalian genomic sequence. TEs are capable of making additional copies of themselves that integrate into new positions in host genomes. This unique property has had an important impact on mammalian genome evolution and on the regulation of gene expression because TE-derived sequences can function as cis-regulatory elements such as enhancers, promoters and silencers. Now, advances in our ability to identify and characterize TEs have revealed that TE-derived sequences also regulate gene expression by both maintaining and shaping 3D genome architecture. Studies are revealing how TEs contribute raw sequence that can give rise to the structures that shape chromatin organization, and thus gene expression, allowing for species-specific genome innovation and evolutionary novelty.

Systematic identification of disease-causing promoter and untranslated region variants in 8,040 undiagnosed individuals with rare disease

Background

Both promoters and untranslated regions (UTRs) have critical regulatory roles, yet variants in these regions are largely excluded from clinical genetic testing due to difficulty in interpreting pathogenicity. The extent to which these regions may harbour diagnoses for individuals with rare disease is currently unknown.

Methods

We present a framework for the identification and annotation of potentially deleterious proximal promoter and UTR variants in known dominant disease genes. We use this framework to annotate de novo variants (DNVs) in 8,040 undiagnosed individuals in the Genomics England 100,000 genomes project, which were subject to strict region-based filtering, clinical review, and validation studies where possible. In addition, we performed region and variant annotation-based burden testing in 7,862 unrelated probands against matched unaffected controls.

Results

We prioritised eleven DNVs and identified an additional variant overlapping one of the eleven. Ten of these twelve variants (82%) are in genes that are a strong match to the individual's phenotype and six had not previously been identified. Through burden testing, we did not observe a significant enrichment of potentially deleterious promoter and/or UTR variants in individuals with rare disease collectively across any of our region or variant annotations.

Conclusions

Overall, we demonstrate the value of screening promoters and UTRs to uncover additional diagnoses for previously undiagnosed individuals with rare disease and provide a framework for doing so without dramatically increasing interpretation burden.

Spinal cord repair is modulated by the neurogenic factor Hb-egf under direction of a regeneration-associated enhancer

Unlike adult mammals, zebrafish regenerate spinal cord tissue and recover locomotor ability after a paralyzing injury. Here, we find that ependymal cells in zebrafish spinal cords produce the neurogenic factor Hb-egfa upon transection injury. Animals with hb-egfa mutations display defective swim capacity, axon crossing, and tissue bridging after spinal cord transection, associated with disrupted indicators of neuron production. Local recombinant human HB-EGF delivery alters ependymal cell cycling and tissue bridging, enhancing functional regeneration. Epigenetic profiling reveals a tissue regeneration enhancer element (TREE) linked to hb-egfa that directs gene expression in spinal cord injuries. Systemically delivered recombinant AAVs containing this zebrafish TREE target gene expression to crush injuries of neonatal, but not adult, murine spinal cords. Moreover, enhancer-based HB-EGF delivery by AAV administration improves axon densities after crush injury in neonatal cords. Our results identify Hb-egf as a neurogenic factor necessary for innate spinal cord regeneration and suggest strategies to improve spinal cord repair in mammals.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.