Characterizing the genome-wide dynamic regulation of gene expression is important and will be of much interest in the future. However, there is currently no established method for identifying differentially expressed genes in a time course study. Here we propose a significance method for analyzing time course microarray studies that can be applied to the typical types of comparisons and sampling schemes. This method is applied to two studies on humans. In one study, genes are identified that show differential expression over time in response to in vivo endotoxin administration. By using our method, 7,409 genes are called significant at a 1% false-discovery rate level, whereas several existing approaches fail to identify any genes. In another study, 417 genes are identified at a 10% false-discovery rate level that show expression changing with age in the kidney cortex. Here it is also shown that as many as 47% of the genes change with age in a manner more complex than simple exponential growth or decay. The methodology proposed here has been implemented in the freely distributed and open-source EDGE software package.T he identification of genes that show changes in expression between varying biological conditions is a frequent goal in microarray experiments (1). Differential expression can be studied from a static or temporal viewpoint. In a static experiment, the arrays are obtained irrespective of time, capturing only a single moment of gene expression. In a temporal experiment the arrays are collected over a time course, allowing one to study the dynamic behavior of gene expression. A large amount of work has been done on the problem of identifying differentially expressed genes in static experiments (2). Because the regulation of gene expression is a dynamic process, it is also important to identify and characterize changes in gene expression over time. Here we present a general statistical method that identifies genes differentially expressed over time.Several clustering methods have been applied to time course microarray data, including hierarchical clustering (3, 4), principal components-based clustering (5), Bayesian model-based clustering (6), and K-means clustering of curves (7,8). None of these clustering methods is directly applicable to identifying genes that show statistically significant changes in expression over time. The K-means clustering method has been modified to compare expression over time between two groups (9), but this method can only be applied to a few hundred genes at a time because of the computational cost of fitting a single model to all genes simultaneously (7,8). This approach also requires that the statistical significance be calculated under the assumption that the clustering model estimated for one of the groups is true, which is nonstandard and potentially problematic.The method that we propose draws on ideas from the extensive statistical literature on time course data analysis (10, 11), particularly spline-based methods (12-16). It is applicable to detecting changes in exp...