Shadow Enhancers Are Pervasive Features of Developmental Regulatory Networks

An Osteoblast-Specific Enhancer and Subenhancer Cooperatively Regulate Runx2 Expression in Chondrocytes

Runx2 is an essential transcription factor for osteoblast differentiation and chondrocyte maturation. The spatiotemporal expression of Runx2 is regulated by enhancers. We previously identified a 1.3 kb osteoblast-specific enhancer; however, mice with this deletion showed no phenotypes. A 0.8 kb conserved region detected near the 1.3 kb enhancer did not exhibit enhancer activity in reporter assays, whereas four tandem repeats of 452 bp (452 × 4) containing the most conserved region of 0.8 kb induced strong reporter activity in chondrocyte cell lines. However, chondrocytes of enhanced green fluorescent protein (EGFP) reporter mice using 452 × 4 did not express EGFP. When 452 × 4 was combined with the 1.3 kb enhancer, hypertrophic chondrocytes highly expressed EGFP. Moreover, the 0.8 kb region combined with the 1.3 kb enhancer induced EGFP expression in prehypertrophic and hypertrophic chondrocytes. The deletion of both the 1.3 kb enhancer and the 0.8 kb conserved region slightly reduced Runx2 expression in the limbs. However, neither homozygous nor heterozygous deletions in the Runx2+/- background showed phenotypes. The 0.8 kb conserved region itself lacked enhancer activity, but when combined with the 1.3 kb enhancer, EGFP expression was induced in chondrocytes with a similar expression pattern to Runx2. Therefore, the 0.8 kb conserved region has a novel function as a subenhancer.

Comprehensive dissection of cis-regulatory elements in a 2.8 Mb topologically associated domain in six human cancers

Cis-regulatory elements (CREs), such as enhancers and promoters, are fundamental regulators of gene expression and, across different cell types, the MYC locus utilizes a diverse regulatory architecture driven by multiple CREs. To better understand differences in CRE function, we perform pooled CRISPR inhibition (CRISPRi) screens to comprehensively probe the 2.8 Mb topologically-associated domain containing MYC in 6 human cancer cell lines with nucleotide resolution. We map 32 CREs where inhibition leads to changes in cell growth, including 8 that overlap previously identified enhancers. Targeting specific CREs decreases MYC expression by as much as 60%, and cell growth by as much as 50%. Using 3-D enhancer contact mapping, we find that these CREs almost always contact MYC but less than 10% of total MYC contacts impact growth when silenced, highlighting the utility of our approach to identify phenotypically-relevant CREs. We also detect an enrichment of lineage-specific transcription factors (TFs) at MYC CREs and, for some of these TFs, find a strong, tumor-specific correlation between TF and MYC expression not found in normal tissue. Taken together, these CREs represent systematically identified, functional regulatory regions and demonstrate how the same region of the human genome can give rise to complex, tissue-specific gene regulation.

Genetic coupling of enhancer activity and connectivity in gene expression control

Gene enhancers often form long-range contacts with promoters, but it remains unclear if the activity of enhancers and their chromosomal contacts are mediated by the same DNA sequences and recruited factors. Here, we study the effects of expression quantitative trait loci (eQTLs) on enhancer activity and promoter contacts in primary monocytes isolated from 34 male individuals. Using eQTL-Capture Hi-C and a Bayesian approach considering both intra- and inter-individual variation, we initially detect 19 eQTLs associated with enhancer-eGene promoter contacts, most of which also associate with enhancer accessibility and activity. Capitalising on these shared effects, we devise a multi-modality Bayesian strategy, identifying 629 "trimodal QTLs" jointly associated with enhancer accessibility, eGene promoter contact, and gene expression. Causal mediation analysis and CRISPR interference reveal causal relationships between these three modalities. Many detected QTLs overlap disease susceptibility loci and influence the predicted binding of myeloid transcription factors, including SPI1, GABPB and STAT3. Additionally, a variant associated with PCK2 promoter contact directly disrupts a CTCF binding motif and impacts promoter insulation from downstream enhancers. Jointly, our findings suggest an inherent genetic coupling of enhancer activity and connectivity in gene expression control relevant to human disease and highlight the regulatory role of genetically determined chromatin boundaries.

Two coacting shadow enhancers regulate twin of eyeless expression during early Drosophila development

The Drosophila PAX6 homolog twin of eyeless (toy) sits at the pinnacle of the genetic pathway controlling eye development, the retinal determination network. Expression of toy in the embryo is first detectable at cellular blastoderm stage 5 in an anterior-dorsal band in the presumptive procephalic neuroectoderm, which gives rise to the primordia of the visual system and brain. Although several maternal and gap transcription factors that generate positional information in the embryo have been implicated in controlling toy, the regulation of toy expression in the early embryo is currently not well characterized. In this study, we adopt an integrated experimental approach utilizing bioinformatics, molecular genetic testing of putative enhancers in transgenic reporter gene assays and quantitative analysis of expression patterns in the early embryo, to identify 2 novel coacting enhancers at the toy gene. In addition, we apply mathematical modeling to dissect the regulatory landscape for toy. We demonstrate that relatively simple thermodynamic-based models, incorporating only 5 TF binding sites, can accurately predict gene expression from the 2 coacting enhancers and that the HUNCHBACK TF plays a critical regulatory role through a dual-modality function as an activator and repressor. Our analysis also reveals that the molecular architecture of the 2 enhancers is very different, indicating that the underlying regulatory logic they employ is distinct.

Decoding the epigenetics and chromatin loop dynamics of androgen receptor-mediated transcription

Androgen receptor (AR)-mediated transcription plays a critical role in 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 investigate the AR-driven epigenetic and chromosomal chromatin looping changes by generating 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 find that AR binding induces sequential changes in the epigenetic features at CREs, independent of gene expression. Further, we show 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 propose and experimentally validate 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 insights into AR-mediated gene expression upon acute androgen simulation and develop a mechanistic framework to investigate nuclear receptor mediated perturbations.

A Foxf1-Wnt-Nr2f1 cascade promotes atrial cardiomyocyte differentiation in zebrafish

Nr2f transcription factors (TFs) are conserved regulators of vertebrate atrial cardiomyocyte (AC) differentiation. However, little is known about the mechanisms directing Nr2f expression in ACs. Here, we identified a conserved enhancer 3' to the nr2f1a locus, which we call 3'reg1-nr2f1a (3'reg1), that can promote Nr2f1a expression in ACs. Sequence analysis of the enhancer identified putative Lef/Tcf and Foxf TF binding sites. Mutation of the Lef/Tcf sites within the 3'reg1 reporter, knockdown of Tcf7l1a, and manipulation of canonical Wnt signaling support that Tcf7l1a is derepressed via Wnt signaling to activate the transgenic enhancer and promote AC differentiation. Similarly, mutation of the Foxf binding sites in the 3'reg1 reporter, coupled with gain- and loss-of-function analysis supported that Foxf1 promotes expression of the enhancer and AC differentiation. Functionally, we find that Wnt signaling acts downstream of Foxf1 to promote expression of the 3'reg1 reporter within ACs and, importantly, both Foxf1 and Wnt signaling require Nr2f1a to promote a surplus of differentiated ACs. CRISPR-mediated deletion of the endogenous 3'reg1 abrogates the ability of Foxf1 and Wnt signaling to produce surplus ACs in zebrafish embryos. Together, our data support that downstream members of a conserved regulatory network involving Wnt signaling and Foxf1 function on a nr2f1a enhancer to promote AC differentiation in the zebrafish heart.

The case against simplistic genetic explanations of evolution

Humans are curious to understand the causes of traits that distinguish us from other animals and that distinguish vastly different species from one another. We also have a proclivity for simple stories and sometimes tend toward seeking and accepting simple genetic explanations for large evolutionary shifts, even to a single gene. Here, I reveal how a biased expectation of mechanistic simplicity threads through the long history of evolutionary and developmental genetics. I argue, however, that expecting a simple mechanism threatens a deeper understanding of evolution, and I define the limitations for interpreting experimental evidence in evolutionary developmental genetics.

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.

Two C. elegans DM domain proteins, DMD-3 and MAB-3, function in late stages of male somatic gonad development

To bring about sexual dimorphism in form, information from the sex determination pathway must trigger sex-specific modifications in developmental programs. DM-domain encoding genes have been found to be involved in sex determination in a multitude of animals, often at the level of male somatic gonad formation. Here we report our findings that the DM-domain transcription factors MAB-3 and DMD-3 function together in multiple steps during the late stages of C. elegans male somatic gonad development. Both mab-3 and dmd-3 are expressed in the linker cell and hindgut of L4 males and dmd-3 is also expressed in presumptive vas deferens cells. Furthermore, dmd-3, but not mab-3, expression in the linker cell is downstream of nhr-67, a nuclear hormone receptor that was previously shown to control late stages of linker cell migration. In mab-3; dmd-3 double mutant males, the last stage of linker cell migration is partially defective, resulting in aberrant linker cell shapes and often a failure of the linker cell to complete its migration to the hindgut. When mab-3; dmd-3 double mutant linker cells do complete their migration, they fail to be engulfed by the hindgut, indicating that dmd-3 and mab-3 activity are essential for this process. Furthermore, linker cell death and clearance are delayed in mab-3; dmd-3 double mutants, resulting in the linker cell persisting into adulthood. Finally, DMD-3 and MAB-3 function to activate expression of the bZIP transcription factor encoding gene zip-5 and downregulate the expression of the zinc metalloprotease ZMP-1 in the linker cell. Taken together, these results demonstrate a requirement for DM-domain transcription factors in controlling C. elegans male gonad formation, supporting the notion that the earliest DM-domain genes were involved in male somatic gonad development in the last common ancestor of the bilaterians.

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.

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Transcriptomic Response to Pyrethroid Treatment in Closely Related Bed Bug Strains Varying in Resistance

The common bed bug, Cimex lectularius, is one of the main human parasites. The world-wide resurgence of this pest is mainly due to globalization, and the spread of insecticide resistance. A few studies have compared the transcriptomes of susceptible and resistant strains; however, these studies usually relied on strains originating from distant locations, possibly explaining their extended candidate gene lists. Here, we compared the transcriptomes of 2 strains originating from the same location and showing low overall genetic differentiation (FST=0.018) but varying in their susceptibility to pyrethroids, before and after insecticide exposure. In sharp contrast with previous studies, only 24 genes showing constitutive differential expression between the strains were identified. Interestingly, most of the genes with increased expression in the resistant strain encoded cuticular proteins. However, those changes were not associated with significant difference in cuticular thickness, suggesting that they might be involved in qualitative changes in the cuticle. In contrast, insecticide exposure induced the expression of a multitude of genes, mostly involved in detoxification. Finally, our set of transcriptome candidate loci showed little overlap with a set of loci strongly genetically differentiated in a previous study using the same strains. Several hypothesis explaining this discrepancy are discussed.

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.

A Foxf1-Wnt-Nr2f1 cascade promotes atrial cardiomyocyte differentiation in zebrafish

Nr2f transcription factors (TFs) are conserved regulators of vertebrate atrial cardiomyocyte (AC) differentiation. However, little is known about the mechanisms directing Nr2f expression in ACs. Here, we identified a conserved enhancer 3' to the nr2f1a locus, which we call 3'reg1-nr2f1a (3'reg1), that can promote Nr2f1a expression in ACs. Sequence analysis of the enhancer identified putative Lef/Tcf and Foxf TF binding sites. Mutation of the Lef/Tcf sites within the 3'reg1 reporter, knockdown of Tcf7l1a, and manipulation of canonical Wnt signaling support that Tcf7l1a is derepressed via Wnt signaling to activate the transgenic enhancer and promote AC differentiation. Similarly, mutation of the Foxf binding sites in the 3'reg1 reporter, coupled with gain- and loss-of-function analysis supported that Foxf1 promotes expression of the enhancer and AC differentiation. Functionally, we find that Wnt signaling acts downstream of Foxf1 to promote expression of the 3'reg1 reporter within ACs and, importantly, both Foxf1 and Wnt signaling require Nr2f1a to promote a surplus of differentiated ACs. CRISPR-mediated deletion of the endogenous 3'reg1 abrogates the ability of Foxf1 and Wnt signaling to produce surplus ACs in zebrafish embryos. Together, our data support that downstream members of a conserved regulatory network involving Wnt signaling and Foxf1 function on a nr2f1a enhancer to promote AC differentiation in the zebrafish heart.

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.

Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs.

Affinity-optimizing enhancer variants disrupt development

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

Deep molecular learning of transcriptional control of a synthetic CRE enhancer and its variants

Massively parallel reporter assay measures transcriptional activities of various cis-regulatory modules (CRMs) in a single experiment. We developed a thermodynamic computational model framework that calculates quantitative levels of gene expression directly from regulatory DNA sequences. Using the framework, we investigated the molecular mechanisms of cis-regulatory mutations of a synthetic enhancer that cause abnormal gene expression. We found that, in a human cell line, competitive binding between family transcription factors (TFs) with slightly different binding preferences significantly increases the accuracy of recapitulating the transcriptional effects of thousands of single- or multi-mutations. We also discovered that even if various harmful mutations occurred in an activator binding site, CRM could stably maintain or even increase gene expression through a certain form of competitive binding between family TFs. These findings enhance understanding the effect of SNPs and indels on CRMs and would help building robust custom-designed CRMs for biologics production and gene therapy.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

Cross-tissue patterns of DNA hypomethylation reveal genetically distinct histories of cell development

BACKGROUND

Establishment of DNA methylation (DNAme) patterns is essential for balanced multi-lineage cellular differentiation, but exactly how these patterns drive cellular phenotypes is unclear. While > 80% of CpG sites are stably methylated, tens of thousands of discrete CpG loci form hypomethylated regions (HMRs). Because they lack DNAme, HMRs are considered transcriptionally permissive, but not all HMRs actively regulate genes. Unlike promoter HMRs, a subset of non-coding HMRs is cell type-specific and enriched for tissue-specific gene regulatory functions. Our data further argues not only that HMR establishment is an important step in enforcing cell identity, but also that cross-cell type and spatial HMR patterns are functionally informative of gene regulation.

RESULTS

To understand the significance of non-coding HMRs, we systematically dissected HMR patterns across diverse human cell types and developmental timepoints, including embryonic, fetal, and adult tissues. Unsupervised clustering of 126,104 distinct HMRs revealed that levels of HMR specificity reflects a developmental hierarchy supported by enrichment of stage-specific transcription factors and gene ontologies. Using a pseudo-time course of development from embryonic stem cells to adult stem and mature hematopoietic cells, we find that most HMRs observed in differentiated cells (~ 60%) are established at early developmental stages and accumulate as development progresses. HMRs that arise during differentiation frequently (~ 35%) establish near existing HMRs (≤ 6 kb away), leading to the formation of HMR clusters associated with stronger enhancer activity. Using SNP-based partitioned heritability from GWAS summary statistics across diverse traits and clinical lab values, we discovered that genetic contribution to trait heritability is enriched within HMRs. Moreover, the contribution of heritability to cell-relevant traits increases with both increasing HMR specificity and HMR clustering, supporting the role of distinct HMR subsets in regulating normal cell function.

CONCLUSIONS

Our results demonstrate that the entire HMR repertoire within a cell-type, rather than just the cell type-specific HMRs, stores information that is key to understanding and predicting cellular phenotypes. Ultimately, these data provide novel insights into how DNA hypo-methylation provides genetically distinct historical records of a cell's journey through development, highlighting HMRs as functionally distinct from other epigenomic annotations.

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.

The interaction between enhancer variants and environmental factors as an overlooked aetiological paradigm in human complex disease

The interactions between genetic and environmental risk factors contribute to the aetiology of complex human diseases. Genome-wide association studies (GWAS) have revealed that most of the genetic variants associated with complex diseases are located in the non-coding part of the genome, preferentially within enhancers. Enhancers are distal cis-regulatory elements composed of clusters of transcription factors binding sites that positively regulate the expression of their target genes. The generation of genome-wide maps for histone marks (e.g., H3K27ac), chromatin accessibility and transcription factor and coactivator (e.g., p300) binding profiles have enabled the identification of enhancers across many human cell types and tissues. Nonetheless, the functional and pathological consequences of the majority of disease-associated genetic variants located within enhancers seem to be rather minor under normal conditions, thus questioning their medical relevance. Here we propose that, due to the prevalence of enhancer redundancy, the pathological effects of many disease-associated non-coding genetic variants might be preferentially (or even only) manifested under environmental stress.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer

BACKGROUND

Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance.

METHODS

Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis.

RESULTS

We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care.

CONCLUSION

In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Accelerated Evolution Analysis Uncovers PKNOX2 as a Key Transcription Factor in the Mammalian Cochlea

The genetic bases underlying the evolution of morphological and functional innovations of the mammalian inner ear are poorly understood. Gene regulatory regions are thought to play an important role in the evolution of form and function. To uncover crucial hearing genes whose regulatory machinery evolved specifically in mammalian lineages, we mapped accelerated noncoding elements (ANCEs) in inner ear transcription factor (TF) genes and found that PKNOX2 harbors the largest number of ANCEs within its transcriptional unit. Using reporter gene expression assays in transgenic zebrafish, we determined that four PKNOX2-ANCEs drive differential expression patterns when compared with ortholog sequences from close outgroup species. Because the functional role of PKNOX2 in cochlear hair cells has not been previously investigated, we decided to study Pknox2 null mice generated by CRISPR/Cas9 technology. We found that Pknox2-/- mice exhibit reduced distortion product otoacoustic emissions (DPOAEs) and auditory brainstem response (ABR) thresholds at high frequencies together with an increase in peak 1 amplitude, consistent with a higher number of inner hair cells (IHCs)-auditory nerve synapsis observed at the cochlear basal region. A comparative cochlear transcriptomic analysis of Pknox2-/- and Pknox2+/+ mice revealed that key auditory genes are under Pknox2 control. Hence, we report that PKNOX2 plays a critical role in cochlear sensitivity at higher frequencies and that its transcriptional regulation underwent lineage-specific evolution in mammals. Our results provide novel insights about the contribution of PKNOX2 to normal auditory function and to the evolution of high-frequency hearing in mammals.

Identification and characterization of transcribed enhancers during cerebellar development through enhancer RNA analysis

BACKGROUND

The development of the brain requires precise coordination of molecular processes across many cell-types. Underpinning these events are gene expression programs which require intricate regulation by non-coding regulatory sequences known as enhancers. In the context of the developing brain, transcribed enhancers (TEs) regulate temporally-specific expression of genes critical for cell identity and differentiation. Transcription of non-coding RNAs at active enhancer sequences, known as enhancer RNAs (eRNAs), is tightly associated with enhancer activity and has been correlated with target gene expression. TEs have been characterized in a multitude of developing tissues, however their regulatory role has yet to be described in the context of embryonic and early postnatal brain development. In this study, eRNA transcription was analyzed to identify TEs active during cerebellar development, as a proxy for the developing brain. Cap Analysis of Gene Expression followed by sequencing (CAGE-seq) was conducted at 12 stages throughout embryonic and early postnatal cerebellar development.

RESULTS

Temporal analysis of eRNA transcription identified clusters of TEs that peak in activity during either embryonic or postnatal times, highlighting their importance for temporally specific developmental events. Functional analysis of putative target genes identified molecular mechanisms under TE regulation revealing that TEs regulate genes involved in biological processes specific to neurons. We validate enhancer activity using in situ hybridization of eRNA expression from TEs predicted to regulate Nfib, a gene critical for cerebellar granule cell differentiation.

CONCLUSIONS

The results of this analysis provide a valuable dataset for the identification of cerebellar enhancers and provide insight into the molecular mechanisms critical for brain development under TE regulation. This dataset is shared with the community through an online resource ( https://goldowitzlab.shinyapps.io/trans-enh-app/ ).

Long non-coding RNAs: definitions, functions, challenges and recommendations

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

Decoding enhancer complexity with machine learning and high-throughput discovery

Enhancers are genomic DNA elements controlling spatiotemporal gene expression. Their flexible organization and functional redundancies make deciphering their sequence-function relationships challenging. This article provides an overview of the current understanding of enhancer organization and evolution, with an emphasis on factors that influence these relationships. Technological advancements, particularly in machine learning and synthetic biology, are discussed in light of how they provide new ways to understand this complexity. Exciting opportunities lie ahead as we continue to unravel the intricacies of enhancer function.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development. 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 Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

Transcriptome sequencing and metabolite analysis reveal the single and combined effects of microplastics and di-(2-ethylhexyl) phthalate on Peneaus vannamei

Due to the rising usage of plastics, plastic debris are present throughout marine ecosystems and detrimentally affects marine biota. Additionally, plastics likely result in elusive toxicity effects due to addition of plasticizers. The aim of the present study was to reveal the potential effects and mechanism of microplastics (MPs), di-(2-ethylhexyl) phthalate (DEHP) and copollution of MPs and DEHP (MPs-DEHP) on Peneaus vannamei (P. vannamei) juveniles regarding oxidative stress, transcriptomics and metabolomics. MPs, DEHP and MPs-DEHP significantly induced the activities of superoxide dismutase (SOD) and catalase (CAT); MPs and DEHP have an antagonistic effect for malondialdehyde (MDA); suggesting that disorders of the antioxidant defence systems. 13, 133 and 58 differentially expressed genes and 21, 82 and 39 differentially expressed metabolites were responsible for the distinction of MPs, DEHP and MPs-DEHP groups, respectively. The combination of transcriptomic and metabolomic analyses showed that MPs, DEHP and MPs-DEHP exposure disturbed amino acid and lipid metabolism, and further induced inflammatory responses and dysfunction of purine metabolism. Furthermore, the presence of MPs might alleviate the biotoxicity of DEHP in P. vannamei. These findings provide new insights into the single and combined toxicological effects of MPs and additives for marine biota.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

Mechanisms of robustness in gene regulatory networks involved in neural development

The functions of living organisms are affected by different kinds of perturbation, both internal and external, which in many cases have functional effects and phenotypic impact. The effects of these perturbations become particularly relevant for multicellular organisms with complex body patterns and cell type heterogeneity, where transcriptional programs controlled by gene regulatory networks determine, for example, the cell fate during embryonic development. Therefore, an essential aspect of development in these organisms is the ability to maintain the functionality of their genetic developmental programs even in the presence of genetic variation, changing environmental conditions and biochemical noise, a property commonly termed robustness. We discuss the implication of different molecular mechanisms of robustness involved in neurodevelopment, which is characterized by the interplay of many developmental programs at a molecular, cellular and systemic level. We specifically focus on processes affecting the function of gene regulatory networks, encompassing transcriptional regulatory elements and post-transcriptional processes such as miRNA-based regulation, but also higher order regulatory organization, such as gene network topology. We also present cases where impairment of robustness mechanisms can be associated with neurodevelopmental disorders, as well as reasons why understanding these mechanisms should represent an important part of the study of gene regulatory networks driving neural development.

Multi-modular structure of the gene regulatory network for specification and commitment of murine T cells

T cells develop from multipotent progenitors by a gradual process dependent on intrathymic Notch signaling and coupled with extensive proliferation. The stages leading them to T-cell lineage commitment are well characterized by single-cell and bulk RNA analyses of sorted populations and by direct measurements of precursor-product relationships. This process depends not only on Notch signaling but also on multiple transcription factors, some associated with stemness and multipotency, some with alternative lineages, and others associated with T-cell fate. These factors interact in opposing or semi-independent T cell gene regulatory network (GRN) subcircuits that are increasingly well defined. A newly comprehensive picture of this network has emerged. Importantly, because key factors in the GRN can bind to markedly different genomic sites at one stage than they do at other stages, the genes they significantly regulate are also stage-specific. Global transcriptome analyses of perturbations have revealed an underlying modular structure to the T-cell commitment GRN, separating decisions to lose "stem-ness" from decisions to block alternative fates. Finally, the updated network sheds light on the intimate relationship between the T-cell program, which depends on the thymus, and the innate lymphoid cell (ILC) program, which does not.

Hooked Up from a Distance: Charting Genome-Wide Long-Range Interaction Maps in Neural Cells Chromatin to Identify Novel Candidate Genes for Neurodevelopmental Disorders

DNA sequence variants (single nucleotide polymorphisms or variants, SNPs/SNVs; copy number variants, CNVs) associated to neurodevelopmental disorders (NDD) and traits often map on putative transcriptional regulatory elements, including, in particular, enhancers. However, the genes controlled by these enhancers remain poorly defined. Traditionally, the activity of a given enhancer, and the effect of its possible alteration associated to the sequence variants, has been thought to influence the nearest gene promoter. However, the obtainment of genome-wide long-range interaction maps in neural cells chromatin challenged this view, showing that a given enhancer is very frequently not connected to the nearest promoter, but to a more distant one, skipping genes in between. In this Perspective, we review some recent papers, who generated long-range interaction maps (by HiC, RNApolII ChIA-PET, Capture-HiC, or PLACseq), and overlapped the identified long-range interacting DNA segments with DNA sequence variants associated to NDD (such as schizophrenia, bipolar disorder and autism) and traits (intelligence). This strategy allowed to attribute the function of enhancers, hosting the NDD-related sequence variants, to a connected gene promoter lying far away on the linear chromosome map. Some of these enhancer-connected genes had indeed been already identified as contributive to the diseases, by the identification of mutations within the gene's protein-coding regions (exons), validating the approach. Significantly, however, the connected genes also include many genes that were not previously found mutated in their exons, pointing to novel candidate contributors to NDD and traits. Thus, long-range interaction maps, in combination with DNA variants detected in association with NDD, can be used as "pointers" to identify novel candidate disease-relevant genes. Functional manipulation of the long-range interaction network involving enhancers and promoters by CRISPR-Cas9-based approaches is beginning to probe for the functional significance of the identified interactions, and the enhancers and the genes involved, improving our understanding of neural development and its pathology.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

Functional characterization of enhancer activity during a long terminal repeat's evolution

Many transposable elements (TEs) contain transcription factor binding sites and are implicated as potential regulatory elements. However, TEs are rarely functionally tested for regulatory activity, which in turn limits our understanding of how TE regulatory activity has evolved. We systematically tested the human LTR18A subfamily for regulatory activity using massively parallel reporter assay (MPRA) and found AP-1- and CEBP-related binding motifs as drivers of enhancer activity. Functional analysis of evolutionarily reconstructed ancestral sequences revealed that LTR18A elements have generally lost regulatory activity over time through sequence changes, with the largest effects occurring owing to mutations in the AP-1 and CEBP motifs. We observed that the two motifs are conserved at higher rates than expected based on neutral evolution. Finally, we identified LTR18A elements as potential enhancers in the human genome, primarily in epithelial cells. Together, our results provide a model for the origin, evolution, and co-option of TE-derived regulatory elements.

A Bcl6 Intronic Element Regulates T Follicular Helper Cell Differentiation

In response to an intracellular infectious agent, the immune system produces a specific cellular response as well as a T cell–dependent Ab response. Precursor T cells differentiate into effector T cells, including Th1 cells, and T follicular helper (TFH) cells. The latter cooperate with B cells to form germinal centers and induce the formation of Ab-forming plasmacytes. One major focal point for control of T cell differentiation is the transcription factor BCL6. In this study, we demonstrated that the Bcl6 gene is regulated by FOXO1-binding, cis-acting sequences located in a highly conserved region of the first Bcl6 intron. In both mouse and human T cells, deletion of the tandem FOXO1 binding sites increased the expression of BCL6 and enhanced the proportion of TFH cells. These results reveal a fundamental control point for cellular versus humoral immunity.

DNA Methylation in Regulatory T Cell Differentiation and Function: Challenges and Opportunities

As a bona fide epigenetic marker, DNA methylation has been linked to the differentiation and function of regulatory T (Treg) cells, a subset of CD4 T cells that play an essential role in maintaining immune homeostasis and suppressing autoimmunity and antitumor immune response. DNA methylation undergoes dynamic regulation involving maintenance of preexisting patterns, passive and active demethylation, and de novo methylation. Scattered evidence suggests that these processes control different stages of Treg cell lifespan ranging from lineage induction to cell fate maintenance, suppression of effector T cells and innate immune cells, and transdifferentiation. Despite significant progress, it remains to be fully explored how differential DNA methylation regulates Treg cell fate and immunological function. Here, we review recent progress and discuss the questions and challenges for further understanding the immunological roles and mechanisms of dynamic DNA methylation in controlling Treg cell differentiation and function. We also explore the opportunities that these processes offer to manipulate Treg cell suppressive function for therapeutic purposes by targeting DNA methylation.

Pancreas agenesis mutations disrupt a lead enhancer controlling a developmental enhancer cluster

Sequence variants in cis-acting enhancers are important for polygenic disease, but their role in Mendelian disease is poorly understood. Redundancy between enhancers that regulate the same gene is thought to mitigate the pathogenic impact of enhancer mutations. Recent findings, however, have shown that loss-of-function mutations in a single enhancer near PTF1A cause pancreas agenesis and neonatal diabetes. Using mouse and human genetic models, we show that this enhancer activates an entire PTF1A enhancer cluster in early pancreatic multipotent progenitors. This leading role, therefore, precludes functional redundancy. We further demonstrate that transient expression of PTF1A in multipotent progenitors sets in motion an epigenetic cascade that is required for duct and endocrine differentiation. These findings shed insights into the genome regulatory mechanisms that drive pancreas differentiation. Furthermore, they reveal an enhancer that acts as a regulatory master key and is thus vulnerable to pathogenic loss-of-function mutations.

•

The pancreas agenesis enhancer (EnhP) activates PTF1A in early pancreatic progenitors

•

EnhP also activates other progenitor PTF1A enhancers

•

This master key function explains why EnhP is vulnerable to loss-of-function mutations

•

Transient PTF1A expression in progenitors controls pancreas growth and endocrinogenesis

Multi-layered transcriptional control of cranial neural crest development

The neural crest (NC) is an emblematic population of embryonic stem-like cells with remarkable migratory ability. These distinctive attributes have inspired the curiosity of developmental biologists for over 150 years, however only recently the regulatory mechanisms controlling the complex features of the NC have started to become elucidated at genomic scales. Regulatory control of NC development is achieved through combinatorial transcription factor binding and recruitment of associated transcriptional complexes to distal cis-regulatory elements. Together, they regulate when, where and to what extent transcriptional programmes are actively deployed, ultimately shaping ontogenetic processes. Here, we discuss how transcriptional networks control NC ontogeny, with a special emphasis on the molecular mechanisms underlying specification of the cephalic NC. We also cover emerging properties of transcriptional regulation revealed in diverse developmental systems, such as the role of three-dimensional conformation of chromatin, and how they are involved in the regulation of NC ontogeny. Finally, we highlight how advances in deciphering the NC transcriptional network have afforded new insights into the molecular basis of human diseases.