skip to main content

Title: Bayesian inference of distributed time delay in transcriptional and translational regulation
Abstract Motivation

Advances in experimental and imaging techniques have allowed for unprecedented insights into the dynamical processes within individual cells. However, many facets of intracellular dynamics remain hidden, or can be measured only indirectly. This makes it challenging to reconstruct the regulatory networks that govern the biochemical processes underlying various cell functions. Current estimation techniques for inferring reaction rates frequently rely on marginalization over unobserved processes and states. Even in simple systems this approach can be computationally challenging, and can lead to large uncertainties and lack of robustness in parameter estimates. Therefore we will require alternative approaches to efficiently uncover the interactions in complex biochemical networks.


We propose a Bayesian inference framework based on replacing uninteresting or unobserved reactions with time delays. Although the resulting models are non-Markovian, recent results on stochastic systems with random delays allow us to rigorously obtain expressions for the likelihoods of model parameters. In turn, this allows us to extend MCMC methods to efficiently estimate reaction rates, and delay distribution parameters, from single-cell assays. We illustrate the advantages, and potential pitfalls, of the approach using a birth–death model with both synthetic and experimental data, and show that we can robustly infer model parameters using a relatively more » small number of measurements. We demonstrate how to do so even when only the relative molecule count within the cell is measured, as in the case of fluorescence microscopy.

Availability and implementation

Accompanying code in R is available at

Supplementary information

Supplementary data are available at Bioinformatics online.

« less
 ;  ;  ;  ;  ;  ;  ;
Award ID(s):
Publication Date:
Journal Name:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Motivation

    Modeling single-cell gene expression trends along cell pseudotime is a crucial analysis for exploring biological processes. Most existing methods rely on nonparametric regression models for their flexibility; however, nonparametric models often provide trends too complex to interpret. Other existing methods use interpretable but restrictive models. Since model interpretability and flexibility are both indispensable for understanding biological processes, the single-cell field needs a model that improves the interpretability and largely maintains the flexibility of nonparametric regression models.


    Here, we propose the single-cell generalized trend model (scGTM) for capturing a gene’s expression trend, which may be monotone, hill-shaped or valley-shaped, along cell pseudotime. The scGTM has three advantages: (i) it can capture non-monotonic trends that are easy to interpret, (ii) its parameters are biologically interpretable and trend informative, and (iii) it can flexibly accommodate common distributions for modeling gene expression counts. To tackle the complex optimization problems, we use the particle swarm optimization algorithm to find the constrained maximum likelihood estimates for the scGTM parameters. As an application, we analyze several single-cell gene expression datasets using the scGTM and show that scGTM can capture interpretable gene expression trends along cell pseudotime and reveal molecular insights underlying biological processes.

    Availability and implementation

    Themore »Python package scGTM is open-access and available at

    Supplementary information

    Supplementary data are available at Bioinformatics online.

    « less
  2. Abstract Background

    The biophysics of an organism span multiple scales from subcellular to organismal and include processes characterized by spatial properties, such as the diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales, through the generation of representative models. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent the spatial and stochastic features of a biological system, limiting their insights and applications. However, spatial models describing biological systems with spatial information are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them and highlights the need for simpler methods able to model the spatial features of biological systems.


    In this work, we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models, the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates ofmore »extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters. We additionally investigate objects and processes implicitly represented by ODE model terms and parameters and improve the reproducibility of spatial, stochastic models.


    We developed and demonstrate a method for generating spatiotemporal, multicellular models from non-spatial population dynamics models of multicellular systems. We envision employing our method to generate new ODE model terms from spatiotemporal and multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.

    « less
  3. Abstract STUDY QUESTION

    Is the combined use of fluorescence lifetime imaging microscopy (FLIM)-based metabolic imaging and second harmonic generation (SHG) spindle imaging a feasible and safe approach for noninvasive embryo assessment?


    Metabolic imaging can sensitively detect meaningful metabolic changes in embryos, SHG produces high-quality images of spindles and the methods do not significantly impair embryo viability.


    Proper metabolism is essential for embryo viability. Metabolic imaging is a well-tested method for measuring metabolism of cells and tissues, but it is unclear if it is sensitive enough and safe enough for use in embryo assessment.


    This study consisted of time-course experiments and control versus treatment experiments. We monitored the metabolism of 25 mouse oocytes with a noninvasive metabolic imaging system while exposing them to oxamate (cytoplasmic lactate dehydrogenase inhibitor) and rotenone (mitochondrial oxidative phosphorylation inhibitor) in series. Mouse embryos (n = 39) were measured every 2 h from the one-cell stage to blastocyst in order to characterize metabolic changes occurring during pre-implantation development. To assess the safety of FLIM illumination, n = 144 illuminated embryos were implanted into n = 12 mice, and n = 108 nonilluminated embryos were implanted into n = 9 mice.


    Experiments were performed in mouse embryos and oocytes. Samples weremore »monitored with noninvasive, FLIM-based metabolic imaging of nicotinamide adenine dinucleotide (NADH) and flavin adenine dinucleotide (FAD) autofluorescence. Between NADH cytoplasm, NADH mitochondria and FAD mitochondria, a single metabolic measurement produces up to 12 quantitative parameters for characterizing the metabolic state of an embryo. For safety experiments, live birth rates and pup weights (mean ± SEM) were used as endpoints. For all test conditions, the level of significance was set at P < 0.05.


    Measured FLIM parameters were highly sensitive to metabolic changes due to both metabolic perturbations and embryo development. For oocytes, metabolic parameter values were compared before and after exposure to oxamate and rotenone. The metabolic measurements provided a basis for complete separation of the data sets. For embryos, metabolic parameter values were compared between the first division and morula stages, morula and blastocyst and first division and blastocyst. The metabolic measurements again completely separated the data sets. Exposure of embryos to excessive illumination dosages (24 measurements) had no significant effect on live birth rate (5.1 ± 0.94 pups/mouse for illuminated group; 5.7 ± 1.74 pups/mouse for control group) or pup weights (1.88 ± 0.10 g for illuminated group; 1.89 ± 0.11 g for control group).


    The study was performed using a mouse model, so conclusions concerning sensitivity and safety may not generalize to human embryos. A limitation of the live birth data is also that although cages were routinely monitored, we could not preclude that some runt pups may have been eaten.


    Promising proof-of-concept results demonstrate that FLIM with SHG provide detailed biological information that may be valuable for the assessment of embryo and oocyte quality. Live birth experiments support the method’s safety, arguing for further studies of the clinical utility of these techniques.


    Supported by the Blavatnik Biomedical Accelerator Grant at Harvard University and by the Harvard Catalyst/The Harvard Clinical and Translational Science Center (National Institutes of Health Award UL1 TR001102), by NSF grants DMR-0820484 and PFI-TT-1827309 and by NIH grant R01HD092550-01. T.S. was supported by a National Science Foundation Postdoctoral Research Fellowship in Biology grant (1308878). S.F. and S.A. were supported by NSF MRSEC DMR-1420382. Becker and Hickl GmbH sponsored the research with the loaning of equipment for FLIM. T.S. and D.N. are cofounders and shareholders of LuminOva, Inc., and co-hold patents (US20150346100A1 and US20170039415A1) for metabolic imaging methods. D.S. is on the scientific advisory board for Cooper Surgical and has stock options with LuminOva, Inc.

    « less
  4. Abstract Background

    Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling.


    We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, ormore »single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric “MW-score” (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut.


    METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at

    « less
  5. Abstract
    Excessive phosphorus (P) applications to croplands can contribute to eutrophication of surface waters through surface runoff and subsurface (leaching) losses. We analyzed leaching losses of total dissolved P (TDP) from no-till corn, hybrid poplar (Populus nigra X P. maximowiczii), switchgrass (Panicum virgatum), miscanthus (Miscanthus giganteus), native grasses, and restored prairie, all planted in 2008 on former cropland in Michigan, USA. All crops except corn (13 kg P ha−1 year−1) were grown without P fertilization. Biomass was harvested at the end of each growing season except for poplar. Soil water at 1.2 m depth was sampled weekly to biweekly for TDP determination during March–November 2009–2016 using tension lysimeters. Soil test P (0–25 cm depth) was measured every autumn. Soil water TDP concentrations were usually below levels where eutrophication of surface waters is frequently observed (&gt; 0.02 mg L−1) but often higher than in deep groundwater or nearby streams and lakes. Rates of P leaching, estimated from measured concentrations and modeled drainage, did not differ statistically among cropping systems across years; 7-year cropping system means ranged from 0.035 to 0.072 kg P ha−1 year−1 with large interannual variation. Leached P was positively related to STP, which decreased over the 7 years in all systems. These results indicate that both P-fertilized and unfertilized cropping systems mayMore>>