Expanded clinical phenotype and untargeted metabolomics analysis in RARS2-related mitochondrial disorder: a case report.

BACKGROUND: RARS2-related mitochondrial disorder is an autosomal recessive mitochondrial encephalopathy caused by biallelic pathogenic variants in the gene encoding the mitochondrial arginyl-transfer RNA synthetase 2 (RARS2, MIM *611524, NM_020320.5). RARS2 catalyzes the transfer of L-arginine to its cognate tRNA during the translation of mitochondrially-encoded proteins. The classical presentation of RARS2-related mitochondrial disorder includes pontocerebellar hypoplasia (PCH), progressive microcephaly, profound developmental delay, feeding difficulties, and hypotonia. Most patients also develop severe epilepsy by three months of age, which consists of focal or generalized seizures that frequently become pharmacoresistant and lead to developmental and epileptic encephalopathy (DEE). CASE PRESENTATION: Here, we describe a six-year-old boy with developmental delay, hypotonia, and failure to thrive who developed an early-onset DEE consistent with Lennox-Gastaut Syndrome (LGS), which has not previously been observed in this disorder. He had dysmorphic features including bilateral macrotia, overriding second toes, a depressed nasal bridge, retrognathia, and downslanting palpebral fissures, and he did not demonstrate progressive microcephaly. Whole genome sequencing identified two variants in RARS2, c.36 + 1G > T, a previously unpublished variant that is predicted to affect splicing and is, therefore, likely pathogenic and c.419 T > G (p.Phe140Cys), a known pathogenic variant. He exhibited significant, progressive generalized brain atrophy and ex vacuo dilation of the supratentorial ventricular system on brain MRI and did not demonstrate PCH. Treatment with a ketogenic diet (KD) reduced seizure frequency and enabled him to make developmental progress. Plasma untargeted metabolomics analysis showed increased levels of lysophospholipid and sphingomyelin-related metabolites. CONCLUSIONS: Our work expands the clinical spectrum of RARS2-related mitochondrial disorder, demonstrating that patients can present with dysmorphic features and an absence of progressive microcephaly, which can help guide the diagnosis of this condition. Our case highlights the importance of appropriate seizure phenotyping in this condition and indicates that patients can develop LGS, for which a KD may be a viable therapeutic option. Our work further suggests that analytes of phospholipid metabolism may serve as biomarkers of mitochondrial dysfunction.

Integrating Genome Sequencing and Untargeted Metabolomics in Monozygotic Twins with a Rare Complex Neurological Disorder.

Multi-omics approaches, which integrate genomics, transcriptomics, proteomics, and metabolomics, have emerged as powerful tools in the diagnosis of rare diseases. We used untargeted metabolomics and whole-genome sequencing (WGS) to gain a more comprehensive understanding of a rare disease with a complex presentation affecting female twins from a consanguineous family. The sisters presented with polymicrogyria, a Dandy-Walker malformation, respiratory distress, and multiorgan dysfunctions. Through WGS, we identified two rare homozygous variants in both subjects, a pathogenic variant in ADGRG1(p.Arg565Trp) and a novel variant in CNTNAP1(p.Glu910Val). These genes have been previously associated with autosomal recessive polymicrogyria and hypomyelinating neuropathy with/without contractures, respectively. The twins exhibited symptoms that overlapped with both of these conditions. The results of the untargeted metabolomics analysis revealed significant metabolic perturbations relating to neurodevelopmental abnormalities, kidney dysfunction, and microbiome. The significant metabolites belong to essential pathways such as lipids and amino acid metabolism. The identification of variants in two genes, combined with the support of metabolic perturbation, demonstrates the rarity and complexity of this phenotype and provides valuable insights into its underlying mechanisms.

Consensus guidelines for the diagnosis and management of succinic semialdehyde dehydrogenase deficiency.

Succinic semialdehyde dehydrogenase deficiency (SSADHD) (OMIM #271980) is a rare autosomal recessive metabolic disorder caused by pathogenic variants of ALDH5A1. Deficiency of SSADH results in accumulation of γ-aminobutyric acid (GABA) and other GABA-related metabolites. The clinical phenotype of SSADHD includes a broad spectrum of non-pathognomonic symptoms such as cognitive disabilities, communication and language deficits, movement disorders, epilepsy, sleep disturbances, attention problems, anxiety, and obsessive-compulsive traits. Current treatment options for SSADHD remain supportive, but there are ongoing attempts to develop targeted genetic therapies. This study aimed to create consensus guidelines for the diagnosis and management of SSADHD. Thirty relevant statements were initially addressed by a systematic literature review, resulting in different evidence levels of strength according to the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) criteria. The highest level of evidence (level A), based on randomized controlled trials, was unavailable for any of the statements. Based on cohort studies, Level B evidence was available for 12 (40%) of the statements. Thereupon, through a process following the Delphi Method and directed by the Appraisal of Guidelines for Research and Evaluation (AGREE II) criteria, expert opinion was sought, and members of an SSADHD Consensus Group evaluated all the statements. The group consisted of neurologists, epileptologists, neuropsychologists, neurophysiologists, metabolic disease specialists, clinical and biochemical geneticists, and laboratory scientists affiliated with 19 institutions from 11 countries who have clinical experience with SSADHD patients and have studied the disorder. Representatives from parent groups were also included in the Consensus Group. An analysis of the survey's results yielded 25 (83%) strong and 5 (17%) weak agreement strengths. These first-of-their-kind consensus guidelines intend to consolidate and unify the optimal care that can be provided to individuals with SSADHD.

Heteroplasmic pathogenic m.12315G>A variant in MT-TL2 presenting with MELAS syndrome and depletion of nitric oxide donors.

The MT-TL2 m.12315G>A pathogenic variant has previously been reported in five individuals with mild clinical phenotypes. Herein we report the case of a 5-year-old child with heteroplasmy for this variant who developed neurological regression and stroke-like episodes similar to those observed in mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS). Biochemical evaluation revealed depletion of arginine on plasma amino acid analysis and low z-scores for citrulline on untargeted plasma metabolomics analysis. These findings suggested that decreased availability of nitric oxide may have contributed to the stroke-like episodes. The use of intravenous arginine during stroke-like episodes and daily enteral L-citrulline supplementation normalized her biochemical values of arginine and citrulline. Untargeted plasma metabolomics showed the absence of nicotinamide and 1-methylnicotinamide, and plasma total glutathione levels were low; thus, nicotinamide riboside and N-acetylcysteine therapies were initiated. This report expands the phenotype associated with the rare mitochondrial variant MT-TL2 m.12315G>A to include neurological regression and a MELAS-like phenotype. Individuals with this variant should undergo in-depth biochemical analysis to include untargeted plasma metabolomics, plasma amino acids, and glutathione levels to help guide a targeted approach to treatment.

Monitoring the treatment of urea cycle disorders using phenylbutyrate metabolite analyses: Still many lessons to learn.

Medications that elicit an alternate pathway for nitrogen excretion such as oral sodium phenylbutyrate (NaPBA) and glycerol phenylbutyrate (GPB) and intravenous sodium phenylacetate (NaPAA) are important for the management of urea cycle disorders (UCDs). Plasma concentrations of their primary metabolite, phenylacetate (PAA), as well as the ratio of PAA to phenylacetylglutamine (PAGN) are useful for guiding dosing and detecting toxicity. However, the frequency of toxic elevations of metabolites and associated clinical covariates is relatively unknown. A retrospective analysis was conducted on 1255 plasma phenylbutyrate metabolite measurements from 387 individuals. An additional analysis was also conducted on a subset of 68 individuals in whom detailed clinical information was available. In the course of these analyses, abnormally elevated plasma PAA and PAA:PAGN were identified in 39 individuals (4.15% of samples) and 42 individuals (4.30% of samples), respectively. Abnormally elevated PAA and PAA:PAGN values were more likely to occur in younger individuals and associate positively with dose of NAPBA and negatively with plasma glutamine and glycine levels. These results demonstrate that during routine clinical management, the majority of patients have PAA levels that are deemed safe. As age is negatively associated with PAA levels however, children undergoing treatment with NaPBA may need close monitoring of their phenylbutyrate metabolite levels.

Prevalence of DDC genotypes in patients with aromatic L-amino acid decarboxylase (AADC) deficiency and in silico prediction of structural protein changes.

Aromatic L-amino acid decarboxylase (AADC) deficiency is a rare autosomal recessive genetic disorder affecting the biosynthesis of dopamine, a precursor of both norepinephrine and epinephrine, and serotonin. Diagnosis is based on the analysis of CSF or plasma metabolites, AADC activity in plasma and genetic testing for variants in the DDC gene. The exact prevalence of AADC deficiency, the number of patients, and the variant and genotype prevalence are not known. Here, we present the DDC variant (n = 143) and genotype (n = 151) prevalence of 348 patients with AADC deficiency, 121 of whom were previously not reported. In addition, we report 26 new DDC variants, classify them according to the ACMG/AMP/ACGS recommendations for pathogenicity and score them based on the predicted structural effect. The splice variant c.714+4A>T, with a founder effect in Taiwan and China, was the most common variant (allele frequency = 32.4%), and c.[714+4A>T];[714+4A>T] was the most common genotype (genotype frequency = 21.3%). Approximately 90% of genotypes had variants classified as pathogenic or likely pathogenic, while 7% had one VUS allele and 3% had two VUS alleles. Only one benign variant was reported. Homozygous and compound heterozygous genotypes were interpreted in terms of AADC protein and categorized as: i) devoid of full-length AADC, ii) bearing one type of AADC homodimeric variant or iii) producing an AADC protein population composed of two homodimeric and one heterodimeric variant. Based on structural features, a score was attributed for all homodimers, and a tentative prediction was advanced for the heterodimer. Almost all AADC protein variants were pathogenic or likely pathogenic.

Clinical Untargeted Metabolomics as a Functional Screen to Improve Variant Classification.

With the rapid increase in clinical exome and genome sequencing, the number of variants of uncertain significance (VUS) that are reported continues to rise, which poses a significant barrier to interpretation of genetic findings. For metabolic disorders, biochemical testing can help alleviate this burden of variant interpretation by providing functional validation of uncertain genetic findings in many cases. However, a major limitation of traditional biochemical testing is the targeted, narrow range of analytes clinically available, resulting in delays in diagnosis if testing is negative. Untargeted metabolomic screening offers higher diagnostic yield and assays for thousands of metabolites across multiple metabolic pathways in a single test, saving time and resources for patients, families, and physicians. When integrated with exome or genome sequencing, untargeted metabolomic screening improves diagnostic outcomes by providing functional validation of genetic findings, particularly for VUS. Here, we present representative cases across the breadth of metabolic pathways as examples of the utility of metabolomics in genomic variant classification. © 2023 Wiley Periodicals LLC.

Association of Glucose-6-Phosphate Dehydrogenase Deficiency With Outcomes in US Veterans With COVID-19.

Importance: The underlying biological risk factors for severe outcome due to SAR-CoV-2 infection are not well defined. Objective: To determine the association between glucose-6-phosphate dehydrogenase (G6PD) deficiency and severity of COVID-19. Design, Setting, and Participants: This retrospective cohort study included analysis of 24 700 veterans with G6PD enzyme testing prior to January 1, 2020, obtained through the US Veterans Health Administration national databases. These veterans were cross-referenced with the Veterans Administration COVID-19 Shared Data Resource for SARS-CoV-2 testing from February 15, 2020, to January 1, 2021. The final study population consisted of 4811 veterans who tested positive for SARS-CoV-2. Statistical analysis was performed from June to December 2021. Exposures: G6PD deficiency. Main Outcomes and Measures: COVID-19 severe illness, as defined by the Centers for Disease Control and Prevention: hospitalization, need for mechanical ventilation and/or intensive care unit admission, or in-hospital mortality after a positive SARS-CoV-2 test. Results: Among 4811 veterans in the Veterans Health Administration who had historical G6PD enzyme activity test results and SARS-CoV-2 positivity included in this study, 3868 (80.4%) were male, 1553 (32.3%) were Black, and 1855 (39%) were White; 1228 (25.5%) were 65 years or older and 3583 (74.5%) were younger than 65 years. There were no significant differences in age, body mass index, or Charlson Comorbidity Index were present between the veterans with G6PD deficiency and without G6PD deficiency. Among these veterans with SARS-CoV-2 infection, G6PD deficiency was more prevalent in Black male veterans (309 of 454 [68.1%]) compared with other racial and ethnic groups. Black male veterans less than 65 years of age with G6PD deficiency had approximately 1.5-fold increased likelihood of developing severe outcomes from SARS-CoV-2 infection compared with Black male veterans without G6PD deficiency (OR, 1.47; 95% CI, 1.03-2.09). In the small subset of White male veterans with G6PD deficiency, we observed an approximately 3.6-fold increased likelihood of developing severe outcomes from SARS-CoV-2 infection compared with White male veterans aged 65 years or older without G6PD deficiency (OR, 3.58; 95% CI, 1.64-7.80). This difference between veterans with and without G6PD deficiency was not observed in younger White male veterans or older Black male veterans, nor in smaller subsets of other male veterans or in female veterans of any age. Conclusions and Relevance: In this cohort study of COVID-19-positive veterans, Black male veterans less than 65 years of age and White male veterans 65 years of age or older with G6PD deficiency had an increased likelihood of developing severe COVID-19 compared with veterans without G6PD deficiency. These data indicate a need to consider the potential for G6PD deficiency prior to treatment of patients with SARS-CoV-2 infection as part of clinical strategies to mitigate severe outcomes.

Novel phenotype of aortic root dilatation and late-onset metabolic decompensation in a patient with TMEM70 deficiency.

TMEM70 deficiency causing mitochondrial complex V deficiency, nuclear type 2 (MIM: 614052) is the most common nuclear encoded defect affecting ATP synthase and has been well described in the literature as being characterized by neonatal or infantile onset of poor feeding, hypotonia, lethargy, respiratory compromise, heart failure, lactic acidosis, hyperammonemia, and 3-methylglutaconic aciduria progressing to a phenotype of developmental delay, failure to thrive, short stature, nonprogressive cardiomyopathy, microcephaly, facial dysmorphisms, hypospadias, persistent pulmonary hypertension of the newborn, and Wolff-Parkinson-White syndrome, as well as metabolic crises followed by developmental regression. The patient with TMEM70 deficiency herein reported has the unique presentation of aortic root dilatation, differing facial dysmorphisms, and no history of neonatal metabolic decompensation or developmental delay, as well as a plasma metabolomics signature, including elevated 3-methylglutaconic acid, 3-methylglutarylcarnitine, alanine, and lactate, in addition to the commonly described increased 3-methylglutaconic acid on urine organic acid analysis that helped aid in the diagnostic interpretation of variants of uncertain significance in TMEM70.

Metabolic individuality: Limitations, challenges, and potential for clinical utility.

In Nature Medicine, Surendran and colleagues recently reported the analysis of human plasma metabolomic data for 913 metabolites in ∼20,000 individuals, identifying 2,599 metabolite-genetic variant associations and >400 metabolite signatures comprised of jointly regulated metabolites. This extensive atlas of variant-metabolite relationships reveals novel genomic mechanisms driving metabolic phenotypes.

Expansion and mechanistic insights into de novo DEAF1 variants in DEAF1-associated neurodevelopmental disorders.

De novo deleterious and heritable biallelic mutations in the DNA binding domain (DBD) of the transcription factor deformed epidermal autoregulatory factor 1 (DEAF1) result in a phenotypic spectrum of disorders termed DEAF1-associated neurodevelopmental disorders (DAND). RNA-sequencing using hippocampal RNA from mice with conditional deletion of Deaf1 in the central nervous system indicate that loss of Deaf1 activity results in the altered expression of genes involved in neuronal function, dendritic spine maintenance, development, and activity, with reduced dendritic spines in hippocampal regions. Since DEAF1 is not a dosage-sensitive gene, we assessed the dominant negative activity of previously identified de novo variants and a heritable recessive DEAF1 variant on selected DEAF1-regulated genes in 2 different cell models. While no altered gene expression was observed in cells over-expressing the recessive heritable variant, the gene expression profiles of cells over-expressing de novo variants resulted in similar gene expression changes as observed in CRISPR-Cas9-mediated DEAF1-deleted cells. Altered expression of DEAF1-regulated genes was rescued by exogenous expression of WT-DEAF1 but not by de novo variants in cells lacking endogenous DEAF1. De novo heterozygous variants within the DBD of DEAF1 were identified in 10 individuals with a phenotypic spectrum including autism spectrum disorder, developmental delays, sleep disturbance, high pain tolerance, and mild dysmorphic features. Functional assays demonstrate these variants alter DEAF1 transcriptional activity. Taken together, this study expands the clinical phenotypic spectrum of individuals with DAND, furthers our understanding of potential roles of DEAF1 on neuronal function, and demonstrates dominant negative activity of identified de novo variants.

Spectrum of DDC variants causing aromatic l-amino acid decarboxylase (AADC) deficiency and pathogenicity interpretation using ACMG-AMP/ACGS recommendations.

Pathogenic variants in dopa decarboxylase (DDC), the gene encoding the aromatic l-amino acid decarboxylase (AADC) enzyme, lead to a severe deficiency of neurotransmitters, resulting in neurological, neuromuscular, and behavioral manifestations clinically characterized by developmental delays, oculogyric crises, dystonia, and severe neurologic dysfunction in infancy. Historically, therapy has been aimed at compensating for neurotransmitter abnormalities, but response to pharmacologic therapy varies, and in most cases, the therapy shows little or no benefit. A novel human DDC gene therapy was recently approved in the European Union that targets the underlying genetic cause of the disorder, providing a new treatment option for patients with AADC deficiency. However, the applicability of human DDC gene therapy depends on the ability of laboratories and clinicians to interpret the results of genetic testing accurately enough to diagnose the patient. An accurate interpretation of genetic variants depends in turn on expert-guided curation of locus-specific databases. The purpose of this research was to identify previously uncharacterized DDC variants that are of pathologic significance in AADC deficiency as well as characterize and curate variants of unknown significance (VUSs) to further advance the diagnostic accuracy of genetic testing for this condition. DDC variants were identified using existing databases and the literature. The pathogenicity of the variants was classified using modified American College of Medical Genetics and Genomics/Association for Molecular Pathology/Association for Clinical Genomic Science (ACMG-AMP/ACGS) criteria. To improve the current variant interpretation recommendations, in silico variant interpretation tools were combined with structural 3D modeling of protein variants and applied comparative analysis to predict the impact of the variant on protein function. A total of 422 variants were identified (http://biopku.org/home/pnddb.asp). Variants were identified on nearly all introns and exons of the DDC gene, as well as the 3' and 5' untranslated regions. The largest percentage of the identified variants (48%) were classified as missense variants. The molecular effects of these missense variants were then predicted, and the pathogenicity of each was classified using a number of variant effect predictors. Using ACMG-AMP/ACGS criteria, 7% of variants were classified as pathogenic, 32% as likely pathogenic, 58% as VUSs of varying subclassifications, 1% as likely benign, and 1% as benign. For 101 out of 108 reported genotypes, at least one allele was classified as pathogenic or likely pathogenic. In silico variant pathogenicity interpretation tools, combined with structural 3D modeling of variant proteins and applied comparative analysis, have improved the current DDC variant interpretation recommendations, particularly of VUSs.

Untargeted metabolomic profiling in a patient with glycogen storage disease Ib receiving empagliflozin treatment.

Glycogen storage disease type Ib (GSD-Ib) is a rare inborn error of glycogen metabolism uniquely associated with neutropenia and neutrophil dysfunction, causing severe infections, inflammatory bowel disease (IBD), and impaired wound healing. Recently, kidney sodium-glucose co-transporter-2 (SGLT2) inhibitors such as empagliflozin known to reduce plasma levels of 1,5-anhydroglucitol (1,5-AG) and its toxic derivatives in neutrophils, have been described as a new treatment option in case reports of patients with GSD-Ib from Europe and Asia. We report our experience with an 11-year-old girl with GSD-Ib presenting with short fasting hypoglycemia, neutropenia with neutrophil dysfunction, recurrent infections, suboptimal growth, iron-deficiency anemia, and IBD. Treatment with daily empagliflozin improved neutrophil counts and function with a significant reduction in G-CSF needs. Significant improvement in IBD has led to weight gain with improved nutritional markers and improved fasting tolerance. Reduction of maximum empagliflozin dose was needed due to arthralgia. No other significant side effects of empagliflozin were observed. This report uniquely highlights the novel use of untargeted metabolomics profiling for monitoring plasma levels of 1,5-AG to assess empagliflozin dose responsiveness and guide dietary management and G-CSF therapy. Clinical improvement correlated to rapid normalization of 1,5-AG levels in plasma sustained after dose reduction. In conclusion, empagliflozin appeared to be a safe treatment option for GSD-Ib-associated neutropenia and neutrophil dysfunction. Global untargeted metabolomics is an efficient method to assess biochemical responsiveness to treatment.

NAXE deficiency: A neurometabolic disorder of NAD(P)HX repair amenable for metabolic correction.

The NAD(P)HX repair system is a metabolite damage repair mechanism responsible for restoration of NADH and NADPH after their inactivation by hydration. Deficiency in either of its two enzymes, NAD(P)HX dehydratase (NAXD) or NAD(P)HX epimerase (NAXE), causes a fatal neurometabolic disorder characterized by decompensations precipitated by inflammatory stress. Clinical findings include rapidly progressive muscle weakness, ataxia, ophthalmoplegia, and motor and cognitive regression, while neuroimaging abnormalities are subtle or nonspecific, making a clinical diagnosis challenging. During stress, nonenzymatic conversion of NAD(P)H to NAD(P)HX increases, and in the absence of repair, NAD(P)H is depleted, and NAD(P)HX accumulates, leading to decompensation; however, the contribution of each to the metabolic derangement is not established. Herein, we summarize the clinical knowledge of NAXE deficiency from 30 cases and lessons learned about disease pathogenesis from cell cultures and model organisms and describe a metabolomics signature obtained by untargeted metabolomics analysis in one case at the time of crisis and after initiation of treatment. Overall, biochemical findings support a model of acute depletion of NAD+, signs of mitochondrial dysfunction, and altered lipidomics. These findings are further substantiated by untargeted metabolomics six months post-crisis showing that niacin supplementation reverses primary metabolomic abnormalities concurrent with improved clinical status.

Relationships between food-related behaviors, obesity, and medication use in individuals with Smith-Magenis syndrome.

BACKGROUND: Smith-Magenis syndrome (SMS) is a complex neurodevelopmental disorder that includes obesity and food-seeking/satiety-related behaviors. AIMS: This study examined associations between food-related/hyperphagic behaviors, weight, and medication use in individuals with SMS. METHODS/PROCEDURES: Caregivers of individuals with SMS in the Parents and Researchers Interested in SMS (PRISMS) Patient Registry completed a demographic/medication questionnaire, the Hyperphagia Questionnaire for Clinical Trials, and the Food Related Problems Questionnaire. OUTCOMES/RESULTS: Among 49 participants (Mage = 16.41 ± 12.73 years, range = 4-69 years, 55% girls/women), individuals with SMS with overweight/obesity (n = 22) had worse overall food-related problems including greater impaired satiety (p < 0.05), maladaptive eating behaviors (p < 0.05), inappropriate response (p < 0.01), and hyperphagia (p < 0.01) compared to individuals of normal/underweight (n = 27). Those taking anti-depressants/anxiolytics (n = 16) had greater maladaptive eating behaviors (p < 0.05), hyperphagic behaviors (p < 0.05), and hyperphagic severity (p < 0.05) than those not taking anti-depressants/anxiolytics (n = 33). Boys/men with SMS had greater maladaptive eating behaviors (p < 0.05), inappropriate response (p < 0.05), and hyperphagic drive (p < 0.01) than girls/women with SMS. CONCLUSIONS/IMPLICATIONS: Maladaptive food-related behaviors were higher in individuals with SMS with overweight/obesity, taking anti-depressants/anxiolytics, or who were male. Medications in this population should be chosen with weight-related side effects in mind.

Untargeted Metabolomics of Slc13a5 Deficiency Reveal Critical Liver-Brain Axis for Lipid Homeostasis.

Though biallelic variants in SLC13A5 are known to cause severe encephalopathy, the mechanism of this disease is poorly understood. SLC13A5 protein deficiency reduces citrate transport into the cell. Downstream abnormalities in fatty acid synthesis and energy generation have been described, though biochemical signs of these perturbations are inconsistent across SLC13A5 deficiency patients. To investigate SLC13A5-related disorders, we performed untargeted metabolic analyses on the liver, brain, and serum from a Slc13a5-deficient mouse model. Metabolomic data were analyzed using the connect-the-dots (CTD) methodology and were compared to plasma and CSF metabolomics from SLC13A5-deficient patients. Mice homozygous for the Slc13a5tm1b/tm1b null allele had perturbations in fatty acids, bile acids, and energy metabolites in all tissues examined. Further analyses demonstrated that for several of these molecules, the ratio of their relative tissue concentrations differed widely in the knockout mouse, suggesting that deficiency of Slc13a5 impacts the biosynthesis and flux of metabolites between tissues. Similar findings were observed in patient biofluids, indicating altered transport and/or flux of molecules involved in energy, fatty acid, nucleotide, and bile acid metabolism. Deficiency of SLC13A5 likely causes a broader state of metabolic dysregulation than previously recognized, particularly regarding lipid synthesis, storage, and metabolism, supporting SLC13A5 deficiency as a lipid disorder.

Clinical diagnosis of metabolic disorders using untargeted metabolomic profiling and disease-specific networks learned from profiling data.

Untargeted metabolomics is a global molecular profiling technology that can be used to screen for inborn errors of metabolism (IEMs). Metabolite perturbations are evaluated based on current knowledge of specific metabolic pathway deficiencies, a manual diagnostic process that is qualitative, has limited scalability, and is not equipped to learn from accumulating clinical data. Our purpose was to improve upon manual diagnosis of IEMs in the clinic by developing novel computational methods for analyzing untargeted metabolomics data. We employed CTD, an automated computational diagnostic method that "connects the dots" between metabolite perturbations observed in individual metabolomics profiling data and modules identified in disease-specific metabolite co-perturbation networks learned from prior profiling data. We also extended CTD to calculate distances between any two individuals (CTDncd) and between an individual and a disease state (CTDdm), to provide additional network-quantified predictors for use in diagnosis. We show that across 539 plasma samples, CTD-based network-quantified measures can reproduce accurate diagnosis of 16 different IEMs, including adenylosuccinase deficiency, argininemia, argininosuccinic aciduria, aromatic L-amino acid decarboxylase deficiency, cerebral creatine deficiency syndrome type 2, citrullinemia, cobalamin biosynthesis defect, GABA-transaminase deficiency, glutaric acidemia type 1, maple syrup urine disease, methylmalonic aciduria, ornithine transcarbamylase deficiency, phenylketonuria, propionic acidemia, rhizomelic chondrodysplasia punctata, and the Zellweger spectrum disorders. Our approach can be used to supplement information from biochemical pathways and has the potential to significantly enhance the interpretation of variants of uncertain significance uncovered by exome sequencing. CTD, CTDdm, and CTDncd can serve as an essential toolset for biological interpretation of untargeted metabolomics data that overcomes limitations associated with manual diagnosis to assist diagnosticians in clinical decision-making. By automating and quantifying the interpretation of perturbation patterns, CTD can improve the speed and confidence by which clinical laboratory directors make diagnostic and treatment decisions, while automatically improving performance with new case data.