Abstract:Disease development and cell differentiation both involve dynamic changes; therefore, the reconstruction of dynamic gene regulatory networks (DGRNs) is an important but difficult problem in systems biology. With recent technical advances in single-cell RNA sequencing (scRNA-seq), large volumes of scRNA-seq data are being obtained for various processes. However, most current methods of inferring DGRNs from bulk samples may not be suitable for scRNA-seq data. In this work, we present scPADGRN, a novel DGRN inference method using "time-series" scRNA-seq data. scPADGRN combines the preconditioned alternating direction method of multipliers with cell clustering for DGRN reconstruction. It exhibits advantages in accuracy, robustness and fast convergence. Moreover, a quantitative index called Differentiation Genes' Interaction Enrichment (DGIE) is presented to quantify the interaction enrichment of genes related to differentiation. From the DGIE scores of relevant subnetworks, we infer that the functions of embryonic stem (ES) cells are most active initially and may gradually fade over time. The communication strength of known contributing genes that facilitate cell differentiation increases from ES cells to terminally differentiated cells. We also identify several genes responsible for the changes in the DGIE scores occurring during cell differentiation based on three real single-cell datasets. Our results demonstrate that single-cell analyses based on network inference coupled with quantitative computations can reveal key transcriptional regulators involved in cell differentiation and disease development.Single-cell RNA sequencing (scRNA-seq) data are gaining popularity for providing access to cell-level measurements. Currently, time-series scRNA-seq data allow researchers to study dynamic changes during biological processes. This work proposes a novel method, scPADGRN, for application to time-series scRNA-seq data to construct dynamic gene regulatory networks, which are informative for investigating dynamic changes during disease development and cell differentiation. The proposed method shows satisfactory performance on both simulated data and three real datasets concerning cell differentiation. To quantify network dynamics, we present a quantitative index, DGIE, to measure the degree of activity of a certain set of genes in a regulatory network. Quantitative computations based on dynamic networks identify key regulators in cell differentiation and reveal the activity states of the identified regulators. Specifically, Bhlhe40, Msx2, Foxa2 and Dnmt3l might be important regulatory genes involved in differentiation from mouse ES cells to primitive endoderm (PrE) cells. For differentiation from mouse embryonic fibroblast cells to myocytes, Scx, Fos and Tcf12 are suggested to be key regulators. Sox5, Meis2, Hoxb3, Tcf7l1 and Plagl1 critically contribute during differentiation from human ES cells to definitive endoderm cells. These results may guide further theoretical and experimental efforts to understand cell differentiation processes and explore cell heterogeneity.

Estimating drivers of cell state transitions using gene regulatory network models

Identification of Regulatory Variables for State Transition of Biological Networks.

Parameterizing Cell-to-cell Regulatory Heterogeneities Via Stochastic Transcriptional Profiles

Inference of Spatiotemporal Effects on Cellular State Transitions from Time-Lapse Microscopy

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

Estimating sample-specific regulatory networks

Functional Dissection of Regulatory Models Using Gene Expression Data of Deletion Mutants.

Detecting phenotype-driven transitions in regulatory network structure

scPADGRN: A preconditioned ADMM approach for reconstructing dynamic gene regulatory network using single-cell RNA sequencing data

State-space modeling of dynamic genetic networks

Network estimation in State Space Model with L1-regularization constraint

Directing cellular transitions on gene graph-enhanced cell state manifold

Flexible modeling of regulatory networks improves transcription factor activity estimation

Toward Modeling Context-Specific EMT Regulatory Networks Using Temporal Single Cell RNA-Seq Data

Uncovering the Cell Fate Decision in Lysis-Lysogeny Transition and Stem Cell Development Via Markov State Modeling

Inferring gene regulatory networks from single-cell data: a mechanistic approach

Synthesising Executable Gene Regulatory Networks from Single-cell Gene Expression Data

LaGrACE: Estimating gene program dysregulation using latent gene regulatory network for biomedical discovery

Generating Ensembles of Gene Regulatory Networks to Assess Robustness of Disease Modules

Dissecting reversible and irreversible single cell state transitions from gene regulatory networks

Inferring gene regulatory networks using transcriptional profiles as dynamical attractors