Maurice Berk1, Timothy Ebbels, Giovanni Montana. 1. Statistics Section, Department of Mathematics, Imperial College London, Huxley Building, South Kensington, London SW7 2AZ, UK.
Abstract
MOTIVATION: Metabolomics is the study of the complement of small molecule metabolites in cells, biofluids and tissues. Many metabolomic experiments are designed to compare changes observed over time under two experimental conditions or groups (e.g. a control and drug-treated group) with the goal of identifying discriminatory metabolites or biomarkers that characterize each condition. A common study design consists of repeated measurements taken on each experimental unit thus producing time courses of all metabolites. We describe a statistical framework for estimating time-varying metabolic profiles and their within-group variability and for detecting between-group differences. Specifically, we propose (i) a smoothing splines mixed effects (SME) model that treats each longitudinal measurement as a smooth function of time and (ii) an associated functional test statistic. Statistical significance is assessed by a non-parametric bootstrap procedure. RESULTS: The methodology has been extensively evaluated using simulated data and has been applied to real nuclear magnetic resonance spectroscopy data collected in a preclinical toxicology study as part of a larger project lead by the COMET (Consortium for Metabonomic Toxicology). Our findings are compatible with the previously published studies. AVAILABILITY: An R script is freely available for download at http://www2.imperial.ac.uk/~gmontana/sme.htm.
MOTIVATION: Metabolomics is the study of the complement of small molecule metabolites in cells, biofluids and tissues. Many metabolomic experiments are designed to compare changes observed over time under two experimental conditions or groups (e.g. a control and drug-treated group) with the goal of identifying discriminatory metabolites or biomarkers that characterize each condition. A common study design consists of repeated measurements taken on each experimental unit thus producing time courses of all metabolites. We describe a statistical framework for estimating time-varying metabolic profiles and their within-group variability and for detecting between-group differences. Specifically, we propose (i) a smoothing splines mixed effects (SME) model that treats each longitudinal measurement as a smooth function of time and (ii) an associated functional test statistic. Statistical significance is assessed by a non-parametric bootstrap procedure. RESULTS: The methodology has been extensively evaluated using simulated data and has been applied to real nuclear magnetic resonance spectroscopy data collected in a preclinical toxicology study as part of a larger project lead by the COMET (Consortium for Metabonomic Toxicology). Our findings are compatible with the previously published studies. AVAILABILITY: An R script is freely available for download at http://www2.imperial.ac.uk/~gmontana/sme.htm.
Authors: Ziv Bar-Joseph; Georg Gerber; Itamar Simon; David K Gifford; Tommi S Jaakkola Journal: Proc Natl Acad Sci U S A Date: 2003-08-21 Impact factor: 11.205
Authors: A K Smilde; J A Westerhuis; H C J Hoefsloot; S Bijlsma; C M Rubingh; D J Vis; R H Jellema; H Pijl; F Roelfsema; J van der Greef Journal: Metabolomics Date: 2009-12-04 Impact factor: 4.290
Authors: John C Lindon; Hector C Keun; Timothy Md Ebbels; Jake Mt Pearce; Elaine Holmes; Jeremy K Nicholson Journal: Pharmacogenomics Date: 2005-10 Impact factor: 2.533
Authors: Joshua Heinemann; Brigit Noon; Mohammad J Mohigmi; Aurélien Mazurie; David L Dickensheets; Brian Bothner Journal: J Am Soc Mass Spectrom Date: 2014-07-08 Impact factor: 3.109