A novel composition of endogenous metabolic modulators improves red blood cell properties in sickle cell disease.

The most common forms of sickle cell disease (SCD) are sickle cell anemia (SCA; HbSS) and HbSC disease. In both, especially the more dense, dehydrated and adherent red blood cells (RBCs) with reduced deformability are prone to hemolysis and sickling, and thereby vaso-occlusion. Based on plasma amino acid profiling in SCD, a composition of 10 amino acids and derivatives (RCitNacQCarLKHVS; Axcella Therapeutics, USA), referred to as endogenous metabolic modulators (EMMs), was designed to target RBC metabolism. The effects of ex vivo treatment with the EMM composition on different RBC properties were studied in SCD (n = 9 SCA, n = 5 HbSC disease). Dose-dependent improvements were observed in RBC hydration assessed by hemocytometry (MCV, MCHC, dense RBCs) and osmotic gradient ektacytometry (Ohyper). Median (interquartile range [IQR]) increase in Ohyper compared to vehicle was 4.9% (4.0%-5.5%), 7.5% (6.9%-9.4%), and 12.8% (11.5%-14.0%) with increasing 20×, 40×, and 80X concentrations, respectively (all p < 0.0001). RBC deformability (EImax using oxygen gradient ektacytometry) increased by 8.1% (2.2%-12.1%; p = 0.0012), 9.6% (2.9%-15.1%; p = 0.0013), and 13.3% (5.7%-25.5%; p = 0.0007), respectively. Besides, RBC adhesion to subendothelial laminin decreased by 43% (6%-68%; p = 0.4324), 58% (48%-72%; p = 0.0185), and 71% (49%-82%; p = 0.0016), respectively. Together, these results provide a rationale for further studies with the EMM composition targeting multiple RBC properties in SCD.

Metabolic Alterations in NADSYN1-Deficient Cells.

NAD synthetase 1 (encoded by the gene NADSYN1) is a cytosolic enzyme that catalyzes the final step in the biosynthesis of nicotinamide adenine dinucleotide (NAD+) from tryptophan and nicotinic acid. NADSYN1 deficiency has recently been added to the spectrum of congenital NAD+ deficiency disorders. To gain insight into the metabolic consequences of NADSYN1 deficiency, the encoding gene was disrupted in A549 and HEK293T cells, and the metabolome was profiled in the presence of different NAD+ precursors, including tryptophan, nicotinamide and nicotinic acid. We demonstrate that when precursors of the NAD+ salvage pathway in the form of nicotinamide become limiting, NADSYN1 deficiency results in a decline in intracellular NAD+ levels even in the presence of other potential NAD+ sources such as tryptophan and nicotinic acid. As a consequence, alterations in 122 and 69 metabolites are observed in NADSYN1-deficient A549 and HEK293T cells compared to the wild-type cell line (FC > 2 and p < 0.05). We thus show that NADSYN1 deficiency results in a metabolic phenotype characterized by alterations in glycolysis, the TCA cycle, the pentose phosphate pathway, and the polyol pathway.

A one-year pilot study comparing direct-infusion high resolution mass spectrometry based untargeted metabolomics to targeted diagnostic screening for inherited metabolic diseases.

Background: Early diagnosis of inherited metabolic diseases (IMDs) is important because treatment may lead to reduced mortality and improved prognosis. Due to their diversity, it is a challenge to diagnose IMDs in time, effecting an emerging need for a comprehensive test to acquire an overview of metabolite status. Untargeted metabolomics has proven its clinical potential in diagnosing IMDs, but is not yet widely used in genetic metabolic laboratories. Methods: We assessed the potential role of plasma untargeted metabolomics in a clinical diagnostic setting by using direct infusion high resolution mass spectrometry (DI-HRMS) in parallel with traditional targeted metabolite assays. We compared quantitative data and qualitative performance of targeted versus untargeted metabolomics in patients suspected of an IMD (n = 793 samples) referred to our laboratory for 1 year. To compare results of both approaches, the untargeted data was limited to polar metabolites that were analyzed in targeted plasma assays. These include amino acid, (acyl)carnitine and creatine metabolites and are suitable for diagnosing IMDs across many of the disease groups described in the international classification of inherited metabolic disorders (ICIMD). Results: For the majority of metabolites, the concentrations as measured in targeted assays correlated strongly with the semi quantitative Z-scores determined with DI-HRMS. For 64/793 patients, targeted assays showed an abnormal metabolite profile possibly indicative of an IMD. In 55 of these patients, similar aberrations were found with DI-HRMS. The remaining 9 patients showed only marginally increased or decreased metabolite concentrations that, in retrospect, were most likely to be clinically irrelevant. Illustrating its potential, DI-HRMS detected additional patients with aberrant metabolites that were indicative of an IMD not detected by targeted plasma analysis, such as purine and pyrimidine disorders and a carnitine synthesis disorder. Conclusion: This one-year pilot study showed that DI-HRMS untargeted metabolomics can be used as a first-tier approach replacing targeted assays of amino acid, acylcarnitine and creatine metabolites with ample opportunities to expand. Using DI-HRMS untargeted metabolomics as a first-tier will open up possibilities to look for new biomarkers.

PPA1 Deficiency Causes a Deranged Galactose Metabolism Recognizable in Neonatal Screening.

Two siblings showed increased galactose and galactose-related metabolites in neonatal screening. Diagnostic workup did not reveal abnormalities in any of the known disease-causing enzymes involved in galactose metabolism. Using whole-exome sequencing, we identified a homozygous missense variant in PPA1 encoding the cytosolic pyrophosphatase 1 (PPA1), c.557C>T (p.Thr186Ile). The enzyme activity of PPA1 was determined using a colorimetric assay, and the protein content was visualized via western blotting in skin fibroblasts from one of the affected individuals. The galactolytic activity of the affected fibroblasts was determined by measuring extracellular acidification with a Seahorse XFe96 analyzer. PPA1 activity decreased to 22% of that of controls in the cytosolic fraction of homogenates from patient fibroblasts. PPA1 protein content decreased by 50% according to western blot analysis, indicating a reduced stability of the variant protein. The extracellular acidification rate was reduced in patient fibroblasts when galactose was used as a substrate. Untargeted metabolomics of blood samples revealed an elevation of other metabolites related to pyrophosphate metabolism. Besides hyperbilirubinemia in the neonatal period in one child, both children were clinically unremarkable at the ages of 3 and 14 years, respectively. We hypothesize that the observed metabolic derangement is a possible mild manifestation of PPA1 deficiency. Unresolved abnormalities in galactosemia screening might result in the identification of more individuals with PPA1 deficiency, a newly discovered inborn metabolic disorder (IMD).

One-year safety and efficacy of mitapivat in sickle cell disease: follow-up results of a phase 2, open-label study.

Targeting the primary pathogenic event of sickle cell disease (SCD), the polymerization of sickle hemoglobin (HbS), may prevent downstream clinical events. Mitapivat, an oral pyruvate kinase (PK) activator, has therapeutic potential by increasing adenosine triphosphate (ATP) and decreasing 2,3-diphosphoglycerate (2,3-DPG), a glycolytic red blood cell (RBC) intermediate. In the previously reported 8-week dose-finding period of this phase 2, investigator-initiated, open-label study, mitapivat was well tolerated and showed efficacy in SCD. Here, the 1-year fixed-dose extension period is reported in which 9 of 10 included patients (90%) aged ≥16 years with SCD (HbSS, HbS/ß0, or HbS/ß+) continued with mitapivat. Mostly mild treatment-emergent adverse events (AEs) (most commonly, transaminase increase and headache) were still reported. Apart from the reported nontreatment-related serious AE (SAE) of a urinary tract infection in the dose-finding period, 1 nontreatment-related SAE occurred in the fixed-dose extension period in a patient who died of massive pulmonary embolism due to COVID-19. Importantly, sustained improvement in Hb level (mean increase, 1.1 ± 0.7 g/dL; P = .0014) was seen, which was accompanied by decreases in markers of hemolysis. In addition, the annualized rate of vaso-occlusive events reduced significantly from a historic baseline of 1.33 ± 1.32 to 0.64 ± 0.87 (P = .0489) when combining the dose-finding period and fixed-dose extension period. Cellularly, the ATP:2,3-DPG ratio and Hb-oxygen affinity significantly increased and RBC sickling (point of sickling) nonsignificantly reduced. Overall, this study demonstrated 1-year safety and efficacy of treatment with mitapivat in SCD, supporting further evaluation in ongoing phase 2/3 study (RISE UP, NCT05031780). This trial was registered at https://www.clinicaltrialsregister.eu/ as NL8517 and EudraCT 2019-003438-18.

Abnormal glucose homeostasis and fasting intolerance in patients with congenital porto-systemic shunts.

In physiological glucose homeostasis, the liver plays a crucial role in the extraction of glucose from the portal circulation and storage as glycogen to enable release through glycogenolysis upon fasting. In addition, insulin secreted by the pancreas is partly eliminated from the systemic circulation by hepatic first-pass. Therefore, patients with a congenital porto-systemic shunt present a unique combination of (a) postabsorptive hyperinsulinemic hypoglycaemia (HH) because of decreased insulin elimination and (b) fasting (ketotic) hypoglycaemia because of decreased glycogenolysis. Patients with porto-systemic shunts therefore provide important insight into the role of the portal circulation and hepatic function in different phases of glucose homeostasis.

Parasitic, bacterial, viral, immune-mediated, metabolic and nutritional factors associated with nodding syndrome.

Nodding syndrome is a neglected, disabling and potentially fatal epileptic disorder of unknown aetiology affecting thousands of individuals mostly confined to Eastern sub-Saharan Africa. Previous studies have identified multiple associations-including Onchocerca volvulus, antileiomodin-1 antibodies, vitamin B6 deficiency and measles virus infection-yet, none is proven causal. We conducted a case-control study of children with early-stage nodding syndrome (symptom onset <1 year). Cases and controls were identified through a household survey in the Greater Mundri area in South Sudan. A wide range of parasitic, bacterial, viral, immune-mediated, metabolic and nutritional risk factors was investigated using conventional and state-of-the-art untargeted assays. Associations were examined by multiple logistic regression analysis, and a hypothetical causal model was constructed using structural equation modelling. Of 607 children with nodding syndrome, 72 with early-stage disease were included as cases and matched to 65 household- and 44 community controls. Mansonella perstans infection (odds ratio 7.04, 95% confidence interval 2.28-21.7), Necator americanus infection (odds ratio 2.33, 95% confidence interval 1.02-5.3), higher antimalarial seroreactivity (odds ratio 1.75, 95% confidence interval 1.20-2.57), higher vitamin E concentration (odds ratio 1.53 per standard deviation increase, 95% confidence interval 1.07-2.19) and lower vitamin B12 concentration (odds ratio 0.56 per standard deviation increase, 95% confidence interval 0.36-0.87) were associated with higher odds of nodding syndrome. In a structural equation model, we hypothesized that Mansonella perstans infection, higher vitamin E concentration and fewer viral exposures increased the risk of nodding syndrome while lower vitamin B12 concentration, Necator americanus and malaria infections resulted from having nodding syndrome. We found no evidence that Onchocerca volvulus, antileiomodin-1 antibodies, vitamin B6 and other factors were associated with nodding syndrome. Our results argue against several previous causal hypotheses including Onchocerca volvulus. Instead, nodding syndrome may be caused by a complex interplay between multiple pathogens and nutrient levels. Further studies need to confirm these associations and determine the direction of effect.

The malate-aspartate shuttle is important for de novo serine biosynthesis.

The malate-aspartate shuttle (MAS) is a redox shuttle that transports reducing equivalents across the inner mitochondrial membrane while recycling cytosolic NADH to NAD+. We genetically disrupted each MAS component to generate a panel of MAS-deficient HEK293 cell lines in which we performed [U-13C]-glucose tracing. MAS-deficient cells have reduced serine biosynthesis, which strongly correlates with the lactate M+3/pyruvate M+3 ratio (reflective of the cytosolic NAD+/NADH ratio), consistent with the NAD+ dependency of phosphoglycerate dehydrogenase in the serine synthesis pathway. Among the MAS-deficient cells, those lacking malate dehydrogenase 1 (MDH1) show the most severe metabolic disruptions, whereas oxoglutarate-malate carrier (OGC)- and MDH2-deficient cells are less affected. Increasing the NAD+-regenerating capacity using pyruvate supplementation resolves most of the metabolic disturbances. Overall, we show that the MAS is important for de novo serine biosynthesis, implying that serine supplementation could be used as a therapeutic strategy for MAS defects and possibly other redox disorders.

A second case of glutaminase hyperactivity: Expanding the phenotype with epilepsy.

Glutaminase (GLS) hyperactivity was first described in 2019 in a patient with profound developmental delay and infantile cataract. Here, we describe a 4-year-old boy with GLS hyperactivity due to a de novo heterozygous missense variant in GLS, detected by trio whole exome sequencing. This boy also exhibits developmental delay without dysmorphic features, but does not have cataract. Additionally, he suffers from epilepsy with tonic clonic seizures. In line with the findings in the previously described patient with GLS hyperactivity, in vivo 3 T magnetic resonance spectroscopy (MRS) of the brain revealed an increased glutamate/glutamine ratio. This increased ratio was also found in urine with UPLC-MS/MS, however, inconsistently. This case indicates that the phenotypic spectrum evoked by GLS hyperactivity may include epilepsy. Clarifying this phenotypic spectrum is of importance for the prognosis and identification of these patients. The combination of phenotyping, genetic testing, and metabolic diagnostics with brain MRS and in urine is essential to identify new patients with GLS hyperactivity and to further extend the phenotypic spectrum of this disease.

Untargeted Metabolomics Identifies Potential Hypertrophic Cardiomyopathy Biomarkers in Carriers of MYBPC3 Founder Variants.

Hypertrophic cardiomyopathy (HCM) is the most prevalent monogenic heart disease, commonly caused by pathogenic MYBPC3 variants, and a significant cause of sudden cardiac death. Severity is highly variable, with incomplete penetrance among genotype-positive family members. Previous studies demonstrated metabolic changes in HCM. We aimed to identify metabolite profiles associated with disease severity in carriers of MYBPC3 founder variants using direct-infusion high-resolution mass spectrometry in plasma of 30 carriers with a severe phenotype (maximum wall thickness ≥20 mm, septal reduction therapy, congestive heart failure, left ventricular ejection fraction <50%, or malignant ventricular arrhythmia) and 30 age- and sex-matched carriers with no or a mild phenotype. Of the top 25 mass spectrometry peaks selected by sparse partial least squares discriminant analysis, XGBoost gradient boosted trees, and Lasso logistic regression (42 total), 36 associated with severe HCM at a p < 0.05, 20 at p < 0.01, and 3 at p < 0.001. These peaks could be clustered to several metabolic pathways, including acylcarnitine, histidine, lysine, purine and steroid hormone metabolism, and proteolysis. In conclusion, this exploratory case-control study identified metabolites associated with severe phenotypes in MYBPC3 founder variant carriers. Future studies should assess whether these biomarkers contribute to HCM pathogenesis and evaluate their contribution to risk stratification.

Heterozygosity for bisphosphoglycerate mutase deficiency expressing clinically as congenital erythrocytosis: A case series and literature review.

Erythrocytosis is associated with increased red blood cell mass and can be either congenital or acquired. Congenital secondary causes are rare and include germline variants increasing haemoglobin (Hb)-oxygen affinity (e.g., Hb or bisphosphoglycerate mutase (BPGM) variants) or affecting oxygen-sensing pathway proteins. Here, we describe five adults from three kindreds with erythrocytosis associated with heterozygosity for BPGM variants, including one novel. Functional analyses showed partial BPGM deficiency, reduced 2,3-bisphosphoglycerate levels and/or increased Hb-oxygen affinity. We also review currently known BPGM variants. This study contributes to raising awareness of BPGM variants, and in particular that heterozygosity for BPGM deficiency may already manifest clinically.

Untargeted metabolic analysis in dried blood spots reveals metabolic signature in 22q11.2 deletion syndrome.

The 22q11.2 deletion syndrome (22q11.2DS) is characterized by a well-defined microdeletion and is associated with increased risk of neurodevelopmental phenotypes including autism spectrum disorders (ASD) and intellectual impairment. The typically deleted region in 22q11.2DS contains multiple genes with the potential of altering metabolism. Deficits in metabolic processes during early brain development may help explain the increased prevalence of neurodevelopmental phenotypes seen in 22q11.2DS. However, relatively little is known about the metabolic impact of the 22q11.2 deletion, while such insight may lead to increased understanding of the etiology. We performed untargeted metabolic analysis in a large sample of dried blood spots derived from 49 22q11.2DS patients and 87 controls, to identify a metabolic signature for 22q11.2DS. We also examined trait-specific metabolomic patterns within 22q11.2DS patients, focusing on intelligence (intelligence quotient, IQ) and ASD. We used the Boruta algorithm to select metabolites distinguishing patients from controls, patients with ASD from patients without, and patients with an IQ score in the lowest range from patients with an IQ score in the highest range. The relevance of the selected metabolites was visualized with principal component score plots, after which random forest analysis and logistic regression were used to measure predictive performance of the selected metabolites. Analysis yielded a distinct metabolic signature for 22q11.2DS as compared to controls, and trait-specific (IQ and ASD) metabolomic patterns within 22q11.2DS patients. The metabolic characteristics of 22q11.2DS provide insights in biological mechanisms underlying the neurodevelopmental phenotype and may ultimately aid in identifying novel therapeutic targets for patients with developmental disorders.

The potential and limitations of intrahepatic cholangiocyte organoids to study inborn errors of metabolism.

Inborn errors of metabolism (IEMs) comprise a diverse group of individually rare monogenic disorders that affect metabolic pathways. Mutations lead to enzymatic deficiency or dysfunction, which results in intermediate metabolite accumulation or deficit leading to disease phenotypes. Currently, treatment options for many IEMs are insufficient. Rarity of individual IEMs hampers therapy development and phenotypic and genetic heterogeneity suggest beneficial effects of personalized approaches. Recently, cultures of patient-own liver-derived intrahepatic cholangiocyte organoids (ICOs) have been established. Since most metabolic genes are expressed in the liver, patient-derived ICOs represent exciting possibilities for in vitro modeling and personalized drug testing for IEMs. However, the exact application range of ICOs remains unclear. To address this, we examined which metabolic pathways can be studied with ICOs and what the potential and limitations of patient-derived ICOs are to model metabolic functions. We present functional assays in patient ICOs with defects in branched-chain amino acid metabolism (methylmalonic acidemia), copper metabolism (Wilson disease), and transporter defects (cystic fibrosis). We discuss the broad range of functional assays that can be applied to ICOs, but also address the limitations of these patient-specific cell models. In doing so, we aim to guide the selection of the appropriate cell model for studies of a specific disease or metabolic process.

Distinct Metabolomic Signatures in Preclinical and Obstructive Hypertrophic Cardiomyopathy.

Hypertrophic Cardiomyopathy (HCM) is a common inherited heart disease with poor risk prediction due to incomplete penetrance and a lack of clear genotype-phenotype correlations. Advanced imaging techniques have shown altered myocardial energetics already in preclinical gene variant carriers. To determine whether disturbed myocardial energetics with the potential to serve as biomarkers are also reflected in the serum metabolome, we analyzed the serum metabolome of asymptomatic carriers in comparison to healthy controls and obstructive HCM patients (HOCM). We performed non-quantitative direct-infusion high-resolution mass spectrometry-based untargeted metabolomics on serum from fasted asymptomatic gene variant carriers, symptomatic HOCM patients and healthy controls (n = 31, 14 and 9, respectively). Biomarker panels that discriminated the groups were identified by performing multivariate modeling with gradient-boosting classifiers. For all three group-wise comparisons we identified a panel of 30 serum metabolites that best discriminated the groups. These metabolite panels performed equally well as advanced cardiac imaging modalities in distinguishing the groups. Seven metabolites were found to be predictive in two different comparisons and may play an important role in defining the disease stage. This study reveals unique metabolic signatures in serum of preclinical carriers and HOCM patients that may potentially be used for HCM risk stratification and precision therapeutics.

Inborn disorders of the malate aspartate shuttle.

Over the last few years, various inborn disorders have been reported in the malate aspartate shuttle (MAS). The MAS consists of four metabolic enzymes and two transporters, one of them having two isoforms that are expressed in different tissues. Together they form a biochemical pathway that shuttles electrons from the cytosol into mitochondria, as the inner mitochondrial membrane is impermeable to the electron carrier NADH. By shuttling NADH across the mitochondrial membrane in the form of a reduced metabolite (malate), the MAS plays an important role in mitochondrial respiration. In addition, the MAS maintains the cytosolic NAD+ /NADH redox balance, by using redox reactions for the transfer of electrons. This explains why the MAS is also important in sustaining cytosolic redox-dependent metabolic pathways, such as glycolysis and serine biosynthesis. The current review provides insights into the clinical and biochemical characteristics of MAS deficiencies. To date, five out of seven potential MAS deficiencies have been reported. Most of them present with a clinical phenotype of infantile epileptic encephalopathy. Although not specific, biochemical characteristics include high lactate, high glycerol 3-phosphate, a disturbed redox balance, TCA abnormalities, high ammonia, and low serine, which may be helpful in reaching a diagnosis in patients with an infantile epileptic encephalopathy. Current implications for treatment include a ketogenic diet, as well as serine and vitamin B6 supplementation.

Metabolomics in Severe Aortic Stenosis Reveals Intermediates of Nitric Oxide Synthesis as Most Distinctive Markers.

BACKGROUND: Calcific aortic valve disease (CAVD) is a rapidly growing global health problem with an estimated 12.6 million cases globally in 2017 and a 112% increase of deaths since 1990 due to aging and population growth. CAVD may develop into aortic stenosis (AS) by progressive narrowing of the aortic valve. AS is underdiagnosed, and if treatment by aortic valve replacement (AVR) is delayed, this leads to poor recovery of cardiac function, absence of symptomatic improvement and marked increase of mortality. Considering the current limitations to define the stage of AS-induced cardiac remodeling, there is need for a novel method to aid in the diagnosis of AS and timing of intervention, which may be found in metabolomics profiling of patients. METHODS: Serum samples of nine healthy controls and 10 AS patients before and after AVR were analyzed by untargeted mass spectrometry. Multivariate modeling was performed to determine a metabolic profile of 30 serum metabolites which distinguishes AS patients from controls. Human cardiac microvascular endothelial cells (CMECs) were incubated with serum of the AS patients and then stained for ICAM-1 with Western Blot to analyze the effect of AS patient serum on endothelial cell activation. RESULTS: The top 30 metabolic profile strongly distinguishes AS patients from healthy controls and includes 17 metabolites related to nitric oxide metabolism and 12 metabolites related to inflammation, in line with the known pathomechanism for calcific aortic valve disease. Nine metabolites correlate strongly with left ventricular mass, of which three show reversal back to control values after AVR. Western blot analysis of CMECs incubated with AS patient sera shows a significant reduction (14%) in ICAM-1 in AS samples taken after AVR compared to AS patient sera before AVR. CONCLUSION: Our study defined a top 30 metabolic profile with biological and clinical relevance, which may be used as blood biomarker to identify AS patients in need of cardiac surgery. Future studies are warranted in patients with mild-to-moderate AS to determine if these metabolites reflect disease severity and can be used to identify AS patients in need of cardiac surgery.

MetaboShiny: interactive analysis and metabolite annotation of mass spectrometry-based metabolomics data.

Direct infusion untargeted mass spectrometry-based metabolomics allows for rapid insight into a sample's metabolic activity. However, analysis is often complicated by the large array of detected m/z values and the difficulty to prioritize important m/z and simultaneously annotate their putative identities. To address this challenge, we developed MetaboShiny, a novel R/RShiny-based metabolomics package featuring data analysis, database- and formula-prediction-based annotation and visualization. To demonstrate this, we reproduce and further explore a MetaboLights metabolomics bioinformatics study on lung cancer patient urine samples. MetaboShiny enables rapid and rigorous analysis and interpretation of direct infusion untargeted mass spectrometry-based metabolomics data.

Cross-Omics: Integrating Genomics with Metabolomics in Clinical Diagnostics.

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.