Abstract:Skeletal muscle regeneration is a complex process orchestrated by diverse cell populations within a dynamic niche. In response to muscle damage and intercellular signaling, these cells undergo cell fate and migration decisions including quiescence, activation, proliferation, differentiation, infiltration, apoptosis, and exfiltration. The emergence of single-cell RNA sequencing (scRNA-seq) studies of muscle regeneration offers a significant opportunity to refine models of regeneration and enhance our understanding of cellular interactions. To better understand how crosstalk between cell types governs cell fate decisions and cell population dynamics, we developed a novel non-linear ordinary differential equation model guided by scRNA-seq data. Our model consists of 9 variables and 19 parameters, capturing the dynamics of key myogenic lineage and immune cell types. We calibrated time-series scRNA-seq data to units of cells per cubic millimeter of tissue and fit our model's parameters to capture the observed dynamics, validating on an independent time series. The model successfully captures regeneration dynamics, particularly after incorporating a novel type of regulatory interaction between M2 macrophages and satellite cells that has been hypothesized in the literature. Our model lays a foundation for future computational explorations of muscle regeneration, modeling of disease conditions, and in silico testing of therapeutic strategies.

What problem does this paper attempt to address?

The problem this paper attempts to address is: **Are the known intercellular communication channels sufficient to explain the dynamic changes in cell populations observed during muscle regeneration from single-cell transcriptome data**. Specifically, the authors developed a new nonlinear ordinary differential equation model aimed at better understanding the interactions between different cell types and their fate decisions during muscle regeneration through single-cell transcriptome data (scRNA-seq). The main research questions include: 1. **Sufficiency of intercellular communication**: Are the known intercellular communication channels sufficient to explain the observed dynamic changes in cell populations? 2. **Calibration of model parameters**: How can single-cell transcriptome data be normalized to absolute cell number estimates, and how can model parameters be calibrated based on this? 3. **Determination of model complexity**: How to construct a mathematical model that matches the complexity of the data? Through these questions, the authors aim to verify whether the known mechanisms of intercellular communication are sufficient to explain the dynamic changes in cells during muscle regeneration and to explore whether new regulatory mechanisms need to be introduced. For example, the authors included a hypothetical regulatory interaction in the model, where M2 macrophages deactivate activated satellite cells, with some cells differentiating into muscle cells while others return to a quiescent state, thus replenishing the satellite cell pool. This hypothetical interaction has been mentioned in the literature but has not been explicitly verified.

Differential Equation Modeling of Cell Population Dynamics in Skeletal Muscle Regeneration from Single-Cell Transcriptomic Data