Independence screening for high dimensional nonlinear additive ODE models with applications to dynamic gene regulatory networks

A statistical framework for recovering pseudo-dynamic networks from static data

MOTIVATION

The collection of temporal or perturbed data is often a prerequisite for reconstructing dynamic networks in most cases. However, these types of data are seldom available for genomic studies in medicine, thus significantly limiting the use of dynamic networks to characterize the biological principles underlying human health and diseases.

RESULTS

We proposed a statistical framework to recover disease risk-associated pseudo-dynamic networks (DRDNet) from steady-state data. We incorporated a varying coefficient model with multiple ordinary differential equations to learn a series of networks. We analyzed the publicly available Genotype-Tissue Expression data to construct networks associated with hypertension risk, and biological findings showed that key genes constituting these networks had pivotal and biologically relevant roles associated with the vascular system. We also provided the selection consistency of the proposed learning procedure and evaluated its utility through extensive simulations.

AVAILABILITY AND IMPLEMENTATION

DRDNet is implemented in the R language, and the source codes are available at https://github.com/chencxxy28/DRDnet/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.