Impulse model-based differential expression analysis of time course sequencing data

The global gene expression trajectories of cellular systems in response to developmental or environmental stimuli often follow the prototypic single-pulse or state-transition patterns which can be modeled with the impulse model. Here we combine the continuous impulse expression model with a sequencing data noise model in ImpulseDE2, a differential expression algorithm for time course sequencing experiments such as RNA-seq, ATAC-seq and ChIP-seq. We show that ImpulseDE2 outperforms currently used differential expression algorithms on data sets with sufficiently many sampled time points. ImpulseDE2 is capable of differentiating between transiently and monotonously changing expression trajectories. This classification separates genes which are responsible for the initial and final cell state phenotypes from genes which drive or are driven by the cell state transition and identifies down-regulation of oxidative-phosphorylation as a molecular signature which can drive human embryonic stem cell differentiation.