The chronODE framework for modelling multi-omic time series with ordinary differential equations and machine learning.

Journal: Nature communications
Published Date:

Abstract

Many genome-wide studies capture isolated moments in cell differentiation or organismal development. Conversely, longitudinal studies provide a more direct way to study these kinetic processes. Here, we present an approach for modeling gene-expression and chromatin kinetics from such studies: chronODE, an interpretable framework based on ordinary differential equations. chronODE incorporates two parameters that capture biophysical constraints governing the initial cooperativity and later saturation in gene expression. These parameters group genes into three major kinetic patterns: accelerators, switchers, and decelerators. Applying chronODE to bulk and single-cell time-series data from mouse brain development reveals that most genes (~87%) follow simple logistic kinetics. Among them, genes with rapid acceleration and high saturation values are rare, highlighting biochemical limitations that prevent cells from attaining both simultaneously. Early- and late-emerging cell types display distinct kinetic patterns, with essential genes ramping up faster. Extending chronODE to chromatin, we find that genes regulated by both enhancer and silencer cis-regulatory elements are enriched in brain-specific functions. Finally, we develop a bidirectional recurrent neural network to predict changes in gene expression from corresponding chromatin changes, successfully capturing the cumulative effect of multiple regulatory elements. Overall, our framework allows investigation of the kinetics of gene regulation in diverse biological systems.

Authors

  • Beatrice Borsari
    Program in Computational Biology and Biomedical Informatics, Yale University, New Haven, CT 06511, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT 06511, USA.
  • Mor Frank
    Program in Computational Biology and Biomedical Informatics, Yale University, New Haven, CT, USA.
  • Eve S Wattenberg
    Program in Computational Biology and Biomedical Informatics, Yale University, New Haven, CT, USA.
  • Ke Xu
    Mechatronics Engineering of University of Electronic Science and Technology of China, Chengdu, 611731, China.
  • Susanna X Liu
    Program in Computational Biology and Biomedical Informatics, Yale University, New Haven, CT 06511, USA; Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT 06511, USA.
  • Xuezhu Yu
    Program in Computational Biology and Biomedical Informatics, Yale University, New Haven, CT, USA.
  • Mark Gerstein
    Program of Computational Biology and Bioinformatics and Department of Molecular Biophysics and Biochemistry and Department of Computer Science, Yale University, New Haven, CT 06511, USA.