A spatio-temporal model to reveal oscillator phenotypes in molecular clocks: Parameter estimation elucidates circadian gene transcription dynamics in single-cells.

We propose a stochastic distributed delay model together with a Markov random field prior and a measurement model for bioluminescence-reporting to analyse spatio-temporal gene expression in intact networks of cells. The model describes the oscillating time evolution of molecular mRNA counts through a negative transcriptional-translational feedback loop encoded in a chemical Langevin equation with a probabilistic delay distribution. The model is extended spatially by means of a multiplicative random effects model with a first order Markov random field prior distribution. Our methodology effectively separates intrinsic molecular noise, measurement noise, and extrinsic noise and phenotypic variation driving cell heterogeneity, while being amenable to parameter identification and inference. Based on the single-cell model we propose a novel computational stability analysis that allows us to infer two key characteristics, namely the robustness of the oscillations, i.e. whether the reaction network exhibits sustained or damped oscillations, and the profile of the regulation, i.e. whether the inhibition occurs over time in a more distributed versus a more direct manner, which affects the cells' ability to phase-shift to new schedules. We show how insight into the spatio-temporal characteristics of the circadian feedback loop in the suprachiasmatic nucleus (SCN) can be gained by applying the methodology to bioluminescence-reported expression of the circadian core clock gene Cry1 across mouse SCN tissue. We find that while (almost) all SCN neurons exhibit robust cell-autonomous oscillations, the parameters that are associated with the regulatory transcription profile give rise to a spatial division of the tissue between the central region whose oscillations are resilient to perturbation in the sense that they maintain a high degree of synchronicity, and the dorsal region which appears to phase shift in a more diversified way as a response to large perturbations and thus could be more amenable to entrainment.

Filtering and inference for stochastic oscillators with distributed delays.

MOTIVATION: The time evolution of molecular species involved in biochemical reaction networks often arises from complex stochastic processes involving many species and reaction events. Inference for such systems is profoundly challenged by the relative sparseness of experimental data, as measurements are often limited to a small subset of the participating species measured at discrete time points. The need for model reduction can be realistically achieved for oscillatory dynamics resulting from negative translational and transcriptional feedback loops by the introduction of probabilistic time-delays. Although this approach yields a simplified model, inference is challenging and subject to ongoing research. The linear noise approximation (LNA) has recently been proposed to address such systems in stochastic form and will be exploited here. RESULTS: We develop a novel filtering approach for the LNA in stochastic systems with distributed delays, which allows the parameter values and unobserved states of a stochastic negative feedback model to be inferred from univariate time-series data. The performance of the methods is tested for simulated data. Results are obtained for real data when the model is fitted to imaging data on Cry1, a key gene involved in the mammalian central circadian clock, observed via a luciferase reporter construct in a mouse suprachiasmatic nucleus. AVAILABILITY AND IMPLEMENTATION: Programmes are written in MATLAB and Statistics Toolbox Release 2016 b, The MathWorks, Inc., Natick, Massachusetts, USA. Sample code and Cry1 data are available on GitHub https://github.com/scalderazzo/FLNADD. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Disentangling juxtacrine from paracrine signalling in dynamic tissue.

Prolactin is a major hormone product of the pituitary gland, the central endocrine regulator. Despite its physiological importance, the cell-level mechanisms of prolactin production are not well understood. Having significantly improved the resolution of real-time-single-cell-GFP-imaging, the authors recently revealed that prolactin gene transcription is highly dynamic and stochastic yet shows space-time coordination in an intact tissue slice. However, it still remains an open question as to what kind of cellular communication mediates the observed space-time organization. To determine the type of interaction between cells we developed a statistical model. The degree of similarity between two expression time series was studied in terms of two distance measures, Euclidean and geodesic, the latter being a network-theoretic distance defined to be the minimal number of edges between nodes, and this was used to discriminate between juxtacrine from paracrine signalling. The analysis presented here suggests that juxtacrine signalling dominates. To further determine whether the coupling is coordinating transcription or post-transcriptional activities we used stochastic switch modelling to infer the transcriptional profiles of cells and estimated their similarity measures to deduce that their spatial cellular coordination involves coupling of transcription via juxtacrine signalling. We developed a computational model that involves an inter-cell juxtacrine coupling, yielding simulation results that show space-time coordination in the transcription level that is in agreement with the above analysis. The developed model is expected to serve as the prototype for the further study of tissue-level organised gene expression for epigenetically regulated genes, such as prolactin.

Bayesian inference on stochastic gene transcription from flow cytometry data.

Motivation: Transcription in single cells is an inherently stochastic process as mRNA levels vary greatly between cells, even for genetically identical cells under the same experimental and environmental conditions. We present a stochastic two-state switch model for the population of mRNA molecules in single cells where genes stochastically alternate between a more active ON state and a less active OFF state. We prove that the stationary solution of such a model can be written as a mixture of a Poisson and a Poisson-beta probability distribution. This finding facilitates inference for single cell expression data, observed at a single time point, from flow cytometry experiments such as FACS or fluorescence in situ hybridization (FISH) as it allows one to sample directly from the equilibrium distribution of the mRNA population. We hence propose a Bayesian inferential methodology using a pseudo-marginal approach and a recent approximation to integrate over unobserved states associated with measurement error. Results: We provide a general inferential framework which can be widely used to study transcription in single cells from the kind of data arising in flow cytometry experiments. The approach allows us to separate between the intrinsic stochasticity of the molecular dynamics and the measurement noise. The methodology is tested in simulation studies and results are obtained for experimental multiple single cell expression data from FISH flow cytometry experiments. Availability and implementation: All analyses were implemented in R. Source code and the experimental data are available at https://github.com/SimoneTiberi/Bayesian-inference-on-stochastic-gene-transcription-from-flow-cytometry-data. Supplementary information: Supplementary data are available at Bioinformatics online.

Time-Series Transcriptomics Reveals That AGAMOUS-LIKE22 Affects Primary Metabolism and Developmental Processes in Drought-Stressed Arabidopsis.

In Arabidopsis thaliana, changes in metabolism and gene expression drive increased drought tolerance and initiate diverse drought avoidance and escape responses. To address regulatory processes that link these responses, we set out to identify genes that govern early responses to drought. To do this, a high-resolution time series transcriptomics data set was produced, coupled with detailed physiological and metabolic analyses of plants subjected to a slow transition from well-watered to drought conditions. A total of 1815 drought-responsive differentially expressed genes were identified. The early changes in gene expression coincided with a drop in carbon assimilation, and only in the late stages with an increase in foliar abscisic acid content. To identify gene regulatory networks (GRNs) mediating the transition between the early and late stages of drought, we used Bayesian network modeling of differentially expressed transcription factor (TF) genes. This approach identified AGAMOUS-LIKE22 (AGL22), as key hub gene in a TF GRN. It has previously been shown that AGL22 is involved in the transition from vegetative state to flowering but here we show that AGL22 expression influences steady state photosynthetic rates and lifetime water use. This suggests that AGL22 uniquely regulates a transcriptional network during drought stress, linking changes in primary metabolism and the initiation of stress responses.

Inferring transcriptional logic from multiple dynamic experiments.

MOTIVATION: The availability of more data of dynamic gene expression under multiple experimental conditions provides new information that makes the key goal of identifying not only the transcriptional regulators of a gene but also the underlying logical structure attainable. RESULTS: We propose a novel method for inferring transcriptional regulation using a simple, yet biologically interpretable, model to find the logic by which a set of candidate genes and their associated transcription factors (TFs) regulate the transcriptional process of a gene of interest. Our dynamic model links the mRNA transcription rate of the target gene to the activation states of the TFs assuming that these interactions are consistent across multiple experiments and over time. A trans-dimensional Markov Chain Monte Carlo (MCMC) algorithm is used to efficiently sample the regulatory logic under different combinations of parents and rank the estimated models by their posterior probabilities. We demonstrate and compare our methodology with other methods using simulation examples and apply it to a study of transcriptional regulation of selected target genes of Arabidopsis Thaliana from microarray time series data obtained under multiple biotic stresses. We show that our method is able to detect complex regulatory interactions that are consistent under multiple experimental conditions. AVAILABILITY AND IMPLEMENTATION: Programs are written in MATLAB and Statistics Toolbox Release 2016b, The MathWorks, Inc., Natick, Massachusetts, United States and are available on GitHub https://github.com/giorgosminas/TRS and at http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software. CONTACT: giorgos.minas@warwick.ac.uk or b.f.finkenstadt@warwick.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Relevance of a Mobile Internet Platform for Capturing Inter- and Intrasubject Variabilities in Circadian Coordination During Daily Routine: Pilot Study.

BACKGROUND: Experimental and epidemiologic studies have shown that circadian clocks' disruption can play an important role in the development of cancer and metabolic diseases. The cellular clocks outside the brain are effectively coordinated by the body temperature rhythm. We hypothesized that concurrent measurements of body temperature and rest-activity rhythms would assess circadian clocks coordination in individual patients, thus enabling the integration of biological rhythms into precision medicine. OBJECTIVE: The objective was to evaluate the circadian clocks' coordination in healthy subjects and patients through simultaneous measurements of rest-activity and body temperature rhythms. METHODS: Noninvasive real-time measurements of rest-activity and chest temperature rhythms were recorded during the subject's daily life, using a dedicated new mobile electronic health platform (PiCADo). It involved a chest sensor that jointly measured accelerations, 3D orientation, and skin surface temperature every 1-5 min and relayed them out to a mobile gateway via Bluetooth Low Energy. The gateway tele-transmitted all stored data to a server via General Packet Radio Service every 24 hours. The technical capabilities of PiCADo were validated in 55 healthy subjects and 12 cancer patients, whose rhythms were e-monitored during their daily routine for 3-30 days. Spectral analyses enabled to compute rhythm parameters values, with their 90% confidence limits, and their dynamics in each subject. RESULTS: All the individuals displayed a dominant circadian rhythm in activity with maxima occurring from 12:09 to 20:25. This was not the case for the dominant temperature period, which clustered around 24 hours for 51 out of 67 subjects (76%), and around 12 hours for 13 others (19%). Statistically significant sex- and age-related differences in circadian coordination were identified in the noncancerous subjects, based upon the range of variations in temperature rhythm amplitudes, maxima (acrophases), and phase relations with rest-activity. The circadian acrophase of chest temperature was located at night for the majority of people, but it occurred at daytime for 26% (14/55) of the noncancerous people and 33% (4/12) of the cancer patients, thus supporting important intersubject differences in circadian coordination. Sex, age, and cancer significantly impacted the circadian coordination of both rhythms, based on their phase relationships. CONCLUSIONS: Complementing rest-activity with chest temperature circadian e-monitoring revealed striking intersubject differences regarding human circadian clocks' coordination and timing during daily routine. To further delineate the clinical importance of such finding, the PiCADo platform is currently applied for both the assessment of health effects resulting from atypical work schedules and the identification of the key determinants of circadian disruption in cancer patients.

ReTrOS: a MATLAB toolbox for reconstructing transcriptional activity from gene and protein expression data.

BACKGROUND: Given the development of high-throughput experimental techniques, an increasing number of whole genome transcription profiling time series data sets, with good temporal resolution, are becoming available to researchers. The ReTrOS toolbox (Reconstructing Transcription Open Software) provides MATLAB-based implementations of two related methods, namely ReTrOS-Smooth and ReTrOS-Switch, for reconstructing the temporal transcriptional activity profile of a gene from given mRNA expression time series or protein reporter time series. The methods are based on fitting a differential equation model incorporating the processes of transcription, translation and degradation. RESULTS: The toolbox provides a framework for model fitting along with statistical analyses of the model with a graphical interface and model visualisation. We highlight several applications of the toolbox, including the reconstruction of the temporal cascade of transcriptional activity inferred from mRNA expression data and protein reporter data in the core circadian clock in Arabidopsis thaliana, and how such reconstructed transcription profiles can be used to study the effects of different cell lines and conditions. CONCLUSIONS: The ReTrOS toolbox allows users to analyse gene and/or protein expression time series where, with appropriate formulation of prior information about a minimum of kinetic parameters, in particular rates of degradation, users are able to infer timings of changes in transcriptional activity. Data from any organism and obtained from a range of technologies can be used as input due to the flexible and generic nature of the model and implementation. The output from this software provides a useful analysis of time series data and can be incorporated into further modelling approaches or in hypothesis generation.

A stochastic transcriptional switch model for single cell imaging data.

Gene expression is made up of inherently stochastic processes within single cells and can be modeled through stochastic reaction networks (SRNs). In particular, SRNs capture the features of intrinsic variability arising from intracellular biochemical processes. We extend current models for gene expression to allow the transcriptional process within an SRN to follow a random step or switch function which may be estimated using reversible jump Markov chain Monte Carlo (MCMC). This stochastic switch model provides a generic framework to capture many different dynamic features observed in single cell gene expression. Inference for such SRNs is challenging due to the intractability of the transition densities. We derive a model-specific birth-death approximation and study its use for inference in comparison with the linear noise approximation where both approximations are considered within the unifying framework of state-space models. The methodology is applied to synthetic as well as experimental single cell imaging data measuring expression of the human prolactin gene in pituitary cells.

Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis.

Transcriptional reprogramming forms a major part of a plant's response to pathogen infection. Many individual components and pathways operating during plant defense have been identified, but our knowledge of how these different components interact is still rudimentary. We generated a high-resolution time series of gene expression profiles from a single Arabidopsis thaliana leaf during infection by the necrotrophic fungal pathogen Botrytis cinerea. Approximately one-third of the Arabidopsis genome is differentially expressed during the first 48 h after infection, with the majority of changes in gene expression occurring before significant lesion development. We used computational tools to obtain a detailed chronology of the defense response against B. cinerea, highlighting the times at which signaling and metabolic processes change, and identify transcription factor families operating at different times after infection. Motif enrichment and network inference predicted regulatory interactions, and testing of one such prediction identified a role for TGA3 in defense against necrotrophic pathogens. These data provide an unprecedented level of detail about transcriptional changes during a defense response and are suited to systems biology analyses to generate predictive models of the gene regulatory networks mediating the Arabidopsis response to B. cinerea.

Inference on periodicity of circadian time series.

Estimation of the period length of time-course data from cyclical biological processes, such as those driven by the circadian pacemaker, is crucial for inferring the properties of the biological clock found in many living organisms. We propose a methodology for period estimation based on spectrum resampling (SR) techniques. Simulation studies show that SR is superior and more robust to non-sinusoidal and noisy cycles than a currently used routine based on Fourier approximations. In addition, a simple fit to the oscillations using linear least squares is available, together with a non-parametric test for detecting changes in period length which allows for period estimates with different variances, as frequently encountered in practice. The proposed methods are motivated by and applied to various data examples from chronobiology.

A temporal switch model for estimating transcriptional activity in gene expression.

MOTIVATION: The analysis and mechanistic modelling of time series gene expression data provided by techniques such as microarrays, NanoString, reverse transcription-polymerase chain reaction and advanced sequencing are invaluable for developing an understanding of the variation in key biological processes. We address this by proposing the estimation of a flexible dynamic model, which decouples temporal synthesis and degradation of mRNA and, hence, allows for transcriptional activity to switch between different states. RESULTS: The model is flexible enough to capture a variety of observed transcriptional dynamics, including oscillatory behaviour, in a way that is compatible with the demands imposed by the quality, time-resolution and quantity of the data. We show that the timing and number of switch events in transcriptional activity can be estimated alongside individual gene mRNA stability with the help of a Bayesian reversible jump Markov chain Monte Carlo algorithm. To demonstrate the methodology, we focus on modelling the wild-type behaviour of a selection of 200 circadian genes of the model plant Arabidopsis thaliana. The results support the idea that using a mechanistic model to identify transcriptional switch points is likely to strongly contribute to efforts in elucidating and understanding key biological processes, such as transcription and degradation.

A hierarchical model of transcriptional dynamics allows robust estimation of transcription rates in populations of single cells with variable gene copy number.

MOTIVATION: cis-regulatory DNA sequence elements, such as enhancers and silencers, function to control the spatial and temporal expression of their target genes. Although the overall levels of gene expression in large cell populations seem to be precisely controlled, transcription of individual genes in single cells is extremely variable in real time. It is, therefore, important to understand how these cis-regulatory elements function to dynamically control transcription at single-cell resolution. Recently, statistical methods have been proposed to back calculate the rates involved in mRNA transcription using parameter estimation of a mathematical model of transcription and translation. However, a major complication in these approaches is that some of the parameters, particularly those corresponding to the gene copy number and transcription rate, cannot be distinguished; therefore, these methods cannot be used when the copy number is unknown. RESULTS: Here, we develop a hierarchical Bayesian model to estimate biokinetic parameters from live cell enhancer-promoter reporter measurements performed on a population of single cells. This allows us to investigate transcriptional dynamics when the copy number is variable across the population. We validate our method using synthetic data and then apply it to quantify the function of two known developmental enhancers in real time and in single cells. AVAILABILITY: Supporting information is submitted with the article.

Network balance via CRY signalling controls the Arabidopsis circadian clock over ambient temperatures.

Circadian clocks exhibit 'temperature compensation', meaning that they show only small changes in period over a broad temperature range. Several clock genes have been implicated in the temperature-dependent control of period in Arabidopsis. We show that blue light is essential for this, suggesting that the effects of light and temperature interact or converge upon common targets in the circadian clock. Our data demonstrate that two cryptochrome photoreceptors differentially control circadian period and sustain rhythmicity across the physiological temperature range. In order to test the hypothesis that the targets of light regulation are sufficient to mediate temperature compensation, we constructed a temperature-compensated clock model by adding passive temperature effects into only the light-sensitive processes in the model. Remarkably, this model was not only capable of full temperature compensation and consistent with mRNA profiles across a temperature range, but also predicted the temperature-dependent change in the level of LATE ELONGATED HYPOCOTYL, a key clock protein. Our analysis provides a systems-level understanding of period control in the plant circadian oscillator.

Dynamic analysis of stochastic transcription cycles.

In individual mammalian cells the expression of some genes such as prolactin is highly variable over time and has been suggested to occur in stochastic pulses. To investigate the origins of this behavior and to understand its functional relevance, we quantitatively analyzed this variability using new mathematical tools that allowed us to reconstruct dynamic transcription rates of different reporter genes controlled by identical promoters in the same living cell. Quantitative microscopic analysis of two reporter genes, firefly luciferase and destabilized EGFP, was used to analyze the dynamics of prolactin promoter-directed gene expression in living individual clonal and primary pituitary cells over periods of up to 25 h. We quantified the time-dependence and cyclicity of the transcription pulses and estimated the length and variation of active and inactive transcription phases. We showed an average cycle period of approximately 11 h and demonstrated that while the measured time distribution of active phases agreed with commonly accepted models of transcription, the inactive phases were differently distributed and showed strong memory, with a refractory period of transcriptional inactivation close to 3 h. Cycles in transcription occurred at two distinct prolactin-promoter controlled reporter genes in the same individual clonal or primary cells. However, the timing of the cycles was independent and out-of-phase. For the first time, we have analyzed transcription dynamics from two equivalent loci in real-time in single cells. In unstimulated conditions, cells showed independent transcription dynamics at each locus. A key result from these analyses was the evidence for a minimum refractory period in the inactive-phase of transcription. The response to acute signals and the result of manipulation of histone acetylation was consistent with the hypothesis that this refractory period corresponded to a phase of chromatin remodeling which significantly increased the cyclicity. Stochastically timed bursts of transcription in an apparently random subset of cells in a tissue may thus produce an overall coordinated but heterogeneous phenotype capable of acute responses to stimuli.

Probing the structure of long DNA molecules in solution using synchrotron radiation linear dichroism.

Linear dichroism (LD), a spectroscopic method for aligned samples, has been used with a synchrotron radiation source to reveal insights into the structure and stability of DNA with increasing salt concentrations (thus stabilizing the base pairing) and increasing temperature while remaining below the melting point (thus destabilizing the base pairing). Measurements have been made from 350 nm to 182 nm, and the spectral changes observed quantified using a Bayesian Markov chain Monte Carlo (MCMC) algorithm, which uses statistical methods to fit to experimental data. Based on literature H-D exchange experiments, we surmise that the cause of the spectral variations is the induction of transient single stranding of tracts in the DNA polymer, particularly those with significant content of the weaker AT base pairs. More detailed analysis of the LD data will require better nucleotide transition polarization assignments.

Digital circadian and sleep health in individual hospital shift workers: A cross sectional telemonitoring study.

BACKGROUND: Telemonitoring of circadian and sleep cycles could identify shift workers at increased risk of poor health, including cancer and cardiovascular diseases, thus supporting personalized prevention. METHODS: The Circadiem cross-sectional study aimed at determining early warning signals of risk of health alteration in hospital nightshifters (NS) versus dayshifters (DS, alternating morning and afternoon shifts). Circadian rhythmicity in activity, sleep, and temperature was telemonitored on work and free days for one week. Participants wore a bluetooth low energy thoracic accelerometry and temperature sensor that was wirelessly connected to a GPRS gateway and a health data hub server. Hidden Markov modelling of activity quantified Rhythm Index, rest quality (probability, p1-1, of remaining at rest), and rest duration. Spectral analyses determined periods in body surface temperature and accelerometry. Parameters were compared and predictors of circadian and sleep disruption were identified by multivariate analyses using information criteria-based model selection. Clusters of individual shift work response profiles were recognized. FINDINGS: Of 140 per-protocol participants (133 females), there were 63 NS and 77 DS. Both groups had similar median rest amount, yet NS had significantly worse median rest-activity Rhythm Index (0·38 [IQR, 0·29-0·47] vs. 0·69 [0·60-0·77], p<0·0001) and rest quality p1-1 (0·94 [0·94-0·95] vs 0·96 [0·94-0·97], p<0·0001) over the whole study week. Only 48% of the NS displayed a circadian period in temperature, as compared to 70% of the DS (p=0·026). Poor p1-1 was associated with nightshift work on both work (p<0·0001) and free days (p=0·0098). The number of years of past night work exposure predicted poor rest-activity Rhythm Index jointly with shift type, age and chronotype on workdays (p= 0·0074), and singly on free days (p=0·0005). INTERPRETATION: A dedicated analysis toolbox of streamed data from a wearable device identified circadian and sleep rhythm markers, that constitute surrogate candidate endpoints of poor health risk in shift-workers. FUNDING: French Agency for Food, Environmental and Occupational Health & Safety (EST-2014/1/064), University of Warwick, Medical Research Council (United Kingdom, MR/M013170), Cancer Research UK(C53561/A19933).

Identifying the Recurrence of Sleep Apnea Using A Harmonic Hidden Markov Model.

We propose to model time-varying periodic and oscillatory processes by means of a hidden Markov model where the states are defined through the spectral properties of a periodic regime. The number of states is unknown along with the relevant periodicities, the role and number of which may vary across states. We address this inference problem by a Bayesian nonparametric hidden Markov model assuming a sticky hierarchical Dirichlet process for the switching dynamics between different states while the periodicities characterizing each state are explored by means of a trans-dimensional Markov chain Monte Carlo sampling step. We develop the full Bayesian inference algorithm and illustrate the use of our proposed methodology for different simulation studies as well as an application related to respiratory research which focuses on the detection of apnea instances in human breathing traces.

Telemonitored Human Circadian Temperature Dynamics During Daily Routine.

BACKGROUND: Circadian rhythms in body temperature coordinate peripheral molecular clocks, hence they could potentially predict optimal treatment timing (chronotherapy) in individual patients. Circadian parameters in chest surface body temperature (Chesttemp) were recorded remotely and in real time through the use of wearable sensors. METHODS: The dynamics of circadian oscillations in Chesttemp and core body temperature (Coretemp) and their moderation by sex and age were analysed in 38 men and 50 women, aged 21-78 years. In two studies (ST1 and ST2), Chesttemp was measured every minute and teletransmitted using a BLE-connected sensor for 3.6-28.3 days. Additionally, in ST2, Coretemp was recorded per minute in 33 age- and sex-stratified subjects using electronic ingestible pills with radio-frequency transmissions. Circadian parameters were computed using spectral analysis and cosinor modelling. The temporal relations between Chesttemp and Coretemp cosinor parameters were summarised with principal component (PC) analysis. The effect of sex and age was analysed through multivariate regression. RESULTS: Using spectral analysis, a dominant period of 24- or 12-h was identified in 93.2% of the Chesttemp and in 100% of the Coretemp time series. The circadian parameters varied largely between-subjects both for Chesttemp (ranges: mesors, 33.2-36.6°C; amplitudes, 0.2-2.5°C; acrophases, 14:05-7:40), and Coretemp (mesors, 36.6-37.5°C; amplitudes, 0.2-0.7°C; bathyphases, 23:50-6:50). Higher PC loadings mainly corresponded to (i) large Chesttemp amplitudes, and phase advance of both temperature rhythms for the first PC (PC1, 27.2% of variance var.), (ii) high mesors in both temperature rhythms for PC2 (22.4% var.), and (iii) large Coretemp amplitudes for PC3 (12.9% var.). Chesttemp and Coretemp mesors and PC2 loadings decreased in females, while remaining quite stable in males as a function of age. In contrast, Coretemp amplitude and PC3 loadings increased with age in females, but decreased in males. Finally, older subjects, both female and male, displayed a reduction in ultradian variabilities, and an increase in both Chesttemp circadian amplitude and PC1 loadings. INTERPRETATION: The dynamics relations between Chesttemp and Coretemp rhythms were largely moderated by age and sex, with results suggesting that treatment timing could be most critical for therapeutic index in women and in order people.

3 '-5 ' crosstalk contributes to transcriptional bursting.

BACKGROUND: Transcription in mammalian cells is a complex stochastic process involving shuttling of polymerase between genes and phase-separated liquid condensates. It occurs in bursts, which results in vastly different numbers of an mRNA species in isogenic cell populations. Several factors contributing to transcriptional bursting have been identified, usually classified as intrinsic, in other words local to single genes, or extrinsic, relating to the macroscopic state of the cell. However, some possible contributors have not been explored yet. Here, we focus on processes at the 3 ' and 5 ' ends of a gene that enable reinitiation of transcription upon termination. RESULTS: Using Bayesian methodology, we measure the transcriptional bursting in inducible transgenes, showing that perturbation of polymerase shuttling typically reduces burst size, increases burst frequency, and thus limits transcriptional noise. Analysis based on paired-end tag sequencing (PolII ChIA-PET) suggests that this effect is genome wide. The observed noise patterns are also reproduced by a generative model that captures major characteristics of the polymerase flux between the ends of a gene and a phase-separated compartment. CONCLUSIONS: Interactions between the 3 ' and 5 ' ends of a gene, which facilitate polymerase recycling, are major contributors to transcriptional noise.