The regulatory landscapes of developmental genes

Two redundant transcription factor binding sites in a single enhancer are essential for mammalian sex determination

Male development in mammals depends on the activity of the two SOX gene: Sry and Sox9, in the embryonic testis. As deletion of Enhancer 13 (Enh13) of the Sox9 gene results in XY male-to-female sex reversal, we explored the critical elements necessary for its function and hence, for testis and male development. Here, we demonstrate that while microdeletions of individual transcription factor binding sites (TFBS) in Enh13 lead to normal testicular development, combined microdeletions of just two SRY/SOX binding motifs can alone fully abolish Enh13 activity leading to XY male-to-female sex reversal. This suggests that for proper male development to occur, these few nucleotides of non-coding DNA must be intact. Interestingly, we show that depending on the nature of these TFBS mutations, dramatically different phenotypic outcomes can occur, providing a molecular explanation for the distinct clinical outcomes observed in patients harboring different variants in the same enhancer.

Genomic context sensitizes regulatory elements to genetic disruption

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Establishment and maintenance of random monoallelic expression

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

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Genomic context sensitizes regulatory elements to genetic disruption

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

Betaglycan promoter activity is differentially regulated during myogenesis in zebrafish embryo somites

BACKGROUND

Betaglycan, also known as the TGFβ type III receptor (Tgfbr3), is a co-receptor that modulates TGFβ family signaling. Tgfbr3 is upregulated during C2C12 myoblast differentiation and expressed in mouse embryos myocytes.

RESULTS

To investigate tgfbr3 transcriptional regulation during zebrafish embryonic myogenesis, we cloned a 3.2 kb promoter fragment that drives reporter transcription during C2C12 myoblasts differentiation and in the Tg(tgfbr3:mCherry) transgenic zebrafish. We detect tgfbr3 protein and mCherry expression in the adaxial cells concomitantly with the onset of their radial migration to become slow-twitch muscle fibers in the Tg(tgfbr3:mCherry). Remarkably, this expression displays a measurable antero-posterior somitic gradient expression.

CONCLUSIONS

tgfbr3 is transcriptionally regulated during somitic muscle development in zebrafish with an antero-posterior gradient expression that preferentially marks the adaxial cells and their descendants.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

Developmental mechanisms underlying the evolution of human cortical circuits

The brain of modern humans has evolved remarkable computational abilities that enable higher cognitive functions. These capacities are tightly linked to an increase in the size and connectivity of the cerebral cortex, which is thought to have resulted from evolutionary changes in the mechanisms of cortical development. Convergent progress in evolutionary genomics, developmental biology and neuroscience has recently enabled the identification of genomic changes that act as human-specific modifiers of cortical development. These modifiers influence most aspects of corticogenesis, from the timing and complexity of cortical neurogenesis to synaptogenesis and the assembly of cortical circuits. Mutations of human-specific genetic modifiers of corticogenesis have started to be linked to neurodevelopmental disorders, providing evidence for their physiological relevance and suggesting potential relationships between the evolution of the human brain and its sensitivity to specific diseases.

Characterisation of the role and regulation of Ultrabithorax in sculpting fine-scale leg morphology

Hox genes are expressed during embryogenesis and determine the regional identity of animal bodies along the antero-posterior axis. However, they also function post-embryonically to sculpt fine-scale morphology. To better understand how Hox genes are integrated into post-embryonic gene regulatory networks, we further analysed the role and regulation of Ultrabithorax (Ubx) during leg development in Drosophila melanogaster. Ubx regulates several aspects of bristle and trichome patterning on the femurs of the second (T2) and third (T3) leg pairs. We found that repression of trichomes in the proximal posterior region of the T2 femur by Ubx is likely mediated by activation of the expression of microRNA-92a and microRNA-92b by this Hox protein. Furthermore, we identified a novel enhancer of Ubx that recapitulates the temporal and regional activity of this gene in T2 and T3 legs. We then used transcription factor (TF) binding motif analysis in regions of accessible chromatin in T2 leg cells to predict and functionally test TFs that may regulate the Ubx leg enhancer. We also tested the role of the Ubx co-factors Homothorax (Hth) and Extradenticle (Exd) in T2 and T3 femurs. We found several TFs that may act upstream or in concert with Ubx to modulate trichome patterning along the proximo-distal axis of developing femurs and that the repression of trichomes also requires Hth and Exd. Taken together our results provide insights into how Ubx is integrated into a post-embryonic gene regulatory network to determine fine-scale leg morphology.

Noncoding Mutations in a Thyroid Hormone Receptor Gene That Impair Cone Photoreceptor Function

The function of a hormone receptor requires mechanisms to control precisely where, when, and at what level the receptor gene is expressed. An intriguing case concerns the selective induction of thyroid hormone receptor β2 (TRβ2), encoded by Thrb, in the pituitary and also in cone photoreceptors, in which it critically regulates expression of the opsin photopigments that mediate color vision. Here, we investigate the physiological significance of a candidate enhancer for induction of TRβ2 by mutagenesis of a conserved intron region in its natural context in the endogenous Thrb gene in mice. Mutation of e-box sites for bHLH (basic-helix-loop-helix) transcription factors preferentially impairs TRβ2 expression in cones whereas mutation of nearby sequences preferentially impairs expression in pituitary. A deletion encompassing all sites impairs expression in both tissues, indicating bifunctional activity. In cones, the e-box mutations disrupt chromatin acetylation, blunt the developmental induction of TRβ2, and ultimately impair cone opsin expression and sensitivity to longer wavelengths of light. These results demonstrate the necessity of studying an enhancer in its natural chromosomal context for defining biological relevance and reveal surprisingly critical nuances of level and timing of enhancer function. Our findings illustrate the influence of noncoding sequences over thyroid hormone functions.

Context-dependent enhancer function revealed by targeted inter-TAD relocation

The expression of some genes depends on large, adjacent regions of the genome that contain multiple enhancers. These regulatory landscapes frequently align with Topologically Associating Domains (TADs), where they integrate the function of multiple similar enhancers to produce a global, TAD-specific regulation. We asked if an individual enhancer could overcome the influence of one of these landscapes, to drive gene transcription. To test this, we transferred an enhancer from its native location, into a nearby TAD with a related yet different functional specificity. We used the biphasic regulation of Hoxd genes during limb development as a paradigm. These genes are first activated in proximal limb cells by enhancers located in one TAD, which is then silenced when the neighboring TAD activates its enhancers in distal limb cells. We transferred a distal limb enhancer into the proximal limb TAD and found that its new context suppresses its normal distal specificity, even though it is bound by HOX13 transcription factors, which are responsible for the distal activity. This activity can be rescued only when a large portion of the surrounding environment is removed. These results indicate that, at least in some cases, the functioning of enhancer elements is subordinated to the host chromatin context, which can exert a dominant control over its activity.

Here the authors show that a strong enhancer sequence can be controlled by the chromatin environment provided by a topologically associated domain (TAD) located nearby. An enhancer relocated by homologous recombination takes all the hallmarks of its new neighboring enhancers located in the recipient TAD.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

The Impact of Space and Time on the Functional Output of the Genome

Over the past two decades, it has become clear that the multiscale spatial and temporal organization of the genome has important implications for nuclear function. This review centers on insights gained from recent advances in light microscopy on our understanding of transcription. We discuss spatial and temporal aspects that shape nuclear order and their consequences on regulatory components, focusing on genomic scales most relevant to function. The emerging picture is that spatiotemporal constraints increase the complexity in transcriptional regulation, highlighting new challenges, such as uncertainty about how information travels from molecular factors through the genome and space to generate a functional output.

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

A fast Myosin super enhancer dictates muscle fiber phenotype through competitive interactions with Myosin genes

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here, we identify by snATAC-seq a 42 kb super-enhancer at the locus regrouping the fast Myosin genes. By 4C-seq we show that active fast Myosin promoters interact with this super-enhancer by DNA looping, leading to the activation of a single promoter per nucleus. A rainbow mouse transgenic model of the locus including the super-enhancer recapitulates the endogenous spatio-temporal expression of adult fast Myosin genes. In situ deletion of the super-enhancer by CRISPR/Cas9 editing demonstrates its major role in the control of associated fast Myosin genes, and deletion of two fast Myosin genes at the locus reveals an active competition of the promoters for the shared super-enhancer. Last, by disrupting the organization of fast Myosin, we uncover positional heterogeneity within limb skeletal muscles that may underlie selective muscle susceptibility to damage in certain myopathies.

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here the authors show that a super enhancer controls the spatiotemporal expression of the genes at the fast myosin heavy chain locus by DNA looping and that this expression profile is recapitulated in a rainbow transgenic mouse model of the locus.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Distal and proximal cis-regulatory elements sense X chromosome dosage and developmental state at the Xist locus

Developmental genes such as Xist, which initiates X chromosome inactivation, are controlled by complex cis-regulatory landscapes, which decode multiple signals to establish specific spatiotemporal expression patterns. Xist integrates information on X chromosome dosage and developmental stage to trigger X inactivation in the epiblast specifically in female embryos. Through a pooled CRISPR screen in differentiating mouse embryonic stem cells, we identify functional enhancer elements of Xist at the onset of random X inactivation. Chromatin profiling reveals that X-dosage controls the promoter-proximal region, while differentiation cues activate several distal enhancers. The strongest distal element lies in an enhancer cluster associated with a previously unannotated Xist-enhancing regulatory transcript, which we named Xert. Developmental cues and X-dosage are thus decoded by distinct regulatory regions, which cooperate to ensure female-specific Xist upregulation at the correct developmental time. With this study, we start to disentangle how multiple, functionally distinct regulatory elements interact to generate complex expression patterns in mammals.

Enhancer-Promoter Communication: It's Not Just About Contact

Cis-regulatory elements such as enhancers can be located even a million base pairs away from their cognate promoter and yet modulate gene transcription. Indeed, the 3D organisation of chromatin enables the establishment of long-range enhancer-promoter communication. The observation of long-range enhancer-promoter chromatin loops at active genes originally led to a model in which enhancers and promoters form physical contacts between each other to control transcription. Yet, recent microscopy data has challenged this prevailing activity-by-contact model of enhancer-promoter communication in transcriptional activation. Live single-cell imaging approaches do not systematically reveal a correlation between enhancer-proximity and transcriptional activation. We therefore discuss the need to move from a static to a dynamic view of enhancer-promoter relationships. We highlight recent studies that not only reveal considerable chromatin movement in specific cell types, but suggest links between chromatin compaction, chromatin movement and transcription. We describe the interplay between enhancer-promoter proximity within the context of biomolecular condensates and the need to understand how condensate microenvironments influence the chromatin binding kinetics of proteins that bind at cis-regulatory elements to activate transcription. Finally, given the complex multi-scale interplay between regulatory proteins, enhancer-promoter proximity and movement, we propose the need to integrate information from complementary single-cell next-generation sequencing and live-cell imaging approaches to derive unified 3D theoretical models of enhancer-promoter communication that are ultimately predictive of transcriptional output and cell fate. In time, improved models will shed light on how tissues grow and diseases emerge.

A multi-enhancer RET regulatory code is disrupted in Hirschsprung disease

The major genetic risk factors for Hirschsprung disease (HSCR) are three common polymorphisms within cis-regulatory elements (CREs) of the receptor tyrosine kinase gene RET, which reduce its expression during enteric nervous system (ENS) development. These risk variants attenuate binding of the transcription factors RARB, GATA2, and SOX10 to their cognate CREs, reduce RET gene expression, and dysregulate other ENS and HSCR genes in the RET-EDNRB gene regulatory network (GRN). Here, we use siRNA, ChIP, and CRISPR-Cas9 deletion analyses in the SK-N-SH cell line to ask how many additional HSCR-associated risk variants reside in RET CREs that affect its gene expression. We identify 22 HSCR-associated variants in candidate RET CREs, of which seven have differential allele-specific in vitro enhancer activity, and four of these seven affect RET gene expression; of these, two enhancers are bound by the transcription factor PAX3. We also show that deleting multiple variant-containing enhancers leads to synergistic effects on RET gene expression. These, coupled with our prior results, show that common sequence variants in at least 10 RET enhancers affect HSCR risk, seven with experimental evidence of affecting RET gene expression, extending the known RET-EDNRB GRN to reveal an extensive regulatory code modulating disease risk at a single gene.

Spatial regulation by multiple Gremlin1 enhancers provides digit development with cis-regulatory robustness and evolutionary plasticity

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

The BMP antagonist Gremlin1 balances BMP and SHH signalling, endowing limb bud development with robustness. Here, the authors identify enhancers controlling Grem1 levels in an additive, and spatial regulation in a synergistic manner, providing digit patterning with cis-regulatory robustness and evolutionary plasticity.

Functional human genes typically exhibit epigenetic conservation

Recent DepMap CRISPR-Cas9 single gene disruptions have identified genes more essential to proliferation in tissue culture. It would be valuable to translate these finding with measurements more practical for human tissues. Here we show that DepMap essential genes and other literature curated functional genes exhibit cell-specific preferential epigenetic conservation when DNA methylation measurements are compared between replicate cell lines and between intestinal crypts from the same individual. Culture experiments indicate that epigenetic drift accumulates through time with smaller differences in more functional genes. In NCI-60 cell lines, greater targeted gene conservation correlated with greater drug sensitivity. These studies indicate that two measurements separated in time allow normal or neoplastic cells to signal through conservation which human genes are more essential to their survival in vitro or in vivo.

The Core Promoter Is a Regulatory Hub for Developmental Gene Expression

The development of multicellular organisms and the uniqueness of each cell are achieved by distinct transcriptional programs. Multiple processes that regulate gene expression converge at the core promoter region, an 80 bp region that directs accurate transcription initiation by RNA polymerase II (Pol II). In recent years, it has become apparent that the core promoter region is not a passive DNA component, but rather an active regulatory module of transcriptional programs. Distinct core promoter compositions were demonstrated to result in different transcriptional outputs. In this mini-review, we focus on the role of the core promoter, particularly its downstream region, as the regulatory hub for developmental genes. The downstream core promoter element (DPE) was implicated in the control of evolutionarily conserved developmental gene regulatory networks (GRNs) governing body plan in both the anterior-posterior and dorsal-ventral axes. Notably, the composition of the basal transcription machinery is not universal, but rather promoter-dependent, highlighting the importance of specialized transcription complexes and their core promoter target sequences as key hubs that drive embryonic development, differentiation and morphogenesis across metazoan species. The extent of transcriptional activation by a specific enhancer is dependent on its compatibility with the relevant core promoter. The core promoter content also regulates transcription burst size. Overall, while for many years it was thought that the specificity of gene expression is primarily determined by enhancers, it is now clear that the core promoter region comprises an important regulatory module in the intricate networks of developmental gene expression.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Orphan CpG islands amplify poised enhancer regulatory activity and determine target gene responsiveness

CpG islands (CGIs) represent a widespread feature of vertebrate genomes, being associated with ~70% of all gene promoters. CGIs control transcription initiation by conferring nearby promoters with unique chromatin properties. In addition, there are thousands of distal or orphan CGIs (oCGIs) whose functional relevance is barely known. Here we show that oCGIs are an essential component of poised enhancers that augment their long-range regulatory activity and control the responsiveness of their target genes. Using a knock-in strategy in mouse embryonic stem cells, we introduced poised enhancers with or without oCGIs within topologically associating domains harboring genes with different types of promoters. Analysis of the resulting cell lines revealed that oCGIs act as tethering elements that promote the physical and functional communication between poised enhancers and distally located genes, particularly those with large CGI clusters in their promoters. Therefore, by acting as genetic determinants of gene-enhancer compatibility, CGIs can contribute to gene expression control under both physiological and potentially pathological conditions.

Allosteric conformational ensembles have unlimited capacity for integrating information

Integration of binding information by macromolecular entities is fundamental to cellular functionality. Recent work has shown that such integration cannot be explained by pairwise cooperativities, in which binding is modulated by binding at another site. Higher-order cooperativities (HOCs), in which binding is collectively modulated by multiple other binding events, appear to be necessary but an appropriate mechanism has been lacking. We show here that HOCs arise through allostery, in which effective cooperativity emerges indirectly from an ensemble of dynamically interchanging conformations. Conformational ensembles play important roles in many cellular processes but their integrative capabilities remain poorly understood. We show that sufficiently complex ensembles can implement any form of information integration achievable without energy expenditure, including all patterns of HOCs. Our results provide a rigorous biophysical foundation for analysing the integration of binding information through allostery. We discuss the implications for eukaryotic gene regulation, where complex conformational dynamics accompanies widespread information integration.

Single-nucleotide-level mapping of DNA regulatory elements that control fetal hemoglobin expression

Pinpointing functional noncoding DNA sequences and defining their contributions to health-related traits is a major challenge for modern genetics. We developed a high-throughput framework to map noncoding DNA functions with single-nucleotide resolution in four loci that control erythroid fetal hemoglobin (HbF) expression, a genetically determined trait that modifies sickle cell disease (SCD) phenotypes. Specifically, we used the adenine base editor ABEmax to introduce 10,156 separate A•T to G•C conversions in 307 predicted regulatory elements and quantified the effects on erythroid HbF expression. We identified numerous regulatory elements, defined their epigenomic structures and linked them to low-frequency variants associated with HbF expression in an SCD cohort. Targeting a newly discovered γ-globin gene repressor element in SCD donor CD34+ hematopoietic progenitors raised HbF levels in the erythroid progeny, inhibiting hypoxia-induced sickling. Our findings reveal previously unappreciated genetic complexities of HbF regulation and provide potentially therapeutic insights into SCD.

The haplolethality paradox of the wupA gene in Drosophila

Haplolethals (HL) are regions of diploid genomes that in one dose are fatal for the organism. Their biological meaning is obscure because heterozygous loss-of-function mutations result in dominant lethality (DL) and, consequently, should be under strong negative selection. We report an in depth study of the HL associated to the gene wings up A (wupA). It encodes 13 transcripts (A-M) that yield 11 protein isoforms (A-K) of Troponin I (TnI). They are functionally diverse in their control of muscle contraction, cell polarity and cell proliferation. Isoform K transfers to the nucleus where it increases transcription of the cell proliferation related genes CDK2, CDK4, Rap and Rab5. The nuclear translocation of isoform K is prevented by the co-expression of A or B isoforms, which illustrates isoform interactions. The corresponding DL mutations are, either DNA rearrangements clustered towards the gene 3’ end, thus affecting the genomic organization of all transcripts, or CRISPR-induced mutations in one of the two ATG sites which eliminate a subset of wupA products. The joint elimination of isoforms C, F, G and H, however, do not cause DL phenotypes. Genetically driven expression of single isoforms rescue neither DL nor any of the mutants known in the gene, suggesting that normal function requires properly regulated expression of specific combinations, rather than single, TnI isoforms. We conclude that the wupA associated HL results from the combined haploinsufficiency of a large set of TnI isoforms. The qualitative and quantitative normal expression of which, requires the chromosomal integrity of the wupA genomic region. Since all fly TnI isoforms are encoded in the same gene, its HL condition becomes unavoidable. These wupA features are comparable to those of dpp, the only other HL studied to some extent, and reveal a scenario of strict dosage dependence with implications for gene expression regulation and splitting.

Most species contain two copies of their genetic endowment, each received from each progenitor. If one of the duplicated genes is non-functional, the remaining copy may supply enough product as to cover the requirements for normal function or, alternatively, may reflect the insufficiency through a visible phenotype. In rare occasions, however, having one copy is so deleterious that causes lethality. These so called “haplolethal regions”, exist across species and represent an evolutionary paradox since they should have been subject to intense negative selection. The inherent difficulties to study haplolethals have precluded their study so far. Here, we analyzed the case of one of the five haplolethal regions of Drosophila, the one associated to the Troponin I encoding gene wupA, by measuring the transcriptional effects of mutations and chromosomal rearrangements affecting this gene. The data show that this haplolethality results from the combined insufficiency of a large number of Troponin I isoforms, which are functionally specialized, show interference and require the integrity of the native chromatin structure for their quantitatively regulated expression. These features unveil novel aspects of gene expression and, possibly, on evolutionary gene splitting. Also, haplolethals underscore the biological significance of protein dosage, in particular for functionally related products.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Alternative Enhancer Usage and Targeted Polycomb Marking Hallmark Promoter Choice during T Cell Differentiation

During thymic development and upon peripheral activation, T cells undergo extensive phenotypic and functional changes coordinated by lineage-specific developmental programs. To characterize the regulatory landscape controlling T cell identity, we perform a wide epigenomic and transcriptional analysis of mouse thymocytes and naive CD4 differentiated T helper cells. Our investigations reveal a dynamic putative enhancer landscape, and we could validate many of the enhancers using the high-throughput CapStarr sequencing (CapStarr-seq) approach. We find that genes using multiple promoters display increased enhancer usage, suggesting that apparent "enhancer redundancy" might relate to isoform selection. Furthermore, we can show that two Runx3 promoters display long-range interactions with specific enhancers. Finally, our analyses suggest a novel function for the PRC2 complex in the control of alternative promoter usage. Altogether, our study has allowed for the mapping of an exhaustive set of active enhancers and provides new insights into their function and that of PRC2 in controlling promoter choice during T cell differentiation.