Clinical, Neuroimaging, and Metabolic Footprint of the Neurodevelopmental Disorder Caused by Monoallelic HK1 Variants.

Background and Objectives: Hexokinase 1 (encoded by HK1) catalyzes the first step of glycolysis, the adenosine triphosphate-dependent phosphorylation of glucose to glucose-6-phosphate. Monoallelic HK1 variants causing a neurodevelopmental disorder (NDD) have been reported in 12 individuals. Methods: We investigated clinical phenotypes, brain MRIs, and the CSF of 15 previously unpublished individuals with monoallelic HK1 variants and an NDD phenotype. Results: All individuals had recurrent variants likely causing gain-of-function, representing mutational hot spots. Eight individuals (c.1370C>T) had a developmental and epileptic encephalopathy with infantile onset and virtually no development. Of the other 7 individuals (n = 6: c.1334C>T; n = 1: c.1240G>A), 3 adults showed a biphasic course of disease with a mild static encephalopathy since early childhood and an unanticipated progressive deterioration with, e.g., movement disorder, psychiatric disease, and stroke-like episodes, epilepsy, starting in adulthood. Individuals who clinically presented in the first months of life had (near)-normal initial neuroimaging and severe cerebral atrophy during follow-up. In older children and adults, we noted progressive involvement of basal ganglia including Leigh-like MRI patterns and cerebellar atrophy, with remarkable intraindividual variability. The CSF glucose and the CSF/blood glucose ratio were below the 5th percentile of normal in almost all CSF samples, while blood glucose was unremarkable. This biomarker profile resembles glucose transporter type 1 deficiency syndrome; however, in HK1-related NDD, CSF lactate was significantly increased in all patients resulting in a substantially different biomarker profile. Discussion: Genotype-phenotype correlations appear to exist for HK1 variants and can aid in counseling. A CSF biomarker profile with low glucose, low CSF/blood glucose, and high CSF lactate may point toward monoallelic HK1 variants causing an NDD. This can help in variant interpretation and may aid in understanding the pathomechanism. We hypothesize that progressive intoxication and/or ongoing energy deficiency lead to the clinical phenotypes and progressive neuroimaging findings.

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.

A spectroscopic test suggests that fragment ion structure annotations in MS/MS libraries are frequently incorrect.

Modern untargeted mass spectrometry (MS) analyses quickly detect and resolve thousands of molecular compounds. Although features are readily annotated with a molecular formula in high-resolution small-molecule MS applications, the large majority of them remains unidentified in terms of their full molecular structure. Collision-induced dissociation tandem mass spectrometry (CID-MS2) provides a diagnostic molecular fingerprint to resolve the molecular structure through a library search. However, for de novo identifications, one must often rely on in silico generated MS2 spectra as reference. The ability of different in silico algorithms to correctly predict MS2 spectra and thus to retrieve correct molecular structures is a topic of lively debate, for instance in the CASMI contest. Underlying the predicted MS2 spectra are the in silico generated product ion structures, which are normally not used in de novo identification, but which can serve to critically assess the fragmentation algorithms. Here we evaluate in silico generated MSn product ion structures by comparison with structures established experimentally by infrared ion spectroscopy (IRIS). For a set of three dozen product ion structures from five precursor molecules, we find that virtually all fragment ion structure annotations in three major in silico MS2 libraries (HMDB, METLIN, mzCloud) are incorrect and caution the reader against their use for structure annotation of MS/MS ions.

Discovery of novel diagnostic biomarkers for Sjögren-Larsson syndrome by untargeted lipidomics.

AIM: Sjögren-Larsson syndrome (SLS) is a rare neurometabolic disorder that mainly affects brain, eye and skin and is caused by deficiency of fatty aldehyde dehydrogenase. Our recent finding of a profoundly disturbed brain tissue lipidome in SLS prompted us to search for similar biomarkers in plasma as no functional test in blood is available for SLS. METHODS AND RESULTS: We performed plasma lipidomics and used a newly developed bioinformatics tool to mine the untargeted part of the SLS plasma and brain lipidome to search for SLS biomarkers. Plasma lipidomics showed disturbed ether lipid metabolism in known lipid classes. Untargeted lipidomics of both plasma and brain (white and grey matter) uncovered two new endogenous lipid classes highly elevated in SLS. The first biomarker group were alkylphosphocholines/ethanolamines containing different lengths of alkyl-chains where some alkylphosphocholines were > 600-fold elevated in SLS plasma. The second group of biomarkers were a set of 5 features of unknown structure. Fragmentation studies suggested that they contain ubiquinol and phosphocholine and one feature was also found as a glucuronide conjugate in plasma. The plasma features were highly distinctive for SLS with levels >100-1000-fold the level in controls, if present at all. We speculate on the origin of the alkylphosphocholines/ethanolamines and the nature of the ubiquinol-containing metabolites. CONCLUSIONS: The metabolites identified in this study represent novel endogenous lipid classes thus far unknown in humans. They represent the first plasma metabolite SLS-biomarkers and may also yield more insight into SLS pathophysiology.

CYP2U1: An emerging treatable neurometabolic disease with cerebral folate deficiency in 2 Chinese brothers.

With the rapid advancement of medical technologies in genomic and molecular medicine, the number of treatable neurometabolic diseases is quickly expanding. Spastic paraplegia 56 (SPG56), one of the severe autosomal recessive forms of neurodegenerative disorders caused by pathogenic variants in the CYP2U1 gene, has no reported specific targeted treatment yet. Here we report 2 Chinese brothers with CYP2U1 bi-allelic pathogenic variants with cerebral folate deficiency who were treated for over a decade with folinic acid supplement. Patients have remained stable under therapy.

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).

A Delphi Survey Study to Formulate Statements on the Treatability of Inherited Metabolic Disorders to Decide on Eligibility for Newborn Screening.

The Wilson and Jungner (W&J) and Andermann criteria are meant to help select diseases eligible for population-based screening. With the introduction of next-generation sequencing (NGS) methods for newborn screening (NBS), more inherited metabolic diseases (IMDs) can technically be included, and a revision of the criteria was attempted. This study aimed to formulate statements and investigate whether those statements could elaborate on the criterion of treatability for IMDs to decide on eligibility for NBS. An online Delphi study was started among a panel of Dutch IMD experts (EPs). EPs evaluated, amended, and approved statements on treatability that were subsequently applied to 10 IMDs. After two rounds of Delphi, consensus was reached on 10 statements. Application of these statements selected 5 out of 10 IMDs proposed for this study as eligible for NBS, including 3 IMDs in the current Dutch NBS. The statement: 'The expected benefit/burden ratio of early treatment is positive and results in a significant health outcome' contributed most to decision-making. Our Delphi study resulted in 10 statements that can help to decide on eligibility for inclusion in NBS based on treatability, also showing that other criteria could be handled in a comparable way. Validation of the statements is required before these can be applied as guidance to authorities.

Exome sequencing links the SUMO protease SENP7 with fatal arthrogryposis multiplex congenita, early respiratory failure and neutropenia.

BACKGROUND: SUMOylation involves the attachment of small ubiquitin-like modifier (SUMO) proteins to specific lysine residues on thousands of substrates with target-specific effects on protein function. Sentrin-specific proteases (SENPs) are proteins involved in the maturation and deconjugation of SUMO. Specifically, SENP7 is responsible for processing polySUMO chains on targeted substrates including the heterochromatin protein 1α (HP1α). METHODS: We performed exome sequencing and segregation studies in a family with several infants presenting with an unidentified syndrome. RNA and protein expression studies were performed in fibroblasts available from one subject. RESULTS: We identified a kindred with four affected subjects presenting with a spectrum of findings including congenital arthrogryposis, no achievement of developmental milestones, early respiratory failure, neutropenia and recurrent infections. All died within four months after birth. Exome sequencing identified a homozygous stop gain variant in SENP7 c.1474C>T; p.(Gln492*) as the probable aetiology. The proband's fibroblasts demonstrated decreased mRNA expression. Protein expression studies showed significant protein dysregulation in total cell lysates and in the chromatin fraction. We found that HP1α levels as well as different histones and H3K9me3 were reduced in patient fibroblasts. These results support previous studies showing interaction between SENP7 and HP1α, and suggest loss of SENP7 leads to reduced heterochromatin condensation and subsequent aberrant gene expression. CONCLUSION: Our results suggest a critical role for SENP7 in nervous system development, haematopoiesis and immune function in humans.

Succinic semialdehyde dehydrogenase deficiency in mice and in humans: An untargeted metabolomics perspective.

Succinic semialdehyde dehydrogenase deficiency (SSADHD) is a rare neurometabolic disorder caused by disruption of the gamma-aminobutyric acid (GABA) pathway. A more detailed understanding of its pathophysiology, beyond the accumulation of GABA and gamma-hydroxybutyric acid (GHB), will increase our understanding of the disease and may support novel therapy development. To this end, we compared biochemical body fluid profiles from SSADHD patients with controls using next-generation metabolic screening (NGMS). Targeted analysis of NGMS data from cerebrospinal fluid (CSF) showed a moderate increase of aspartic acid, glutaric acid, glycolic acid, 4-guanidinobutanoic acid, and 2-hydroxyglutaric acid, and prominent elevations of GHB and 4,5-dihydroxyhexanoic acid (4,5-DHHA) in SSADHD samples. Remarkably, the intensities of 4,5-DHHA and GHB showed a significant positive correlation in control CSF, but not in patient CSF. In an established zebrafish epilepsy model, 4,5-DHHA showed increased mobility that may reflect limited epileptogenesis. Using untargeted metabolomics, we identified 12 features in CSF with high biomarker potential. These had comparable increased fold changes as GHB and 4,5-DHHA. For 10 of these features, a similar increase was found in plasma, urine and/or mouse brain tissue for SSADHD compared to controls. One of these was identified as the novel biomarker 4,5-dihydroxyheptanoic acid. The intensities of selected features in plasma and urine of SSADHD patients positively correlated with the clinical severity score of epilepsy and psychiatric symptoms of those patients, and also showed a high mutual correlation. Our findings provide new insights into the (neuro)metabolic disturbances in SSADHD and give leads for further research concerning SSADHD pathophysiology.

An In Silico Infrared Spectral Library of Molecular Ions for Metabolite Identification.

Infrared ion spectroscopy (IRIS) continues to see increasing use as an analytical tool for small-molecule identification in conjunction with mass spectrometry (MS). The IR spectrum of an m/z selected population of ions constitutes a unique fingerprint that is specific to the molecular structure. However, direct translation of an IR spectrum to a molecular structure remains challenging, as reference libraries of IR spectra of molecular ions largely do not exist. Quantum-chemically computed spectra can reliably be used as reference, but the challenge of selecting the candidate structures remains. Here, we introduce an in silico library of vibrational spectra of common MS adducts of over 4500 compounds found in the human metabolome database. In total, the library currently contains more than 75,000 spectra computed at the DFT level that can be queried with an experimental IR spectrum. Moreover, we introduce a database of 189 experimental IRIS spectra, which is employed to validate the automated spectral matching routines. This demonstrates that 75% of the metabolites in the experimental data set are correctly identified, based solely on their exact m/z and IRIS spectrum. Additionally, we demonstrate an approach for specifically identifying substructures by performing a search without m/z constraints to find structural analogues. Such an unsupervised search paves the way toward the de novo identification of unknowns that are absent in spectral libraries. We apply the in silico spectral library to identify an unknown in a plasma sample as 3-hydroxyhexanoic acid, highlighting the potential of the method.

Distinguishing Oligosaccharide Isomers Using Far-Infrared Ion Spectroscopy: Identification of Biomarkers for Inborn Errors of Metabolism.

Distinguishing isomeric saccharides poses a major challenge for analytical workflows based on (liquid chromatography) mass spectrometry (LC-MS). In recent years, many studies have proposed infrared ion spectroscopy as a possible solution as the orthogonal, spectroscopic characterization of mass-selected ions can often distinguish isomeric species that remain unresolved using conventional MS. However, the high conformational flexibility and extensive hydrogen bonding in saccharides cause their room-temperature fingerprint infrared spectra to have broad features that often lack diagnostic value. Here, we show that room-temperature infrared spectra of ion-complexed saccharides recorded in the previously unexplored far-infrared wavelength range (300-1000 cm-1) provide well-resolved and highly diagnostic features. We show that this enables distinction of isomeric saccharides that differ either by their composition of monosaccharide units and/or the orientation of their glycosidic linkages. We demonstrate the utility of this approach from single monosaccharides up to isomeric tetrasaccharides differing only by the configuration of a single glycosidic linkage. Furthermore, through hyphenation with hydrophilic interaction liquid chromatography, we identify oligosaccharide biomarkers in patient body fluid samples, demonstrating a generalized and highly sensitive MS-based method for the identification of saccharides found in complex sample matrices.

Targeted Small-Molecule Identification Using Heartcutting Liquid Chromatography-Infrared Ion Spectroscopy.

Infrared ion spectroscopy (IRIS) can be used to identify molecular structures detected in mass spectrometry (MS) experiments and has potential applications in a wide range of analytical fields. However, MS-based approaches are often combined with orthogonal separation techniques, in many cases liquid chromatography (LC). The direct coupling of LC and IRIS is challenging due to the mismatching timescales of the two technologies: an IRIS experiment typically takes several minutes, whereas an LC fraction typically elutes in several seconds. To resolve this discrepancy, we present a heartcutting LC-IRIS approach using a setup consisting of two switching valves and two sample loops as an alternative to direct online LC-IRIS coupling. We show that this automated setup enables us to record multiple IR spectra for two LC-features from a single injection without degrading the LC-separation performance. We demonstrate the setup for application in drug metabolism research by recording six m/z-selective IR spectra for two drug metabolites from a single 2 µL sample of cell incubation extract. Additionally, we measure the IR spectra of two closely eluting diastereomeric biomarkers for the inborn error of metabolism pyridoxine-dependent epilepsy (PDE-ALDH7A1), which shows that the heartcutting LC-IRIS setup has good sensitivity (requiring ∼µL injections of ∼µM samples) and that the separation between closely eluting isomers is maintained. We envision applications in a range of research fields, where the identification of molecular structures detected by LC-MS is required.

Newborn screening for Cerebrotendinous Xanthomatosis: A retrospective biomarker study using both flow-injection and UPLC-MS/MS analysis in 20,000 newborns.

BACKGROUND AND AIMS: Cerebrotendinous Xanthomatosis (CTX) is a treatable disorder of bile acid synthesis caused by deficiency of 27-sterol hydroxylase -encoded by CYP27A1- leading to gastrointestinal and progressive neuropsychiatric symptoms. Biochemically, CTX is characterized by accumulation of the bile alcohol cholestanetetrol glucuronide (GlcA-tetrol) and the deficiency of tauro-chenodeoxycholic acid (t-CDCA) and tauro-trihydroxycholestanoic acid (t-THCA). MATERIALS AND METHODS: To ascertain the feasibility of CTX newborn screening (NBS) we performed a study with deidentified Dutch dried blood spots using reagents and equipment that is frequently used in NBS laboratories. 20,076 deidentified newborn blood spots were subjected to flow-injection (FIA)-MS/MS and UPLC-MS/MS analysis to determine the concentration of GlcA-tetrol and calculate the GlcA-tetrol/t-CDCA and t-THCA/GlcA-tetrol ratios. RESULTS: Using UPLC-MS/MS analysis both GlcA-tetrol concentration and/or metabolite ratios GlcA-tetrol/t-CDCA proved to be informative biomarkers; newborn DBS results did not overlap with those of the CTX patients. For FIA-MS/MS, GlcA-tetrol also was an excellent marker but when using the combination of the GlcA-tetrol/t-CDCA and t-THCA/GlcA-tetrol ratios also did not yield any screen positives. CONCLUSION: Newborn screening for CTX using only metabolite ratios following the measurement of three CTX biomarkers is possible using either FIA-MS/MS or UPLC-MS/MS, which paves the way for introduction of CTX NBS.

Novel cerebrospinal fluid biomarkers of glucose transporter type 1 deficiency syndrome: Implications beyond the brain's energy deficit.

We used next-generation metabolic screening to identify new biomarkers for improved diagnosis and pathophysiological understanding of glucose transporter type 1 deficiency syndrome (GLUT1DS), comparing metabolic cerebrospinal fluid (CSF) profiles from 12 patients to those of 116 controls. This confirmed decreased CSF glucose and lactate levels in patients with GLUT1DS and increased glutamine at group level. We identified three novel biomarkers significantly decreased in patients, namely gluconic + galactonic acid, xylose-α1-3-glucose, and xylose-α1-3-xylose-α1-3-glucose, of which the latter two have not previously been identified in body fluids. CSF concentrations of gluconic + galactonic acid may be reduced as these metabolites could serve as alternative substrates for the pentose phosphate pathway. Xylose-α1-3-glucose and xylose-α1-3-xylose-α1-3-glucose may originate from glycosylated proteins; their decreased levels are hypothetically the consequence of insufficient glucose, one of two substrates for O-glucosylation. Since many proteins are O-glucosylated, this deficiency may affect cellular processes and thus contribute to GLUT1DS pathophysiology. The novel CSF biomarkers have the potential to improve the biochemical diagnosis of GLUT1DS. Our findings imply that brain glucose deficiency in GLUT1DS may cause disruptions at the cellular level that go beyond energy metabolism, underlining the importance of developing treatment strategies that directly target cerebral glucose uptake.

Identification of Δ-1-pyrroline-5-carboxylate derived biomarkers for hyperprolinemia type II.

Hyperprolinemia type II (HPII) is an inborn error of metabolism due to genetic variants in ALDH4A1, leading to a deficiency in Δ-1-pyrroline-5-carboxylate (P5C) dehydrogenase. This leads to an accumulation of toxic levels of P5C, an intermediate in proline catabolism. The accumulating P5C spontaneously reacts with, and inactivates, pyridoxal 5'-phosphate, a crucial cofactor for many enzymatic processes, which is thought to be the pathophysiological mechanism for HPII. Here, we describe the use of a combination of LC-QTOF untargeted metabolomics, NMR spectroscopy and infrared ion spectroscopy (IRIS) to identify and characterize biomarkers for HPII that result of the spontaneous reaction of P5C with malonic acid and acetoacetic acid. We show that these biomarkers can differentiate between HPI, caused by a deficiency of proline oxidase activity, and HPII. The elucidation of their molecular structures yields insights into the disease pathophysiology of HPII.

Abnormal concentrations of acetylated amino acids in cerebrospinal fluid in acetyl-CoA transporter deficiency.

Acetyl-CoA transporter 1 (AT-1) is a transmembrane protein which regulates influx of acetyl-CoA from the cytosol to the lumen of the endoplasmic reticulum and is therefore important for the posttranslational modification of numerous proteins. Pathological variants in the SLC33A1 gene coding for AT-1 have been linked to a disorder called Huppke-Brendel syndrome, which is characterized by congenital cataracts, hearing loss, severe developmental delay and early death. It has been described in eight patients so far, who all had the abovementioned symptoms together with low serum copper and ceruloplasmin concentrations. The link between AT-1 and low ceruloplasmin concentrations is not clear, nor is the complex pathogenesis of the disease. Here we describe a further case of Huppke-Brendel syndrome with a novel and truncating homozygous gene variant and provide novel biochemical data on N-acetylated amino acids in cerebrospinal fluid (CSF) and plasma. Our results indicate that decreased levels of many N-acetylated amino acids in CSF are a typical metabolic fingerprint for AT-1 deficiency and are potential biomarkers for the defect. As acetyl-CoA is an important substrate for protein acetylation, we performed N-terminal proteomics, but found only minor effects on this particular protein modification. The acetyl-CoA content in patient's fibroblasts was insignificantly decreased. Our data may help to better understand the mechanisms underlying the metabolic disturbances, the pathophysiology and the clinical phenotype of the disease.

Long-term MRI Findings in Patients With Cerebrotendinous Xanthomatosis Treated With Chenodeoxycholic Acid.

OBJECTIVES: To describe long-term follow-up brain MRI findings in patients with cerebrotendinous xanthomatosis (CTX) treated with chenodeoxycholic acid (CDCA). METHODS: Of a cohort of 79 Dutch patients with CTX, we retrospectively reviewed brain MRI findings of patients at diagnosis (before the start of treatment) and after long-term follow-up (7-27 years) in 12 patients. In addition, we report on 2 families with remarkable brain MRI findings. RESULTS: MRI abnormalities showed progression in all 7 patients diagnosed at 24 years or older and only in 1 of 5 patients diagnosed younger than 24 years. MRI findings in the other patients diagnosed younger than 24 years were normal at baseline and remained normal even after follow-up of more than 25 years. The total MRI scores at baseline were 2 and 19 and at follow-up 4 and 37, respectively, for patients diagnosed before or after the age of 24 years, despite a comparable number of treatment years. DISCUSSION: MRI findings are fully in line with our long-term treatment effect article, emphasizing the importance of early diagnosis and treatment in CTX. Expanding the spectrum of brain MRI findings (including the finding of a posterior leukoencephalopathy) leads to a better understanding of the heterogeneity of this treatable disease.

Preanalytical Pitfalls in Untargeted Plasma Nuclear Magnetic Resonance Metabolomics of Endocrine Hypertension.

Despite considerable morbidity and mortality, numerous cases of endocrine hypertension (EHT) forms, including primary aldosteronism (PA), pheochromocytoma and functional paraganglioma (PPGL), and Cushing's syndrome (CS), remain undetected. We aimed to establish signatures for the different forms of EHT, investigate potentially confounding effects and establish unbiased disease biomarkers. Plasma samples were obtained from 13 biobanks across seven countries and analyzed using untargeted NMR metabolomics. We compared unstratified samples of 106 PHT patients to 231 EHT patients, including 104 PA, 94 PPGL and 33 CS patients. Spectra were subjected to a multivariate statistical comparison of PHT to EHT forms and the associated signatures were obtained. Three approaches were applied to investigate and correct confounding effects. Though we found signatures that could separate PHT from EHT forms, there were also key similarities with the signatures of sample center of origin and sample age. The study design restricted the applicability of the corrections employed. With the samples that were available, no biomarkers for PHT vs. EHT could be identified. The complexity of the confounding effects, evidenced by their robustness to correction approaches, highlighted the need for a consensus on how to deal with variabilities probably attributed to preanalytical factors in retrospective, multicenter metabolomics studies.