Exploiting the full power of temporal gene expression profiling through a new statistical test: Application to the analysis of muscular dystrophy data

BackgroundThe identification of biologically interesting genes in a temporal expression profiling dataset is challenging and complicated by high levels of experimental noise. Most statistical methods used in the literature do not fully exploit the temporal ordering in the dataset and are not suited to the case where temporal profiles are measured for a number of different biological conditions. We present a statistical test that makes explicit use of the temporal order in the data by fitting polynomial functions to the temporal profile of each gene and for each biological condition. A Hotelling T2-statistic is derived to detect the genes for which the parameters of these polynomials are significantly different from each other.ResultsWe validate the temporal Hotelling T2-test on muscular gene expression data from four mouse strains which were profiled at different ages: dystrophin-, beta-sarcoglycan and gamma-sarcoglycan deficient mice, and wild-type mice. The first three are animal models for different muscular dystrophies. Extensive biological validation shows that the method is capable of finding genes with temporal profiles significantly different across the four strains, as well as identifying potential biomarkers for each form of the disease. The added value of the temporal test compared to an identical test which does not make use of temporal ordering is demonstrated via a simulation study, and through confirmation of the expression profiles from selected genes by quantitative PCR experiments. The proposed method maximises the detection of the biologically interesting genes, whilst minimising false detections.ConclusionThe temporal Hotelling T2-test is capable of finding relatively small and robust sets of genes that display different temporal profiles between the conditions of interest. The test is simple, it can be used on gene expression data generated from any experimental design and for any number of conditions, and it allows fast interpretation of the temporal behaviour of genes. The R code is available from V.V. The microarray data have been submitted to GEO under series GSE1574 and GSE3523.

[1]  J. D. Porter,et al.  A chronic inflammatory response dominates the skeletal muscle molecular signature in dystrophin-deficient mdx mice. , 2002, Human molecular genetics.

[2]  D. Botstein,et al.  Cluster analysis and display of genome-wide expression patterns. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[3]  K. Moore,et al.  X chromosome-linked muscular dystrophy (mdx) in the mouse. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[4]  G. Serratrice,et al.  Muscle hypertrophy in Duchenne muscular dystrophy A pathological and morphometric study , 2022 .

[5]  C. Angelini,et al.  Decorin and biglycan expression is differentially altered in several muscular dystrophies. , 2005, Brain : a journal of neurology.

[6]  Michel Georges,et al.  The callipyge mutation enhances the expression of coregulated imprinted genes in cis without affecting their imprinting status , 2001, Nature Genetics.

[7]  J. Olson,et al.  A regression-based method to identify differentially expressed genes in microarray time course studies and its application in an inducible Huntington's disease transgenic model. , 2002, Human molecular genetics.

[8]  G. Valle,et al.  Limb-girdle muscular dystrophy type 2G is caused by mutations in the gene encoding the sarcomeric protein telethonin , 2000, Nature Genetics.

[9]  L. Kunkel,et al.  Gene expression comparison of biopsies from Duchenne muscular dystrophy (DMD) and normal skeletal muscle , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[10]  G. Churchill,et al.  Experimental design for gene expression microarrays. , 2001, Biostatistics.

[11]  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.

[12]  B. Petrof.  Molecular Pathophysiology of Myofiber Injury in Deficiencies of the Dystrophin-Glycoprotein Complex , 2002, American journal of physical medicine & rehabilitation.

[13]  Satoru Miyano,et al.  Statistical analysis of a small set of time-ordered gene expression data using linear splines , 2002, Bioinform..

[14]  Dean Y. Li,et al.  Altered vascular remodeling in fibulin-5-deficient mice reveals a role of fibulin-5 in smooth muscle cell proliferation and migration. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Korbinian Strimmer,et al.  Identifying periodically expressed transcripts in microarray time series data , 2008, Bioinform..

[16]  Paola Sebastiani,et al.  Cluster analysis of gene expression dynamics , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[17]  P. Hoen,et al.  Muscle regeneration in dystrophin-deficient mdx mice studied by gene expression profiling , 2005, BMC Genomics.

[18]  S. Noguchi,et al.  Common pathological mechanisms in mouse models for muscular dystrophies , 2006, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[19]  Susan K. Murphy,et al.  Abnormal postnatal maintenance of elevated DLK1 transcript levels in callipyge sheep , 2005, Mammalian Genome.

[20]  G. Serratrice,et al.  Muscle hypertrophy in Duchenne muscular dystrophy , 2004, Journal of Neurology.

[21]  Taesung Park,et al.  Statistical tests for identifying differentially expressed genes in time-course microarray experiments , 2003, Bioinform..

[22]  H. Hama,et al.  Loss of the sarcoglycan complex and sarcospan leads to muscular dystrophy in beta-sarcoglycan-deficient mice. , 1999, Human molecular genetics.

[23]  S. Yamaguchi,et al.  Role of DBP in the Circadian Oscillatory Mechanism , 2000, Molecular and Cellular Biology.

[24]  Kiyotoshi Kaneko,et al.  Pathological analysis of muscle hypertrophy and degeneration in muscular dystrophy in γ-sarcoglycan-deficient mice , 2003, Neuromuscular Disorders.

[25]  T. Rando.  The dystrophin–glycoprotein complex, cellular signaling, and the regulation of cell survival in the muscular dystrophies , 2001, Muscle & nerve.

[26]  Xiaohui Liu,et al.  An experimental evaluation of a loop versus a reference design for two-channel microarrays , 2005, Bioinform..

[27]  S. Iwata,et al.  CD26 up-regulates expression of CD86 on antigen-presenting cells by means of caveolin-1. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Eric P. Hoffman,et al.  Expression Profiling in the Muscular Dystrophies Identification of Novel Aspects of Molecular Pathophysiology , 2000 .

[29]  Ernst Wit,et al.  Statistics for Microarrays : Design, Analysis and Inference , 2004 .

[30]  Hongzhe Li,et al.  Model-based methods for identifying periodically expressed genes based on time course microarray gene expression data , 2004, Bioinform..