A single-cell gene regulatory network inference method for identifying complex regulatory dynamics across cell phenotypes.

Single-cell RNA-Sequencing has made it possible to infer high-resolution gene regulatory networks (GRNs), providing deep biological insights by revealing regulatory interactions at single-cell resolution. However, current single-cell GRN analysis methods produce only a single GRN per input dataset, potentially missing relationships between cells from different phenotypes. To address this issue, we present SimiC, a single-cell GRN inference method that produces a GRN per phenotype while imposing a similarity constraint that forces a smooth transition between GRNs, allowing for a direct comparison between different states, treatments, or conditions. We show that jointly inferring GRNs can uncover variation in regulatory relationships across phenotypes that would have otherwise been missed. Moreover, SimiC can recapitulate complex regulatory dynamics across a range of systems, both model and non-model alike. Taken together, we establish a new approach to quantitating regulatory architectures between the GRNs of distinct cellular phenotypes, with far-reaching implications for systems biology.