Abstract:Multiple biological processes are driven by oscillatory gene expression at different time scales. Pulsatile dynamics are thought to be widespread, and single-cell live imaging of gene expression has lead to a surge of dynamic, possibly oscillatory, data for different gene networks. However, the regulation of gene expression at the level of an individual cell involves reactions between finite numbers of molecules, and this can result in inherent randomness in expression dynamics, which blurs the boundaries between aperiodic fluctuations and noisy oscillators. This underlies a new challenge to the experimentalist because neither intuition nor pre-existing methods work well for identifying oscillatory activity in noisy biological time series. Thus, there is an acute need for an objective statistical method for classifying whether an experimentally derived noisy time series is periodic. Here, we present a new data analysis method that combines mechanistic stochastic modelling with the powerful methods of non-parametric regression with Gaussian processes. Our method can distinguish oscillatory gene expression from random fluctuations of non-oscillatory expression in single-cell time series, despite peak-to-peak variability in period and amplitude of single-cell oscillations. We show that our method outperforms the Lomb-Scargle periodogram in successfully classifying cells as oscillatory or non-oscillatory in data simulated from a simple genetic oscillator model and in experimental data. Analysis of bioluminescent live-cell imaging shows a significantly greater number of oscillatory cells when luciferase is driven by a Hes1 promoter (10/19), which has previously been reported to oscillate, than the constitutive MoMuLV 5' LTR (MMLV) promoter (0/25). The method can be applied to data from any gene network to both quantify the proportion of oscillating cells within a population and to measure the period and quality of oscillations. It is publicly available as a MATLAB package.

Accurate reconstruction of dynamic gene expression and growth rate profiles from noisy measurements

Accurate characterization of dynamic microbial gene expression and growth rate profiles

Pathway dynamics can delineate the sources of transcriptional noise in gene expression

Dynamic Characterization of Growth and Gene Expression Using High-Throughput Automated Flow Cytometry

Parameterizing Cell-to-cell Regulatory Heterogeneities Via Stochastic Transcriptional Profiles

Deconvolution of gene expression noise into physical dynamics of cognate promoters

Simulated maximum likelihood method for estimating kinetic rates in gene expression

Multi-Experiment Nonlinear Mixed Effect Modeling of Single-Cell Translation Kinetics after Transfection

Probing noise in gene expression and protein production

Inferring gene regulatory networks from single-cell data: a mechanistic approach

Predicting Growth Rate from Gene Expression

Reconstruction of Transcriptional Dynamics from Gene Reporter Data Using Differential Equations.

Cgcorrect: a Method to Correct for Confounding Cell-Cell Variation Due to Cell Growth in Single-Cell Transcriptomics

Noise in gene expression is coupled to growth rate

Single-cell measurement of plasmid copy number and promoter activity

Disentangling intrinsic and extrinsic gene expression noise in growing cells

Stochastic Simulation And Statistical Inference Platform For Visualization And Estimation Of Transcriptional Kinetics

Quantifying how post-transcriptional noise and gene copy number variation bias transcriptional parameter inference from mRNA distributions

Estimating RNA dynamics using one time point for one sample in a single-pulse metabolic labeling experiment

Identifying stochastic oscillations in single-cell live imaging time series using Gaussian processes

Advanced methods for gene network identification and noise decomposition from single-cell data