Visualization and data exploration of chromosome conformation capture data using Voronoi diagrams with v3c-viz.

Chromosome conformation capture (3C) sequencing approaches, like Hi-C or micro-C, allow for an unbiased view of chromatin interactions. Most analysis methods rely on so-called interaction matrices, which are derived from counting read pairs in bins of fixed size. Here, we propose the Voronoi diagram, as implemented in Voronoi for chromosome conformation capture data visualization (v3c-viz) to visualize 3C data. The Voronoi diagram corresponds to an adaptive-binning strategy that adapts to the local densities of points. In this way, visualization of data obtained by moderate sequencing depth pinpoint many, if not most, interesting features such as high frequency contacts. The favorable visualization properties of the Voronoi diagram indicate that the Voronoi diagram as density estimator can be used to identify high frequency contacts at a resolution approaching the typical size of enhancers and promoters. v3c-viz is available at https://github.com/imbbLab/v3c-viz .

Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq.

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.

Androgen and glucocorticoid receptor direct distinct transcriptional programs by receptor-specific and shared DNA binding sites.

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

CRUP: a comprehensive framework to predict condition-specific regulatory units.

We present the software Condition-specific Regulatory Units Prediction (CRUP) to infer from epigenetic marks a list of regulatory units consisting of dynamically changing enhancers with their target genes. The workflow consists of a novel pre-trained enhancer predictor that can be reliably applied across cell types and species, solely based on histone modification ChIP-seq data. Enhancers are subsequently assigned to different conditions and correlated with gene expression to derive regulatory units. We thoroughly test and then apply CRUP to a rheumatoid arthritis model, identifying enhancer-gene pairs comprising known disease genes as well as new candidate genes.

Expanding the repertoire of glucocorticoid receptor target genes by engineering genomic response elements.

The glucocorticoid receptor (GR), a hormone-activated transcription factor, binds to a myriad of genomic binding sites yet seems to regulate a much smaller number of genes. Genome-wide analysis of GR binding and gene regulation has shown that the likelihood of GR-dependent regulation increases with decreased distance of its binding to the transcriptional start site of a gene. To test if we can adopt this knowledge to expand the repertoire of GR target genes, we used CRISPR/Cas-mediated homology-directed repair to add a single GR-binding site directly upstream of the transcriptional start site of each of four genes. To our surprise, we found that the addition of a single GR-binding site can be enough to convert a gene into a GR target. The gain of GR-dependent regulation was observed for two of four genes analyzed and coincided with acquired GR binding at the introduced binding site. However, the gene-specific gain of GR-dependent regulation could not be explained by obvious differences in chromatin accessibility between converted genes and their non-converted counterparts. Furthermore, by introducing GR-binding sequences with different nucleotide compositions, we show that activation can be facilitated by distinct sequences without obvious differences in activity between the GR-binding sequence variants we tested. The approach to use genome engineering to build genomic response elements facilitates the generation of cell lines with tailored repertoires of GR-responsive genes and a framework to test and refine our understanding of the cis-regulatory logic of gene regulation by testing if engineered response elements behave as predicted.

Retrograde BMP signaling activates neuronal gene expression through widespread deployment of a conserved BMP-responsive cis-regulatory activation element.

Retrograde Bone Morphogenetic Protein (BMP) signaling in neurons is essential for the differentiation and synaptic function of many neuronal subtypes. BMP signaling regulates these processes via Smad transcription factor activity, yet the scope and nature of Smad-dependent gene regulation in neurons are mostly unknown. Here, we applied a computational approach to predict Smad-binding cis-regulatory BMP-Activating Elements (BMP-AEs) in Drosophila, followed by transgenic in vivo reporter analysis to test their neuronal subtype enhancer activity in the larval central nervous system (CNS). We identified 34 BMP-AE-containing genomic fragments that are responsive to BMP signaling in neurons, and showed that the embedded BMP-AEs are required for this activity. RNA-seq analysis identified BMP-responsive genes in the CNS and revealed that BMP-AEs selectively enrich near BMP-activated genes. These data suggest that functional BMP-AEs control nearby BMP-activated genes, which we validated experimentally. Finally, we demonstrated that the BMP-AE motif mediates a conserved Smad-responsive function in the Drosophila and vertebrate CNS. Our results provide evidence that BMP signaling controls neuronal function by directly coordinating the expression of a battery of genes through widespread deployment of a conserved Smad-responsive cis-regulatory motif.

Fgfr1 conditional-knockout in neural crest cells induces heterotopic chondrogenesis and osteogenesis in mouse frontal bones.

Most facial bones, including frontal bones, are derived from neural crest cells through intramembranous ossification. Fibroblast growth factor receptor 1 (Fgfr1) plays a pivotal role in craniofacial bone development, and loss of Fgfr1 leads to cleft palate and facial cleft defects in newborn mice. However, the potential role of the Fgfr1 gene in neural crest cell-mediated craniofacial development remains unclear. To investigate the role of Fgfr1 in neural crest cells, we analyzed Wnt1-Cre;Fgfr1flox/flox mice. Our results show that specific knockout of Fgfr1 in neural crest cells induced heterotopic chondrogenesis and osteogenesis at the interface of the anterior portions of frontal bones. We observed that heterotopic bone formation continued through postnatal day 28, whereas heterotopic chondrogenesis lasted only through the embryonic period. In summary, our results indicate that loss of Fgfr1 in neural crest cells leads to heterotopic chondrogenesis and osteogenesis.

PHF13: A new player involved in RNA polymerase II transcriptional regulation and co-transcriptional splicing.

We recently identified PHF13 as an H3K4me2/3 chromatin reader and transcriptional co-regulator. We found that PHF13 interacts with RNAPIIS5P and PRC2 stabilizing their association with active and bivalent promoters. Furthermore, mass spectrometry analysis identified ∼50 spliceosomal proteins in PHF13s interactome. Here, we will discuss the potential role of PHF13 in RNAPII pausing and co-transcriptional splicing.

PHF13 is a molecular reader and transcriptional co-regulator of H3K4me2/3.

PHF13 is a chromatin affiliated protein with a functional role in differentiation, cell division, DNA damage response and higher chromatin order. To gain insight into PHF13's ability to modulate these processes, we elucidate the mechanisms targeting PHF13 to chromatin, its genome wide localization and its molecular chromatin context. Size exclusion chromatography, mass spectrometry, X-ray crystallography and ChIP sequencing demonstrate that PHF13 binds chromatin in a multivalent fashion via direct interactions with H3K4me2/3 and DNA, and indirectly via interactions with PRC2 and RNA PolII. Furthermore, PHF13 depletion disrupted the interactions between PRC2, RNA PolII S5P, H3K4me3 and H3K27me3 and resulted in the up and down regulation of genes functionally enriched in transcriptional regulation, DNA binding, cell cycle, differentiation and chromatin organization. Together our findings argue that PHF13 is an H3K4me2/3 molecular reader and transcriptional co-regulator, affording it the ability to impact different chromatin processes.

BMP-dependent gene repression cascade in Drosophila eggshell patterning.

Bone Morphogenetic Proteins (BMPs) signal by activating Smad transcription factors to control a number of decisions during animal development. In Drosophila, signaling by the BMP ligand Decapentaplegic (Dpp) involves the activity of brinker (brk) which, in most contexts, is repressed by Dpp. Brk encodes a transcription factor which represses BMP signaling output by antagonizing Smad-dependent target gene activation. Here, we study BMP-dependent gene regulation during Drosophila oogenesis by following the signal transmission from Dpp to its target broad (br), a gene with a crucial function in eggshell patterning. We identify regulatory sequences that account for expression of both brk and br, and connect these to the transcription factors of the pathway. We show that Dpp directly regulates brk transcription through Smad- and Schnurri (Shn)-dependent repression. Brk is epistatic to Dpp in br expression and activates br indirectly, through removal of a repressor, which is yet to be identified. Our work provides first cis-regulatory insights into transcriptional interpretation of BMP signaling in eggshell morphogenesis and defines a transcriptional cascade that connects Dpp to target gene regulation.

Dynamic model for the coordination of two enhancers of broad by EGFR signaling.

Although it is widely appreciated that a typical developmental control gene is regulated by multiple enhancers, coordination of enhancer activities remains poorly understood. We propose a mechanism for such coordination in Drosophila oogenesis, when the expression of the transcription factor Broad (BR) evolves from a uniform to a two-domain pattern that prefigures the formation of two respiratory eggshell appendages. This change reflects sequential activities of two enhancers of the br gene, early and late, both of which are controlled by the epidermal growth factor receptor (EGFR) pathway. The late enhancer controls br in the appendage-producing cells, but the function of the early enhancer remained unclear. We found that the early enhancer is essential for the activity of the late enhancer and induction of eggshell appendages. This requirement can be explained by a mechanism whereby the BR protein produced by the early enhancer protects the late enhancer from EGFR-dependent repression. We illustrate this complex mechanism using a computational model that correctly predicts the wild-type dynamics of BR expression and its response to genetic perturbations.

Transcriptional interpretation of the EGF receptor signaling gradient.

Epidermal growth factor receptor (EGFR) controls a wide range of developmental events, from body axes specification in insects to cardiac development in humans. During Drosophila oogenesis, a gradient of EGFR activation patterns the follicular epithelium. Multiple transcriptional targets of EGFR in this tissue have been identified, but their regulatory elements are essentially unknown. We report the regulatory elements of broad (br) and pipe (pip), two important targets of EGFR signaling in Drosophila oogenesis. br is expressed in a complex pattern that prefigures the formation of respiratory eggshell appendages. We found that this pattern is generated by dynamic activities of two regulatory elements, which display different responses to Pointed, Capicua, and Mirror, transcription factors involved in the EGFR-mediated gene expression. One of these elements is active in a pattern similar to pip, a gene repressed by EGFR and essential for establishing the dorsoventral polarity of the embryo. We demonstrate that this similarity of expression depends on a common sequence motif that binds Mirror in vitro and is essential for transcriptional repression in vivo.

Regulation of Tbx22 during facial and palatal development.

Mutations in the gene encoding the T-box transcription factor TBX22 cause X-linked cleft palate and ankyloglossia in humans. Here we show that Tbx22 expression during facial and palatal development is regulated by FGF and BMP signaling. Our results demonstrate that FGF8 induces Tbx22 in the early face while BMP4 represses and thus restricts its expression. This regulation is conserved between chicken and mouse, although the Tbx22-expression patterns differ considerably between these two species. We suggest that these species-specific differences may result at least in part from differences in the spatiotemporal patterns of BMP activity, but we exclude a direct repression of Tbx22 by the BMP-inducible transcriptional repressor MSX1. Together these findings help to integrate Tbx22 into the molecular network of factors regulating facial development.

Nuclear localization of long-VEGF is associated with hypoxia and tumor angiogenesis.

Vascular endothelial growth factor (VEGF) is a potent angiogenic factor that has a pivotal role in normal and pathological angiogenesis. VEGF has a long 5' untranslated region harboring an open reading frame (ORF) initiated by a CUG codon that is in-frame with the VEGF coding region. The ORF translation leads to the expression of a long isoform termed L-VEGF that is extended by an additional 180 amino acids. In this communication, we provide evidence that L-VEGF is subjected to proteolytic cleavage leading to the detachment of the 180 aa extension from the VEGF moiety. Using immunofluorescence staining, we show that upon hypoxia this 180 aa extension translocates to the nuclei of expressing cells. Accordingly, immunohistochemical staining of both normal and tumor tissue samples demonstrated restricted nuclear localization of the ORF, which was correlated with cytoplasmic localization of VEGF. This suggests that the 180 aa ORF is involved in VEGF-mediated angiogenic processes.