In metabolic diagnostics, there is an emerging need for a comprehensive test to acquire a complete view of metabolite status. Here, we describe a non-quantitative direct-infusion high-resolution mass spectrometry (DI-HRMS) based metabolomics method and evaluate the method for both dried blood spots (DBS) and plasma. 110 DBS of 42 patients harboring 23 different inborn errors of metabolism (IEM) and 86 plasma samples of 38 patients harboring 21 different IEM were analyzed using DI-HRMS. A peak calling pipeline developed in R programming language provided Z-scores for ~1875 mass peaks corresponding to ~3835 metabolite annotations (including isomers) per sample. Based on metabolite Z-scores, patients were assigned a ‘most probable diagnosis’ by an investigator blinded for the known diagnoses of the patients. Based on DBS sample analysis, 37/42 of the patients, corresponding to 22/23 IEM, could be correctly assigned a ‘most probable diagnosis’. Plasma sample analysis, resulted in a correct ‘most probable diagnosis’ in 32/38 of the patients, corresponding to 19/21 IEM. The added clinical value of the method was illustrated by a case wherein DI-HRMS metabolomics aided interpretation of a variant of unknown significance (VUS) identified by whole-exome sequencing. In summary, non-quantitative DI-HRMS metabolomics in DBS and plasma is a very consistent, high-throughput and nonselective method for investigating the metabolome in genetic disease.

Background: Pathophysiology of life-threatening acute metabolic decompensations (AMD) in propionic acidemia (PA) and isolated methylmalonic acidemia (MMA) is insufficiently understood. Here, we study the metabolomes of PA and MMA patients over time, to improve insight in which biochemical processes are at play during AMD.Methods: Longitudinal data from clinical chemistry analyses and metabolic assays over the life-course of 11 PA and 13 MMA patients were studied retrospectively. Direct-infusion high-resolution mass spectrometry was performed on 234 and 154 remnant dried blood spot and plasma samples of PA and MMA patients, respectively. In addition, a systematic literature search was performed on reported biomarkers. All results were integrated in an assessment of biochemical processes at play during AMD.Results: We confirmed many of the metabolite alterations reported in literature, including increases of plasma valine and isoleucine during AMD in PA patients. We revealed that plasma leucine and phenylalanine, and urinary pyruvic acid were increased during AMD in PA patients. 3-hydroxyisovaleric acid correlated positively with plasma ammonia. We found that known diagnostic biomarkers were not significantly further increased, while intermediates of the branched-chain amino acid (BCAA) degradation pathway were significantly increased during AMD.Conclusions: We revealed that during AMD in PA and MMA, BCAA and BCAA intermediates accumulate, while known diagnostic biomarkers remain essentially unaltered. This implies that these acidic BCAA intermediates are responsible for metabolic acidosis. Based on this, we suggest to measure plasma 3-hydroxyisovaleric acid and urinary ketones or 3-hydroxybutyric acid for the biochemical follow-up of a patient’s metabolic stability.

Next-generation sequencing and next-generation metabolic screening are, independently, increasingly applied in clinical diagnostics of inborn errors of metabolism (IEM). Integrated into a single bioinformatic method, these two –omics technologies can potentially further improve the diagnostic yield for IEM. Here, we present cross-omics: a method that uses untargeted metabolomics results of patient’s dried blood spots (DBSs), indicated by Z-scores and mapped onto human metabolic pathways, to prioritize potentially affected genes. We demonstrate the optimization of three parameters: (1) maximum distance to the primary reaction of the affected protein, (2) an extension stringency threshold reflecting in how many reactions a metabolite can participate, to be able to extend the metabolite set associated with a certain gene, and (3) a biochemical stringency threshold reflecting paired Z-score thresholds for untargeted metabolomics results. Patients with known IEMs were included. We performed untargeted metabolomics on 168 DBSs of 97 patients with 46 different disease-causing genes, and we simulated their whole-exome sequencing results in silico. We showed that for accurate prioritization of disease-causing genes in IEM, it is essential to take into account not only the primary reaction of the affected protein but a larger network of potentially affected metabolites, multiple steps away from the primary reaction.

Untargeted metabolomics may become a standard approach to address diagnostic requests, but, at present, data interpretation is very labor-intensive. To facilitate its implementation in metabolic diagnostic screening, we developed a method for automated data interpretation that preselects the most likely inborn errors of metabolism (IEM). The input parameters of the knowledge-based algorithm were (1) weight scores assigned to 268 unique metabolites for 119 different IEM based on literature and expert opinion, and (2) metabolite Z-scores and ranks based on direct-infusion high resolution mass spectrometry. The output was a ranked list of differential diagnoses (DD) per sample. The algorithm was first optimized using a training set of 110 dried blood spots (DBS) comprising 23 different IEM and 86 plasma samples comprising 21 different IEM. Further optimization was performed using a set of 96 DBS consisting of 53 different IEM. The diagnostic value was validated in a set of 115 plasma samples, which included 58 different IEM and resulted in the correct diagnosis being included in the DD of 72% of the samples, comprising 44 different IEM. The median length of the DD was 10 IEM, and the correct diagnosis ranked first in 37% of the samples. Here, we demonstrate the accuracy of the diagnostic algorithm in preselecting the most likely IEM, based on the untargeted metabolomics of a single sample. We show, as a proof of principle, that automated data interpretation has the potential to facilitate the implementation of untargeted metabolomics for metabolic diagnostic screening, and we provide suggestions for further optimization of the algorithm to improve diagnostic accuracy.