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.

A Bayesian Change Point Model for Dynamic Alternative Transcription Start Site Usage During Cellular Differentiation.

An Empirical Bayes Change-Point Model for Identifying 3′ and 5′ Alternative Splicing by Next-Generation RNA Sequencing

A Statistical Change-Point Analysis Approach for Modeling the Ratio of Next Generation Sequencing Reads

Identification of differential alternative splicing events with an adjusted beta-distribution model

A single‐cell genomic strategy for alternative transcript start sites identification

A branch point on differentiation trajectory is the bifurcating event revealed by dynamical network biomarker analysis of single-cell data.

Regsnps-Asb: A Computational Framework for Identifying Allele-Specific Transcription Factor Binding from ATAC-seq Data

Accounting for technical noise in differential expression analysis of single-cell RNA sequencing data

A Bayesian model selection approach for identifying differentially expressed transcripts from RNA-Seq data

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

Detecting Change Points in fMRI Data via Bayesian Inference and Genetic Algorithm Model.

An Optimized Method for Bayesian Connectivity Change Point Model.

Identification of changes in gene expression

Mathematical modeling of GATA-switching for regulating the differentiation of hematopoietic stem cell

Bayesian Change-Point Joint Models for Multivariate Longitudinal and Time-to-Event Data

Temporal dynamics and transcriptional control using single-cell gene expression analysis

A Novel Computational Framework for Genome-Scale Alternative Transcription Units Prediction.

PCB: A Pseudotemporal Causality-Based Bayesian Approach to Identify EMT-associated Regulatory Relationships of AS Events and RBPs During Breast Cancer Progression.

Studying the evolution of transcription factor binding events using multi-species ChIP-Seq data.

Revealing the Mechanism of Lymphoid and Myeloid Cell Differentiation and Transdifferentiation Through Landscape Quantification

Detection of Transcription Factor Co-Binding Patterns in Human Cells Via Point Process Models