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.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development.

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

The little skate genome and the evolutionary emergence of wing-like fins.

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

The Transcriptional Regulator Prdm1 Is Essential for the Early Development of the Sensory Whisker Follicle and Is Linked to the Beta-Catenin First Dermal Signal.

Prdm1 mutant mice are one of the rare mutant strains that do not develop whisker hair follicles while still displaying a pelage. Here, we show that Prdm1 is expressed at the earliest stage of whisker development in clusters of mesenchymal cells before placode formation. Its conditional knockout in the murine soma leads to the loss of expression of Bmp2, Shh, Bmp4, Krt17, Edar, and Gli1, though leaving the ß-catenin-driven first dermal signal intact. Furthermore, we show that Prdm1 expressing cells not only act as a signaling center but also as a multipotent progenitor population contributing to the several lineages of the adult whisker. We confirm by genetic ablation experiments that the absence of macro vibrissae reverberates on the organization of nerve wiring in the mystacial pads and leads to the reorganization of the barrel cortex. We demonstrate that Lef1 acts upstream of Prdm1 and identify a primate-specific deletion of a Lef1 enhancer named Leaf. This loss may have been significant in the evolutionary process, leading to the progressive defunctionalization and disappearance of vibrissae in primates.

Mammalian-specific ectodermal enhancers control the expression of Hoxc genes in developing nails and hair follicles.

Vertebrate Hox genes are critical for the establishment of structures during the development of the main body axis. Subsequently, they play important roles either in organizing secondary axial structures such as the appendages, or during homeostasis in postnatal stages and adulthood. Here, we set up to analyze their elusive function in the ectodermal compartment, using the mouse limb bud as a model. We report that the HoxC gene cluster was co-opted to be transcribed in the distal limb ectoderm, where it is activated following the rule of temporal colinearity. These ectodermal cells subsequently produce various keratinized organs such as nails or claws. Accordingly, deletion of the HoxC cluster led to mice lacking nails (anonychia), a condition stronger than the previously reported loss of function of Hoxc13, which is the causative gene of the ectodermal dysplasia 9 (ECTD9) in human patients. We further identified two mammalian-specific ectodermal enhancers located upstream of the HoxC gene cluster, which together regulate Hoxc gene expression in the hair and nail ectodermal organs. Deletion of these regulatory elements alone or in combination revealed a strong quantitative component in the regulation of Hoxc genes in the ectoderm, suggesting that these two enhancers may have evolved along with the mammalian taxon to provide the level of HOXC proteins necessary for the full development of hair and nail.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18.

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

Comparative Transcriptomics Analyses across Species, Organs, and Developmental Stages Reveal Functionally Constrained lncRNAs.

The functionality of long noncoding RNAs (lncRNAs) is disputed. In general, lncRNAs are under weak selective pressures, suggesting that the majority of lncRNAs may be nonfunctional. However, although some surveys showed negligible phenotypic effects upon lncRNA perturbation, key biological roles were demonstrated for individual lncRNAs. Most lncRNAs with proven functions were implicated in gene expression regulation, in pathways related to cellular pluripotency, differentiation, and organ morphogenesis, suggesting that functional lncRNAs may be more abundant in embryonic development, rather than in adult organs. To test this hypothesis, we perform a multidimensional comparative transcriptomics analysis, across five developmental time points (two embryonic stages, newborn, adult, and aged individuals), four organs (brain, kidney, liver, and testes), and three species (mouse, rat, and chicken). We find that, overwhelmingly, lncRNAs are preferentially expressed in adult and aged testes, consistent with the presence of permissive transcription during spermatogenesis. LncRNAs are often differentially expressed among developmental stages and are less abundant in embryos and newborns compared with adult individuals, in agreement with a requirement for tighter expression control and less tolerance for noisy transcription early in development. For differentially expressed lncRNAs, we find that the patterns of expression variation among developmental stages are generally conserved between mouse and rat. Moreover, lncRNAs expressed above noise levels in somatic organs and during development show higher evolutionary conservation, in particular, at their promoter regions. Thus, we show that functionally constrained lncRNA loci are enriched in developing organs, and we suggest that many of these loci may function in an RNA-independent manner.

The constrained architecture of mammalian Hox gene clusters.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

Control of growth and gut maturation by HoxD genes and the associated lncRNA Haglr.

During embryonic development, Hox genes participate in the building of a functional digestive system in metazoans, and genetic conditions involving these genes lead to important, sometimes lethal, growth retardation. Recently, this phenotype was obtained after deletion of Haglr, the Hoxd antisense growth-associated long noncoding RNA (lncRNA) located between Hoxd1 and Hoxd3 In this study, we have analyzed the function of Hoxd genes in delayed growth trajectories by looking at several nested targeted deficiencies of the mouse HoxD cluster. Mutant pups were severely stunted during the suckling period, but many recovered after weaning. After comparing seven distinct HoxD alleles, including CRISPR/Cas9 deletions involving Haglr, we identified Hoxd3 as the critical component for the gut to maintain milk-digestive competence. This essential function could be abrogated by the dominant-negative effect of HOXD10 as shown by a genetic rescue approach, thus further illustrating the importance of posterior prevalence in Hox gene function. A role for the lncRNA Haglr in the control of postnatal growth could not be corroborated.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci.

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Convergent evolution of complex regulatory landscapes and pleiotropy at Hox loci.

Hox genes are required during the morphogenesis of both vertebrate digits and external genitals. We investigated whether transcription in such distinct contexts involves a shared enhancer-containing landscape. We show that the same regulatory topology is used, yet with some tissue-specific enhancer-promoter interactions, suggesting the hijacking of a regulatory backbone from one context to the other. In addition, comparable organizations are observed at both HoxA and HoxD clusters, which separated through genome duplication in an ancestral invertebrate animal. We propose that this convergent regulatory evolution was triggered by the preexistence of some chromatin architecture, thus facilitating the subsequent recruitment of the appropriate transcription factors. Such regulatory topologies may have both favored and constrained the evolution of pleiotropic developmental loci in vertebrates.

A genetic approach to the recruitment of PRC2 at the HoxD locus.

Polycomb group (PcG) proteins are essential for the repression of key factors during early development. In Drosophila, the polycomb repressive complexes (PRC) associate with defined polycomb response DNA elements (PREs). In mammals, however, the mechanisms underlying polycomb recruitment at targeted loci are poorly understood. We have used an in vivo approach to identify DNA sequences of importance for the proper recruitment of polycomb proteins at the HoxD locus. We report that various genomic re-arrangements of the gene cluster do not strongly affect PRC2 recruitment and that relatively small polycomb interacting sequences appear necessary and sufficient to confer polycomb recognition and targeting to ectopic loci. In addition, a high GC content, while not sufficient to recruit PRC2, may help its local spreading. We discuss the importance of PRC2 recruitment over Hox gene clusters in embryonic stem cells, for their subsequent coordinated transcriptional activation during development.