Statistical inference with a manifold-constrained RNA velocity model uncovers cell cycle speed modulations.

Across biological systems, cells undergo coordinated changes in gene expression, resulting in transcriptome dynamics that unfold within a low-dimensional manifold. While low-dimensional dynamics can be extracted using RNA velocity, these algorithms can be fragile and rely on heuristics lacking statistical control. Moreover, the estimated vector field is not dynamically consistent with the traversed gene expression manifold. To address these challenges, we introduce a Bayesian model of RNA velocity that couples velocity field and manifold estimation in a reformulated, unified framework, identifying the parameters of an explicit dynamical system. Focusing on the cell cycle, we implement VeloCycle to study gene regulation dynamics on one-dimensional periodic manifolds and validate its ability to infer cell cycle periods using live imaging. We also apply VeloCycle to reveal speed differences in regionally defined progenitors and Perturb-seq gene knockdowns. Overall, VeloCycle expands the single-cell RNA sequencing analysis toolkit with a modular and statistically consistent RNA velocity inference framework.

Single Live Cell Monitoring of Protein Turnover Reveals Intercellular Variability and Cell-Cycle Dependence of Degradation Rates.

Cells need to reliably control their proteome composition to maintain homeostasis and regulate growth. How protein synthesis and degradation interplay to control protein expression levels remains unclear. Here, we combined a tandem fluorescent timer and pulse-chase protein labeling to disentangle how protein synthesis and degradation control protein homeostasis in single live mouse embryonic stem cells. We discovered substantial cell-cycle dependence in protein synthesis rates and stabilization of a large number of proteins around cytokinesis. Protein degradation rates were highly variable between cells, co-varied within individual cells for different proteins, and were positively correlated with synthesis rates. This suggests variability in proteasome activity as an important source of global extrinsic noise in gene expression. Our approach paves the way toward understanding the complex interplay of synthesis and degradation processes in determining protein levels of individual mammalian cells.

Cross-regulatory circuits linking inflammation, high-fat diet, and the circadian clock.

Mammalian physiology resonates with the daily changes in the external environment, allowing processes such as rest-activity cycles, metabolism, and body temperature to synchronize with daily changes in the surroundings. Studies have identified the molecular underpinnings of robust oscillations in gene expression occurring over the 24-h day, but how acute or chronic perturbations modulate gene expression rhythms, physiology, and behavior is still relatively unknown. In this issue of Genes & Development, Hong and colleagues (pp. 1367-1379) studied how acute and chronic inflammation interacts with the circadian clock. They found that NF-κB signaling can modify chromatin states and modulate expression of genes in the core clock network as well as circadian locomotor behavior. Interestingly, a high-fat diet (HFD) fed to mice also triggers this inflammation pathway, suggesting that cross-regulatory circuits link inflammation, HFD, and the circadian clock.

Clock-dependent chromatin topology modulates circadian transcription and behavior.

The circadian clock in animals orchestrates widespread oscillatory gene expression programs, which underlie 24-h rhythms in behavior and physiology. Several studies have shown the possible roles of transcription factors and chromatin marks in controlling cyclic gene expression. However, how daily active enhancers modulate rhythmic gene transcription in mammalian tissues is not known. Using circular chromosome conformation capture (4C) combined with sequencing (4C-seq), we discovered oscillatory promoter-enhancer interactions along the 24-h cycle in the mouse liver and kidney. Rhythms in chromatin interactions were abolished in arrhythmic Bmal1 knockout mice. Deleting a contacted intronic enhancer element in the Cryptochrome 1 (Cry1) gene was sufficient to compromise the rhythmic chromatin contacts in tissues. Moreover, the deletion reduced the daily dynamics of Cry1 transcriptional burst frequency and, remarkably, shortened the circadian period of locomotor activity rhythms. Our results establish oscillating and clock-controlled promoter-enhancer looping as a regulatory layer underlying circadian transcription and behavior.

Comprehensive analysis of the circadian nuclear and cytoplasmic transcriptome in mouse liver.

In eukaryotes, RNA is synthesised in the nucleus, spliced, and exported to the cytoplasm where it is translated and finally degraded. Any of these steps could be subject to temporal regulation during the circadian cycle, resulting in daily fluctuations of RNA accumulation and affecting the distribution of transcripts in different subcellular compartments. Our study analysed the nuclear and cytoplasmic, poly(A) and total transcriptomes of mouse livers collected over the course of a day. These data provide a genome-wide temporal inventory of enrichment in subcellular RNA, and revealed specific signatures of splicing, nuclear export and cytoplasmic mRNA stability related to transcript and gene lengths. Combined with a mathematical model describing rhythmic RNA profiles, we could test the rhythmicity of export rates and cytoplasmic degradation rates of approximately 1400 genes. With nuclear export times usually much shorter than cytoplasmic half-lives, we found that nuclear export contributes to the modulation and generation of rhythmic profiles of 10% of the cycling nuclear mRNAs. This study contributes to a better understanding of the dynamic regulation of the transcriptome during the day-night cycle.

Oscillating and stable genome topologies underlie hepatic physiological rhythms during the circadian cycle.

The circadian clock drives extensive temporal gene expression programs controlling daily changes in behavior and physiology. In mouse liver, transcription factors dynamics, chromatin modifications, and RNA Polymerase II (PolII) activity oscillate throughout the 24-hour (24h) day, regulating the rhythmic synthesis of thousands of transcripts. Also, 24h rhythms in gene promoter-enhancer chromatin looping accompany rhythmic mRNA synthesis. However, how chromatin organization impinges on temporal transcription and liver physiology remains unclear. Here, we applied time-resolved chromosome conformation capture (4C-seq) in livers of WT and arrhythmic Bmal1 knockout mice. In WT, we observed 24h oscillations in promoter-enhancer loops at multiple loci including the core-clock genes Period1, Period2 and Bmal1. In addition, we detected rhythmic PolII activity, chromatin modifications and transcription involving stable chromatin loops at clock-output gene promoters representing key liver function such as glucose metabolism and detoxification. Intriguingly, these contacts persisted in clock-impaired mice in which both PolII activity and chromatin marks no longer oscillated. Finally, we observed chromatin interaction hubs connecting neighbouring genes showing coherent transcription regulation across genotypes. Thus, both clock-controlled and clock-independent chromatin topology underlie rhythmic regulation of liver physiology.

Systematic analysis of differential rhythmic liver gene expression mediated by the circadian clock and feeding rhythms.

The circadian clock and feeding rhythms are both important regulators of rhythmic gene expression in the liver. To further dissect the respective contributions of feeding and the clock, we analyzed differential rhythmicity of liver tissue samples across several conditions. We developed a statistical method tailored to compare rhythmic liver messenger RNA (mRNA) expression in mouse knockout models of multiple clock genes, as well as PARbZip output transcription factors (Hlf/Dbp/Tef). Mice were exposed to ad libitum or night-restricted feeding under regular light-dark cycles. During ad libitum feeding, genetic ablation of the core clock attenuated rhythmic-feeding patterns, which could be restored by the night-restricted feeding regimen. High-amplitude mRNA expression rhythms in wild-type livers were driven by the circadian clock, but rhythmic feeding also contributed to rhythmic gene expression, albeit with significantly lower amplitudes. We observed that Bmal1 and Cry1/2 knockouts differed in their residual rhythmic gene expression. Differences in mean expression levels between wild types and knockouts correlated with rhythmic gene expression in wild type. Surprisingly, in PARbZip knockout mice, the mean expression levels of PARbZip targets were more strongly impacted than their rhythms, potentially due to the rhythmic activity of the D-box-repressor NFIL3. Genes that lost rhythmicity in PARbZip knockouts were identified to be indirect targets. Our findings provide insights into the diurnal transcriptome in mouse liver as we identified the differential contributions of several core clock regulators. In addition, we gained more insights on the specific effects of the feeding-fasting cycle.

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.

Ribo-DT: An automated pipeline for inferring codon dwell times from ribosome profiling data.

Protein synthesis is an energy consuming process characterised as a pivotal and highly regulated step in gene expression. The net protein output is dictated by a combination of translation initiation, elongation and termination rates that have remained difficult to measure. Recently, the development of ribosome profiling has enabled the inference of translation parameters through modelling, as this method informs on the ribosome position along the mRNA. Here, we present an automated, reproducible and portable computational pipeline to infer relative single-codon and codon-pair dwell times as well as gene flux from raw ribosome profiling sequencing data. As a case study, we applied our workflow to a publicly available yeast ribosome profiling dataset consisting of 57 independent gene knockouts related to RNA and tRNA modifications. We uncovered the effects of those modifications on translation elongation and codon selection during decoding. In particular, knocking out mod5 and trm7 increases codon-specific dwell times which indicates their potential tRNA targets, and highlights effects of nucleotide modifications on ribosome decoding rate.

Robust landscapes of ribosome dwell times and aminoacyl-tRNAs in response to nutrient stress in liver.

Translation depends on messenger RNA (mRNA)-specific initiation, elongation, and termination rates. While translation elongation is well studied in bacteria and yeast, less is known in higher eukaryotes. Here we combined ribosome and transfer RNA (tRNA) profiling to investigate the relations between translation elongation rates, (aminoacyl-) tRNA levels, and codon usage in mammals. We modeled codon-specific ribosome dwell times from ribosome profiling, considering codon pair interactions between ribosome sites. In mouse liver, the model revealed site- and codon-specific dwell times that differed from those in yeast, as well as pairs of adjacent codons in the P and A site that markedly slow down or speed up elongation. While translation efficiencies vary across diurnal time and feeding regimen, codon dwell times were highly stable and conserved in human. Measured tRNA levels correlated with codon usage and several tRNAs showed reduced aminoacylation, which was conserved in fasted mice. Finally, we uncovered that the longest codon dwell times could be explained by aminoacylation levels or high codon usage relative to tRNA abundance.

Engineered signaling centers for the spatially controlled patterning of human pluripotent stem cells.

Signaling centers, localized groups of cells that secrete morphogens, play a key role in early development and organogenesis by orchestrating spatial cell fate patterning. Here we present a microfluidic approach that exposes human pluripotent stem cell (hPSC) colonies to spatiotemporally controlled morphogen gradients generated from artificial signaling centers. In response to a localized source of bone morphogenetic protein 4 (BMP4), hPSC colonies reproducibly break their intrinsic radial symmetry to produce distinct, axially arranged differentiation domains. Counteracting sources of the BMP antagonist NOGGIN enhance this spatial control of cell fate patterning. We also show how morphogen concentration and cell density affect the BMP response and germ layer patterning. These results demonstrate that the intrinsic capacity of stem cells for self-organization can be extrinsically controlled through the use of engineered signaling centers.

The circadian oscillator analysed at the single-transcript level.

The circadian clock is an endogenous and self-sustained oscillator that anticipates daily environmental cycles. While rhythmic gene expression of circadian genes is well-described in populations of cells, the single-cell mRNA dynamics of multiple core clock genes remain largely unknown. Here we use single-molecule fluorescence in situ hybridisation (smFISH) at multiple time points to measure pairs of core clock transcripts, Rev-erbα (Nr1d1), Cry1 and Bmal1, in mouse fibroblasts. The mean mRNA level oscillates over 24 h for all three genes, but mRNA numbers show considerable spread between cells. We develop a probabilistic model for multivariate mRNA counts using mixtures of negative binomials, which accounts for transcriptional bursting, circadian time and cell-to-cell heterogeneity, notably in cell size. Decomposing the mRNA variability into distinct noise sources shows that clock time contributes a small fraction of the total variability in mRNA number between cells. Thus, our results highlight the intrinsic biological challenges in estimating circadian phase from single-cell mRNA counts and suggest that circadian phase in single cells is encoded post-transcriptionally.

Sleep-wake-driven and circadian contributions to daily rhythms in gene expression and chromatin accessibility in the murine cortex.

The timing and duration of sleep results from the interaction between a homeostatic sleep-wake-driven process and a periodic circadian process, and involves changes in gene regulation and expression. Unraveling the contributions of both processes and their interaction to transcriptional and epigenomic regulatory dynamics requires sampling over time under conditions of unperturbed and perturbed sleep. We profiled mRNA expression and chromatin accessibility in the cerebral cortex of mice over a 3-d period, including a 6-h sleep deprivation (SD) on day 2. We used mathematical modeling to integrate time series of mRNA expression data with sleep-wake history, which established that a large proportion of rhythmic genes are governed by the homeostatic process with varying degrees of interaction with the circadian process, sometimes working in opposition. Remarkably, SD caused long-term effects on gene-expression dynamics, outlasting phenotypic recovery, most strikingly illustrated by a damped oscillation of most core clock genes, including Arntl/Bmal1, suggesting that enforced wakefulness directly impacts the molecular clock machinery. Chromatin accessibility proved highly plastic and dynamically affected by SD. Dynamics in distal regions, rather than promoters, correlated with mRNA expression, implying that changes in expression result from constitutively accessible promoters under the influence of enhancers or repressors. Serum response factor (SRF) was predicted as a transcriptional regulator driving immediate response, suggesting that SRF activity mirrors the build-up and release of sleep pressure. Our results demonstrate that a single, short SD has long-term aftereffects at the genomic regulatory level and highlights the importance of the sleep-wake distribution to diurnal rhythmicity and circadian processes.

Two distinct promoter architectures centered on dynamic nucleosomes control ribosomal protein gene transcription.

In yeast, ribosome production is controlled transcriptionally by tight coregulation of the 138 ribosomal protein genes (RPGs). RPG promoters display limited sequence homology, and the molecular basis for their coregulation remains largely unknown. Here we identify two prevalent RPG promoter types, both characterized by upstream binding of the general transcription factor (TF) Rap1 followed by the RPG-specific Fhl1/Ifh1 pair, with one type also binding the HMG-B protein Hmo1. We show that the regulatory properties of the two promoter types are remarkably similar, suggesting that they are determined to a large extent by Rap1 and the Fhl1/Ifh1 pair. Rapid depletion experiments allowed us to define a hierarchy of TF binding in which Rap1 acts as a pioneer factor required for binding of all other TFs. We also uncovered unexpected features underlying recruitment of Fhl1, whose forkhead DNA-binding domain is not required for binding at most promoters, and Hmo1, whose binding is supported by repeated motifs. Finally, we describe unusually micrococcal nuclease (MNase)-sensitive nucleosomes at all RPG promoters, located between the canonical +1 and -1 nucleosomes, which coincide with sites of Fhl1/Ifh1 and Hmo1 binding. We speculate that these "fragile" nucleosomes play an important role in regulating RPG transcriptional output.

Rhythms of the Genome: Circadian Dynamics from Chromatin Topology, Tissue-Specific Gene Expression, to Behavior.

Circadian rhythms in physiology and behavior evolved to resonate with daily cycles in the external environment. In mammals, organs orchestrate temporal physiology over the 24-h day, which requires extensive gene expression rhythms targeted to the right tissue. Although a core set of gene products oscillates across virtually all cell types, gene expression profiling across tissues over the 24-h day showed that rhythmic gene expression programs are tissue specific. We highlight recent progress in uncovering how the circadian clock interweaves with tissue-specific gene regulatory networks involving functions such as xenobiotic metabolism, glucose homeostasis, and sleep. This progress hinges on not only comprehensive experimental approaches but also computational methods for multivariate analysis of periodic functional genomics data. We emphasize dynamic chromatin interactions as a novel regulatory layer underlying circadian gene transcription, core clock functions, and ultimately behavior. Finally, we discuss perspectives on extending the knowledge of the circadian clock in animals to human chronobiology.

Transcription factor activity rhythms and tissue-specific chromatin interactions explain circadian gene expression across organs.

Temporal control of physiology requires the interplay between gene networks involved in daily timekeeping and tissue function across different organs. How the circadian clock interweaves with tissue-specific transcriptional programs is poorly understood. Here, we dissected temporal and tissue-specific regulation at multiple gene regulatory layers by examining mouse tissues with an intact or disrupted clock over time. Integrated analysis uncovered two distinct regulatory modes underlying tissue-specific rhythms: tissue-specific oscillations in transcription factor (TF) activity, which were linked to feeding-fasting cycles in liver and sodium homeostasis in kidney; and colocalized binding of clock and tissue-specific transcription factors at distal enhancers. Chromosome conformation capture (4C-seq) in liver and kidney identified liver-specific chromatin loops that recruited clock-bound enhancers to promoters to regulate liver-specific transcriptional rhythms. Furthermore, this looping was remarkably promoter-specific on the scale of less than 10 kilobases (kb). Enhancers can contact a rhythmic promoter while looping out nearby nonrhythmic alternative promoters, confining rhythmic enhancer activity to specific promoters. These findings suggest that chromatin folding enables the clock to regulate rhythmic transcription of specific promoters to output temporal transcriptional programs tailored to different tissues.

Modulation of transcriptional burst frequency by histone acetylation.

Many mammalian genes are transcribed during short bursts of variable frequencies and sizes that substantially contribute to cell-to-cell variability. However, which molecular mechanisms determine bursting properties remains unclear. To probe putative mechanisms, we combined temporal analysis of transcription along the circadian cycle with multiple genomic reporter integrations, using both short-lived luciferase live microscopy and single-molecule RNA-FISH. Using the Bmal1 circadian promoter as our model, we observed that rhythmic transcription resulted predominantly from variations in burst frequency, while the genomic position changed the burst size. Thus, burst frequency and size independently modulated Bmal1 transcription. We then found that promoter histone-acetylation level covaried with burst frequency, being greatest at peak expression and lowest at trough expression, while remaining unaffected by the genomic location. In addition, specific deletions of ROR-responsive elements led to constitutively elevated histone acetylation and burst frequency. We then investigated the suggested link between histone acetylation and burst frequency by dCas9p300-targeted modulation of histone acetylation, revealing that acetylation levels influence burst frequency more than burst size. The correlation between acetylation levels at the promoter and burst frequency was also observed in endogenous circadian genes and in embryonic stem cell fate genes. Thus, our data suggest that histone acetylation-mediated control of transcription burst frequency is a common mechanism to control mammalian gene expression.

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.

Transcriptional regulatory logic of the diurnal cycle in the mouse liver.

Many organisms exhibit temporal rhythms in gene expression that propel diurnal cycles in physiology. In the liver of mammals, these rhythms are controlled by transcription-translation feedback loops of the core circadian clock and by feeding-fasting cycles. To better understand the regulatory interplay between the circadian clock and feeding rhythms, we mapped DNase I hypersensitive sites (DHSs) in the mouse liver during a diurnal cycle. The intensity of DNase I cleavages cycled at a substantial fraction of all DHSs, suggesting that DHSs harbor regulatory elements that control rhythmic transcription. Using chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq), we found that hypersensitivity cycled in phase with RNA polymerase II (Pol II) loading and H3K27ac histone marks. We then combined the DHSs with temporal Pol II profiles in wild-type (WT) and Bmal1-/- livers to computationally identify transcription factors through which the core clock and feeding-fasting cycles control diurnal rhythms in transcription. While a similar number of mRNAs accumulated rhythmically in Bmal1-/- compared to WT livers, the amplitudes in Bmal1-/- were generally lower. The residual rhythms in Bmal1-/- reflected transcriptional regulators mediating feeding-fasting responses as well as responses to rhythmic systemic signals. Finally, the analysis of DNase I cuts at nucleotide resolution showed dynamically changing footprints consistent with dynamic binding of CLOCK:BMAL1 complexes. Structural modeling suggested that these footprints are driven by a transient heterotetramer binding configuration at peak activity. Together, our temporal DNase I mappings allowed us to decipher the global regulation of diurnal transcription rhythms in the 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.