Characterizing human postprandial metabolic response using multiway data analysis

Analysis of time-resolved postprandial metabolomics data can enhance our knowledge about human metabolism by providing a better understanding of similarities and differences in postprandial responses of individuals, with the potential to advance precision nutrition and medicine. Traditional data analysis methods focus on clustering methods relying on summaries of data across individuals or use univariate methods analyzing one metabolite at a time. However, they fail to provide a compact summary revealing the underlying patterns, i.e., groups of subjects, clusters of metabolites, and their temporal profiles. In this study, we analyze NMR (Nuclear Magnetic Resonance) spectroscopy measurements of plasma samples collected at multiple time points during a meal challenge test from 299 individuals from the COPSAC2000 cohort. We arrange the data as a three-way array: subjects by metabolites by time, and use the CAN-DECOMP/PARAFAC (CP) tensor factorization model to capture the underlying patterns. We analyze the fasting state data to reveal static patterns of subject group differences, and the fasting state-corrected postprandial data to reveal dynamic markers of group differences. Our analysis demonstrates that the CP model reveals replicable and biologically meaningful patterns capturing certain metabolite groups and their temporal profiles, and showing differences among males according to their body mass index (BMI). Furthermore, we observe that certain lipoproteins relate to the group difference differently in the fasting vs. dynamic state in males. While similar dynamic patterns are observed in response to the challenge test in males and females, the BMI-related group difference is only observed in males in the dynamic state.