More than the clock: distinct regulation of muscle function and metabolism by PER2 and RORα.

Circadian rhythms, governed by the dominant central clock, in addition to various peripheral clocks, regulate almost all biological processes, including sleep-wake cycles, hormone secretion and metabolism. In certain contexts, the regulation and function of the peripheral oscillations can be decoupled from the central clock. However, the specific mechanisms underlying muscle-intrinsic clock-dependent modulation of muscle function and metabolism remain unclear. We investigated the outcome of perturbations of the primary and secondary feedback loops of the molecular clock in skeletal muscle by specific gene ablation of Period circadian regulator 2 (Per2) and RAR-related orphan receptor alpha (Rorα), respectively. In both models, a dampening of core clock gene oscillation was observed, while the phase was preserved. Moreover, both loops seem to be involved in the homeostasis of amine groups. Highly divergent outcomes were seen for overall muscle gene expression, primarily affecting circadian rhythmicity in the PER2 knockouts and non-oscillating genes in the RORα knockouts, leading to distinct outcomes in terms of metabolome and phenotype. These results highlight the entanglement of the molecular clock and muscle plasticity and allude to specific functions of different clock components, i.e. the primary and secondary feedback loops, in this context. The reciprocal interaction between muscle contractility and circadian clocks might therefore be instrumental to determining a finely tuned adaptation of muscle tissue to perturbations in health and disease. KEY POINTS: Specific perturbations of the primary and secondary feedback loop of the molecular clock result in specific outcomes on muscle metabolism and function. Ablation of Per2 (primary loop) or Rorα (secondary loop) blunts the amplitude of core clock genes, in absence of a shift in phase. Perturbation of the primary feedback loop by deletion of PER2 primarily affects muscle gene oscillation. Knockout of RORα and the ensuing modulation of the secondary loop results in the aberrant expression of a large number of non-clock genes and proteins. The deletion of PER2 and RORα affects muscle metabolism and contractile function in a circadian manner, highlighting the central role of the molecular clock in modulating muscle plasticity.

Transcriptomic, proteomic and phosphoproteomic underpinnings of daily exercise performance and zeitgeber activity of training in mouse muscle.

KEY POINTS: Maximal endurance performance is greater in the early daytime. Timed exercise differentially alters the muscle transcriptome and (phospho)-proteome. Early daytime exercise triggers energy provisioning and tissue regeneration. Early night-time exercise activates stress-related and catabolic pathways. Scheduled training has limited effects on the muscle and liver circadian clocks. ABSTRACT: Timed physical activity might potentiate the health benefits of training. The underlying signalling events triggered by exercise at different times of day are, however, poorly understood. Here, we found that time-dependent variations in maximal treadmill exercise capacity of naïve mice were associated with energy stores, mostly hepatic glycogen levels. Importantly, running at different times of day resulted in a vastly different activation of signalling pathways, e.g. related to stress response, vesicular trafficking, repair and regeneration. Second, voluntary wheel running at the opposite phase of the dark, feeding period surprisingly revealed a minimal zeitgeber (i.e. phase-shifting) effect of training on the muscle clock. This integrated study provides important insights into the circadian regulation of endurance performance and the control of the circadian clock by exercise. In future studies, these results could contribute to better understanding circadian aspects of training design in athletes and the application of chrono-exercise-based interventions in patients.

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.

Principles for circadian orchestration of metabolic pathways.

Circadian rhythms govern multiple aspects of animal metabolism. Transcriptome-, proteome- and metabolome-wide measurements have revealed widespread circadian rhythms in metabolism governed by a cellular genetic oscillator, the circadian core clock. However, it remains unclear if and under which conditions transcriptional rhythms cause rhythms in particular metabolites and metabolic fluxes. Here, we analyzed the circadian orchestration of metabolic pathways by direct measurement of enzyme activities, analysis of transcriptome data, and developing a theoretical method called circadian response analysis. Contrary to a common assumption, we found that pronounced rhythms in metabolic pathways are often favored by separation rather than alignment in the times of peak activity of key enzymes. This property holds true for a set of metabolic pathway motifs (e.g., linear chains and branching points) and also under the conditions of fast kinetics typical for metabolic reactions. By circadian response analysis of pathway motifs, we determined exact timing separation constraints on rhythmic enzyme activities that allow for substantial rhythms in pathway flux and metabolite concentrations. Direct measurements of circadian enzyme activities in mouse skeletal muscle confirmed that such timing separation occurs in vivo.

Diurnal variation in egg excretion by Heterakis gallinarum.

Periodicity in nematode egg excretion may be of evolutionary origin as it can favour dispersal of the eggs in the environment. We investigated whether egg excretion by Heterakis gallinarum shows a repeatable pattern of periodicity. The faecal egg concentration and total number of eggs excreted within 4-h intervals were significantly affected by the sampling time within 1 day, but remained unaffected by the sampling day or interaction effects. By contrast, the total number of eggs excreted within 24 h did not differ among the 4 days of the study, collectively indicating repeatable egg excretion patterns. Both host feces and parasite egg excretion increased from night to late afternoon, followed by a decrease in the evening, resulting in higher egg excretion during daytime than the dark period. Feces excretion and worm fecundity showed overlapping diurnal rhythms with similarly timed phases, suggesting the existence of synchronicity between the host feces and nematode egg excretion patterns. We conclude that egg excretion by H. gallinarum is synchronized with host feces excretion and is higher during the daytime than during the dark period. This overlaps with the maximum activity of the day-active host and allows a maximal dispersal of the eggs in the environment.

Linking Core Promoter Classes to Circadian Transcription.

Circadian rhythms in transcription are generated by rhythmic abundances and DNA binding activities of transcription factors. Propagation of rhythms to transcriptional initiation involves the core promoter, its chromatin state, and the basal transcription machinery. Here, I characterize core promoters and chromatin states of genes transcribed in a circadian manner in mouse liver and in Drosophila. It is shown that the core promoter is a critical determinant of circadian mRNA expression in both species. A distinct core promoter class, strong circadian promoters (SCPs), is identified in mouse liver but not Drosophila. SCPs are defined by specific core promoter features, and are shown to drive circadian transcriptional activities with both high averages and high amplitudes. Data analysis and mathematical modeling further provided evidence for rhythmic regulation of both polymerase II recruitment and pause release at SCPs. The analysis provides a comprehensive and systematic view of core promoters and their link to circadian mRNA expression in mouse and Drosophila, and thus reveals a crucial role for the core promoter in regulated, dynamic transcription.

Differential rhythmicity: detecting altered rhythmicity in biological data.

MOTIVATION: Biological rhythms, such as rhythms in gene expression controlled by the cell cycle or the circadian clock, are important in cell physiology. A common type of experiment compares rhythmicity in tissues or cells either kept under different conditions or having different genotypes. Such investigations provide insights into underlying mechanisms as well as functions of rhythms. RESULTS: We present and benchmark a set of statistical and computational methods for this type of analysis, here termed differential rhythmicity analysis. The methods detect alterations in rhythm amplitude, phase and signal to noise ratio in one set of measurements compared to another. Using these methods, we compared circadian rhythms in liver mRNA expression in mice held under two different lighting conditions: constant darkness and light-dark cycles, respectively. This analysis revealed widespread and reproducible amplitude increases in mice kept in light-dark cycles. Further analysis of the subset of differentially rhythmic transcripts implied the immune system in mediating ambient light-dark cycles to rhythmic transcriptional activities. The methods are suitable for genome- or proteome-wide studies, and provide rigorous P values against well-defined null hypotheses. AVAILABILITY AND IMPLEMENTATION: The methods were implemented as the accompanying R software package DODR, available on CRAN. CONTACT: pal-olof.westermark@charite.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Circadian mRNA expression: insights from modeling and transcriptomics.

Circadian clocks synchronize organisms to the 24 h rhythms of the environment. These clocks persist under constant conditions, have their origin at the cellular level, and produce an output of rhythmic mRNA expression affecting thousands of transcripts in many mammalian cell types. Here, we review the charting of circadian output rhythms in mRNA expression, focusing on mammals. We emphasize the challenges in statistics, interpretation, and quantitative descriptions that such investigations have faced and continue to face, and outline remaining outstanding questions.

Dynamics of fat cell turnover in humans.

Obesity is increasing in an epidemic manner in most countries and constitutes a public health problem by enhancing the risk for cardiovascular disease and metabolic disorders such as type 2 diabetes. Owing to the increase in obesity, life expectancy may start to decrease in developed countries for the first time in recent history. The factors determining fat mass in adult humans are not fully understood, but increased lipid storage in already developed fat cells (adipocytes) is thought to be most important. Here we show that adipocyte number is a major determinant for the fat mass in adults. However, the number of fat cells stays constant in adulthood in lean and obese individuals, even after marked weight loss, indicating that the number of adipocytes is set during childhood and adolescence. To establish the dynamics within the stable population of adipocytes in adults, we have measured adipocyte turnover by analysing the integration of 14C derived from nuclear bomb tests in genomic DNA. Approximately 10% of fat cells are renewed annually at all adult ages and levels of body mass index. Neither adipocyte death nor generation rate is altered in early onset obesity, suggesting a tight regulation of fat cell number in this condition during adulthood. The high turnover of adipocytes establishes a new therapeutic target for pharmacological intervention in obesity.

Tuning the mammalian circadian clock: robust synergy of two loops.

The circadian clock is accountable for the regulation of internal rhythms in most living organisms. It allows the anticipation of environmental changes during the day and a better adaptation of physiological processes. In mammals the main clock is located in the suprachiasmatic nucleus (SCN) and synchronizes secondary clocks throughout the body. Its molecular constituents form an intracellular network which dictates circadian time and regulates clock-controlled genes. These clock-controlled genes are involved in crucial biological processes including metabolism and cell cycle regulation. Its malfunction can lead to disruption of biological rhythms and cause severe damage to the organism. The detailed mechanisms that govern the circadian system are not yet completely understood. Mathematical models can be of great help to exploit the mechanism of the circadian circuitry. We built a mathematical model for the core clock system using available data on phases and amplitudes of clock components obtained from an extensive literature search. This model was used to answer complex questions for example: how does the degradation rate of Per affect the period of the system and what is the role of the ROR/Bmal/REV-ERB (RBR) loop? Our findings indicate that an increase in the RNA degradation rate of the clock gene Period (Per) can contribute to increase or decrease of the period--a consequence of a non-monotonic effect of Per transcript stability on the circadian period identified by our model. Furthermore, we provide theoretical evidence for a potential role of the RBR loop as an independent oscillator. We carried out overexpression experiments on members of the RBR loop which lead to loss of oscillations consistent with our predictions. These findings challenge the role of the RBR loop as a merely auxiliary loop and might change our view of the clock molecular circuitry and of the function of the nuclear receptors (REV-ERB and ROR) as a putative driving force of molecular oscillations.

Daily running enhances molecular and physiological circadian rhythms in skeletal muscle.

OBJECTIVE: Exercise is a critical component of a healthy lifestyle and a key strategy for the prevention and management of metabolic disease. Identifying molecular mechanisms underlying adaptation in response to chronic physical activity is of critical interest in metabolic physiology. Circadian rhythms broadly modulate metabolism, including muscle substrate utilization and exercise capacity. Here, we define the molecular and physiological changes induced across the daily cycle by voluntary low intensity daily exercise. METHODS: Wildtype C57BL6/J male and female mice were housed with or without access to a running wheel for six weeks. Maximum running speed was measured at four different zeitgeber times (ZTs, hours after lights on) using either electrical or manual stimulation to motivate continued running on a motorized treadmill. RNA isolated from plantaris muscles at six ZTs was sequenced to establish the impact of daily activity on genome-wide transcription. Patterns of gene expression were analyzed using Gene Set Enrichment Analysis (GSEA) and Detection of Differential Rhythmicity (DODR). Blood glucose, lactate, and ketones, and muscle and liver glycogen were measured before and after exercise. RESULTS: We demonstrate that the use of mild electrical shocks to motivate running negatively impacts maximum running speed in mice, and describe a manual method to motivate running in rodent exercise studies. Using this method, we show that time of day influences the increase in exercise capacity afforded by six weeks of voluntary wheel running: when maximum running speed is measured at the beginning of the nighttime active period in mice, there is no measurable benefit from a history of daily voluntary running, while maximum increase in performance occurs at the end of the night. We show that daily voluntary exercise dramatically remodels the murine muscle circadian transcriptome. Finally, we describe daily rhythms in carbohydrate metabolism associated with the time-dependent response to moderate daily exercise in mice. CONCLUSIONS: Collectively, these data indicate that chronic nighttime physical activity dramatically remodels daily rhythms of murine muscle gene expression, which in turn support daily fluctuations in exercise performance.

Coupling governs entrainment range of circadian clocks.

Circadian clocks are endogenous oscillators driving daily rhythms in physiology and behavior. Synchronization of these timers to environmental light-dark cycles ('entrainment') is crucial for an organism's fitness. Little is known about which oscillator qualities determine entrainment, i.e., entrainment range, phase and amplitude. In a systematic theoretical and experimental study, we uncovered these qualities for circadian oscillators in the suprachiasmatic nucleus (SCN-the master clock in mammals) and the lung (a peripheral clock): (i) the ratio between stimulus (zeitgeber) strength and oscillator amplitude and (ii) the rigidity of the oscillatory system (relaxation rate upon perturbation) determine entrainment properties. Coupling among oscillators affects both qualities resulting in increased amplitude and rigidity. These principles explain our experimental findings that lung clocks entrain to extreme zeitgeber cycles, whereas SCN clocks do not. We confirmed our theoretical predictions by showing that pharmacological inhibition of coupling in the SCN leads to larger ranges of entrainment. These differences between master and the peripheral clocks suggest that coupling-induced rigidity in the SCN filters environmental noise to create a robust circadian system.

Molecular insights into human daily behavior.

Human beings exhibit wide variation in their timing of daily behavior. We and others have suggested previously that such differences might arise because of alterations in the period length of the endogenous human circadian oscillator. Using dermal fibroblast cells from skin biopsies of 28 subjects of early and late chronotype (11 "larks" and 17 "owls"), we have studied the circadian period lengths of these two groups, as well as their ability to phase-shift and entrain to environmental and chemical signals. We find not only period length differences between the two classes, but also significant changes in the amplitude and phase-shifting properties of the circadian oscillator among individuals with identical "normal" period lengths. Mathematical modeling shows that these alterations could also account for the extreme behavioral phenotypes of these subjects. We conclude that human chronotype may be influenced not only by the period length of the circadian oscillator, but also by cellular components that affect its amplitude and phase. In many instances, these changes can be studied at the molecular level in primary dermal cells.

Simple Kinetic Models in Molecular Chronobiology.

Circadian rhythms are constituted by a complex dynamical system with intertwined feedback loops, molecular switches, and self-sustained oscillations. Mathematical modeling supports understanding available heterogeneous kinetic data, highlights basic mechanisms, and can guide experimental research. Here, we introduce the basic steps from a biological question to simple models providing insight into gene-regulatory mechanisms. We illustrate the general approach by three examples: modeling decay processes, clock-controlled genes, and self-sustained oscillations.

Quantification of circadian rhythms in single cells.

Bioluminescence techniques allow accurate monitoring of the circadian clock in single cells. We have analyzed bioluminescence data of Per gene expression in mouse SCN neurons and fibroblasts. From these data, we extracted parameters such as damping rate and noise intensity using two simple mathematical models, one describing a damped oscillator driven by noise, and one describing a self-sustained noisy oscillator. Both models describe the data well and enabled us to quantitatively characterize both wild-type cells and several mutants. It has been suggested that the circadian clock is self-sustained at the single cell level, but we conclude that present data are not sufficient to determine whether the circadian clock of single SCN neurons and fibroblasts is a damped or a self-sustained oscillator. We show how to settle this question, however, by testing the models' predictions of different phases and amplitudes in response to a periodic entrainment signal (zeitgeber).

A mesoscale model of G1/S phase transition in liver regeneration.

The liver regenerates and maintains its function and size after injury by counterbalancing cell death with compensatory cell division. During liver regeneration, injured sites release cytokines, which stimulate normally quiescent hepatocytes to re-enter cell division cycle. Using a mesoscale approach, we have implemented the first mathematical model that describes cytokine-induced dedifferentiation of hepatocytes and the subsequent initiation of DNA synthesis (G0/G1 and G1/S phase transitions of the cell cycle). The model accurately reproduces experimentally measured kinetics of various signaling intermediates and DNA synthesis in hepatocytes for varying degrees of liver damage, in both wild type and knockout backgrounds. Liver regeneration is known to be a robust process, as liver mass reconstitution still occurs in various knockout mice (albeit with different kinetics). We analyze the robustness of the model using methods of control analysis. Moreover, we discuss the system's bandpass filtering properties and delays, which arise from feedbacks and nested feed-forward loops.

Beta-TrCP1-mediated degradation of PERIOD2 is essential for circadian dynamics.

Regulated degradation of circadian clock proteins is a crucial step for rhythm generation per se but also for establishing a normal circadian period. Here, the authors show that the F-box protein beta-transducin repeat containing protein 1 (beta-TrCP1) as part of the E3 ubiquitin ligase complex is an essential component of the mammalian circadian oscillator. Down-regulation of endogenous beta-TrCP1 as well as expression of a dominant-negative form both result in lengthening of the circadian period in oscillating fibroblasts. These phenotypes are due to an impaired degradation of PERIOD (PER) proteins, since expression of beta-TrCP interaction-deficient PER2 variants--but not wild-type PER2--results in a dramatic stabilization of PER2 protein as well as in the disruption of circadian rhythmicity. Mathematical modeling conceptualizes the authors' findings and suggests that loss of sustained rhythmicity in cells with eliminated beta-TrCP-mediated PER2 degradation is due to excessive nuclear repression, a prediction they verified experimentally.

A mathematical model of the Pyrosequencing reaction system.

The Pyrosequencing technology is a newly developed DNA sequencing method that monitors DNA nucleotide incorporation in real-time. A set of coupled enzymatic reactions, together with bioluminescence, detects incorporated nucleotides in the form of light pulses, yielding a characteristic light profile. In this study, a biochemical model of the Pyrosequencing reaction system is suggested and implemented. The model is constructed utilizing an assumption of irreversible Michaelis-Menten rate equations and a constant incorporation efficiency. The kinetic parameters are studied and values are chosen to obtain as reliable simulation results as possible. The results presented here show strong resemblance with real experiments. The model is able to capture the dynamics of a single light pulse with great accuracy, as well as the overall characteristics of a whole pyrogram trade mark. The plus- and minus-shift effects observed in experiments are successfully reconstructed by two constant efficiency factors. Furthermore, pulse broadening can partly be explained by apyrase inhibition and successive dilution.

Detecting rhythms in time series with RAIN.

A fundamental problem in research on biological rhythms is that of detecting and assessing the significance of rhythms in large sets of data. Classic methods based on Fourier theory are often hampered by the complex and unpredictable characteristics of experimental and biological noise. Robust nonparametric methods are available but are limited to specific wave forms. We present RAIN, a robust nonparametric method for the detection of rhythms of prespecified periods in biological data that can detect arbitrary wave forms. When applied to measurements of the circadian transcriptome and proteome of mouse liver, the sets of transcripts and proteins with rhythmic abundances were significantly expanded due to the increased detection power, when we controlled for false discovery. Validation against independent data confirmed the quality of these results. The large expansion of the circadian mouse liver transcriptomes and proteomes reflected the prevalence of nonsymmetric wave forms and led to new conclusions about function. RAIN was implemented as a freely available software package for R/Bioconductor and is presently also available as a web interface.