Multiple latent clustering model for the inference of RNA life-cycle kinetic rates from sequencing data

Summary We here propose a hierarchical Bayesian model to infer RNA synthesis, processing, and degradation rates from sequencing data, based on an ordinary differential equation system that models the RNA life cycle. We parametrize the latent kinetic rates, that rule the system, with a novel functional form, and estimate their parameters through 6 Dirichlet process mixture models. Owing to the complexity of this approach, we are able to simultaneously perform inference, clustering and model selection. We apply our method to investigate transcriptional and post-transcriptional responses of murine fibroblasts to the activation of proto-oncogene Myc. Our approach uncovers simultaneous regulations of the rates, which had not previously been observed in this biological system.