Rapid and pervasive changes in genome-wide enhancer usage during mammalian development.

Enhancers are distal regulatory elements that can activate tissue-specific gene expression and are abundant throughout mammalian genomes. Although substantial progress has been made toward 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 midgestation across all tissues examined. The dynamic enhancer activities uncovered in this study illuminate rapid and pervasive temporal in vivo changes in enhancer usage that underlie processes central to development and disease.

Differential alternative splicing analysis links variation in ZRSR2 to a novel type of oral-facial-digital syndrome.

PURPOSE: Oral-facial-digital (OFD) syndromes are genetically heterogeneous developmental disorders, caused by pathogenic variants in genes involved in primary cilia formation and function. We identified a previously undescribed type of OFD with brain anomalies, ranging from alobar holoprosencephaly to pituitary anomalies, in 6 unrelated families. METHODS: Exome sequencing of affected probands was supplemented with alternative splicing analysis in patient and control lymphoblastoid and fibroblast cell lines, and primary cilia structure analysis in patient fibroblasts. RESULTS: In 1 family with 2 affected males, we identified a germline variant in the last exon of ZRSR2, NM_005089.4:c.1211_1212del NP_005080.1:p.(Gly404GlufsTer23), whereas 7 affected males from 5 unrelated families were hemizygous for the ZRSR2 variant NM_005089.4:c.1207_1208del NP_005080.1:p.(Arg403GlyfsTer24), either occurring de novo or inherited in an X-linked recessive pattern. ZRSR2, located on chromosome Xp22.2, encodes a splicing factor of the minor spliceosome complex, which recognizes minor introns, representing 0.35% of human introns. Patient samples showed significant enrichment of minor intron retention. Among differentially spliced targets are ciliopathy-related genes, such as TMEM107 and CIBAR1. Primary fibroblasts containing the NM_005089.4:c.1207_1208del ZRSR2 variant had abnormally elongated cilia, confirming an association between defective U12-type intron splicing, OFD and abnormal primary cilia formation. CONCLUSION: We introduce a novel type of OFD associated with elongated cilia and differential splicing of minor intron-containing genes due to germline variation in ZRSR2.

Periodontal Ligament Reactions to Orthodontic Force: A Transcriptomic Study on Maxillary and Mandibular Human Premolars.

AIMS: Orthodontic force (OF) induces a variety of reactions in the periodontal ligament (PDL) that could potentially account for individual variability regarding orthodontic tooth movement (OTM). This study investigates the transcriptomic profile of human PDL tissue subjected to OF in vivo for 7 and 28 days, additionally comparing the differences between maxillary and mandibular PDL. METHODS: Healthy patients requiring orthodontic premolar extractions were randomly assigned to one of three groups: control (CG) where no OF was applied, 7 days and 28 days, where premolars were extracted either 7 or 28 days after the application of a 50-100 g OF. Total RNA was extracted from the PDL tissue and analyzed via RNA-seq. Differentially expressed genes (DEGs) were identified using a false discovery rate and fold change threshold of < 0.05 and ≥ 1.5 respectively. Functional and Protein-Protein Interaction analysis were performed. RESULTS: After 7 days of OF, the reaction of PDL to OF is characterized by cell responses to stress, increased bone resorption, inflammation and immune response, and decreased bone formation. In contrast, after 28 days, bone regeneration is more prominent, and processes of bone homeostasis, immune response, and cell migration are present. The response of maxillary and mandibular PDL was different. Bone resorption was observed in the maxilla at 7 and 28 days, while in the mandible expression of cell proliferation and transcriptional activity were predominant after 28 days of OF. CONCLUSIONS: The early reaction of the PDL to OF corresponds with increased bone resorption and decreased bone formation. After 28 days, bone formation became more prominent. The maxillary and mandibular PDL present asynchronous responses during OTM. These findings enhance our comprehension of the mechanisms underlying the origin-specific responses of PDL to different lengths of OF, which is potentially relevant in the development of personalized therapeutic strategies.

How does orthodontic tooth movement influence the dental pulp? RNA-sequencing on human premolars.

OBJECTIVES: The objective of this study is to analyse the gene expression profile of the dental pulp (DP) of human premolars subjected to 7 and 28 days of orthodontic force (OF) in vivo by using RNA sequencing. The maxillary and mandibular DP were additionally compared. METHODS: Healthy patients requiring orthodontic premolar extractions were randomly assigned to one of the three groups: control (CG) where no OF was applied, 7 and 28 days, where premolars were extracted either 7 or 28 days after the application of a 50-100 g OF. Total RNA was extracted from the DP and analysed via RNA-seq. Differentially expressed genes (DEGs) were identified using a false discovery rate and fold change threshold of <0.05 and ≥1.5, respectively. Functional analysis was performed. RESULTS: After 7 days of OF, pulp reaction indicates immune response, hypoxia, DNA damage and epigenetic regulation. After 28 days, cell adhesion, migration, organization and tissue repair are evident. The maxillary and mandibular pulp tissues react differently to OF. The maxilla exhibits minimal alterations, mostly related to immune response at 7 days and tissue repair at 28 days, whereas the mandible shows mostly DNA damage and epigenetic regulation at 7 days and return to the original state at 28 days. CONCLUSIONS: This study demonstrates that the early reaction of the DP to OF is marked by immune response, hypoxia and DNA damage. In contrast, after 28 days, cell adhesion, migration, organization, tissue repair and dentine formation are observed. Maxillary and mandibular premolars react differently to OF: although the maxilla exhibits minimal alterations at both time points, the mandible mostly shows DNA damage, epigenetic regulation, and immune response at 7 days. These disparities could stem from different blood supplies or the lower maxillary bone density, potentially triggering faster biological changes. Our findings provide insights into the gene regulatory networks modulating DP response to OF.

Tissue-specific SMARCA4 binding at active and repressed regulatory elements during embryogenesis.

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 mid-gestation (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.

An etiologic regulatory mutation in IRF6 with loss- and gain-of-function effects.

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.

Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice.

Sequence polymorphisms in a 58-kilobase (kb) interval on chromosome 9p21 confer a markedly increased risk of coronary artery disease (CAD), the leading cause of death worldwide. The variants have a substantial effect on the epidemiology of CAD and other life-threatening vascular conditions because nearly one-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 70-kb non-coding interval on mouse chromosome 4 affects cardiac expression of neighbouring genes, as well as proliferation properties of vascular cells. Chr4(Delta70kb/Delta70kb) mice are viable, but show increased mortality both during development and as adults. Cardiac expression of two genes near the non-coding interval, Cdkn2a and Cdkn2b, is severely reduced in chr4(Delta70kb/Delta70kb) mice, indicating that distant-acting gene regulatory functions are located in the non-coding CAD risk interval. Allele-specific expression of Cdkn2b transcripts in heterozygous mice showed that the deletion affects expression through a cis-acting mechanism. Primary cultures of chr4(Delta70kb/Delta70kb) 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 has 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.

Identification of novel craniofacial regulatory domains located far upstream of SOX9 and disrupted in Pierre Robin sequence.

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 harboring 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 noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

Pleiotropic effects of trisomy and pharmacologic modulation on structural, functional, molecular, and genetic systems in a Down syndrome mouse model.

Down syndrome (DS) is characterized by skeletal and brain structural malformations, cognitive impairment, altered hippocampal metabolite concentration and gene expression imbalance. These alterations were usually investigated separately, and the potential rescuing effects of green tea extracts enriched in epigallocatechin-3-gallate (GTE-EGCG) provided disparate results due to different experimental conditions. We overcame these limitations by conducting the first longitudinal controlled experiment evaluating genotype and GTE-EGCG prenatal chronic treatment effects before and after treatment discontinuation. Our findings revealed that the Ts65Dn mouse model reflected the pleiotropic nature of DS, exhibiting brachycephalic skull, ventriculomegaly, neurodevelopmental delay, hyperactivity, and impaired memory robustness with altered hippocampal metabolite concentration and gene expression. GTE-EGCG treatment modulated most systems simultaneously but did not rescue DS phenotypes. On the contrary, the treatment exacerbated trisomic phenotypes including body weight, tibia microarchitecture, neurodevelopment, adult cognition, and metabolite concentration, not supporting the therapeutic use of GTE-EGCG as a prenatal chronic treatment. Our results highlight the importance of longitudinal experiments assessing the co-modulation of multiple systems throughout development when characterizing preclinical models in complex disorders and evaluating the pleiotropic effects and general safety of pharmacological treatments.

A liver enhancer in the fibrinogen gene cluster.

The plasma concentration of fibrinogen varies in the healthy human population between 1.5 and 3.5 g/L. Understanding the basis of this variability has clinical importance because elevated fibrinogen levels are associated with increased cardiovascular disease risk. To identify novel regulatory elements involved in the control of fibrinogen expression, we used sequence conservation and in silico-predicted regulatory potential to select 14 conserved noncoding sequences (CNCs) within the conserved block of synteny containing the fibrinogen locus. The regulatory potential of each CNC was tested in vitro using a luciferase reporter gene assay in fibrinogen-expressing hepatoma cell lines (HuH7 and HepG2). 4 potential enhancers were tested for their ability to direct enhanced green fluorescent protein expression in zebrafish embryos. CNC12, a sequence equidistant from the human fibrinogen alpha and beta chain genes, activates strong liver enhanced green fluorescent protein expression in injected embryos and their transgenic progeny. A transgenic assay in embryonic day 14.5 mouse embryos confirmed the ability of CNC12 to activate transcription in the liver. While additional experiments are necessary to prove the role of CNC12 in the regulation of fibrinogen, our study reveals a novel regulatory element in the fibrinogen locus that is active in the liver and may contribute to variable fibrinogen expression in humans.

Comparison of human dental tissue RNA extraction methods for RNA sequencing.

OBJECTIVE: The purpose of this study was to identify an efficient RNA extraction method for periodontal ligament (PDL) and dental pulp (DP) tissues to be used in RNA sequencing studies, given the increased use of these techniques in dental research and the lack of standard protocols. DESIGN: PDL and DP were harvested from extracted third molars. Total RNA was extracted with four RNA extraction kits. RNA concentration, purity and integrity were assessed by means of NanoDrop and Bioanalyzer and statistically compared. RESULTS: RNA from PDL was more likely to be degraded than that of DP. The TRIzol method yielded the highest RNA concentration from both tissues. All methods harvested RNA with A260/A280 close to 2.0 and with A260/A230 above 1.5, except for the A260/A230 from PDL obtained with the RNeasy Mini kit. For RNA integrity, the RNeasy Fibrous Tissue Mini kit yielded the highest RIN values and 28 S/18 S from PDL, while the RNeasy Mini kit obtained relatively high RIN values with an appropriate 28 S/18 S for DP. CONCLUSION: Significantly different results were obtained for PDL and DP when using the RNeasy Mini kit. The RNeasy Mini kit provided the highest RNA yields and quality for DP, while the RNeasy Fibrous Tissue Mini kit obtained the highest quality RNA from PDL.

Disease-causing 7.4 kb cis-regulatory deletion disrupting conserved non-coding sequences and their interaction with the FOXL2 promotor: implications for mutation screening.

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.

Fine tuning of craniofacial morphology by distant-acting enhancers.

The shape of the human face and skull is largely genetically determined. However, the genomic basis of craniofacial morphology is incompletely understood and hypothesized to involve protein-coding genes, as well as gene regulatory sequences. 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 spatially complex developmental expression patterns. Analysis of mouse lines in which individual craniofacial enhancers had been deleted revealed significant alterations of craniofacial shape, demonstrating the functional importance of enhancers in defining face and skull morphology. These results demonstrate that enhancers are involved in craniofacial development and suggest that enhancer sequence variation contributes to the diversity of human facial morphology.

A systematic enhancer screen using lentivector transgenesis identifies conserved and non-conserved functional elements at the Olig1 and Olig2 locus.

Finding sequences that control expression of genes is central to understanding genome function. Previous studies have used evolutionary conservation as an indicator of regulatory potential. Here, we present a method for the unbiased in vivo screen of putative enhancers in large DNA regions, using the mouse as a model. We cloned a library of 142 overlapping fragments from a 200 kb-long murine BAC in a lentiviral vector expressing LacZ from a minimal promoter, and used the resulting vectors to infect fertilized murine oocytes. LacZ staining of E11 embryos obtained by first using the vectors in pools and then testing individual candidates led to the identification of 3 enhancers, only one of which shows significant evolutionary conservation. In situ hybridization and 3C/4C experiments suggest that this enhancer, which is active in the neural tube and posterior diencephalon, influences the expression of the Olig1 and/or Olig2 genes. This work provides a new approach for the large-scale in vivo screening of transcriptional regulatory sequences, and further demonstrates that evolutionary conservation alone seems too limiting a criterion for the identification of enhancers.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Transcriptional activation by bidirectional RNA polymerase II elongation over a silent promoter.

Transcriptional interference denotes negative cis effects between promoters. Here, we show that promoters can also interact positively. Bidirectional RNA polymerase II (Pol II) elongation over the silent human endogenous retrovirus (HERV)-K 18 promoter (representative of 2.5 x 10(3) similar promoters genomewide) activates transcription. In tandem constructs, an upstream promoter activates HERV-K 18 transcription. This is abolished by inversion of the upstream promoter, or by insertion of a poly(A) signal between the promoters; transcription is restored by poly(A) signal mutants. TATA-box mutants in the upstream promoter reduce HERV-K 18 transcription. Experiments with the same promoters in a convergent orientation produce similar effects. A small promoter deletion partially restores HERV-K 18 activity, consistent with activation resulting from repressor repulsion by the elongating Pol II. Transcriptional elongation over this class of intragenic promoters will generate co-regulated sense-antisense transcripts, or, alternatively initiating transcripts, thus expanding the diversity and complexity of the human transcriptome.

Outcome of donor splice site mutations accounting for congenital afibrinogenemia reflects order of intron removal in the fibrinogen alpha gene (FGA).

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.