Inference of gene regulatory networks using time-series single-cell RNA-seq data with CRISPR perturbations

Single-cell RNA-seq analysis coupled with CRISPR-based perturbation (scCRISPR) has enabled the inference of gene regulatory networks (GRNs) with causal relationships. However, a snapshot of scCRISPR data may not lead to an accurate inference, since a gene knockout can influence multi-layered downstream over time. Here, we developed RENGE, a new computational method that infers GRNs using a time-series scCRISPR dataset. RENGE models the propagation process of the effects elicited by a gene knockout on its GRN. It can distinguish between direct and indirect regulations, which allows for the inference of regulations by genes that are not knocked out. RENGE therefore outperforms current methods in the accuracy of inferring GRNs. When used on a dataset we derived from human-induced pluripotent stem cells, RENGE yielded a GRN consistent with multiple databases and literature. Accurate inference of GRNs by RENGE would enable the identification of key factors for various biological systems.