Summary Enhancers are distal regulatory elements that can activate tissue-specific gene expression and are abundant throughout mammalian genomes. While substantial progress has been made towards genome-wide annotation of mammalian enhancers, their temporal activity patterns and global contributions in the context of developmental in vivo processes remain poorly explored. Here we used epigenomic profiling for H3K27ac, a mark of active enhancers, coupled to transgenic mouse assays to examine the genome-wide utilization of enhancers in three different mouse tissues across seven developmental stages. The majority of the ~90,000 enhancers identified exhibited tightly temporally restricted predicted activity windows and were associated with stage-specific biological functions and regulatory pathways in individual tissues. Comparative genomic analysis revealed that evolutionary conservation of enhancers decreases following mid-gestation across all tissues examined. The dynamic enhancer activities uncovered in this study illuminate rapid and pervasive temporal in vivo changes in enhancer usage underlying processes central to development and disease.

Sequence polymorphisms in a 58kb interval on chromosome 9p21 confer a markedly increased risk for coronary artery disease (CAD), the leading cause of death worldwide 1,2. The variants have a substantial impact on the epidemiology of CAD and other life-threatening vascular conditions since nearly a quarter of Caucasians are homozygous for risk alleles. However, the risk interval is devoid of protein-coding genes and the mechanism linking the region to CAD risk has remained enigmatic. Here we show that deletion of the orthologous 70kb noncoding interval on mouse chromosome4 affects cardiac expression of neighboring genes, as well as proliferation properties of vascular cells. Chr4Δ70kb/Δ70kb mice are viable, but show increased mortality both during development and as adults. Cardiac expression of two genes near the noncoding interval, Cdkn2a and Cdkn2b, is severely reduced in chr4Δ70kb/Δ70kb mice, indicating that distant-acting gene regulatory functions are located in the noncoding CAD risk interval. Allele-specific expression of Cdkn2b transcripts in heterozygous mice revealed that the deletion affects expression through a cis-acting mechanism. Primary cultures of chr4Δ70kb/Δ70kb aortic smooth muscle cells exhibited excessive proliferation and diminished senescence, a cellular phenotype consistent with accelerated CAD pathogenesis. Taken together, our results provide direct evidence that the CAD risk interval plays a pivotal role in regulation of cardiac Cdkn2a/b expression and suggest that this region affects CAD progression by altering the dynamics of vascular cell proliferation.

The shape of the human face and skull is largely genetically determined, but the genetic drivers of craniofacial morphology remain poorly understood. Here we used a combination of epigenomic profiling, in vivo characterization of candidate enhancer sequences in transgenic mice, and targeted deletion experiments to examine the role of distant-acting enhancers in craniofacial development. We identified complex regulatory landscapes, consisting of enhancers that drive a remarkable spatial complexity of developmental expression patterns. Deletion of individual craniofacial enhancers from the mouse genome resulted in significant alterations of craniofacial shape, demonstrating their functional importance in defining face and skull morphology. These results demonstrate that enhancers play a pervasive role in mammalian craniofacial development and suggest that enhancer sequence variation contributes to human facial morphology.

The SMARCA4 (also known as BRG1 in humans) chromatin remodeling factor is critical for establishing lineage-specific chromatin states during early mammalian development. However, the role of SMARCA4 in tissue-specific gene regulation during embryogenesis remains poorly defined. To investigate the genome-wide binding landscape of SMARCA4 in differentiating tissues, we engineered a Smarca4 FLAG knock-in mouse line. Using ChIP-seq, we identified~51,000 SMARCA4-associated regions across six embryonic mouse tissues (forebrain, hindbrain, neural tube, heart, limb, and face) at midgestation (E11.5). The majority of these regions was distal from promoters and showed dynamic occupancy, with most distal SMARCA4 sites (73%) confined to a single or limited subset of tissues. To further characterize these regions, we profiled active and repressive histone marks in the same tissues and examined the intersection of informative chromatin states and SMARCA4 binding. This revealed distinct classes of distal SMARCA4-associated elements characterized by activating and repressive chromatin signatures that were associated with tissue-specific up-or down-regulation of gene expression and relevant active/repressed biological pathways. We further demonstrate the predicted active regulatory properties of SMARCA4-associated elements by retrospective analysis of tissue-specific enhancers and direct testing of SMARCA4-bound regions in transgenic mouse assays. Our results indicate a dual active/repressive function of SMARCA4 at distal regulatory sequences in vivo and support its role in tissue-specific gene regulation during embryonic development.

To date, the contribution of disrupted potentially cis-regulatory conserved non-coding sequences (CNCs) to human disease is most likely underestimated, as no systematic screens for putative deleterious variations in CNCs have been conducted. As a model for monogenic disease we studied the involvement of genetic changes of CNCs in the cis-regulatory domain of FOXL2 in blepharophimosis syndrome (BPES). Fifty-seven molecularly unsolved BPES patients underwent high-resolution copy number screening and targeted sequencing of CNCs. Apart from three larger distant deletions, a de novo deletion as small as 7.4 kb was found at 283 kb 5′ to FOXL2. The deletion appeared to be triggered by an H-DNA-induced double-stranded break (DSB). In addition, it disrupts a novel long non-coding RNA (ncRNA) PISRT1 and 8 CNCs. The regulatory potential of the deleted CNCs was substantiated by in vitro luciferase assays. Interestingly, Chromosome Conformation Capture (3C) of a 625 kb region surrounding FOXL2 in expressing cellular systems revealed physical interactions of three upstream fragments and the FOXL2 core promoter. Importantly, one of these contains the 7.4 kb deleted fragment. Overall, this study revealed the smallest distant deletion causing monogenic disease and impacts upon the concept of mutation screening in human disease and developmental disorders in particular.

DNA variation in Interferon Regulatory Factor 6 (IRF6) causes Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate (CLP). However, an etiologic variant in IRF6 has been found in only 70% of VWS families. To test whether DNA variants in regulatory elements cause VWS, we sequenced three conserved elements near IRF6 in 70 VWS families that lack an etiologic mutation within IRF6 exons. A rare mutation (350dupA) was found in a conserved IRF6 enhancer element (MCS9.7) in a Brazilian family. The 350dupA mutation abrogated the binding of p63 and E47 transcription factors to cis-overlapping motifs, and significantly disrupted enhancer activity in human cell cultures. Moreover, using a transgenic assay in mice, the 350dupA mutation disrupted the activation of MCS9.7 enhancer element and led to failure of lacZ expression in all head and neck pharyngeal arches. Interestingly, disruption of the p63 Motif1 and/or E47 binding sites by nucleotide substitution did not fully recapitulate the effect of the 350dupA mutation. Rather, we recognized that the 350dupA created a CAAAGT motif, a binding site for Lef1 protein. We showed that Lef1 binds to the mutated site and that overexpression of Lef1/β-Catenin chimeric protein repressed MCS9.7-350dupA enhancer activity. In conclusion, our data strongly suggest that 350dupA variant is an etiologic mutation in VWS patients and disrupts enhancer activity by a loss- and gain-of-function mechanism, and thus support the rationale for additional screening for regulatory mutations in patients with CLP.

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions and duplications within a ~2 Mb region upstream of SOX9 can recapitulate the CD-DSD phenotype fully or partially, suggesting the existence of an unusually large cis-regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ~1.2-1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harbouring candidate craniofacial regulatory elements. In parallel, ChIP-Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple non-coding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

Congenital afibrinogenemia (Mendelian Inheritance in Man #202400) is a rare, autosomal recessive disorder characterized by the complete absence of circulating fibrinogen. Our recent studies on the molecular basis of the disease showed that the most common genetic defect is a donor splice mutation in fibrinogen alpha gene (FGA) intron 4, IVS4؉1G>T. Two other FGA donor splice mutations, in intron 1 (IVS1؉3A>G) and intron 3 (IVS3؉1_؉4delGTAA), were identified in afibrinogenemia patients. Because it was impossible to directly study the effect of these mutations on mRNA splicing in patient hepatocytes, we used a transfected cell approach, which previously allowed us to show that the common IVS4 mutation causes afibrinogenemia due to the activation of multiple cryptic donor splice sites. In this study, analysis of the IVS3delGTAA mutation showed exon 3 skipping in 99% of transcripts and exons 2 and 3 skipping in 1% of transcripts. The different outcomes of these donor splice mutations appear to follow the model proposed in a study of fibrillar collagen genes, where donor splice mutations occurring in a rapidly spliced intron with respect to upstream introns lead in most cases to exon skipping, while mutations in later-spliced introns lead to intron inclusion or cryptic splice-site utilization. Indeed, we found that in FGA intron 3 was preferentially spliced first, followed by intron 2, intron 4, and intron 1. (Blood. 2003;