The BTB-domain transcription factor ZBTB2 recruits chromatin remodelers and a histone chaperone during the exit from pluripotency.

Transcription factors (TFs) harboring broad-complex, tramtrack, and bric-a-brac (BTB) domains play important roles in development and disease. These BTB domains are thought to recruit transcriptional modulators to target DNA regions. However, a systematic molecular understanding of the mechanism of action of this TF family is lacking. Here, we identify the zinc finger BTB-TF Zbtb2 from a genetic screen for regulators of exit from pluripotency and demonstrate that its absence perturbs embryonic stem cell differentiation and the gene expression dynamics underlying peri-implantation development. We show that ZBTB2 binds the chromatin remodeler Ep400 to mediate downstream transcription. Independently, the BTB domain directly interacts with nucleosome remodeling and deacetylase and histone chaperone histone regulator A. Nucleosome remodeling and deacetylase recruitment is a common feature of BTB TFs, and based on phylogenetic analysis, we propose that this is a conserved evolutionary property. Binding to UBN2, in contrast, is specific to ZBTB2 and requires a C-terminal extension of the BTB domain. Taken together, this study identifies a BTB-domain TF that recruits chromatin modifiers and a histone chaperone during a developmental cell state transition and defines unique and shared molecular functions of the BTB-domain TF family.

A histone deacetylase adjusts transcription kinetics at coding sequences during Candida albicans morphogenesis.

Despite their classical role as transcriptional repressors, several histone deacetylases, including the baker's yeast Set3/Hos2 complex (Set3C), facilitate gene expression. In the dimorphic human pathogen Candida albicans, the homologue of the Set3C inhibits the yeast-to-filament transition, but the precise molecular details of this function have remained elusive. Here, we use a combination of ChIP-Seq and RNA-Seq to show that the Set3C acts as a transcriptional co-factor of metabolic and morphogenesis-related genes in C. albicans. Binding of the Set3C correlates with gene expression during fungal morphogenesis; yet, surprisingly, deletion of SET3 leaves the steady-state expression level of most genes unchanged, both during exponential yeast-phase growth and during the yeast-filament transition. Fine temporal resolution of transcription in cells undergoing this transition revealed that the Set3C modulates transient expression changes of key morphogenesis-related genes. These include a transcription factor cluster comprising of NRG1, EFG1, BRG1, and TEC1, which form a regulatory circuit controlling hyphal differentiation. Set3C appears to restrict the factors by modulating their transcription kinetics, and the hyperfilamentous phenotype of SET3-deficient cells can be reverted by mutating the circuit factors. These results indicate that the chromatin status at coding regions represents a dynamic platform influencing transcription kinetics. Moreover, we suggest that transcription at the coding sequence can be transiently decoupled from potentially conflicting promoter information in dynamic environments.

Identification of transcription factor binding sites from ChIP-seq data at high resolution.

MOTIVATION: Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) is widely used to study the in vivo binding sites of transcription factors (TFs) and their regulatory targets. Recent improvements to ChIP-seq, such as increased resolution, promise deeper insights into transcriptional regulation, yet require novel computational tools to fully leverage their advantages. RESULTS: To this aim, we have developed peakzilla, which can identify closely spaced TF binding sites at high resolution (i.e. resolves individual binding sites even if spaced closely), as we demonstrate using semisynthetic datasets, performing ChIP-seq for the TF Twist in Drosophila embryos with different experimental fragment sizes, and analyzing ChIP-exo datasets. We show that the increased resolution reached by peakzilla is highly relevant, as closely spaced Twist binding sites are strongly enriched in transcriptional enhancers, suggesting a signature to discriminate functional from abundant non-functional or neutral TF binding. Peakzilla is easy to use, as it estimates all the necessary parameters from the data and is freely available. AVAILABILITY AND IMPLEMENTATION: The peakzilla program is available from https://github.com/steinmann/peakzilla or http://www.starklab.org/data/peakzilla/. CONTACT: stark@starklab.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals.

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

E2F6 initiates stable epigenetic silencing of germline genes during embryonic development.

In mouse development, long-term silencing by CpG island DNA methylation is specifically targeted to germline genes; however, the molecular mechanisms of this specificity remain unclear. Here, we demonstrate that the transcription factor E2F6, a member of the polycomb repressive complex 1.6 (PRC1.6), is critical to target and initiate epigenetic silencing at germline genes in early embryogenesis. Genome-wide, E2F6 binds preferentially to CpG islands in embryonic cells. E2F6 cooperates with MGA to silence a subgroup of germline genes in mouse embryonic stem cells and in embryos, a function that critically depends on the E2F6 marked box domain. Inactivation of E2f6 leads to a failure to deposit CpG island DNA methylation at these genes during implantation. Furthermore, E2F6 is required to initiate epigenetic silencing in early embryonic cells but becomes dispensable for the maintenance in differentiated cells. Our findings elucidate the mechanisms of epigenetic targeting of germline genes and provide a paradigm for how transient repression signals by DNA-binding factors in early embryonic cells are translated into long-term epigenetic silencing during mouse development.

B7-H3 is a potent inhibitor of human T-cell activation: No evidence for B7-H3 and TREML2 interaction.

B7-H3 belongs to the B7 superfamily, a group of molecules that costimulate or down-modulate T-cell responses. Although it was shown that B7-H3 could inhibit T-cell responses, several studies - most of them performed in murine systems - found B7-H3 to act in a costimulatory manner. In this study, we have specifically addressed a potential functional dualism of human B7-H3 by assessing the effect of this molecule under varying experimental conditions as well as on different T-cell subsets. We show that B7-H3 does not costimulate human T cells. In the presence of strong activating signals, B7-H3 potently and consistently down-modulated human T-cell responses. This inhibitory effect was evident when analysing proliferation and cytokine production and affected naïve as well as pre-activated T cells. Furthermore, we demonstrate that B7-H3-T-cell interaction is characterised by an early suppression of IL-2 and that T-cell inhibition can be reverted by exogenous IL-2. Since the triggering receptor expressed on myeloid cells like transcript 2 (TREML2/TLT-2) has been recently described as costimulatory receptor of murine B7-H3 we have extensively analysed interaction of human B7-H3 with TREML2/TLT-2. In these experiments we found no evidence for such an interaction. Furthermore, our data do not point to a role for murine TREML2 as a receptor for murine B7-H3.

Genome-wide analysis in the mouse embryo reveals the importance of DNA methylation for transcription integrity.

Mouse embryos acquire global DNA methylation of their genome during implantation. However the exact roles of DNA methyltransferases (DNMTs) in embryos have not been studied comprehensively. Here we systematically analyze the consequences of genetic inactivation of Dnmt1, Dnmt3a and Dnmt3b on the methylome and transcriptome of mouse embryos. We find a strict division of function between DNMT1, responsible for maintenance methylation, and DNMT3A/B, solely responsible for methylation acquisition in development. By analyzing severely hypomethylated embryos, we uncover multiple functions of DNA methylation that is used as a mechanism of repression for a panel of genes including not only imprinted and germline genes, but also lineage-committed genes and 2-cell genes. DNA methylation also suppresses multiple retrotransposons and illegitimate transcripts from cryptic promoters in transposons and gene bodies. Our work provides a thorough analysis of the roles of DNA methyltransferases and the importance of DNA methylation for transcriptome integrity in mammalian embryos.

Loss of Apc Rapidly Impairs DNA Methylation Programs and Cell Fate Decisions in Lgr5+ Intestinal Stem Cells.

Colorectal cancer initiation and progression result from the accumulation of genetic and epigenetic alterations. Although aberrant gene expression and DNA methylation profiles are considered hallmarks of colorectal cancer development, the precise timing at which these are produced during tumor establishment remains elusive. Here we investigated the early transcriptional and epigenetic changes induced by adenomatous polyposis coli (Apc) inactivation in intestinal crypts. Hyperactivation of the Wnt pathway via Apc inactivation in crypt base columnar intestinal stem cells (ISC) led to their rapid accumulation driven by an impaired molecular commitment to differentiation, which was associated with discrete alterations in DNA methylation. Importantly, inhibiting the enzymes responsible for de novo DNA methylation restored the responsiveness of Apc-deficient intestinal organoids to stimuli regulating the proliferation-to-differentiation transition in ISC. This work reveals that early DNA methylation changes play critical roles in the establishment of the impaired fate decision program consecutive to Apc loss of function. SIGNIFICANCE: This study demonstrates the functional impact of changes in DNA methylation to determine the colorectal cancer cell phenotype following loss of Apc function.

cis-regulatory requirements for tissue-specific programs of the circadian clock.

BACKGROUND: Broadly expressed transcriptions factors (TFs) control tissue-specific programs of gene expression through interactions with local TF networks. A prime example is the circadian clock: although the conserved TFs CLOCK (CLK) and CYCLE (CYC) control a transcriptional circuit throughout animal bodies, rhythms in behavior and physiology are generated tissue specifically. Yet, how CLK and CYC determine tissue-specific clock programs has remained unclear. RESULTS: Here, we use a functional genomics approach to determine the cis-regulatory requirements for clock specificity. We first determine CLK and CYC genome-wide binding targets in heads and bodies by ChIP-seq and show that they have distinct DNA targets in the two tissue contexts. Computational dissection of CLK/CYC context-specific binding sites reveals sequence motifs for putative partner factors, which are predictive for individual binding sites. Among them, we show that the opa and GATA motifs, differentially enriched in head and body binding sites respectively, can be bound by OPA and SERPENT (SRP). They act synergistically with CLK/CYC in the Drosophila feedback loop, suggesting that they help to determine their direct targets and therefore orchestrate tissue-specific clock outputs. In addition, using in vivo transgenic assays, we validate that GATA motifs are required for proper tissue-specific gene expression in the adult fat body, midgut, and Malpighian tubules, revealing a cis-regulatory signature for enhancers of the peripheral circadian clock. CONCLUSIONS: Our results reveal how universal clock circuits can regulate tissue-specific rhythms and, more generally, provide insights into the mechanism by which universal TFs can be modulated to drive tissue-specific programs of gene expression.

A computational pipeline for comparative ChIP-seq analyses.

Chromatin immunoprecipitation (ChIP) followed by deep sequencing can now easily be performed across different conditions, time points and even species. However, analyzing such data is not trivial and standard methods are as yet unavailable. Here we present a protocol to systematically compare ChIP-sequencing (ChIP-seq) data across conditions. We first describe technical guidelines for data preprocessing, read mapping, read-density visualization and peak calling. We then describe methods and provide code with specific examples to compare different data sets across species and across conditions, including a threshold-free approach to measure global similarity, a strategy to assess the binary conservation of binding events and measurements for quantitative changes of binding. We discuss how differences in binding can be related to gene functions, gene expression and sequence changes. Once established, this protocol should take about 2 d to complete and be generally applicable to many data sets.

High conservation of transcription factor binding and evidence for combinatorial regulation across six Drosophila species.

The binding of some transcription factors has been shown to diverge substantially between closely related species. Here we show that the binding of the developmental transcription factor Twist is highly conserved across six Drosophila species, revealing strong functional constraints at its enhancers. Conserved binding correlates with sequence motifs for Twist and its partners, permitting the de novo discovery of their combinatorial binding. It also includes over 10,000 low-occupancy sites near the detection limit, which tend to mark enhancers of later developmental stages. These results suggest that developmental enhancers can be highly evolutionarily constrained, presumably because of their complex combinatorial nature.

Single amino acid repeats in signal peptides.

There has been an increasing interest in single amino acid repeats ever since it was shown that these are the cause of a variety of diseases. Although a systematic study of single amino acid repeats is challenging, they have subsequently been implicated in a number of functional roles. In general surveys, leucine runs were among the most frequent. In the present study, we present a detailed investigation of repeats in signal peptides of secreted and type I membrane proteins in comparison with their mature parts. We focus on eukaryotic species because single amino acid repeats are generally rather rare in archaea and bacteria. Our analysis of over 100 species shows that repeats of leucine (but not of other hydrophobic amino acids) are over-represented in signal peptides. This trend is most pronounced in higher eukaryotes, particularly in mammals. In the human proteome, although less than one-fifth of all proteins have a signal peptide, approximately two-thirds of all leucine repeats are located in these transient regions. Signal peptides are cleaved early from the growing polypeptide chain and then degraded rapidly. This may explain why leucine repeats, which can be toxic, are tolerated at such high frequencies. The substantial fraction of proteins affected by the strong enrichment of repeats in these transient segments highlights the bias that they can introduce for systematic analyses of protein sequences. In contrast to a general lack of conservation of single amino acid repeats, leucine repeats were found to be more conserved than the remaining signal peptide regions, indicating that they may have an as yet unknown functional role.