Workflow for Integrated Processing of Multicohort Untargeted 1H NMR Metabolomics Data in Large-Scale Metabolic Epidemiology

• Published in: Journal of proteome research Vol. 15; no. 12; pp. 4188 - 4194
• Main Authors: Karaman, Ibrahim, Ferreira, Diana L. S, Boulangé, Claire L, Kaluarachchi, Manuja R, Herrington, David, Dona, Anthony C, Castagné, Raphaële, Moayyeri, Alireza, Lehne, Benjamin, Loh, Marie, de Vries, Paul S, Dehghan, Abbas, Franco, Oscar H, Hofman, Albert, Evangelou, Evangelos, Tzoulaki, Ioanna, Elliott, Paul, Lindon, John C, Ebbels, Timothy M. D
• Format: Journal Article
• Language: English
• Published: American Chemical Society 02.12.2016
• Subjects: large scale; epidemiology; NMR; quality control; metabolomics; preprocessing; normalization; alignment; multicohort
• ISSN: 1535-3893; 1535-3907
• DOI: 10.1021/acs.jproteome.6b00125

Large-scale metabolomics studies involving thousands of samples present multiple challenges in data analysis, particularly when an untargeted platform is used. Studies with multiple cohorts and analysis platforms exacerbate existing problems such as peak alignment and normalization. Therefore, there is a need for robust processing pipelines that can ensure reliable data for statistical analysis. The COMBI-BIO project incorporates serum from ∼8000 individuals, in three cohorts, profiled by six assays in two phases using both 1H NMR and UPLC–MS. Here we present the COMBI-BIO NMR analysis pipeline and demonstrate its fitness for purpose using representative quality control (QC) samples. NMR spectra were first aligned and normalized. After eliminating interfering signals, outliers identified using Hotelling’s T 2 were removed and a cohort/phase adjustment was applied, resulting in two NMR data sets (CPMG and NOESY). Alignment of the NMR data was shown to increase the correlation-based alignment quality measure from 0.319 to 0.391 for CPMG and from 0.536 to 0.586 for NOESY, showing that the improvement was present across both large and small peaks. End-to-end quality assessment of the pipeline was achieved using Hotelling’s T 2 distributions. For CPMG spectra, the interquartile range decreased from 1.425 in raw QC data to 0.679 in processed spectra, while the corresponding change for NOESY spectra was from 0.795 to 0.636, indicating an improvement in precision following processing. PCA indicated that gross phase and cohort differences were no longer present. These results illustrate that the pipeline produces robust and reproducible data, successfully addressing the methodological challenges of this large multifaceted study.