Molecular Co-occupancy Identifies Transcription Factor Binding Cooperativity In Vivo.

Gene activation requires the cooperative activity of multiple transcription factors at cis-regulatory elements (CREs). Yet, most transcription factors have short residence time, questioning the requirement of their physical co-occupancy on DNA to achieve cooperativity. Here, we present a DNA footprinting method that detects individual molecular interactions of transcription factors and nucleosomes with DNA in vivo. We apply this strategy to quantify the simultaneous binding of multiple transcription factors on single DNA molecules at mouse CREs. Analysis of the binary occupancy patterns at thousands of motif combinations reveals that high DNA co-occupancy occurs for most types of transcription factors, in the absence of direct physical interaction, at sites of competition with nucleosomes. Perturbation of pairwise interactions demonstrates the function of molecular co-occupancy in binding cooperativity. Our results reveal the interactions regulating CREs at molecular resolution and identify DNA co-occupancy as a widespread cooperativity mechanism used by transcription factors to remodel chromatin.

Temporal specificity and heterogeneity of Drosophila immune cells.

Immune cells provide defense against non-self and have recently been shown to also play key roles in diverse processes such as development, metabolism, and tumor progression. The heterogeneity of Drosophila immune cells (hemocytes) remains an open question. Using bulk RNA sequencing, we find that the hemocytes display distinct features in the embryo, a closed and rapidly developing system, compared to the larva, which is exposed to environmental and metabolic challenges. Through single-cell RNA sequencing, we identify fourteen hemocyte clusters present in unchallenged larvae and associated with distinct processes, e.g., proliferation, phagocytosis, metabolic homeostasis, and humoral response. Finally, we characterize the changes occurring in the hemocyte clusters upon wasp infestation, which triggers the differentiation of a novel hemocyte type, the lamellocyte. This first molecular atlas of hemocytes provides insights and paves the way to study the biology of the Drosophila immune cells in physiological and pathological conditions.

Temperature regulates splicing efficiency of the cold-inducible RNA-binding protein gene Cirbp.

In mammals, body temperature fluctuates diurnally around a mean value of 36°C-37°C. Despite the small differences between minimal and maximal values, body temperature rhythms can drive robust cycles in gene expression in cultured cells and, likely, animals. Here we studied the mechanisms responsible for the temperature-dependent expression of cold-inducible RNA-binding protein (CIRBP). In NIH3T3 fibroblasts exposed to simulated mouse body temperature cycles, Cirbp mRNA oscillates about threefold in abundance, as it does in mouse livers. This daily mRNA accumulation cycle is directly controlled by temperature oscillations and does not depend on the cells' circadian clocks. Here we show that the temperature-dependent accumulation of Cirbp mRNA is controlled primarily by the regulation of splicing efficiency, defined as the fraction of Cirbp pre-mRNA processed into mature mRNA. As revealed by genome-wide "approach to steady-state" kinetics, this post-transcriptional mechanism is widespread in the temperature-dependent control of gene expression.

Regulation of transcription reactivation dynamics exiting mitosis.

Proliferating cells experience a global reduction of transcription during mitosis, yet their cell identity is maintained and regulatory information is propagated from mother to daughter cells. Mitotic bookmarking by transcription factors has been proposed as a potential mechanism to ensure the reactivation of transcription at the proper set of genes exiting mitosis. Recently, mitotic transcription and waves of transcription reactivation have been observed in synchronized populations of human hepatoma cells. However, the study did not consider that mitotic-arrested cell populations progressively desynchronize leading to measurements of gene expression on a mixture of cells at different internal cell-cycle times. Moreover, it is not well understood yet what is the precise role of mitotic bookmarking on mitotic transcription as well as on the transcription reactivation waves. Ultimately, the core gene regulatory network driving the precise transcription reactivation dynamics remains to be identified. To address these questions, we developed a mathematical model to correct for the progressive desynchronization of cells and estimate gene expression dynamics with respect to a cell-cycle pseudotime. Furthermore, we used a multiple linear regression model to infer transcription factor activity dynamics. Our analysis allows us to characterize waves of transcription factor activities exiting mitosis and predict a core gene regulatory network responsible of the transcription reactivation dynamics. Moreover, we identified more than 60 transcription factors that are highly active during mitosis and represent new candidates of mitotic bookmarking factors which could be relevant therapeutic targets to control cell proliferation.

Circadian clock-dependent and -independent posttranscriptional regulation underlies temporal mRNA accumulation in mouse liver.

The mammalian circadian clock coordinates physiology with environmental cycles through the regulation of daily oscillations of gene expression. Thousands of transcripts exhibit rhythmic accumulations across mouse tissues, as determined by the balance of their synthesis and degradation. While diurnally rhythmic transcription regulation is well studied and often thought to be the main factor generating rhythmic mRNA accumulation, the extent of rhythmic posttranscriptional regulation is debated, and the kinetic parameters (e.g., half-lives), as well as the underlying regulators (e.g., mRNA-binding proteins) are relatively unexplored. Here, we developed a quantitative model for cyclic accumulations of pre-mRNA and mRNA from total RNA-seq data, and applied it to mouse liver. This allowed us to identify that about 20% of mRNA rhythms were driven by rhythmic mRNA degradation, and another 15% of mRNAs regulated by both rhythmic transcription and mRNA degradation. The method could also estimate mRNA half-lives and processing times in intact mouse liver. We then showed that, depending on mRNA half-life, rhythmic mRNA degradation can either amplify or tune phases of mRNA rhythms. By comparing mRNA rhythms in wild-type and Bmal1-/- animals, we found that the rhythmic degradation of many transcripts did not depend on a functional BMAL1. Interestingly clock-dependent and -independent degradation rhythms peaked at distinct times of day. We further predicted mRNA-binding proteins (mRBPs) that were implicated in the posttranscriptional regulation of mRNAs, either through stabilizing or destabilizing activities. Together, our results demonstrate how posttranscriptional regulation temporally shapes rhythmic mRNA accumulation in mouse liver.

Circadian Dbp transcription relies on highly dynamic BMAL1-CLOCK interaction with E boxes and requires the proteasome.

The transcription factors BMAL1 and CLOCK drive the circadian transcription of clock and clock-controlled genes, such as Dbp. To investigate the kinetics of BMAL1 binding to target genes in real time, we generated a cell line harboring tandem arrays of Dbp repeats and monitored the binding of a fluorescent BMAL1 fusion protein to these arrays by time-lapse microscopy. BMAL1 occupancy at the Dbp locus was highly circadian and strictly dependent on CLOCK. Moreover, BMAL1-CLOCK associations with Dbp were extremely unstable and displayed stochastic, proteasome-dependent fluctuations. Proteasome inhibition prolonged the residence time of BMAL1-CLOCK but resulted in an immediate attenuation of Dbp transcription. In cells harboring a single Dbp-luciferase reporter gene copy, this silencing was shown to be caused by a decrease in both the frequencies and sizes of transcriptional bursts. Thus, BMAL1 and CLOCK may act as Kamikaze activators, in that they are rapidly degraded once bound to Dbp chromatin.

Revealing Assembly of a Pore-Forming Complex Using Single-Cell Kinetic Analysis and Modeling.

Many biological processes depend on the sequential assembly of protein complexes. However, studying the kinetics of such processes by direct methods is often not feasible. As an important class of such protein complexes, pore-forming toxins start their journey as soluble monomeric proteins, and oligomerize into transmembrane complexes to eventually form pores in the target cell membrane. Here, we monitored pore formation kinetics for the well-characterized bacterial pore-forming toxin aerolysin in single cells in real time to determine the lag times leading to the formation of the first functional pores per cell. Probabilistic modeling of these lag times revealed that one slow and seven equally fast rate-limiting reactions best explain the overall pore formation kinetics. The model predicted that monomer activation is the rate-limiting step for the entire pore formation process. We hypothesized that this could be through release of a propeptide and indeed found that peptide removal abolished these steps. This study illustrates how stochasticity in the kinetics of a complex process can be exploited to identify rate-limiting mechanisms underlying multistep biomolecular assembly pathways.

Structure of silent transcription intervals and noise characteristics of mammalian genes.

Mammalian transcription occurs stochastically in short bursts interspersed by silent intervals showing a refractory period. However, the underlying processes and consequences on fluctuations in gene products are poorly understood. Here, we use single allele time-lapse recordings in mouse cells to identify minimal models of promoter cycles, which inform on the number and durations of rate-limiting steps responsible for refractory periods. The structure of promoter cycles is gene specific and independent of genomic location. Typically, five rate-limiting steps underlie the silent periods of endogenous promoters, while minimal synthetic promoters exhibit only one. Strikingly, endogenous or synthetic promoters with TATA boxes show simplified two-state promoter cycles. Since transcriptional bursting constrains intrinsic noise depending on the number of promoter steps, this explains why TATA box genes display increased intrinsic noise genome-wide in mammals, as revealed by single-cell RNA-seq. These findings have implications for basic transcription biology and shed light on interpreting single-cell RNA-counting experiments.

Stimulus-induced modulation of transcriptional bursting in a single mammalian gene.

Mammalian genes are often transcribed discontinuously as short bursts of RNA synthesis followed by longer silent periods. However, how these "on" and "off" transitions, together with the burst sizes, are modulated in single cells to increase gene expression upon stimulation is poorly characterized. By combining single-cell time-lapse luminescence imaging with stochastic modeling of the time traces, we quantified the transcriptional responses of the endogenous connective tissue growth factor gene to different physiological stimuli: serum and TGF-ß1. Both stimuli caused a rapid and acute increase in burst sizes. Whereas TGF-ß1 showed prolonged transcriptional activation mediated by an increase of transcription rate, serum stimulation resulted in a large and temporally tight first transcriptional burst, followed by a refractory period in the range of hours. Our study thus reveals how different physiological stimuli can trigger kinetically distinct transcriptional responses of the same gene.

Gene-specific RNA homeostasis revealed by perturbation of coactivator complexes.

Transcript buffering entails the reciprocal modulation of mRNA synthesis and degradation rates to maintain stable RNA levels under varying cellular conditions. Current research supports a global, non-sequence-specific connection between mRNA synthesis and degradation, but the underlying mechanisms are still unclear. In this study, we investigated changes in RNA metabolism following acute depletion of TIP60/KAT5, the acetyltransferase subunit of the NuA4 transcriptional coactivator complex, in mouse embryonic stem cells. By combining RNA sequencing of nuclear, cytoplasmic, and newly synthesised transcript fractions with biophysical modelling, we demonstrate that TIP60 predominantly enhances transcription of numerous genes, while a smaller set of genes undergoes TIP60-dependent transcriptional repression. Surprisingly, transcription changes caused by TIP60 depletion were offset by corresponding changes in RNA nuclear export and cytoplasmic stability, indicating gene-specific buffering mechanisms. Similarly, disruption of the unrelated ATAC coactivator complex also resulted in gene-specific transcript buffering. These findings reveal that transcript buffering functions at a gene-specific level and suggest that cells dynamically adjust RNA splicing, export, and degradation in response to individual RNA synthesis alterations, thereby sustaining cellular homeostasis.

MotEvo: integrated Bayesian probabilistic methods for inferring regulatory sites and motifs on multiple alignments of DNA sequences.

MOTIVATION: Probabilistic approaches for inferring transcription factor binding sites (TFBSs) and regulatory motifs from DNA sequences have been developed for over two decades. Previous work has shown that prediction accuracy can be significantly improved by incorporating features such as the competition of multiple transcription factors (TFs) for binding to nearby sites, the tendency of TFBSs for co-regulated TFs to cluster and form cis-regulatory modules and explicit evolutionary modeling of conservation of TFBSs across orthologous sequences. However, currently available tools only incorporate some of these features, and significant methodological hurdles hampered their synthesis into a single consistent probabilistic framework. RESULTS: We present MotEvo, a integrated suite of Bayesian probabilistic methods for the prediction of TFBSs and inference of regulatory motifs from multiple alignments of phylogenetically related DNA sequences, which incorporates all features just mentioned. In addition, MotEvo incorporates a novel model for detecting unknown functional elements that are under evolutionary constraint, and a new robust model for treating gain and loss of TFBSs along a phylogeny. Rigorous benchmarking tests on ChIP-seq datasets show that MotEvo's novel features significantly improve the accuracy of TFBS prediction, motif inference and enhancer prediction. AVAILABILITY: Source code, a user manual and files with several example applications are available at www.swissregulon.unibas.ch.

Characterization of cell-fate decision landscapes by estimating transcription factor dynamics.

Time-specific modulation of gene expression during differentiation by transcription factors promotes cell diversity. However, estimating their dynamic regulatory activity at the single-cell level and in a high-throughput manner remains challenging. We present FateCompass, an integrative approach that utilizes single-cell transcriptomics data to identify lineage-specific transcription factors throughout differentiation. By combining a probabilistic framework with RNA velocities or differentiation potential, we estimate transition probabilities, while a linear model of gene regulation is employed to compute transcription factor activities. Considering dynamic changes and correlations of expression and activities, FateCompass identifies lineage-specific regulators. Our validation using in silico data and application to pancreatic endocrine cell differentiation datasets highlight both known and potentially novel lineage-specific regulators. Notably, we uncovered undescribed transcription factors of an enterochromaffin-like population during in vitro differentiation toward ß-like cells. FateCompass provides a valuable framework for hypothesis generation, advancing our understanding of the gene regulatory networks driving cell-fate decisions.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics.

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the spatiotemporal arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find a close link between chromatin mobility and transcriptional status: active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

Competition between transcription and loop extrusion modulates promoter and enhancer dynamics.

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the 4D arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find that alterations in chromatin mobility, not promoter-enhancer distance, is more informative about transcriptional status. Active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

Scaling laws in functional genome content across prokaryotic clades and lifestyles.

For high-level functional categories that are represented in almost all prokaryotic genomes, the numbers of genes in these categories scale as power-laws in the total number of genes. We present a comprehensive analysis of the variation in these scaling laws across prokaryotic clades and lifestyles. For the large majority of functional categories, including transcription regulators, the inferred scaling laws are statistically indistinguishable across clades and lifestyles, supporting the simple hypothesis that these scaling laws are universally shared by all prokaryotes.

Cell cycle gene regulation dynamics revealed by RNA velocity and deep-learning.

Despite the fact that the cell cycle is a fundamental process of life, a detailed quantitative understanding of gene regulation dynamics throughout the cell cycle is far from complete. Single-cell RNA-sequencing (scRNA-seq) technology gives access to these dynamics without externally perturbing the cell. Here, by generating scRNA-seq libraries in different cell systems, we observe cycling patterns in the unspliced-spliced RNA space of cell cycle-related genes. Since existing methods to analyze scRNA-seq are not efficient to measure cycling gene dynamics, we propose a deep learning approach (DeepCycle) to fit these patterns and build a high-resolution map of the entire cell cycle transcriptome. Characterizing the cell cycle in embryonic and somatic cells, we identify major waves of transcription during the G1 phase and systematically study the stages of the cell cycle. Our work will facilitate the study of the cell cycle in multiple cellular models and different biological contexts.

Precise measurements of chromatin diffusion dynamics by modeling using Gaussian processes.

The spatiotemporal organization of chromatin influences many nuclear processes: from chromosome segregation to transcriptional regulation. To get a deeper understanding of these processes, it is essential to go beyond static viewpoints of chromosome structures, to accurately characterize chromatin's diffusion properties. We present GP-FBM: a computational framework based on Gaussian processes and fractional Brownian motion to extract diffusion properties from stochastic trajectories of labeled chromatin loci. GP-FBM uses higher-order temporal correlations present in the data, therefore, outperforming existing methods. Furthermore, GP-FBM allows to interpolate incomplete trajectories and account for substrate movement when two or more particles are present. Using our method, we show that average chromatin diffusion properties are surprisingly similar in interphase and mitosis in mouse embryonic stem cells. We observe surprising heterogeneity in local chromatin dynamics, correlating with potential regulatory activity. We also present GP-Tool, a user-friendly graphical interface to facilitate usage of GP-FBM by the research community.

First Responders Shape a Prompt and Sharp NF-κB-Mediated Transcriptional Response to TNF-α.

Nuclear factor (NF)-κB controls the transcriptional response to inflammatory signals by translocating into the nucleus, but we lack a single-cell characterization of the resulting transcription dynamics. Here we show that upon tumor necrosis factor (TNF)-α transcription of NF-κB target genes is heterogeneous in individual cells but results in an average nascent transcription profile that is prompt (i.e., occurs almost immediately) and sharp (i.e., increases and decreases rapidly) compared with NF-κB nuclear localization. Using an NF-κB-controlled MS2 reporter we show that the single-cell nascent transcription is more heterogeneous than NF-κB translocation dynamics, with a fraction of synchronized "first responders" that shape the average transcriptional profile and are more prone to respond to multiple TNF-α stimulations. A mathematical model combining NF-κB-mediated gene activation and a gene refractory state is able to reproduce these features. Our work shows how the expression of target genes induced by transcriptional activators can be heterogeneous across single cells and yet time resolved on average.

SwissRegulon: a database of genome-wide annotations of regulatory sites.

SwissRegulon (http://www.swissregulon.unibas.ch) is a database containing genome-wide annotations of regulatory sites in the intergenic regions of genomes. The regulatory site annotations are produced using a number of recently developed algorithms that operate on multiple alignments of orthologous intergenic regions from related genomes in combination with, whenever available, known sites from the literature, and ChIP-on-chip binding data. Currently SwissRegulon contains annotations for yeast and 17 prokaryotic genomes. The database provides information about the sequence, location, orientation, posterior probability and, whenever available, binding factor of each annotated site. To enable easy viewing of the regulatory site annotations in the context of other features annotated on the genomes, the sites are displayed using the GBrowse genome browser interface and can be queried based on any annotated genomic feature. The database can also be queried for regulons, i.e. sites bound by a common factor.

Visualization of Endogenous Transcription Factors in Single Cells Using an Antibody Electroporation-Based Imaging Approach.

In this chapter, we describe an antibody electroporation-based imaging approach that allows for precise imaging and quantification of endogenous transcription factor (i.e., RNA Polymerase II) distributions in single cells using 3D structured illumination microscopy (3D-SIM). The labeling is achieved by the efficient and harmless delivery of fluorescent dye-conjugated antibodies into living cells and the specific binding of these antibodies to the targeted factors. Our step-by-step protocol describes the procedure of the labeling of the specific antibodies, their electroporation into living cells, the sample preparation and 3D-SIM imaging as well as the postimaging analyses of the labeled endogenous transcription factors to obtain information about their nuclear distribution as well as their function. This protocol can be applied to a plethora of endogenous nuclear factors by using target specific noninhibiting antibodies.