Mapping the Shh long-range regulatory domain

Genetic factors mediating long-range enhancer-promoter communication in mammalian development

Enhancers are remotely located noncoding DNA sequences that regulate gene expression in response to developmental, homeostatic, and environmental cues. Canonical short-range enhancers located <50 kb from their cognate promoters function by binding transcription factors, coactivators, and chromatin modifiers. In this review, we discuss recent evidence that medium-range (50-400 kb) and long-range (>400 kb) enhancers rely on additional mechanisms, including cohesin, CCCTC-binding factor, and high-affinity protein-protein interactions. These mechanisms are crucial for establishing the physical proximity and interaction between enhancers and their target promoters over extended genomic distances and ensuring robust gene activation during mammalian development. Future studies will be critical to unravel their prevalence and evolutionary significance across various genomic loci, cell types, and species.

Multiple allelic configurations govern long-range Shh enhancer-promoter communication in the embryonic forebrain

Developmental gene regulation requires input from enhancers spread over large genomic distances. Our understanding of long-range enhancer-promoter (E-P) communication, characterized as loops, remains incomplete without addressing the role of intervening chromatin. Here, we examine the topology of the entire Sonic hedgehog (Shh) regulatory domain in individual alleles from the mouse embryonic forebrain. Through sequential Oligopaint labeling and super-resolution microscopy, we find that the Shh locus maintains a compact structure that adopts several diverse configurations independent of Shh expression. The most frequent configuration contained distal E-P contacts at the expense of those more proximal to Shh, consistent with an interconnected loop. Genetic perturbations demonstrate that this long-range E-P communication operates by Shh-expression-independent and dependent mechanisms, involving CTCF binding sites and active enhancers, respectively. We propose a model whereby gene regulatory elements secure long-range E-P interactions amid an inherent architectural framework to coordinate spatiotemporal patterns of gene expression.

Common variants increase risk for congenital diaphragmatic hernia within the context of de novo variants

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly often accompanied by other structural anomalies and/or neurobehavioral manifestations. Rare de novo protein-coding variants and copy-number variations contribute to CDH in the population. However, most individuals with CDH remain genetically undiagnosed. Here, we perform integrated de novo and common-variant analyses using 1,469 CDH individuals, including 1,064 child-parent trios and 6,133 ancestry-matched, unaffected controls for the genome-wide association study. We identify candidate CDH variants in 15 genes, including eight novel genes, through deleterious de novo variants. We further identify two genomic loci contributing to CDH risk through common variants with similar effect sizes among Europeans and Latinx. Both loci are in putative transcriptional regulatory regions of developmental patterning genes. Estimated heritability in common variants is ∼19%. Strikingly, there is no significant difference in estimated polygenic risk scores between isolated and complex CDH or between individuals harboring deleterious de novo variants and individuals without these variants. The data support a polygenic model as part of the CDH genetic architecture.

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.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Whole Genome Resequencing Revealed the Genetic Relationship and Selected Regions among Baicheng-You, Beijing-You, and European-Origin Broilers

As the only two You-chicken breeds in China, Baicheng-You (BCY) and Beijing-You (BJY) chickens are famous for their good meat quality. However, so far, the molecular basis of germplasm of the two You-chicken breeds is not yet clear. The genetic relationship among BCY, BJY, and European-origin broilers (BRs) was analyzed using whole genome resequencing data to contribute to this issue. A total of 18,852,372 single nucleotide polymorphisms (SNPs) were obtained in this study. After quality control, 8,207,242 SNPs were applied to subsequent analysis. The data indicated that BJY chickens possessed distant distance with BRs (genetic differentiation coefficient (FST) = 0.1681) and BCY (FST = 0.1231), respectively, while BCY and BRs had a closer relationship (FST = 0.0946). In addition, by using FST, cross-population extended haplotype homozygosity (XP-EHH), and cross-population composite likelihood ratio (XP-CLR) methods, we found 374 selected genes between BJY and BRs chickens and 279 selected genes between BCY and BJY chickens, respectively, which contained a number of important candidates or genetic variations associated with feather growth and fat deposition of BJY chickens and potential disease resistance of BCY chickens. Our study demonstrates a genome-wide view of genetic diversity and differentiation among BCY, BJY, and BRs. These results may provide useful information on a molecular basis related to the special characteristics of these broiler breeds, thus enabling us to better understand the formation mechanism of Chinese-You chickens.

Hnrnpk is essential for embryonic limb bud development as a transcription activator and a collaborator of insulator protein Ctcf

Proper development of the limb bud relies on the concordance of various signals, but its molecular mechanisms have not yet been fully illustrated. Here we report that heterogeneous nuclear ribonucleoprotein K (hnRNPK) is essential for limb bud development. Its ablation in the limb bud results in limbless forelimbs and severe deformities of the hindlimbs. In terms of mechanism, hnRNPK functions as a transcription activator for the vital genes involved in the three regulatory axes of limb bud development. Simultaneously, for the first time we elucidate that hnRNPK binds to and coordinates with the insulator protein CCCTC binding factor (CTCF) to maintain a three-dimensional chromatin architecture. Ablation of hnRNPK weakens the binding strength of CTCF to topologically associating domain (TAD) boundaries, then leading to the loose TADs, and decreased interactions between promoters and enhancers, and further decreased transcription of developmental genes. Our study establishes a fundamental and novel role of hnRNPK in regulating limb bud development.

What is an enhancer?

Tight control of the transcription process is essential for the correct spatial and temporal gene expression pattern during development and in homeostasis. Enhancers are at the core of correct transcriptional activation. The original definition of an enhancer is straightforward: a DNA sequence that activates transcription independent of orientation and direction. Dissection of numerous enhancer loci has shown that many enhancer-like elements might not conform to the original definition, suggesting that enhancers and enhancer-like elements might use multiple different mechanisms to contribute to transcriptional activation. Here, we review methodologies to identify enhancers and enhancer-like elements and discuss pitfalls and consequences for our understanding of transcriptional regulation.

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.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

Structural elements promote architectural stripe formation and facilitate ultra-long-range gene regulation at a human disease locus

Enhancer clusters overlapping disease-associated mutations in Pierre Robin sequence (PRS) patients regulate SOX9 expression at genomic distances over 1.25 Mb. We applied optical reconstruction of chromatin architecture (ORCA) imaging to trace 3D locus topology during PRS-enhancer activation. We observed pronounced changes in locus topology between cell types. Subsequent analysis of single-chromatin fiber traces revealed that these ensemble-average differences arise through changes in the frequency of commonly sampled topologies. We further identified two CTCF-bound elements, internal to the SOX9 topologically associating domain, which promote stripe formation, are positioned near the domain's 3D geometric center, and bridge enhancer-promoter contacts in a series of chromatin loops. Ablation of these elements results in diminished SOX9 expression and altered domain-wide contacts. Polymer models with uniform loading across the domain and frequent cohesin collisions recapitulate this multi-loop, centrally clustered geometry. Together, we provide mechanistic insights into architectural stripe formation and gene regulation over ultra-long genomic ranges.

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.

Chemotherapy suppresses SHH gene expression via a specific enhancer

Sonic hedgehog (SHH) signaling is a key regulator of embryonic development and tissue homeostasis that is involved in gastrointestinal (GI) cancer progression. Regulation of SHH gene expression is a paradigm of long-range enhancer function. Using the classical chemotherapy drug 5-fluorouracil (5FU) as an example, here we show that SHH gene expression is suppressed by chemotherapy. SHH is downstream of immediate early genes (IEGs), including Early growth response 1 (Egr1). A specific 139 kb upstream enhancer is responsible for its down-regulation. Knocking down EGR1 expression or blocking its binding to this enhancer renders SHH unresponsive to chemotherapy. We further demonstrate that down-regulation of SHH expression does not depend on 5FU's impact on nucleotide metabolism or DNA damage; rather, a sustained oxidative stress response mediates this rapid suppression. This enhancer is present in a wide range of tumors and normal tissues, thus providing a target for cancer chemotherapy and its adverse effects on normal tissues. We propose that SHH is a stress-responsive gene downstream of IEGs, and that traditional chemotherapy targets a specific enhancer to suppress its expression.

The role of timing in the development and evolution of the limb

The term heterochrony was coined to describe changes in the timing of developmental processes relative to an ancestral state. Limb development is a well-suited system to address the contribution of heterochrony to morphological evolution. We illustrate how timing mechanisms have been used to establish the correct pattern of the limb and provide cases where natural variations in timing have led to changes in limb morphology.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Cohesin is required for long-range enhancer action at the Shh locus

The regulatory landscapes of developmental genes in mammals can be complex, with enhancers spread over many hundreds of kilobases. It has been suggested that three-dimensional genome organization, particularly topologically associating domains formed by cohesin-mediated loop extrusion, is important for enhancers to act over such large genomic distances. By coupling acute protein degradation with synthetic activation by targeted transcription factor recruitment, here we show that cohesin, but not CTCF, is required for activation of the target gene Shh by distant enhancers in mouse embryonic stem cells. Cohesin is not required for activation directly at the promoter or by an enhancer located closer to the Shh gene. Our findings support the hypothesis that chromatin compaction via cohesin-mediated loop extrusion allows for genes to be activated by enhancers that are located many hundreds of kilobases away in the linear genome and suggests that cohesin is dispensable for enhancers located more proximally.

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.

Nonlinear control of transcription through enhancer-promoter interactions

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

The transcription factor activity gradient (TAG) model: contemplating a contact-independent mechanism for enhancer-promoter communication

How distal cis-regulatory elements (e.g., enhancers) communicate with promoters remains an unresolved question of fundamental importance. Although transcription factors and cofactors are known to mediate this communication, the mechanism by which diffusible molecules relay regulatory information from one position to another along the chromosome is a biophysical puzzle-one that needs to be revisited in light of recent data that cannot easily fit into previous solutions. Here we propose a new model that diverges from the textbook enhancer-promoter looping paradigm and offer a synthesis of the literature to make a case for its plausibility, focusing on the coactivator p300.

Cis-Regulatory Control of Mammalian Sex Determination

Sex determination is the process by which an initial bipotential gonad adopts either a testicular or ovarian cell fate. The inability to properly complete this process leads to a group of developmental disorders classified as disorders of sex development (DSD). To date, dozens of genes were shown to play roles in mammalian sex determination, and mutations in these genes can cause DSD in humans or gonadal sex reversal/dysfunction in mice. However, exome sequencing currently provides genetic diagnosis for only less than half of DSD patients. This points towards a major role for the non-coding genome during sex determination. In this review, we highlight recent advances in our understanding of non-coding, cis-acting gene regulatory elements and discuss how they may control transcriptional programmes that underpin sex determination in the context of the 3-dimensional folding of chromatin. As a paradigm, we focus on the Sox9 gene, a prominent pro-male factor and one of the most extensively studied genes in gonadal cell fate determination.

Complete lung agenesis caused by complex genomic rearrangements with neo-TAD formation at the SHH locus

During human organogenesis, lung development is a timely and tightly regulated developmental process under the control of a large number of signaling molecules. Understanding how genetic variants can disturb normal lung development causing different lung malformations is a major goal for dissecting molecular mechanisms during embryogenesis. Here, through exome sequencing (ES), array CGH, genome sequencing (GS) and Hi-C, we aimed at elucidating the molecular basis of bilateral isolated lung agenesis in three fetuses born to a non-consanguineous family. We detected a complex genomic rearrangement containing duplicated, triplicated and deleted fragments involving the SHH locus in fetuses presenting complete agenesis of both lungs and near-complete agenesis of the trachea, diagnosed by ultrasound screening and confirmed at autopsy following termination. The rearrangement did not include SHH itself, but several regulatory elements for lung development, such as MACS1, a major SHH lung enhancer, and the neighboring genes MNX1 and NOM1. The rearrangement incorporated parts of two topologically associating domains (TADs) including their boundaries. Hi-C of cells from one of the affected fetuses showed the formation of two novel TADs each containing SHH enhancers and the MNX1 and NOM1 genes. Hi-C together with GS indicate that the new 3D conformation is likely causative for this condition by an inappropriate activation of MNX1 included in the neo-TADs by MACS1 enhancer, further highlighting the importance of the 3D chromatin conformation in human disease.

A Highly Conserved Shh Enhancer Coordinates Hypothalamic and Craniofacial Development

Enhancers that are conserved deep in evolutionary time regulate characteristics held in common across taxonomic classes. Here, deletion of the highly conserved Shh enhancer SBE2 (Shh brain enhancer 2) in mouse markedly reduced Shh expression within the embryonic brain specifically in the rostral diencephalon; however, no abnormal anatomical phenotype was observed. Secondary enhancer activity was subsequently identified which likely mediates low levels of expression. In contrast, when crossing the SBE2 deletion with the Shh null allele, brain and craniofacial development were disrupted; thus, linking SBE2 regulated Shh expression to multiple defects and further enabling the study of the effects of differing levels of Shh on embryogenesis. Development of the hypothalamus, derived from the rostral diencephalon, was disrupted along both the anterior-posterior (AP) and the dorsal-ventral (DV) axes. Expression of DV patterning genes and subsequent neuronal population induction were particularly sensitive to Shh expression levels, demonstrating a novel morphogenic context for Shh. The role of SBE2, which is highlighted by DV gene expression, is to step-up expression of Shh above the minimal activity of the second enhancer, ensuring the necessary levels of Shh in a regional-specific manner. We also show that low Shh levels in the diencephalon disrupted neighbouring craniofacial development, including mediolateral patterning of the bones along the cranial floor and viscerocranium. Thus, SBE2 contributes to hypothalamic morphogenesis and ensures there is coordination with the formation of the adjacent midline cranial bones that subsequently protect the neural tissue.

Long-range cis-regulatory elements controlling GDF6 expression are essential for ear development

Molecular mechanisms governing the development of the mammalian cochlea, the hearing organ, remain largely unknown. Through genome sequencing in 3 subjects from 2 families with nonsyndromic cochlear aplasia, we identified homozygous 221-kb and 338-kb deletions in a noncoding region on chromosome 8 with an approximately 200-kb overlapping section. Genomic location of the overlapping deleted region started from approximately 350 kb downstream of GDF6, which codes for growth and differentiation factor 6. Otic lineage cells differentiated from induced pluripotent stem cells derived from an affected individual showed reduced expression of GDF6 compared with control cells. Knockout of Gdf6 in a mouse model resulted in cochlear aplasia, closely resembling the human phenotype. We conclude that GDF6 plays a necessary role in early cochlear development controlled by cis-regulatory elements located within an approximately 500-kb region of the genome in humans and that its disruption leads to deafness due to cochlear aplasia.

A 300-kb microduplication of 7q36.3 in a patient with triphalangeal thumb-polysyndactyly syndrome combined with congenital heart disease and optic disc coloboma: a case report

Background

Triphalangeal thumb-polysyndactyly syndrome (TPT-PS) is a rare well-defined autosomal dominant disorder characterized by long thumbs with three phalanges combined with pre- and postaxial polydactyly/syndactyly of limbs. By now, the syndrome has been reported in several large families from different ethnic backgrounds, with a high degree of inter- and intrafamilial variability. The genome locus responsible for TPT-PS has been mapped to the 7q36.3 region harboring a long-range sonic hedgehog (SHH) regulatory sequence (ZRS). Both single-nucleotide variants and complete duplications of ZRS were shown to cause TPT-PS and similar limb phenotypes. TPT-PS usually forms as isolated limb pathology not associated with additional malformations, in particular, with cardiovascular abnormalities.

Case presentation

Here we report on a rare Russian neonatal case of TPT-PS combined with severe congenital heart disease, namely double outlet right ventricle, and microphthalmia with optic disc coloboma. Pedigree analysis revealed TPT-PS of various expressivity in 10 family members throughout five generations, while the cardiac defect and the eye pathology were detected only in the proband. To extend the knowledge on genotype–phenotype spectrum of TPT-PS, the careful clinical and genomic analysis of the family was performed. High-resolution array-based comparative genomic hybridization (array-CGH) revealed a ~ 300 kb microduplication of 7q36.3 locus (arr[GRCh37] 7q36.3(156385810_156684811) × 3) that co-segregated with TPT-PS in the proband and her mother. The duplication encompassed three genes including LMBR1, the intron 5 of which is known to harbor ZRS. Based on whole-exome sequencing data, no additional pathogenic mutations or variants of uncertain clinical significance were found in morbid cardiac genes or genes associated with a microphthalmia/anophthalmia/coloboma spectrum of ocular malformations.

Conclusions

The results support the previous data, indicating that complete ZRS duplication underlies TPT-PS, and suggest a broader phenotypic impact of the 7q36.3 microduplication. Potential involvement of the 7q36.3 microduplication in the patient’s cardiac and eye malformations is discussed. However, the contribution of some additional genetic/epigenetic factors to the complex patient`s phenotype cannot be excluded entirely. Further comprehensive functional studies are needed to prove the possible involvement of the 7q36.3 locus in congenital heart disease and eye pathology.

Transcriptional enhancers: from prediction to functional assessment on a genome-wide scale

Enhancers are cis-regulatory sequences located distally to target genes. These sequences consolidate developmental and environmental cues to coordinate gene expression in a tissue-specific manner. Enhancer function and tissue specificity depend on the expressed set of transcription factors, which recognize binding sites and recruit cofactors that regulate local chromatin organization and gene transcription. Unlike other genomic elements, enhancers are challenging to identify because they function independently of orientation, are often distant from their promoters, have poorly defined boundaries, and display no reading frame. In addition, there are no defined genetic or epigenetic features that are unambiguously associated with enhancer activity. Over recent years there have been developments in both empirical assays and computational methods for enhancer prediction. We review genome-wide tools, CRISPR advancements, and high-throughput screening approaches that have improved our ability to both observe and manipulate enhancers in vitro at the level of primary genetic sequences, chromatin states, and spatial interactions. We also highlight contemporary animal models and their importance to enhancer validation. Together, these experimental systems and techniques complement one another and broaden our understanding of enhancer function in development, evolution, and disease.

Gene regulatory landscape of the sonic hedgehog locus in embryonic development

The organs of vertebrate species display a wide variety of morphology. A remaining challenge in evolutionary developmental biology is to elucidate how vertebrate lineages acquire distinct morphological features. Developmental programs are driven by spatiotemporal regulation of gene expression controlled by hundreds of thousands of cis-regulatory elements. Changes in the regulatory elements caused by the introduction of genetic variants can confer regulatory innovation that may underlie morphological novelties. Recent advances in sequencing technology have revealed a number of potential regulatory variants that can alter gene expression patterns. However, a limited number of studies demonstrate causal dependence between genetic and morphological changes. Regulation of Shh expression is a good model to understand how multiple regulatory elements organize tissue-specific gene expression patterns. This model also provides insights into how evolution of molecular traits, such as gene regulatory networks, lead to phenotypic novelty.

Regulatory Landscaping: How Enhancer-Promoter Communication Is Sculpted in 3D

During embryogenesis, precise gene transcription in space and time requires that distal enhancers and promoters communicate by physical proximity within gene regulatory landscapes. To achieve this, regulatory landscapes fold in nuclear space, creating complex 3D structures that influence enhancer-promoter communication and gene expression and that, when disrupted, can cause disease. Here, we provide an overview of how enhancers and promoters construct regulatory landscapes and how multiple scales of 3D chromatin structure sculpt their communication. We focus on emerging views of what enhancer-promoter contacts and chromatin domains physically represent and how two antagonistic fundamental forces-loop extrusion and homotypic attraction-likely form them. We also examine how these same forces spatially separate regulatory landscapes by functional state, thereby creating higher-order compartments that reconfigure during development to enable proper enhancer-promoter communication.

Temporospatial sonic hedgehog signalling is essential for neural crest-dependent patterning of the intrinsic tongue musculature

The tongue is a highly specialised muscular organ with a complex anatomy required for normal function. We have utilised multiple genetic approaches to investigate local temporospatial requirements for sonic hedgehog (SHH) signalling during tongue development. Mice lacking a Shh cis-enhancer, MFCS4 (ShhMFCS4/-), with reduced SHH in dorsal tongue epithelium have perturbed lingual septum tendon formation and disrupted intrinsic muscle patterning, with these defects reproduced following global Shh deletion from E10.5 in pCag-CreERTM; Shhflox/flox embryos. SHH responsiveness was diminished in local cranial neural crest cell (CNCC) populations in both mutants, with SHH targeting these cells through the primary cilium. CNCC-specific deletion of orofaciodigital syndrome 1 (Ofd1), which encodes a ciliary protein, in Wnt1-Cre; Ofdfl/Y mice led to a complete loss of normal myotube arrangement and hypoglossia. In contrast, mesoderm-specific deletion of Ofd1 in Mesp1-Cre; Ofdfl/Y embryos resulted in normal intrinsic muscle arrangement. Collectively, these findings suggest key temporospatial requirements for local SHH signalling in tongue development (specifically, lingual tendon differentiation and intrinsic muscle patterning through signalling to CNCCs) and provide further mechanistic insight into the tongue anomalies seen in patients with disrupted hedgehog signalling.

Developmentally regulated Shh expression is robust to TAD perturbations

Topologically associating domains (TADs) have been proposed to both guide and constrain enhancer activity. Shh is located within a TAD known to contain all its enhancers. To investigate the importance of chromatin conformation and TAD integrity on developmental gene regulation, we have manipulated the Shh TAD - creating internal deletions, deleting CTCF sites, and deleting and inverting sequences at TAD boundaries. Chromosome conformation capture and fluorescence in situ hybridisation assays were used to investigate the changes in chromatin conformation that result from these manipulations. Our data suggest that these substantial alterations in TAD structure have no readily detectable effect on Shh expression patterns or levels of Shh expression during development - except where enhancers are deleted - and result in no detectable phenotypes. Only in the case of a larger deletion at one TAD boundary could ectopic influence of the Shh limb enhancer be detected on a gene (Mnx1) in the neighbouring TAD. Our data suggests that, contrary to expectations, the developmental regulation of Shh expression is remarkably robust to TAD perturbations.

Chromatin topology, condensates and gene regulation: shifting paradigms or just a phase?

In the past decade, two major advances in our understanding of nuclear organization have taken the field of gene regulation by storm. First, technologies that can analyze the three-dimensional conformation of chromatin have revealed how the genome is organized and have provided novel insights into how regulatory regions in the genome interact. Second, the recognition that many proteins can form membraneless compartments through liquid-liquid phase separation (LLPS) has challenged long-standing notions of how proteins within the nucleus are organized and has offered a tantalizing general mechanism by which many aspects of nuclear function may be regulated. However, the functional roles of chromatin topology and LLPS in regulating gene expression remain poorly understood. These topics were discussed with great fervor during an open discussion held at a recent workshop titled 'Chromatin-based regulation of development' organized by The Company of Biologists. Here, we summarize the major points covered during this debate and discuss how they tie into current thinking in the field of gene regulation.

A multidisciplinary review of triphalangeal thumb

Despite being a rare congenital limb anomaly, triphalangeal thumb is a subject of research in various scientific fields, providing new insights in clinical research and evolutionary biology. The findings of triphalangeal thumb can be predictive for other congenital anomalies as part of an underlying syndrome. Furthermore, triphalangeal thumb is still being used as a model in molecular genetics to study gene regulation by long-range regulatory elements. We present a review that summarizes a number of scientifically relevant topics that involve the triphalangeal thumb phenotype. Future initiatives involving multidisciplinary teams collaborating in the field of triphalangeal thumb research can lead to a better understanding of the pathogenesis and molecular mechanisms of this condition as well as other congenital upper limb anomalies.

CNEr: A toolkit for exploring extreme noncoding conservation

Conserved Noncoding Elements (CNEs) are elements exhibiting extreme noncoding conservation in Metazoan genomes. They cluster around developmental genes and act as long-range enhancers, yet nothing that we know about their function explains the observed conservation levels. Clusters of CNEs coincide with topologically associating domains (TADs), indicating ancient origins and stability of TAD locations. This has suggested further hypotheses about the still elusive origin of CNEs, and has provided a comparative genomics-based method of estimating the position of TADs around developmentally regulated genes in genomes where chromatin conformation capture data is missing. To enable researchers in gene regulation and chromatin biology to start deciphering this phenomenon, we developed CNEr, a R/Bioconductor toolkit for large-scale identification of CNEs and for studying their genomic properties. We apply CNEr to two novel genome comparisons—fruit fly vs tsetse fly, and two sea urchin genomes—and report novel insights gained from their analysis. We also show how to reveal interesting characteristics of CNEs by coupling CNEr with existing Bioconductor packages. CNEr is available at Bioconductor (https://bioconductor.org/packages/CNEr/) and maintained at github (https://github.com/ge11232002/CNEr).

Preformed chromatin topology assists transcriptional robustness of Shh during limb development

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS, in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh-ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh-ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

Conserved non-coding elements: developmental gene regulation meets genome organization

Comparative genomics has revealed a class of non-protein-coding genomic sequences that display an extraordinary degree of conservation between two or more organisms, regularly exceeding that found within protein-coding exons. These elements, collectively referred to as conserved non-coding elements (CNEs), are non-randomly distributed across chromosomes and tend to cluster in the vicinity of genes with regulatory roles in multicellular development and differentiation. CNEs are organized into functional ensembles called genomic regulatory blocks–dense clusters of elements that collectively coordinate the expression of shared target genes, and whose span in many cases coincides with topologically associated domains. CNEs display sequence properties that set them apart from other sequences under constraint, and have recently been proposed as useful markers for the reconstruction of the evolutionary history of organisms. Disruption of several of these elements is known to contribute to diseases linked with development, and cancer. The emergence, evolutionary dynamics and functions of CNEs still remain poorly understood, and new approaches are required to enable comprehensive CNE identification and characterization. Here, we review current knowledge and identify challenges that need to be tackled to resolve the impasse in understanding extreme non-coding conservation.

Fibroblast growth factors (FGFs) prime the limb specific Shh enhancer for chromatin changes that balance histone acetylation mediated by E26 transformation-specific (ETS) factors

Sonic hedgehog (Shh) expression in the limb bud organizing centre called the zone of polarizing activity is regulated by the ZRS enhancer. Here, we examine in mouse and in a mouse limb-derived cell line the dynamic events that activate and restrict the spatial activity of the ZRS. Fibroblast growth factor (FGF) signalling in the distal limb primes the ZRS at early embryonic stages maintaining a poised, but inactive state broadly across the distal limb mesenchyme. The E26 transformation-specific transcription factor, ETV4, which is induced by FGF signalling and acts as a repressor of ZRS activity, interacts with the histone deacetylase HDAC2 and ensures that the poised ZRS remains transcriptionally inactive. Conversely, GABPα, an activator of the ZRS, recruits p300, which is associated with histone acetylation (H3K27ac) indicative of an active enhancer. Hence, the primed but inactive state of the ZRS is induced by FGF signalling and in combination with balanced histone modification events establishes the restricted, active enhancer responsible for patterning the limb bud during development.

As an animal embryo develops, specific genes need to be switched on and off at the right time and place to ensure that the embryo’s tissues and organs form properly. Proteins called transcription factors control the activity of individual genes by binding to regions of DNA known as enhancers. Changes in the way DNA is packaged inside cells can affect the ability of transcription factors to access the enhancers, and therefore also influence when particular genes are switched on or off.

Sonic hedgehog (or Shh for short) is a gene that helps to control various aspects of development including the formation of the limbs and brain. The limb forms from a structure in the embryo referred to as the limb bud. An enhancer called ZRS regulates the precise position within the limb bud where the Shh gene is active in a region designated as the “zone of polarizing activity”. Yet, it was not known how the enhancer is controlled to ensure this pattern is achieved. Peluso et al. investigated the events that lead to ZRS becoming active in mice embryos.

The experiments show that the ZRS enhancer exists in three different states in cells across the limb bud: poised, active and inactive. The enhancer is poised in a broad region of the limb bud in cells that are potentially able to switch on the Shh gene. Proteins called fibroblast growth factors drive the enhancer to enter this poised state by altering the way the DNA containing the enhancer is packaged in the cell. Specific transcription factors are able to bind to the poised enhancer and it is the balance between these different transcription factors that activates the enhancer in the zone of polarizing activity. Furthermore in the region of the limb bud where the fibroblast growth factors are not present the ZRS is inactive.

These findings show that fibroblast growth factors, in combination with other changes to the ZRS enhancer, restrict the area in which the enhancer is active to a particular region of the limb bud. Differences in enhancer elements are known to underlie a range of inherited characteristics and may influence whether an individual develops many common diseases. In the future, investigating how cells control the activity of enhancers may provide clues to identifying new targets for drugs to treat some of these diseases.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

Crossed wires: 3D genome misfolding in human disease

Norton and Phillips-Cremins review the 3D architecture of the genome and discuss links between chromatin misfolding and human diseases.

Mammalian genomes are folded into unique topological structures that undergo precise spatiotemporal restructuring during healthy development. Here, we highlight recent advances in our understanding of how the genome folds inside the 3D nucleus and how these folding patterns are miswired during the onset and progression of mammalian disease states. We discuss potential mechanisms underlying the link among genome misfolding, genome dysregulation, and aberrant cellular phenotypes. We also discuss cases in which the endogenous 3D genome configurations in healthy cells might be particularly susceptible to mutation or translocation. Together, these data support an emerging model in which genome folding and misfolding is critically linked to the onset and progression of a broad range of human diseases.

Prenatal identification of two discontinuous maternally inherited chromosome 7q36.3 microduplications totaling 507 kb including the sonic hedgehog gene in a fetus with multiple congenital anomalies

Duplications of the SHH gene, an important developmental gene, are rare. Disruption of this gene produces a variable phenotype in humans from major anomalies to isolated facial defects. This is the first reported case of a maternally inherited 507 kb discontinuous chromosome 7q36.3 microduplication resulting in duplication of SHH and nearby enhancer sequences.

Smoothened antagonist GDC-0449 (Vismodegib) inhibits proliferation and triggers apoptosis in colon cancer cell lines

The sonic hedgehog (Shh) pathway has been proven to be involved in embryonic development and cancer growth. GDC-0449, an antagonist of the hedgehog signaling receptor Smoothened (Smo), was recently approved by the US Food and Drug Administration as a prescription for skin basal cell carcinoma. However, the efficacy of GDC-0449 in the treatment of colon cancer and other malignancies, such as basal cell carcinoma and pancreatic cancer, has remained to be proven. The present study assessed the effect of GDC-0449 on the colon cancer cell lines Caco-2 and Ht-29. A Cell Counting Kit-8 assay was applied to assess the cell proliferation rate and apoptosis was tested by flow cytometry. Reverse-transcription quantitative PCR and western blot analysis were used for analyzing expression levels of target genes. Cell proliferation was inhibited, while apoptosis was increased by GDC-0449, whereas the expression of B-cell lymphoma 2 (Bcl-2), a downstream target of Shh signaling, was decreased. Consistent with the inhibition of Gli1 expression, the cancer stem cell markers CD44 and ALDH were decreased in the presence of GDC-0449. In conclusion, GDC-0449 was shown to inhibit the replication of colon cancer cells and trigger apoptosis through downregulating Bcl-2. This may also influence the stemness of cancer stem cells as indicated by the decreased stem cell surface markers.

Evolution of Shh endoderm enhancers during morphological transition from ventral lungs to dorsal gas bladder

Shh signalling plays a crucial role for endoderm development. A Shh endoderm enhancer, MACS1, is well conserved across terrestrial animals with lungs. Here, we first show that eliminating mouse MACS1 causes severe defects in laryngeal development, indicating that MACS1-directed Shh signalling is indispensable for respiratory organogenesis. Extensive phylogenetic analyses revealed that MACS1 emerged prior to the divergence of cartilaginous and bony fishes, and even euteleost fishes have a MACS1 orthologue. Meanwhile, ray-finned fishes evolved a novel conserved non-coding sequence in the neighbouring region. Transgenic assays showed that MACS1 drives reporter expression ventrally in laryngeal epithelium. This activity has been lost in the euteleost lineage, and instead, the conserved non-coding sequence of euteleosts acquired an enhancer activity to elicit dorsal epithelial expression in the posterior pharynx and oesophagus. These results implicate that evolution of these two enhancers is relevant to the morphological transition from ventral lungs to dorsal gas bladder.

Endoderm enhancer MACS1 of Sonic Hedgehog is conserved in animals with lungs. Here, the authors show that mouse without MACS1 has defective laryngeal development, and use phylogenetic analyses to show association of evolutionary lung-gas bladder transition with change of the enhancer.

Pan-cancer analysis of somatic copy-number alterations implicates IRS4 and IGF2 in enhancer hijacking

Extensive prior research focused on somatic copy-number alterations (SCNAs) affecting cancer genes, yet the extent to which recurrent SCNAs exert their influence through rearrangement of cis-regulatory elements (CREs) remains unclear. Here we present a framework for inferring cancer-related gene overexpression resulting from CRE reorganization (e.g., enhancer hijacking) by integrating SCNAs, gene expression data and information on topologically associating domains (TADs). Analysis of 7,416 cancer genomes uncovered several pan-cancer candidate genes, including IRS4, SMARCA1 and TERT. We demonstrate that IRS4 overexpression in lung cancer is associated with recurrent deletions in cis, and we present evidence supporting a tumor-promoting role. We additionally pursued cancer-type-specific analyses and uncovered IGF2 as a target for enhancer hijacking in colorectal cancer. Recurrent tandem duplications intersecting with a TAD boundary mediate de novo formation of a 3D contact domain comprising IGF2 and a lineage-specific super-enhancer, resulting in high-level gene activation. Our framework enables systematic inference of CRE rearrangements mediating dysregulation in cancer.

The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances

Gene expression often requires interaction between promoters and distant enhancers, which occur within the context of highly organized topologically associating domains (TADs). Using a series of engineered chromosomal rearrangements at the Shh locus, we carried out an extensive fine-scale characterization of the factors that govern the long-range regulatory interactions controlling Shh expression. We show that Shh enhancers act pervasively, yet not uniformly, throughout the TAD. Importantly, changing intra-TAD distances had no impact on Shh expression. In contrast, inversions disrupting the TAD altered global folding of the region and prevented regulatory contacts in a distance-dependent manner. Our data indicate that the Shh TAD promotes distance-independent contacts between distant regions that would otherwise interact only sporadically, enabling functional communication between them. In large genomes where genomic distances per se can limit regulatory interactions, this function of TADs could be as essential for gene expression as the formation of insulated neighborhoods.

Structural Chromosomal Rearrangements Require Nucleotide-Level Resolution: Lessons from Next-Generation Sequencing in Prenatal Diagnosis

In this exciting era of "next-gen cytogenetics," integrating genomic sequencing into the prenatal diagnostic setting is possible within an actionable time frame and can provide precise delineation of balanced chromosomal rearrangements at the nucleotide level. Given the increased risk of congenital abnormalities in newborns with de novo balanced chromosomal rearrangements, comprehensive interpretation of breakpoints could substantially improve prediction of phenotypic outcomes and support perinatal medical care. Herein, we present and evaluate sequencing results of balanced chromosomal rearrangements in ten prenatal subjects with respect to the location of regulatory chromatin domains (topologically associated domains [TADs]). The genomic material from all subjects was interpreted to be "normal" by microarray analyses, and their rearrangements would not have been detected by cell-free DNA (cfDNA) screening. The findings of our systematic approach correlate with phenotypes of both pregnancies with untoward outcomes (5/10) and with healthy newborns (3/10). Two pregnancies, one with a chromosomal aberration predicted to be of unknown clinical significance and another one predicted to be likely benign, were terminated prior to phenotype-genotype correlation (2/10). We demonstrate that the clinical interpretation of structural rearrangements should not be limited to interruption, deletion, or duplication of specific genes and should also incorporate regulatory domains of the human genome with critical ramifications for the control of gene expression. As detailed in this study, our molecular approach to both detecting and interpreting the breakpoints of structural rearrangements yields unparalleled information in comparison to other commonly used first-tier diagnostic methods, such as non-invasive cfDNA screening and microarray analysis, to provide improved genetic counseling for phenotypic outcome in the prenatal setting.

GATA4 and GATA6 regulate pancreatic endoderm identity through inhibition of hedgehog signaling

GATA4 and GATA6 are zinc finger transcription factors that have important functions in several mesodermal and endodermal organs, including heart, liver and pancreas. In humans, heterozygous mutations of either factor are associated with pancreatic agenesis; however, homozygous deletion of both Gata4 and Gata6 is necessary to disrupt pancreas development in mice. In this study, we demonstrate that arrested pancreatic development in Gata4(fl/fl); Gata6(fl/fl); Pdx1:Cre (pDKO) embryos is accompanied by the transition of ventral and dorsal pancreatic fates into intestinal or stomach lineages, respectively. These results indicate that GATA4 and GATA6 play essential roles in maintaining pancreas identity by regulating foregut endodermal fates. Remarkably, pancreatic anlagen derived from pDKO embryos also display a dramatic upregulation of hedgehog pathway components, which are normally absent from the presumptive pancreatic endoderm. Consistent with the erroneous activation of hedgehog signaling, we demonstrate that GATA4 and GATA6 are able to repress transcription through the sonic hedgehog (Shh) endoderm-specific enhancer MACS1 and that GATA-binding sites within this enhancer are necessary for this repressive activity. These studies establish the importance of GATA4/6-mediated inhibition of hedgehog signaling as a major mechanism regulating pancreatic endoderm specification during patterning of the gut tube.