Extracting transcriptional events from temporal gene expression patterns during Dictyostelium development

Motivation: The DNA microarray technology can generate a large amount of data describing the time-course of gene expression. These data, when properly interpreted, can yield a great deal of information concerning differential gene expression during development. Much current effort in bioinformatics has been devoted to the analysis of gene expression data, usually via some 'clustering analysis' on the raw data in some abstract high dimensional space. Here, we describe a method where we first 'process' the raw time-course data using a simple biologically based kinetic model of gene expression. This allows us to reduce the vast data to a few vital attributes characterizing each expression profile, e.g. the times of the onset and cessation of the expression of the developmentally regulated genes. These vital attributes can then be trivially clustered by visual inspection to reveal biologically significant effects. Results: We have applied this approach to microarray expression data from samples isolated every 2 h throughout the 24 h developmental program of Dictyostelium discoideum. mRNA accumulation patterns for 50 developmental genes were found to fit the kinetic model with a p-value of 0.05 or better. Transcription of these genes appears to be initiated in bursts at well-defined periods during development, in a manner suggestive of a dependent sequence. This approach can be applied to analyses of other temporal gene expression patterns, including those of the cell cycle.