System identification of dynamic patterns of genome-wide gene expression

High-throughput methods systematically measure the internal state of the entire cell, but powerful computational tools are needed to infer dynamics from their raw data. Therefore, we have developed a new computational method, Eigen-genomic System Dynamic-pattern Analysis (ESDA), which uses systems theory to infer dynamic parameters from a time series of gene expression measurements. As many genes are measured at a modest number of time points, estimation of the system matrix is underdetermined and traditional approaches for estimating dynamic parameters are ineffective; thus, ESDA uses the principle of dimensionality reduction to overcome the data imbalance. We identify degradation dynamic patterns of a genomic system using ESDA. We also combine ESDA and Principal-oscillation-pattern (POP) analysis, which has been widely used in geosciences, to identify oscillation patterns. We demonstrate the first application of POP analysis to genome-wide time-series gene-expression data. Both simulation data and real-world data are used in this study to demonstrate the applicability of ESDA to genomic data. The biological interpretations of dynamic patterns are provided. We also show that ESDA not only compares favorably with previous experimental methods and existing computational methods, but that it also provides complementary information relative to other approaches.