Ultraconserved Enhancers Are Required for Normal Development.

Non-coding "ultraconserved" regions containing hundreds of consecutive bases of perfect sequence conservation across mammalian genomes can function as distant-acting enhancers. However, initial deletion studies in mice revealed that loss of such extraordinarily constrained sequences had no immediate impact on viability. Here, we show that ultraconserved enhancers are required for normal development. Focusing on some of the longest ultraconserved sites genome wide, located near the essential neuronal transcription factor Arx, we used genome editing to create an expanded series of knockout mice lacking individual or combinations of ultraconserved enhancers. Mice with single or pairwise deletions of ultraconserved enhancers were viable and fertile but in nearly all cases showed neurological or growth abnormalities, including substantial alterations of neuron populations and structural brain defects. Our results demonstrate the functional importance of ultraconserved enhancers and indicate that remarkably strong sequence conservation likely results from fitness deficits that appear subtle in a laboratory setting.

Progressive Loss of Function in a Limb Enhancer during Snake Evolution.

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.

Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions.

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.

Single cell evaluation of endocardial Hand2 gene regulatory networks reveals HAND2-dependent pathways that impact cardiac morphogenesis.

The transcription factor HAND2 plays essential roles during cardiogenesis. Hand2 endocardial deletion (H2CKO) results in tricuspid atresia or double inlet left ventricle with accompanying intraventricular septum defects, hypo-trabeculated ventricles and an increased density of coronary lumens. To understand the regulatory mechanisms of these phenotypes, single cell transcriptome analysis of mouse E11.5 H2CKO hearts was performed revealing a number of disrupted endocardial regulatory pathways. Using HAND2 DNA occupancy data, we identify several HAND2-dependent enhancers, including two endothelial enhancers for the shear-stress master regulator KLF2. A 1.8 kb enhancer located 50 kb upstream of the Klf2 TSS imparts specific endothelial/endocardial expression within the vasculature and endocardium. This enhancer is HAND2-dependent for ventricular endocardium expression but HAND2-independent for Klf2 vascular and valve expression. Deletion of this Klf2 enhancer results in reduced Klf2 expression within ventricular endocardium. These data reveal that HAND2 functions within endocardial gene regulatory networks including shear-stress response.

Induction of kidney-related gene programs through co-option of SALL1 in mole ovotestes.

Changes in gene expression represent an important source of phenotypic innovation. Yet how such changes emerge and impact the evolution of traits remains elusive. Here, we explore the molecular mechanisms associated with the development of masculinizing ovotestes in female moles. By performing integrative analyses of epigenetic and transcriptional data in mole and mouse, we identified the co-option of SALL1 expression in mole ovotestes formation. Chromosome conformation capture analyses highlight a striking conservation of the 3D organization at the SALL1 locus, but an evolutionary divergence of enhancer activity. Interspecies reporter assays support the capability of mole-specific enhancers to activate transcription in urogenital tissues. Through overexpression experiments in transgenic mice, we further demonstrate the capability of SALL1 to induce kidney-related gene programs, which are a signature of mole ovotestes. Our results highlight the co-option of gene expression, through changes in enhancer activity, as a plausible mechanism for the evolution of traits.

Wt1 transcription factor impairs cardiomyocyte specification and drives a phenotypic switch from myocardium to epicardium.

During development, the heart grows by addition of progenitor cells to the poles of the primordial heart tube. In the zebrafish, Wilms tumor 1 transcription factor a (wt1a) and b (wt1b) genes are expressed in the pericardium, at the venous pole of the heart. From this pericardial layer, the proepicardium emerges. Proepicardial cells are subsequently transferred to the myocardial surface and form the epicardium, covering the myocardium. We found that while wt1a and wt1b expression is maintained in proepicardial cells, it is downregulated in pericardial cells that contributes cardiomyocytes to the developing heart. Sustained wt1b expression in cardiomyocytes reduced chromatin accessibility of specific genomic loci. Strikingly, a subset of wt1a- and wt1b-expressing cardiomyocytes changed their cell-adhesion properties, delaminated from the myocardium and upregulated epicardial gene expression. Thus, wt1a and wt1b act as a break for cardiomyocyte differentiation, and ectopic wt1a and wt1b expression in cardiomyocytes can lead to their transdifferentiation into epicardial-like cells.

Noncoding deletions reveal a gene that is critical for intestinal function.

Large-scale genome sequencing is poised to provide a substantial increase in the rate of discovery of disease-associated mutations, but the functional interpretation of such mutations remains challenging. Here we show that deletions of a sequence on human chromosome 16 that we term the intestine-critical region (ICR) cause intractable congenital diarrhoea in infants1,2. Reporter assays in transgenic mice show that the ICR contains a regulatory sequence that activates transcription during the development of the gastrointestinal system. Targeted deletion of the ICR in mice caused symptoms that recapitulated the human condition. Transcriptome analysis revealed that an unannotated open reading frame (Percc1) flanks the regulatory sequence, and the expression of this gene was lost in the developing gut of mice that lacked the ICR. Percc1-knockout mice displayed phenotypes similar to those observed upon ICR deletion in mice and patients, whereas an ICR-driven Percc1 transgene was sufficient to rescue the phenotypes found in mice that lacked the ICR. Together, our results identify a gene that is critical for intestinal function and underscore the need for targeted in vivo studies to interpret the growing number of clinical genetic findings that do not affect known protein-coding genes.

Identification and characterization of Hand2 upstream genomic enhancers active in developing stomach and limbs.

BACKGROUND: The bHLH transcription factor HAND2 plays important roles in the development of the embryonic heart, face, limbs, and sympathetic and enteric nervous systems. To define how and when HAND2 regulates these developmental systems, requires understanding the transcriptional regulation of Hand2. RESULTS: Remarkably, Hand2 is flanked by an extensive upstream gene desert containing a potentially diverse enhancer landscape. Here, we screened the regulatory interval 200 kb proximal to Hand2 for putative enhancers using evolutionary conservation and histone marks in Hand2-expressing tissues. H3K27ac signatures across embryonic tissues pointed to only two putative enhancer regions showing deep sequence conservation. Assessment of the transcriptional enhancer potential of these elements using transgenic reporter lines uncovered distinct in vivo enhancer activities in embryonic stomach and limb mesenchyme, respectively. Activity of the identified stomach enhancer was restricted to the developing antrum and showed expression within the smooth muscle and enteric neurons. Surprisingly, the activity pattern of the limb enhancer did not overlap Hand2 mRNA but consistently yielded a defined subectodermal anterior expression pattern within multiple transgenic lines. CONCLUSIONS: Together, these results start to uncover the diverse regulatory potential inherent to the Hand2 upstream regulatory interval.

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals.

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

Enhancer redundancy provides phenotypic robustness in mammalian development.

Distant-acting tissue-specific enhancers, which regulate gene expression, vastly outnumber protein-coding genes in mammalian genomes, but the functional importance of this regulatory complexity remains unclear. Here we show that the pervasive presence of multiple enhancers with similar activities near the same gene confers phenotypic robustness to loss-of-function mutations in individual enhancers. We used genome editing to create 23 mouse deletion lines and inter-crosses, including both single and combinatorial enhancer deletions at seven distinct loci required for limb development. Unexpectedly, none of the ten deletions of individual enhancers caused noticeable changes in limb morphology. By contrast, the removal of pairs of limb enhancers near the same gene resulted in discernible phenotypes, indicating that enhancers function redundantly in establishing normal morphology. In a genetic background sensitized by reduced baseline expression of the target gene, even single enhancer deletions caused limb abnormalities, suggesting that functional redundancy is conferred by additive effects of enhancers on gene expression levels. A genome-wide analysis integrating epigenomic and transcriptomic data from 29 developmental mouse tissues revealed that mammalian genes are very commonly associated with multiple enhancers that have similar spatiotemporal activity. Systematic exploration of three representative developmental structures (limb, brain and heart) uncovered more than one thousand cases in which five or more enhancers with redundant activity patterns were found near the same gene. Together, our data indicate that enhancer redundancy is a remarkably widespread feature of mammalian genomes that provides an effective regulatory buffer to prevent deleterious phenotypic consequences upon the loss of individual enhancers.

HAND transcription factors cooperatively specify the aorta and pulmonary trunk.

Congenital heart defects (CHDs) affecting the cardiac outflow tract (OFT) constitute a significant cause of morbidity and mortality. The OFT develops from migratory cell populations which include the cardiac neural crest cells (cNCCs) and secondary heart field (SHF) derived myocardium and endocardium. The related transcription factors HAND1 and HAND2 have been implicated in human CHDs involving the OFT. Although Hand1 is expressed within the OFT, Hand1 NCC-specific conditional knockout mice (H1CKOs) are viable. Here we show that these H1CKOs present a low penetrance of OFT phenotypes, whereas SHF-specific Hand1 ablation does not reveal any cardiac phenotypes. Further, HAND1 and HAND2 appear functionally redundant within the cNCCs, as a reduction/ablation of Hand2 on an NCC-specific H1CKO background causes pronounced OFT defects. Double conditional Hand1 and Hand2 NCC knockouts exhibit persistent truncus arteriosus (PTA) with 100% penetrance. NCC lineage-tracing and Sema3c in situ mRNA expression reveal that Sema3c-expressing cells are mis-localized, resulting in a malformed septal bridge within the OFTs of H1CKO;H2CKO embryos. Interestingly, Hand1 and Hand2 also genetically interact within the SHF, as SHF H1CKOs on a heterozygous Hand2 background exhibit Ventricular Septal Defects (VSDs) with incomplete penetrance. Previously, we identified a BMP, HAND2, and GATA-dependent Hand1 OFT enhancer sufficient to drive reporter gene expression within the nascent OFT and aorta. Using these transcription inputs as a probe, we identify a novel Hand2 OFT enhancer, suggesting that a conserved BMP-GATA dependent mechanism transcriptionally regulates both HAND factors. These findings support the hypothesis that HAND factors interpret BMP signaling within the cNCCs to cooperatively coordinate OFT morphogenesis.

Interrogating the Grainyhead-like 2 (Grhl2) genomic locus identifies an enhancer element that regulates palatogenesis in mouse.

The highly-conserved Grainyhead-like (Grhl) transcription factors are critical regulators of embryogenesis that regulate cellular survival, proliferation, migration and epithelial integrity, especially during the formation of the craniofacial skeleton. Family member Grhl2 is expressed throughout epithelial tissues during development, and loss of Grhl2 function leads to significant defects in neurulation, abdominal wall closure, formation of the face and fusion of the maxilla/palate. Whereas numerous downstream target genes of Grhl2 have been identified, very little is known about how this crucial developmental transcription factor itself is regulated. Here, using in silico and in utero expression analyses and functional deletion in mice, we have identified a novel 2.4 â€‹kb enhancer element (mm1286) that drives reporter gene expression in a pattern that strongly recapitulates endogenous Grhl2 in the craniofacial primordia, modulates Grhl2 expression in these tissues, and augments Grhl2-mediated closure of the secondary palate. Deletion of this genomic element, in the context of inactivation of one allele of Grhl2 (through generation of double heterozygous Grhl2+/-;mm1286+/- mice), results in a significant predisposition to palatal clefting at birth. Moreover, we found that a highly conserved 325 bp region of mm1286 is both necessary and sufficient for mediating the craniofacial-specific enhancer activity of this region, and that an extremely well-conserved 12-bp sequence within this element (CTGTCAAACAGGT) substantially determines full enhancer function. Together, these data provide valuable new insights into the upstream genomic regulatory landscape responsible for transcriptional control of Grhl2 during palatal closure.

A unique stylopod patterning mechanism by Shox2-controlled osteogenesis.

Vertebrate appendage patterning is programmed by Hox-TALE factor-bound regulatory elements. However, it remains unclear which cell lineages are commissioned by Hox-TALE factors to generate regional specific patterns and whether other Hox-TALE co-factors exist. In this study, we investigated the transcriptional mechanisms controlled by the Shox2 transcriptional regulator in limb patterning. Harnessing an osteogenic lineage-specific Shox2 inactivation approach we show that despite widespread Shox2 expression in multiple cell lineages, lack of the stylopod observed upon Shox2 deficiency is a specific result of Shox2 loss of function in the osteogenic lineage. ChIP-Seq revealed robust interaction of Shox2 with cis-regulatory enhancers clustering around skeletogenic genes that are also bound by Hox-TALE factors, supporting a lineage autonomous function of Shox2 in osteogenic lineage fate determination and skeleton patterning. Pbx ChIP-Seq further allowed the genome-wide identification of cis-regulatory modules exhibiting co-occupancy of Pbx, Meis and Shox2 transcriptional regulators. Integrative analysis of ChIP-Seq and RNA-Seq data and transgenic enhancer assays indicate that Shox2 patterns the stylopod as a repressor via interaction with enhancers active in the proximal limb mesenchyme and antagonizes the repressive function of TALE factors in osteogenesis.

Limb-Enhancer Genie: An accessible resource of accurate enhancer predictions in the developing limb.

Epigenomic mapping of enhancer-associated chromatin modifications facilitates the genome-wide discovery of tissue-specific enhancers in vivo. However, reliance on single chromatin marks leads to high rates of false-positive predictions. More sophisticated, integrative methods have been described, but commonly suffer from limited accessibility to the resulting predictions and reduced biological interpretability. Here we present the Limb-Enhancer Genie (LEG), a collection of highly accurate, genome-wide predictions of enhancers in the developing limb, available through a user-friendly online interface. We predict limb enhancers using a combination of >50 published limb-specific datasets and clusters of evolutionarily conserved transcription factor binding sites, taking advantage of the patterns observed at previously in vivo validated elements. By combining different statistical models, our approach outperforms current state-of-the-art methods and provides interpretable measures of feature importance. Our results indicate that including a previously unappreciated score that quantifies tissue-specific nuclease accessibility significantly improves prediction performance. We demonstrate the utility of our approach through in vivo validation of newly predicted elements. Moreover, we describe general features that can guide the type of datasets to include when predicting tissue-specific enhancers genome-wide, while providing an accessible resource to the general biological community and facilitating the functional interpretation of genetic studies of limb malformations.

Islet1 regulates establishment of the posterior hindlimb field upstream of the Hand2-Shh morphoregulatory gene network in mouse embryos.

How divergent genetic systems regulate a common pathway during the development of two serial structures, forelimbs and hindlimbs, is not well understood. Specifically, HAND2 has been shown to regulate Shh directly to initiate its expression in the posterior margin of the limb mesenchyme. Although the Hand2-Shh morphoregulatory system operates in both the forelimb and hindlimb bud, a recent analysis suggested that its upstream regulation is different in the forelimb and hindlimb bud. A combination of all four Hox9 genes is required for Hand2 expression in the forelimb-forming region; however, it remains elusive what genetic system regulates the Hand2-Shh pathway in the hindlimb-forming region. By conditional inactivation of Islet1 in the hindlimb-forming region using the Hoxb6Cre transgene, we show that Islet1 is required for establishing the posterior hindlimb field, but not the forelimb field, upstream of the Hand2-Shh pathway. Inactivation of Islet1 caused the loss of posterior structures in the distal and proximal regions, specifically in the hindlimb. We found that Hand2 expression was downregulated in the hindlimb field and that Shh expression was severely impaired in the hindlimb bud. In the Hoxb6Cre; Islet1 mutant pelvis, the proximal element that is formed in a Shh-independent manner, displayed complementary defects in comparison with Pitx1(-/-) hindlimbs. This suggests that Islet1 and Pitx1 function in parallel during girdle development in hindlimbs, which is in contrast with the known requirement for Tbx5 in girdle development in forelimbs. Our studies have identified a role for Islet1 in hindlimb-specific development and have revealed Islet1 functions in two distinct processes: regulation upstream of the Hand2-Shh pathway and contributions to girdle development.

Dynamic enhancer landscapes in human craniofacial development.

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

Dual RMCE for efficient re-engineering of mouse mutant alleles.

We have developed dual recombinase-mediated cassette exchange (dRMCE) to efficiently re-engineer the thousands of available conditional alleles in mouse embryonic stem cells. dRMCE takes advantage of the wild-type loxP and FRT sites present in these conditional alleles and in many gene-trap lines. dRMCE is a scalable, flexible tool to introduce tags, reporters and mutant coding regions into an endogenous locus of interest in an easy and highly efficient manner.

Distinct roles of Hand2 in initiating polarity and posterior Shh expression during the onset of mouse limb bud development.

The polarization of nascent embryonic fields and the endowment of cells with organizer properties are key to initiation of vertebrate organogenesis. One such event is antero-posterior (AP) polarization of early limb buds and activation of morphogenetic Sonic Hedgehog (SHH) signaling in the posterior mesenchyme, which in turn promotes outgrowth and specifies the pentadactylous autopod. Inactivation of the Hand2 transcriptional regulator from the onset of mouse forelimb bud development disrupts establishment of posterior identity and Shh expression, which results in a skeletal phenotype identical to Shh deficient limb buds. In wild-type limb buds, Hand2 is part of the protein complexes containing Hoxd13, another essential regulator of Shh activation in limb buds. Chromatin immunoprecipitation shows that Hand2-containing chromatin complexes are bound to the far upstream cis-regulatory region (ZRS), which is specifically required for Shh expression in the limb bud. Cell-biochemical studies indicate that Hand2 and Hoxd13 can efficiently transactivate gene expression via the ZRS, while the Gli3 repressor isoform interferes with this positive transcriptional regulation. Indeed, analysis of mouse forelimb buds lacking both Hand2 and Gli3 reveals the complete absence of antero-posterior (AP) polarity along the entire proximo-distal axis and extreme digit polydactyly without AP identities. Our study uncovers essential components of the transcriptional machinery and key interactions that set-up limb bud asymmetry upstream of establishing the SHH signaling limb bud organizer.

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.