A novel Bayesian approach to quantify clinical variables and to determine their spectroscopic counterparts in ¹H NMR metabonomic data

• Main Authors: Kaski Kimmo, Hannuksela Minna L, Savolainen Markku J, Mäkelä Sanna M, Ingman Petri, Soininen Pasi, Mäkinen Ville-Petteri, Vehtari Aki, Ala-Korpela Mika
• Format: Article
• Language: English
• Published: BMC 2007-05-01
• Series: BMC Bioinformatics

<p>Abstract</p> <p>Background</p> <p>A key challenge in metabonomics is to uncover quantitative associations between multidimensional spectroscopic data and biochemical measures used for disease risk assessment and diagnostics. Here we focus on clinically relevant estimation of lipoprotein lipids by ¹H NMR spectroscopy of serum.</p> <p>Results</p> <p>A Bayesian methodology, with a biochemical motivation, is presented for a real ¹H NMR metabonomics data set of 75 serum samples. Lipoprotein lipid concentrations were independently obtained for these samples via ultracentrifugation and specific biochemical assays. The Bayesian models were constructed by Markov chain Monte Carlo (MCMC) and they showed remarkably good quantitative performance, the predictive R-values being 0.985 for the very low density lipoprotein triglycerides (VLDL-TG), 0.787 for the intermediate, 0.943 for the low, and 0.933 for the high density lipoprotein cholesterol (IDL-C, LDL-C and HDL-C, respectively). The modelling produced a kernel-based reformulation of the data, the parameters of which coincided with the well-known biochemical characteristics of the ¹H NMR spectra; particularly for VLDL-TG and HDL-C the Bayesian methodology was able to clearly identify the most characteristic resonances within the heavily overlapping information in the spectra. For IDL-C and LDL-C the resulting model kernels were more complex than those for VLDL-TG and HDL-C, probably reflecting the severe overlap of the IDL and LDL resonances in the ¹H NMR spectra.</p> <p>Conclusion</p> <p>The systematic use of Bayesian MCMC analysis is computationally demanding. Nevertheless, the combination of high-quality quantification and the biochemical rationale of the resulting models is expected to be useful in the field of metabonomics.</p>