Craniosynostosis, defined by premature fusion of one or multiple cranial sutures, is a common congenital defect affecting more than 1/2000 infants and results in restricted brain expansion. Single gene mutations account for 15%-20% of cases, largely as part of a syndrome, but the majority are nonsyndromic with complex underlying genetics. We hypothesized that the two noncoding genomic regions identified by a GWAS for craniosynostosis contain distal regulatory elements for the risk genes BMPER and BMP2. To identify such regulatory elements, we surveyed conserved noncoding sequences from both risk loci for enhancer activity in transgenic Danio rerio. We identified enhancers from both regions that direct expression to skeletal tissues, consistent with the endogenous expression of bmper and bmp2. For each locus, we also found a skeletal enhancer that also contains a sequence variant associated with craniosynostosis risk. We examined the activity of each enhancer during craniofacial development and found that the BMPER-associated enhancer is active in the restricted region of cartilage closely associated with frontal bone initiation. The same enhancer is active in mouse skeletal tissues, demonstrating evolutionarily conserved activity. Using enhanced yeast one-hybrid assays, we identified transcription factors that bind each enhancer and observed differential binding between alleles, implicating multiple signaling pathways. Our findings help unveil the genetic mechanism of the two craniosynostosis risk loci. More broadly, our combined in vivo approach is applicable to many complex genetic diseases to build a link between association studies and specific genetic mechanisms.

Regulatory elements (enhancers) are major drivers of gene expression in mammals and harbor many genetic variants associated with human diseases. Here, we present an updated VISTA Enhancer Browser (https://enhancer.lbl.gov), a database of transgenic enhancer assays conducted in developing mouse embryos in vivo. Since the original publication in 2007, the database grew nearly 20-fold from 250 to over 4500 experiments and currently harbors over 23 500 images. The updated database provides structured information on experiments conducted at different stages of embryonic development, including enhancer activities of human pathogenic and synthetic variants and sequences derived from a variety of species. In addition to manually curated results of thousands of individual experiments, the new database also features hundreds of manually curated comparisons between alleles. The VISTA Enhancer Browser provides a crucial resource for study of human genetic variation, gene regulation and developmental biology.

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.