Analysis of factorial time-course microarrays with application to a clinical study of burn injury

Time-course microarray experiments are capable of capturing dynamic gene expression profiles. It is important to study how these dynamic profiles depend on the multiple factors that characterize the experimental condition under which the time course is observed. Analytic methods are needed to simultaneously handle the time course and factorial structure in the data. We developed a method to evaluate factor effects by pooling information across the time course while accounting for multiple testing and nonnormality of the microarray data. The method effectively extracts gene-specific response features and models their dependency on the experimental factors. Both longitudinal and cross-sectional time-course data can be handled by our approach. The method was used to analyze the impact of age on the temporal gene response to burn injury in a large-scale clinical study. Our analysis reveals that 21% of the genes responsive to burn are age-specific, among which expressions of mitochondria and immunoglobulin genes are differentially perturbed in pediatric and adult patients by burn injury. These new findings in the body’s response to burn injury between children and adults support further investigations of therapeutic options targeting specific age groups. The methodology proposed here has been implemented in R package “TANOVA” and submitted to the Comprehensive R Archive Network at http://www.r-project.org/. It is also available for download at http://gluegrant1.stanford.edu/TANOVA/.

[1]  Tommi S. Jaakkola,et al.  Continuous Representations of Time-Series Gene Expression Data , 2003, J. Comput. Biol..

[2]  Wing Hung Wong,et al.  A BOOTSTRAP-BASED NON-PARAMETRIC ANOVA METHOD WITH APPLICATIONS TO FACTORIAL MICROARRAY DATA , 2011 .

[3]  V. Carey,et al.  Mixed-Effects Models in S and S-Plus , 2001 .

[4]  M. Crawley.  Mixed‐Effects Models , 2007 .

[5]  R. Wolfe,et al.  Insulin Sensitivity and Mitochondrial Function Are Improved in Children With Burn Injury During a Randomized Controlled Trial of Fenofibrate , 2007, Annals of surgery.

[6]  Hongzhe Li,et al.  Clustering of time-course gene expression data using a mixed-effects model with B-splines , 2003, Bioinform..

[7]  Hongzhe Li,et al.  A hidden spatial-temporal Markov random field model for network-based analysis of time course gene expression data , 2008, 0803.3942.

[8]  B. Molinari,et al.  Age‐dependent change in reactive oxygen species and nitric oxide generation by rat alveolar macrophages * , 2003, Aging cell.

[9]  R. Gamelli,et al.  Temporal Cytokine Profiles in Severely Burned Patients: A Comparison of Adults and Children , 2008, Molecular medicine.

[10]  G Arturson,et al.  Changes in immunoglobulin levels in severely burned patients. , 1969, Lancet.

[11]  Christina Kendziorski,et al.  Hidden Markov Models for Microarray Time Course Data in Multiple Biological Conditions , 2006 .

[12]  T. Speed,et al.  A multivariate empirical Bayes statistic for replicated microarray time course data , 2006, math/0702685.

[13]  L. Astrakas,et al.  Burn injury causes mitochondrial dysfunction in skeletal muscle. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[14]  Wenxuan Zhong,et al.  A data-driven clustering method for time course gene expression data , 2006, Nucleic acids research.

[15]  C. Ryan,et al.  Objective estimates of the probability of death from burn injuries. , 1998, The New England journal of medicine.

[16]  Ralph Weissleder,et al.  Identification of the target self-antigens in reperfusion injury , 2006, The Journal of experimental medicine.

[17]  R. Tompkins,et al.  Support of the metabolic response to burn injury , 2004, The Lancet.

[18]  R. Tibshirani,et al.  Diagnosis of multiple cancer types by shrunken centroids of gene expression , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[19]  H. Hechtman,et al.  Identification of a specific self-reactive IgM antibody that initiates intestinal ischemia/reperfusion injury. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Chuan Zhou,et al.  Modelling Gene Expression Data over Time: Curve Clustering with Informative Prior Distributions , 2003 .

[21]  John D. Storey,et al.  Significance analysis of time course microarray experiments. , 2005, Proceedings of the National Academy of Sciences of the United States of America.