Search

Untargeted metabolomic profiling through liquid chromatography-mass spectrometry (LC-MS) measures a vast array of metabolites within biospecimens, advancing drug development, disease diagnosis, and risk prediction. However, the low throughput of LC-MS poses a major challenge for biomarker discovery, annotation, and experimental comparison, necessitating the merging of multiple datasets. Current data pooling methods encounter practical limitations due to their vulnerability to data variations and hyperparameter dependence. Here, we introduce GromovMatcher, a flexible and user-friendly algorithm that automatically combines LC-MS datasets using optimal transport. By capitalizing on feature intensity correlation structures, GromovMatcher delivers superior alignment accuracy and robustness compared to existing approaches. This algorithm scales to thousands of features requiring minimal hyperparameter tuning. Manually curated datasets for validating alignment algorithms are limited in the field of untargeted metabolomics, and hence we develop a dataset split procedure to generate pairs of validation datasets to test the alignments produced by GromovMatcher and other methods. Applying our method to experimental patient studies of liver and pancreatic cancer, we discover shared metabolic features related to patient alcohol intake, demonstrating how GromovMatcher facilitates the search for biomarkers associated with lifestyle risk factors linked to several cancer types.

Abstract Background Amino acid metabolism is dysregulated in colorectal cancer patients; however, it is not clear whether pre-diagnostic levels of amino acids are associated with subsequent risk of colorectal cancer. We investigated circulating levels of amino acids in relation to colorectal cancer risk in the European Prospective Investigation into Cancer and Nutrition (EPIC) and UK Biobank cohorts. Methods Concentrations of 13-21 amino acids were determined in baseline fasting plasma or serum samples in 654 incident colorectal cancer cases and 654 matched controls in EPIC. Amino acids associated with colorectal cancer risk following adjustment for the false discovery rate (FDR) were then tested for associations in the UK Biobank, for which measurements of 9 amino acids were available in 111,323 participants, of which 1221 were incident colorectal cancer cases. Results Histidine levels were inversely associated with colorectal cancer risk in EPIC (odds ratio [OR] 0.80 per standard deviation [SD], 95% confidence interval [CI] 0.69–0.92, FDR P-value=0.03) and in UK Biobank (HR 0.93 per SD, 95% CI 0.87–0.99, P-value=0.03). Glutamine levels were borderline inversely associated with colorectal cancer risk in EPIC (OR 0.85 per SD, 95% CI 0.75–0.97, FDR P-value=0.08) and similarly in UK Biobank (HR 0.95, 95% CI 0.89–1.01, P=0.09) In both cohorts, associations changed only minimally when cases diagnosed within 2 or 5 years of follow-up were excluded. Conclusions Higher circulating levels of histidine were associated with a lower risk of colorectal cancer in two large prospective cohorts. Further research to ascertain the role of histidine metabolism and potentially that of glutamine in colorectal cancer development is warranted.

Abstract Background Epidemiological studies of associations between metabolites and cancer risk have typically focused on specific cancer types separately. Here, we designed a multivariate pan-cancer analysis to identify metabolites potentially associated with multiple cancer types, while also allowing the investigation of cancer type-specific associations. Methods We analysed targeted metabolomics data available for 5828 matched case-control pairs from cancer-specific case-control studies on breast, colorectal, endometrial, gallbladder, kidney, localized and advanced prostate cancer, and hepatocellular carcinoma nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. From pre-diagnostic blood levels of an initial set of 117 metabolites, 33 cluster representatives of strongly correlated metabolites and 17 single metabolites were derived by hierarchical clustering. The mutually adjusted associations of the resulting 50 metabolites with cancer risk were examined in penalized conditional logistic regression models adjusted for body mass index, using the data-shared lasso penalty. Results Out of the 50 studied metabolites, (i) six were inversely associated with the risk of most cancer types: glutamine, butyrylcarnitine, lysophosphatidylcholine a C18:2, and three clusters of phosphatidylcholines (PCs); (ii) three were positively associated with most cancer types: proline, decanoylcarnitine, and one cluster of PCs; and (iii) 10 were specifically associated with particular cancer types, including histidine that was inversely associated with colorectal cancer risk and one cluster of sphingomyelins that was inversely associated with risk of hepatocellular carcinoma and positively with endometrial cancer risk. Conclusions These results could provide novel insights for the identification of pathways for cancer development, in particular those shared across different cancer types.

Pooling metabolomics data across studies is often desirable to increase the statistical power of the analysis. However, this can raise methodological challenges as several preanalytical and analytical factors could introduce differences in measured concentrations and variability between datasets. Specifically, different studies may use variable sample types (e.g., serum versus plasma) collected, treated, and stored according to different protocols, and assayed in different laboratories using different instruments. To address these issues, a new pipeline was developed to normalize and pool metabolomics data through a set of sequential steps: (i) exclusions of the least informative observations and metabolites and removal of outliers; imputation of missing data; (ii) identification of the main sources of variability through principal component partial R-square (PC-PR2) analysis; (iii) application of linear mixed models to remove unwanted variability, including samples’ originating study and batch, and preserve biological variations while accounting for potential differences in the residual variances across studies. This pipeline was applied to targeted metabolomics data acquired using Biocrates AbsoluteIDQ kits in eight case-control studies nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. Comprehensive examination of metabolomics measurements indicated that the pipeline improved the comparability of data across the studies. Our pipeline can be adapted to normalize other molecular data, including biomarkers as well as proteomics data, and could be used for pooling molecular datasets, for example in international consortia, to limit biases introduced by inter-study variability. This versatility of the pipeline makes our work of potential interest to molecular epidemiologists.