Search

1H nuclear magnetic resonance (NMR) spectroscopy of biofluids generates reproducible data allowing detection and quantification of small molecules in large population cohorts.Statistical models to analyse such data are now well-established, and the use of univariate metabolome wide association studies (MWAS) investigating the spectral features separately has emerged as a computationally efficient and interpretable alternative to multivariate models. The MWAS relies on the accurate estimation of a metabolome wide significance level (MWSL) to be applied to control the family-wise error rate. Subsequent interpretation requires efficient visualization and formal feature annotation, which, in-turn, calls for efficient prioritisation of spectral variables of interest.Using human serum 1H NMR spectroscopic profiles from 3,948 participants from the Multi-Ethnic Study of Atherosclerosis (MESA), we have performed a series of MWAS for serum levels of glucose. We first propose an extension of the conventional MWSL that yields stable estimates of the MWSL across the different model parametrisations and distributional features of the outcome. We propose both efficient visualisation methods, and a strategy based on subsampling and internal validation to prioritize the associations. Our work proposes and illustrates practical and scalable solutions to facilitate the implementation of the MWAS approach and improve interpretation in large cohort studies.

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 aneed for robust processing pipelines which can ensure reliable data for statistical analysis. The COMBIBIO project incorporates serum from approximately 8000 individuals, in 3 cohorts, profiled by 6 assays in 2 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 T2 were removed and a cohort/phase adjustment was applied, resulting in two NMR datasets (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 T2 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 multi-faceted study.