Next Article in Journal
Multi-Omics Integrative Analyses Identified Two Endotypes of Hip Osteoarthritis
Previous Article in Journal
Sitting Interruption Modalities during Prolonged Sitting Acutely Improve Postprandial Metabolome in a Crossover Pilot Trial among Postmenopausal Women
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Phenylketonuria on the Serum Metabolome and Plasma Lipidome: A Study in Early-Treated Patients

by
Jorine C. van der Weerd
1,
Annemiek M. J. van Wegberg
2,
Theo S. Boer
1,
Udo F. H. Engelke
3,
Karlien L. M. Coene
3,4,
Ron A. Wevers
3,
Stephan J. L. Bakker
5,
Pim de Blaauw
1,
Joost Groen
1,
Francjan J. van Spronsen
2 and
M. Rebecca Heiner-Fokkema
1,*
1
Department of Laboratory Medicine, Laboratory of Metabolic Disease, University Medical Center Groningen, University of Groningen, 9700 RB Groningen, The Netherlands
2
Division of Metabolic Diseases, Beatrix Children’s Hospital, University Medical Center Groningen, University of Groningen, 9700 RB Groningen, The Netherlands
3
Department of Human Genetics, Translational Metabolic Laboratory (TML), Radboud University Medical Center, 6525 GA Nijmegen, The Netherlands
4
Laboratory of Clinical Chemistry and Hematology, Máxima Medical Centre, 5504 DB Veldhoven, The Netherlands
5
Division of Nephrology, Department of Internal Medicine, University Medical Center Groningen, University of Groningen, 9700 RB Groningen, The Netherlands
*
Author to whom correspondence should be addressed.
Metabolites 2024, 14(9), 479; https://doi.org/10.3390/metabo14090479
Submission received: 18 July 2024 / Revised: 22 August 2024 / Accepted: 26 August 2024 / Published: 30 August 2024
(This article belongs to the Section Endocrinology and Clinical Metabolic Research)

Abstract

:
Background: Data suggest that metabolites, other than blood phenylalanine (Phe), better and independently predict clinical outcomes in patients with phenylketonuria (PKU). Methods: To find new biomarkers, we compared the results of untargeted lipidomics and metabolomics in treated adult PKU patients to those of matched controls. Samples (lipidomics in EDTA-plasma (22 PKU and 22 controls) and metabolomics in serum (35 PKU and 20 controls)) were analyzed using ultra-high-performance liquid chromatography and high-resolution mass spectrometry. Data were subjected to multivariate (PCA, OPLS-DA) and univariate (Mann–Whitney U test, p < 0.05) analyses. Results: Levels of 33 (of 20,443) lipid features and 56 (of 5885) metabolite features differed statistically between PKU patients and controls. For lipidomics, findings include higher glycerolipids, glycerophospholipids, and sphingolipids species. Significantly lower values were found for sterols and glycerophospholipids species. Seven features had unknown identities. Total triglyceride content was higher. Higher Phe and Phe catabolites, tryptophan derivatives, pantothenic acid, and dipeptides were observed for metabolomics. Ornithine levels were lower. Twenty-six metabolite features were not annotated. Conclusions: This study provides insight into the metabolic phenotype of PKU patients. Additional studies are required to establish whether the observed changes result from PKU itself, diet, and/or an unknown reason.

1. Introduction

Phenylketonuria (PKU, McKusick/OMIM# 261600) is a rare autosomal recessive disorder affecting aromatic amino acid metabolism. Worldwide, PKU has an estimated incidence of 1 in 23,930 live births [1]. The disease is characterized by mutations in the gene encoding phenylalanine hydroxylase (PAH; EC 1.14.16.1; OMIM* 612349), an enzyme responsible for converting L-phenylalanine (Phe) into L-tyrosine (Tyr). In individuals with PKU, the activity of the enzyme PAH is significantly reduced or absent, leading to hyperphenylalaninemia (HPA), which is the accumulation of Phe in blood and brain tissue. HPA is not specific for PKU because it may also be a consequence of other hereditary causes, like defects in the pterin metabolism, resulting in a deficiency of the PAH cosubstrate tetrahydrobiopterin (BH4) and DNAJC12 [2,3]. Depending on the residual enzyme activity, patients have a widely varying Phe tolerance. If late-treated, PKU can lead to progressive and irreversible damage to the brain, resulting in intellectual disability [4,5,6,7,8,9,10,11,12].
Systematic measurement of Phe in dried blood spots or plasma is crucial to the diagnosis and monitoring of treatment of PKU patients [13]. Management of PKU primarily consists of a natural protein-restricted diet enriched with minerals, vitamins, and other micronutrients and/or low Phe casein glycomacropeptide to keep blood Phe concentrations within an acceptable range [2,13,14,15]. Supplementation is essential for meeting nutritional needs, supporting normal growth, and preventing developmental deficiencies [2,13,14,15]. A subset of patients may benefit from sapropterin dihydrochloride, an orally active synthetic form of tetrahydrobiopterin that enhances the residual activity of PAH. All patients may benefit from pegvaliase, which converts Phe into harmless compounds that can be excreted from the body. Both drugs lower Phe concentrations and enable patients to relax their protein restriction. However, pegvaliase is only reimbursed in a minority of countries around the world [16,17,18]. Except perhaps for treatment with large neutral amino acids (LNAAs), intensive Phe monitoring is crucial but is not able to predict the clinical outcome fully [19]. For this reason, there is a need to expand the biomarker signature beyond blood Phe to better reflect neurocognitive and psychosocial function.
Over the years, high-resolution mass spectrometry (HRMS) has become a powerful tool for discovering novel disease biomarkers for diagnosis and monitoring various inborn errors of metabolism, including PKU [20,21]. This technique features high resolution and, thus, high mass accuracy, enabling comprehensive profiling of known and unknown metabolites in complex biological matrices. In the last decade, the impact of PAH deficiency on the metabolome has become evident [21,22]. The levels of Tyr and downstream metabolites are lower, while the alternative Phe degradation pathway results in higher levels of phenylpyruvate and other (neuro)toxic metabolic products. In contrast, the non-polar part of the metabolome remains rather unexplored. Guerra et al. [23] emphasized the need to study the lipidome in this group of patients. Therefore, we compared the plasma lipid and serum metabolite signature of adult PKU patients with non-PKU controls.

2. Materials and Methods

Chemicals. The following chemicals were used for lipidomics analyses: methyl tert-butyl ether UPLC-grade (MTBE; Biosolve BV, Valkenswaard, The Netherlands), methanol absolute UPLC-grade (MeOH; Biosolve BV, Valkenswaard, The Netherlands), Milli-Q (MQ; Millipore, MA, USA), chloroform (CHCl3; Merck KGaA, Darmstadt, Germany), acetonitrile (ACN; Biosolve BV, Valkenswaard, The Netherlands), 2-propanol (IPA; Biosolve BV, Valkenswaard, The Netherlands), and ammonium formate UPLC-grade (AmF; Sigma-Aldrich, Zwijndrecht, The Netherlands). The following non-physiologic or stable isotope-labeled lipid standards were purchased from Avanti Polar Lipids, Inc (Alabaster, AL, USA): LPC(17:1/0:0), PC(17:0/17:0), LPE(17:1/0:0), PE(17:0/17:0), LPG(17:0/0:0), PG(17:0/17:0), LPS(17:0/0:0), PS(17:0/17:0), CL(14:0.14:0)(14:0/14:0), CL(16:0/16:0)/(16:0/16:0), DG(16:0/16:0/0:0), and TG(17:0/17:1/17:0)d5. All lipids were dissolved in CHCl3/MeOH/MQ (60:30:4.5, v/v/v) or chloroform to a final concentration of 20 µM and stored at −20 °C (compound details are provided in Table S1). Chemicals used for metabolomics analyses are described elsewhere [21].
PKU samples. The PKU-COBESO study, a multicenter research project conducted in The Netherlands, focused on exploring the cognitive, behavioral, and social outcomes in relation to metabolic control in early-diagnosed and treated PKU patients [24,25,26,27]. The study included children and adults (age [7–42 years]). Participants were included between 2012 and 2015. Venous blood samples were taken after an overnight fast on the day of testing before the neuropsychological assessment. The EDTA- and serum samples were stored for up to 11 years at −80 °C without thawing before analysis. In the current study, individuals under the age of 18 were excluded from participation.
Control samples. No blood samples were taken from the control subjects in the COBESO study. For this reason, blood samples of the control group for the current study were sourced from pre-transplant kidney donors through the TransplantLines Biobank [28]. Transplantlines is a longitudinal cohort study investigating the outcomes after transplantation in patients in the UMCG. We consider this control group to be a good representative of the general population, with the sole selection criteria as kidney donors that they have two kidneys and a normal kidney function. We opted to use control samples from this cohort study rather than obtaining new samples from healthy individuals to partially account for any potential impacts of extended frozen storage. The EDTA- and serum samples from the TransplantLines Biobank, which were gathered between 2016 and 2020, have been preserved at −80 °C before their analysis.
For the statistical analysis, samples were matched by age and gender. In lipidomics, 22 PKU samples (median age: 33, age range: 21–43, gender ratio: 14M/8F) were compared with 22 control samples (median age: 34.5, age range: 24-43, gender ratio: 14M/8F). In metabolomics, 20 PKU serum samples (median age: 33, age range: 21–43, gender ratio: 12M/8F) were compared with 20 control samples (median age: 34.5, age range: 24–43, gender ratio: 12M/8F). In total, 15 metabolomics samples could not be matched with a suitable control; these were measured but were not included in the statistical analysis.

Lipidomics LC-MS Workflow and Data Processing

Lipid extraction. Lipid extraction was performed with the one-phase MMC extraction method, as previously described by Gil et al. [29], with slight modifications. In short, a volume of 40 μL of an internal standard mix containing 12 non-physiologic or stable isotope-labeled lipid standards was added to an Eppendorf tube and vacuum dried at 40 °C. Afterward, 10 μL of EDTA-anticoagulated plasma and 200 μL of MMC extraction solution containing MeOH/MTBE/CHCl3 (1.33:1:1, v/v/v) was added. The samples were vortexed (~10 s) and incubated on an Eppendorf Thermomixer (Eppendorf, Hamburg, Germany) at a rate of 800 rpm for 1 h at room temperature. Next, the samples were vortexed again and centrifuged at a speed of 1000 g for 10 min at 4 °C. Subsequently, the supernatant was collected and transferred into a clean Eppendorf tube with a spin filter. The extraction procedure was repeated with 40 μL MMC extraction solution, and the supernatant was collected and combined with the first fraction. The samples were centrifugated for 5 min at 1020 RCF. All samples were vacuum dried at 40 °C. In the final step, the lipids were dissolved in a final volume of 100 μL composed of 25 μL CHCl3/MeOH/MQ (60:30:4.5, v/v/v) and 75 μL ACN/IPA/MQ (2:1:1, v/v/v). All samples were transferred into a vial with an insert shortly before LC-MS analysis.
Lipidomic analysis. Plasma total triglyceride (TTG) and cholesterol (TC) concentrations were analyzed using the Roche Triglycerides (08058687190) and Cholesterol Gen.2 (08057443190) kits on a Cobas Pro analyzer (Roche Diagnostics, Risch-Rotkreuz, Switzerland). In 2021, untargeted lipidomic analysis was performed using ultra-high-performance liquid chromatography-electrospray ionization mass spectrometry (UPLC-ESI MS) on a Synapt G2-Si high-resolution quadrupole time-of-flight (QTOF) mass spectrometer equipped with a Jetstream ESI source (Waters, Manchester, UK), as described by Gil et al. [29], with some minor adjustments. Liquid chromatography was performed on an Acquity UPLC CSH C18 column (2.1 × 100 mm 1.7 μm) (Waters, Manchester, UK) column, using a binary solvent system with mobile phase A (10 mM AmF in Milli-Q) and B (10 mM AmF in MeOH). The column was pre-heated to 80 °C, and the autosampler was set to 15 °C. Lipid extracts (2 μL) were injected and separated with a flow rate of 0.5 mL/min. Linear gradient elution proceeded as follows: 0 min 40% eluent A, 7.5–15 min from 40 to 10% eluent A, followed by an isocratic elution of 0% eluent A and 100% eluent B for 17.6 min.
The QTOF settings were as follows: the capillary voltage was 0.50 kV for ESI(+) and 0.70 kV for ESI(−); the cone voltage of the samples was set at 30 V; sampling cone 20 °C, the desolvation gas temperature was set at 600 °C, the desolvation gas flow was set at 1000 L/h, source offset 20, source temperature 120 °C, cone gas 100 L/h and nebulizer gas 6.5 bar. Accurate mass determination was ensured using 0.1 ng mL−1 leucine-enkephalin solution ([M + H]+ = 556.2771 and [M − H] = 554.2615) as the lock mass solution. The instrument was operated in data-independent acquisition mode. Fragmentation was achieved with argon gas at a low energy of 1 kV and a high energy of 30–60 kV. Data were collected over the m/z range from 50 to 2000 Da in high-resolution mode for ESI(+) and resolution mode for ESI(−) with an acquisition rate of 1 spectrum/0.3 s. Detailed information about UPLC and MS parameters is provided in Table S2.
Batch design and quality control. To reduce and monitor analytical variation and bias, each batch included (i) randomization of the sequence order and (ii) stabilization of the UHPLC-QTOF before the actual sequence of samples was measured by injecting 10 times the quality control (QC) sample before each batch, (iii) measuring a QC sample at regular intervals (every 10–15 samples) in both series to evaluate the stability and repeatability of the system, (iv) injection of blank samples with and without IS, and (v) the use and analysis of internal non-physiological and/or isotopically labeled standards in the QC samples (see also Figure S1). All samples were measured in singlicate.
Metabolomic analysis. Serum samples from PKU patients and controls were analyzed by the UHPLC-ESI-QTOF-MS method previously described in detail for heparin-anticoagulated plasma, with some minor adjustments [21]. All serum samples were measured in duplicate and antiparallel order. The measurements were conducted in 2021 for ESI(+) and in 2023 for ESI(−), respectively.
Data processing. Both lipidomics and metabolomics data are available upon reasonable request but only in anonymized form. Both datasets were processed in the same manner. A schematic description of the data processing workflow is shown in Figure 1. Waters or Agilent raw data files were analyzed using Progenesis QI software Version 2.4 (Nonlinear Dynamics, Newcastle, UK) for alignment, peak picking, and deconvolution. Data were aligned on a reference run that was automatically selected from all QC measurements before peak picking. The run used as a reference file was most representative of all QC runs within or between batch(es). Alignment was considered to be successful when the score was above 95%. Runs with a score below the threshold were excluded. Afterward, the data were subjected to automatic peak picking. The following list of adducts was used for peak picking in ESI(+) mode: [M + H]+, [M + NH4]+, [M + Na]+, [M + K]+, [M + 2Na − H]+, [M+H − H2O]+, [M + H-2H2O]+, [M + CH3OH + H]+, [2M + H]+, [2M + NH4]+, [2M + Na]+, [2M + K]+, [M + 2H]2+, [M + H + Na]2+, [M + 2Na]2+, [M + 2H + Na]3+, [M + 2Na + H]3+. The following adducts were used in ESI(−) mode: [M − H], [M − H2O − H], [M + Na − 2H], [M + K − 2H], [M + FA − H], [2M − H], and [2M + FA − H].
Peak picking limits were set at ‘default’ sensitivity mode. Features were excluded if the predominant adduct detected had a charge (z) of 2 or 3 (i) and/or when the isotopic pattern was absent (Figure S2) (ii). Subsequently, the detected compound ions were deconvoluted to generate the final list with features. Raw feature abundances were exported from Progenesis QI. Data processing was performed with MetaboAnalystR [30]. Features with 80% missing values were excluded and missing values were imputed with 1/5 of the minimum positive value of each variable. Data were filtered on the pooled quality control samples that were included in each run (coefficient of variation (CV) < 30%).
Multivariate analyses. Multivariate analyses (MVA) were performed with SIMCA P version 17.0 (Umetrics AB, Umea, Sweden) by using principal component analysis (PCA) and orthogonal projections to latent structures (OPLS) analysis (PKU versus Control). Metabolomics data were normalized on L-phenyl-D5-alanine before MVA. Lipidomics data were normalized on the lipid standard that was most representative of the feature (e.g., LPC(18:0) was normalized on the standard LPC(17:1)); this was possible only after annotation following MVA. If there was not a suitable standard, the average intensity of the standards detected in the sample was taken. Before analysis, the lipidomics and metabolomics datasets were log-transformed and Pareto-scaled. For each MVA model, the quality was controlled by computing Hotelling’s T2 statistics. Model performance is reported as cumulative correlation coefficients for the model (R2), predictive performance based on K-fold cross-validation (Q2), as well as cross-validated ANOVA (CV-ANOVA) p-values for OPLS-based group separation. Feature selection was performed using variable importance in projection (VIP) scores, applying a cutoff for the VIP of greater than 0.5.
Compound reviewing. To gain confidence in the validity of the discriminating compounds from MVA, each feature was extensively reviewed. From the list of discriminating features, we aimed to determine the number of unique features. The number of features detected does not necessarily reflect the number of unique features present in the dataset [31]. Therefore, we aimed to remove redundant features before proceeding to feature annotation. Detailed information about compound reviewing is reported in Figures S3–S6. Reasons to reject features for further analysis were: (A) Deconvolution failure: compound ions are grouped by the algorithm but do not belong to the same feature; (B) Noise: features that were considered noise based on manual inspection of the peak shape of the raw extracted ion chromatograms; (C) Chromatographic resolution failure: features that were integrated with other compounds that were not separated chromatographically or not recognized by the algorithm as different features; (D) In-source fragmentation (ISF): features that were suspected of being a result of ISF. This type of feature was recognized based on the identical abundance pattern, peak shape, and retention time.
Univariate analysis. Univariate statistical analysis was employed to evaluate whether the age difference, total triglyceride content (TTG), total cholesterol (TC), and normalized feature abundances were significant. For this, the Wilcoxon rank sum test (p < 0.05; false discovery rate (FDR) adjusted, Benjamini–Hochberg) was applied. We used the Benjamini–Hochberg procedure to control the false discovery rate, i.e., the proportion of “discoveries” (significant results) that are actually false positives [32]. Chi-squared test was conducted to evaluate whether the sex difference was significant (p < 0.05) between matched PKU and control samples.
Feature annotation and nomenclature. Features were putatively annotated with mass information present in public databases. In both lipidomics and metabolomics data, the accurate mass could not deviate more than |5| parts per million (ppm) from the theoretical mass. The lipidomics features were annotated based on information from the LipidMaps structure database and ChemSpider databases [33,34]. The mass spectra were checked for the presence of qualitative fragments that were previously described in 2015 by T’Kindt et al. [35] using a comparable MS method (Tables S3 and S4). All lipids described in this article were annotated according to the nomenclature proposed by the International Lipid Classification and Nomenclature Committee [36,37]. The majority of the lipid features were annotated at the subclass level, e.g., PG(34:1). Only lysophospholipids annotations were made at the fatty acyl level, e.g., PG(16:0_0:0).
For metabolomic features, the Human Metabolome Database (HMDB) was used in combination with an in-house IEM panel that consisted of 340 IEM-related metabolites (version 0.8) [38,39,40,41,42]. Both lipidomic and metabolomic feature annotations were systematically classified by the guidelines pointed out by the Metabolomics Standards Initiative [43]. Features were annotated across four distinct levels: identified metabolites or lipids (level 1; confirmed by an authentic chemical standard analyzed under the same conditions, matching at least two orthogonal criteria here, accurate mass and retention time), putatively annotated features (level 2), putatively characterized feature classes (level 3), and features of unknown identity (level 4). In the last case, the identity was reported as the m/z value, or if possible, the neutral mass calculated by the software after deconvolution of the compound ions.
Correlation with concurrent dried blood spot Phe levels. Spearman’s rank-order correlation analysis was performed to explore the relationship between concurrent dried blood spot Phe levels and the statistically different features. All results are given as two-tailed p-values, and p < 0.05 was considered statistically significant.
Pathway enrichment analysis. Over-representation analysis (ORA) was performed using the pathway analysis module from the web-based platform MetaboAnalyst 6.0. Annotated metabolites with a VIP > 0.5 were included. Metabolites whose identification was uncertain (e.g., due to isomerism) or that lacked an available KEGG compound identification number were excluded from the analysis. The analysis was facilitated using the hypergeometric test. Human pathways were collected from the KEGG (KEGG pathway information was obtained in December 2023). To assess the significance and control for false discoveries in the pathway analysis, both the raw p-value and the false discovery rate (FDR) p-value were considered. The established thresholds for statistical significance were set at raw p-value < 0.05 and FDR p-value < 0.1.

3. Results

3.1. Population Characteristics

For all PKU patients, compliance to diet was controlled by regular monitoring of blood Phe levels. Therefore, for the majority of the patients, concurrent Phe levels were available. All patients from this study were diagnosed early by neonatal screening and treated continuously at clinical centers in The Netherlands. Treatment consisted of diet, sapropterin therapy, or both. All characteristics are displayed in Table 1 and Table S5. None of the selected controls was diagnosed with HPA or PKU, and Phe levels were, therefore, considered to be normal (local reference range for adults: 35–85 µmol/L). For the (sub)sets used for MVA of the metabolomics and lipidomics datasets, the age and sex of PKU and controls did not differ significantly.

3.2. Data Processing

For lipidomic analysis, all samples were measured in a single batch in both ESI(+) and ESI(−) modes. Two samples for which alignment failed were excluded from the analysis. Processing of the raw lipidomics data in Progenesis QI resulted in the detection of 15,220 and 5223 features in ESI(+) and ESI(−) ion mode, respectively. After the exclusion of features without an isotopic pattern, with a charge of 2+ or 3+, and with a high (>30%) CV in the pooled sample (Table 2), 30% (4438 features, ESI(+)) and 84% (4416 features, ESI(−)) of the features met our criteria. Metabolomics was performed in two separate batches, ESI(+) and ESI(−) modes, resulting in four datasets in total. The first batch contained 20 control samples and 20 PKU samples, the second batch included also the 20 control samples and the remaining 15 PKU samples. Processing of the raw metabolomics data in Progenesis QI resulted in the detection of 3627 and 2258 features in ESI(+) and ESI(−) ion modes, respectively. Both batches were processed in a single Progenesis QI experiment; as a result, the outcomes of some processing steps (e.g., peak picking) are consistent for batches A and B. After applying the primary exclusion criteria (Table 2), 20% (725 features, ESI(+)) and 22% (503 features, ESI(−)) of the features passed the criteria in both batches and were subjected to further statistical analyses.

3.3. Lipidomics

3.3.1. Multivariate Analyses of Plasma Lipidomic Data

Total TG content was statistically significantly higher in PKU patients compared to controls (p = 0.04), while there was no difference in plasma total cholesterol concentrations. Unsupervised PCA was applied to the lipidomic datasets as an exploratory data analysis tool to gain insight into quality, general trends, and intrinsic clustering of samples (Figures S7 and S8). The PCA score plot revealed clustering according to the experimental group, indicating different lipid phenotypes for PKU and controls (Figure 2).
To assess which features are most discriminating between both groups, supervised MVA was employed. For this, an OPLS-DA model was constructed including the features that remained after filtering. For both ESI(−) and ESI(+), significant models were obtained comparing the plasma lipid profiles between PKU and control groups. The OPLS-DA score plots illustrate a clear discrimination between the groups (Figure 3A,B). Because the plasma lipidome is significantly different in treated PKU, we further examined which features were most discriminating.
The VIP scores—a metric summarizing the importance of each variable in driving the observed group separation—were used for variable selection. The number of discriminating features (VIP > 0.5) was 1819 in ESI(+) mode and 1816 in ESI(−) mode. A total of 272 unique differential features were found between both groups, among which 196 were elevated and 76 were lower in the plasma of patients with PKU. The top 50 of the most discriminating features/lipids according to VIP score in ESI(+) and ESI(−)are displayed in Figure 3C,D.

3.3.2. Annotation and Univariate Significance of Discriminating Lipid Features

For annotation of lipidomics data, the queries of accurate mass values in compound databases provided several matches with a mass error below 5 ppm. In total, 205 out of 274 features (ESI(+) 154; ESI(−) 51) were annotated, and the identity of 69 features (ESI(+) 57 features; ESI(−) 12 features) remained unknown. Lipid features were annotated at level 2 due to a lack of standards. Tables S6 and S7 present a detailed summary of lipid features, including their retention times, m/z ratio, detected adducts, annotation when available, and level of identification. Examples of obtained fragmentation patterns are provided in the Supplementary Materials (Figure S9).
Alterations were predominantly noted in two key classes of lipids: glycerolipids and glycerophospholipids. Ranking the top 50 lipid species with the highest VIP scores for ESI(+), 15 positions were assigned to glycerolipid species (Figure 2). Remarkably, within this class, triglycerides (TG) were particularly discriminating (range VIP score [1.2–1.8]). The total number of carbon atoms ranged from 45 to 58. The degree of saturation varied from fully saturated up to species with 10 double bonds. Furthermore, the list included two diacylglycerol (DG) species: 34:2 [iso 3] and 40:6.
Regarding glycerophospholipids, the top 50 VIP values for ESI(+) and ESI(−) demonstrated that phosphatidylcholines (PC) and lysoPC (LPC) exhibit prominent influence within the lipidomic profile (range VIP score [0.87–1.9]). Specifically, 16 (ESI(+) and 12 (ESI(−) PCs and LPCs were in the top 50. The VIP score of multiple species exceeded 1.5, indicating a significant contribution to the predictive power of the models, including LPC: 20:3, 22:6 [iso 1], 22:6 [iso 2], P-22:6 and PC: 32:1, 36:5, 38:7 [iso 2], 42:6, O-44:4, O-44:5, O-34:2 or its isomer P-43:1. Conversely, LPC: 17:0, 18:2, 20:4, 20:5 [iso 2], 20:5 [iso 1], 22:5, LPC O-18:0, and PC: 32:2, 34:1, 36:1, 38:3, 38:7 [iso 1], 40:6, O-44:5 had a VIP score below 1.5.
Additional discriminating PL subgroups comprised lysophosphatidylethanolamines (LPE): 17:0, 17:1, 18:2, 22:0, 22:6, lysophosphatidylglycerols (LPG): 17:0, phosphatidylethanolamines (PE): 34:0, O-38:5 or P-38:4, phosphatidylinositols (PI): 34:1, 34:2, 36:1, 36:4, 39:1, 38:2, 38:3, 40:3 [iso 1], 40:3 [iso 2], O-42:6 and phosphatidylserines (PS): PS 34:0, O-38:6. Other notable lipid subgroups/species within the top 50 encompassed Ceramides (Cer): 34:1; O2 and 38:1; O2, Ceramide Phosphates (CerP): 40:1; O2, Ceramide Phosphoethanolamines (CerPE): 36:1; O2, 38:2; O2 [iso 2], Hexosylceramides (Hex2Cer): 32:0; O2, 42:2; O2, 40:1; O2, Inositol Phosphoryl Ceramides (IPC): 36:0; O2, Sphingomyelins (SM): 42:3; O2, Sulfatides (ST): 27:1; O; S Fatty Acids (FA): 16:0, 17:0, 18:3, 22:5, 22:6 and Cholesteryl Esters (CE): 18:2.
Univariate statistical testing of features with a VIP above 0.5 revealed that 33 lipid features differed statistically (adjusted p-value < 0.05) (Tables S8 and S9). Results are displayed in Figure 4. In total, 30 features were found to be higher, and 3 features were lower in PKU patients. More precisely, the following lipids were significantly higher in PKU patients: LPC 20:3, PI 36:4, 38:3, 40:3, PC: 32:0, 32:1, 34:1, 36:1, 36:5, 38:3, 38:6, 38:7 [iso 2], 40:6, 42:6, O-44:4, CerPE 38:2; O2 [iso 2], Hex2Cer 32:0; O2 and TG 45:5, 50:1, 50:2, 56:7, 58:7. In contrast, the following lipid species were significantly lower: CE 18:2, LPC 17:0, O-34:2/P-34:1. Furthermore, seven features without annotation were significantly higher in PKU patients compared to healthy controls. In ESI(+), the neutral masses (adducts ≥ 2) of the features were 375.2526 u, 384.1917 u, 360.1796 u, and 646.4526 u, and for ESI(−): 833.6073 u and m/z 752.5923. Boxplots of all significantly altered lipid features are shown in Figures S10 and S11.
Because of the relatively high proportion of females in the lipidomics datasets, we conducted the non-parametric Wilcoxon rank sum test within the PKU population (n = 21) to determine if sex is related to the potential biomarkers. The results show that there was no statistically significant difference in the intensity of lipidomic features between the male and female groups.

3.4. Metabolomics

3.4.1. Multivariate Analyses of Serum Metabolomics Data

The datasets were processed according to the procedure applied to the lipidomics data. After normalization, the majority of the quality control samples clustered together, indicating the repeatability and stability of the analytical platform. It should be noted that in ESI(+), the batch effect was not fully corrected by normalization on the IS. The serum metabolome showed a separation trend in the PCA scatter plots (Figures S12 and S13), denoting a difference in metabolic profiles. The pattern identified in the score plot was further analyzed with OPLSDA to understand the relationships and differences between samples, facilitating the identification of discriminating features (Figure S14). Among them, 160 features (ESI(+) 51 features, ESI(−) 67 features) were found to be reliable for statistical interpretation based on the secondary exclusion criteria.

3.4.2. Annotation and Univariate Significance of Discriminating Metabolite Features

For metabolomics, 60 out of 118 features were successfully annotated (ESI(+) 29 features and ESI(−) 31 features) (Tables S10 and S11). The identities of the remaining features are yet to be determined, comprising 22 features in ESI(+) and 35 features in ESI(−). Univariate testing revealed that 56 metabolite features (ESI(+) 13 features, ESI(−) 43 features) were significantly different (adjusted p-value < 0.05) compared to controls (Tables S12 and S13). Results are presented in the volcano plot (Figure 5). In total, 48 features were found to be lower (FC [0.08–0.75]), and 8 features (FC [1.15–1267]) were higher in PKU patients compared to controls. Significantly altered features were visualized using boxplots (Figures S15 and S16).
Phe and Phe catabolites were significantly higher in PKU patients, including phenyl pyruvic acid, p-hydroxyphenyl acetic acid, N-lactoyl-Phe, phenyllactic acid, glutamyl-Phe (Glu-Phe), hexanoyl-Phe (Hex-Phe), leucyl-Phe (Leu-Phe), Phe, glycyl-Phe (Gly-Phe), and phenylacetylglutamine. Furthermore, the abundance of tryptophan (Trp) derivatives indole-3-lactic acid, indole-3-carboxaldehyde, and N1-methyl-2-pyridone-5-carboxamide were also significantly higher (Figure 6).
The level of three dipeptides composed of large neutral amino acids was significantly higher, including isoleucyl-valine/leucyl-valine, tryptophan-histidine, and isoleucyl-isoleucine/leucyl-isoleucine/leucyl-leucine. Additionally, the dipeptides alanylglycine (Ala-Gly) and prolylhydroxyproline (Pro-Hyp) were significantly higher. Also, pantothenic acid, better known as vitamin B5, was higher. One annotated metabolite was significantly lower: ornithine. Metabolomics data partly confirmed alterations found with lipidomics. More specifically, there was a significantly higher level of FA 18:3; LPE 20:0 [iso 1], 20:0 [iso 2], 22:1; and LPC 20:3. In contrast, LPE 22:6 and 20:4 were significantly lower
Furthermore, there was a notable significant difference in 26 metabolite features without annotation. In more detail, one feature in ESI(+): m/z 186.1124, and three features in ESI(−): m/z 602.3458, 199.0061, and 209.0818 were lower. Likewise, a higher abundance was observed for a single feature in ESI(+): m/z 282.1195 and 21 features in ESI(−) m/z 144.0452, 143.1075, 201.0299, 339.2332, 174.0592, 151.0389, 150.0018, 329.1505, 163.0762, 147.0448, 712.2200, 91.0551, 511.2909, 385.1688, 369.1740, 383.1533, 604.3615, 187.1336, 539.2494, 165.9789, and a single feature for which a neutral mass (≥2 adducts present) was determined: 204.1360 u. With this method and the established criteria (mass error < 5 ppm), no or multiple (>3) matches were found for the identity of these features.
To determine if sex is a determinant of the potential biomarkers, we conducted the non-parametric Wilcoxon rank sum test within the PKU population (n = 35). Only the abundance of the feature P_Ile-Val/Leu-Val in the female group was statistically significantly higher than the male group (p-adjusted < 0.01); no other significant differences were found.

3.5. Correlation between Concurrent Dried Blood Spot Phe and Potential Biomarkers

In the PKU patients, concurrent dried blood spot Phe levels previously analyzed during the COBESO study (μmol/L) showed a strong significant correlation with Phe abundance of 0.88 (p  <  0.01). Additionally, several features demonstrated a strong significant correlation (Spearman ρ 0.50–1.00, p < 0.01) with concurrent Phe levels, including Glu-Phe, phenylpyruvic acid, 3.25_147.0448 m/z, phenyllactic acid, hydroxyphenylacetic acid, 4.49_151.0398 m/z, 3.25_329.1505 m/z, 7.93_91.0551 m/z, indolelactic acid, 7.40_144.0452 m/z, indole-3-carboxaldehyde, N-lactoyl Phe, phenylacetylglutamine, Ala-Gly, and 3.26_712.2200 m/z. Moderate significant correlations (0.30–0.49, p < 0.05) were observed with 10.71_163.0762 m/z, Trp-His, 1.94_282.1195 m/z, 4.68_174.0592 m/z, 12.05_143.1075 m/z, and 13.20_369.1740 m/z. No significant correlation was found with any of the other biomarkers. The detailed results of these correlations are presented in Spearman rank order in Table S14.

3.6. Metabolic Pathway Enrichment Analysis

Annotated metabolomic features were subjected to pathway analysis. Results of the analysis are shown in Figure 7 and Table 3. Results demonstrate that the identified metabolites important for PKU were mainly representative (raw p-value> 0.05 and FDR p-value < 0.1) for the following metabolism pathways: phenylalanine metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, and nicotinate and nicotinamide metabolism.

4. Discussion

Early-diagnosed and treated patients suffering from PKU can exhibit cognitive and/or behavioral abnormalities despite intensive treatment monitoring of Phe levels. Although blood Phe concentrations correlate with neurocognitive and behavioral outcomes, blood Phe only partly explains these outcomes [49]. This reinforces the hypothesis that other metabolites beyond Phe are involved in the pathogenesis of PKU. Identification of novel metabolites that better explain the pathophysiology underlying PKU disease could provide markers that are of diagnostic, prognostic, and/or therapeutic value. In this study, we conducted untargeted lipidomic and metabolomic analysis using high-resolution mass spectrometry platforms in a population of treated PKU patients and non-PKU controls. The PKU patients were part of the COBESO study, whose neurocognitive functioning was assessed. In this work, we were able to show that the plasma lipid and serum metabolite signature of PKU patients is different from healthy controls.
In the analysis of the metabolomics data, Phe, Phe catabolites (phenylpyruvate, phenyllactate, phenylacetate, hydroxyphenylacetic acid), Phe-conjugates (Phe-Hex, Lac-Phe, Phenylacetylglutamine) and Phe-containing dipeptides (Glu-Phe, Leu-Phe, Gly-Phe) were significantly elevated in the PKU population. The results are consistent with our previously conducted targeted analysis of this metabolomics dataset using a different bioinformatics workflow, which searched specifically for 11 previously described PKU biomarkers, including Phe, N-lactoyl-Phe, Glu-Phe, Phe-Hexose, Glu-Glu-Phe, Phe-Phe, Phe-Leu, N-acetyl-Phe, phenyllactate, phenylacetate, and phenylpyruvate [49]. Phenylacetylglutamine and Gly-Phe were not included previously, and also have not been described as altered in plasma before. However, a higher level of phenylacetylglutamine has been reported as a urinary marker for non-treated PKU patients [50,51,52].
Most Phe and Phe derivatives show a significant strong positive correlation with concurrent Phe levels. These components may seem less interesting as new biomarkers because the increases in the Phe-containing metabolites are inherent to elevated Phe levels in the blood and probably the intake of balanced amino acid supplements. However, Wegberg et al. [49] demonstrated that N-lactoyl-Phe outperformed Phe as a predictor of working memory and mental health outcomes. These findings suggest that N-lactoyl-Phe may be a better biomarker to reflect neurocognitive outcomes than Phe. Additionally, phenylacetylglutamine is produced from phenylacetic acid and Phe, which may competitively inhibit L-amino acid decarboxylase, thereby decreasing neurotransmitter levels [53]. Regarding Phe-containing dipeptides, limited information exists about their general physiological significance or their clinical and pathophysiological relevance for PKU.
The level of pantothenic acid, also referred to as vitamin B5, was reported to be higher in the PKU population. This micronutrient is a precursor of coenzyme-A (CoA) and the acyl carrier protein. It is involved in cellular energy production and the synthesis of carbohydrates, proteins, and fats. Further, it is involved in the biosynthesis of essential lipids, steroids, hormones, neurotransmitters, and porphyrin, and the formation of red blood cells, sex, and stress-related hormones. Stroup et al. [54] reported that without supplementation, PKU patients are at risk of inadequate intake of pantothenic acid. Recently, Bokayeva et al. [55] pointed out the lack of publications describing vitamin status, including vitamin B5, in PKU patients and the effect on clinical outcomes. We speculate that the higher level of vitamin B5 in our study is most likely a result of its presence in supplements.
Interestingly, multiple Trp derivatives were found to be significantly higher in the serum of PKU patients, including indole-3-lactic acid (I3LA), indole-3-carboxaldehyde (I3A), and N1-methyl-2-pyridone-5-carboxamide (2PY). Moreover, N1-methyl-4-pyridone-3-carboxamide (4PY) and indole-3-propionic acid (IPA) levels were higher in PKU patients and almost reached statistical significance (adjusted p-value = 0.05). Trp is an essential amino acid, meaning that it is obtained from the diet or protein catabolism but not synthesized endogenously. Trp itself is a precursor for various biologically active metabolites, including the neurotransmitter serotonin. The gut mucosa’s enterochromaffin cells (ECs) are primarily responsible for the body’s serotonin synthesis. Serotonin is released by activated endothelial cells (ECs) into the interstitium of the lamina propria, where it binds intrinsic sensory neurons’ nerve endings [53]. Being part of the microbiota–gut–brain axis, Trp and its metabolites could be of interest for clinical outcomes in PKU.
Trp is metabolized along four pathways, including serotonin (serotonin and melatonin), tryptamine (tryptamine), indole (indolepyruvic acid), and kynurenine (niacin) pathways. The majority (~90%) of Trp is metabolized through the kynurenine pathway. The metabolites that were found to be elevated in our PKU patients are downstream catabolites of the kynurenine (2PY, 4PY) and indole (I3LA, I3A, IPA) pathways (see also Figure 6). Higher levels of 2PY and 4PY have been described before in urine samples of early-treated PKU patients [56]. The reason for this relation and the clinical relevance is uncertain. Both compounds are considered to be uremic toxins based on the fact that their concentrations increase with decreased kidney function and based on their ability to inhibit poly(ADP-ribose) polymerase-1 (PARP1). This enzyme is involved in cellular response to DNA damage [57,58,59,60]. Both compounds are, however, end products of Trp catabolism (kynurenine pathway) via the degradation of nicotinamide-adenine dinucleotide (NAD). In theory, higher levels might merely reflect an increased Trp and/or niacin status or a rerouting of Trp metabolism in the kynurenine pathway. The niacin status of PKU patients is not well described. The implications of this finding need further research.
Several studies suggest changes in the gut bacterial microbes of PKU patients who have undergone treatment compared to individuals without PKU. In more detail, PKU patients exhibited a higher number of Trp-consuming bacteria in the colon [61,62]. Alterations in bacteria consuming Trp might affect serotonin synthesis in the brain since it modulates the Trp availability in the circulation. Only I3A has been detected before in the urine of patients with untreated PKU. While more research is needed, these findings support those described by Parolisi et al. [63,64], highlighting the relevance of exploring the potential effects of Trp availability through gut microbiota on cognitive and behavioral functions in PKU.
We observed that PKU patients exhibit a different lipid signature compared to healthy individuals. Differences were predominantly observed between lipids belonging to the glycerophospholipid and glycerolipid classes. This observation confirms that closer evaluation of the lipidome in PKU patients might be interesting, as previously emphasized by Guerra et al. [65,66,67,68]. We identified a significant increase in various glycerophospholipids, including (L)PE, (L)PC, and PI. Most species incorporated at least one fatty acid with an unsaturated aliphatic chain, as indicated by the total number of carbons and double bonds. None of the features were correlated with concurrent Phe, which makes them potential biomarkers of interest for PKU, reflecting either pathophysiological processes or consequences of the dietary treatment.
Although we were unable to elucidate the full molecular structure of the fatty acyl chain(s), previous studies and data provided herein provide a level of insight. To date, alterations in plasma-free FA status, particularly in n-3 FA were observed in PKU patients treated with a Phe-restricted diet [69,70,71,72,73]. The n-3 long-chain PUFA species DHA (22:6 n-3) and eicosapentaenoic acid (EPA, 20:5n-3) were consistently demonstrated to be reduced [70,72,74,75]. These findings were attributed to dietary restrictions. PKU patients are advised to consume no or only a limited amount of animal products because of the high protein content in these products. As a result, foods rich in saturated fatty acids and long-chain PUFAs are excluded from their diet [74,76]. Some protein substitutes are enriched with DHA, but not all of them. Long-chain PUFAs, particularly arachidonic acid (AA, 20:4n-6) and DHA, are crucial for normal neurological development [77,78,79,80]. Based on these findings, it is recommended to supplement EPA and DHA if these are not already added to the protein substitute [17]. Unfortunately, dietary data were not available for this cohort. Studies show that suppletion increases the free and PL plasma levels of these FAs in PKU [81,82,83,84,85]. Although there is a debate whether the difference is clinically relevant it is known to improve fine motor skills and short-term visual response in children with PKU [69,70,74,81,82,83,84,85,86,87].
Furthermore, the FA composition of PL-enriched plasma extracts has been reported in several studies. For instance, Drzymała-Czyż et al. [88] reported a reduction in FA: 20:3n-9, 18:2n-6 (LA, linoleic acid), and 22:6n-3 (DHA). In 2021, Guerra et al. [89] investigated the plasma PL and sphingolipid profile in children treated for PKU. In short, they found a significant increase in FA 18:1n-9 (oleic acid), 14:0 (myristic acid), 18:0 (stearic acid), EPA, and DHA, while no significant changes in statistical difference were observed for FA 16:0, LA, and AA. Consistent with these findings, our detailed analysis of PL species confirmed a significant increase in PC species (32:0, 32:1, 34:1, 36:5, 38:7, 40:6, and 42:6) and PI species (36:4 and 38:3). Altogether, changes in FA status due to the PKU diet may explain the observed alterations in the PL profile of PKU patients. Additional measurements to fully elucidate the structure of the FA chains are essential to give further meaning to these findings.
The total triglyceride (TTG) concentration was significantly elevated in PKU patients. This only captures a fraction of the total lipid content and provides limited detailed information about the level of individual species. To the best of our knowledge, this is the first study that described detailed information about TG alterations. In addition, we observed a significant increase in various TG (45:5, 50:1, 50:2, 56:7, and 58:7) containing saturated and unsaturated fatty acyl chains. To date, some perturbations in lipoprotein components of PKU individuals have been described in plasma and serum. However, the results are contradictory. To demonstrate, four studies report a reduction in TC, LDL-C, and/or HDL-C [90,91,92,93]. However, Azabdaftari et al. [94] reported an increase in TC content in serum. Similarly, VLDL-C has been found to be higher [90,95] but also reduced [92,94]. Studies in both plasma and serum show that total TG levels increase in children and older PKU patients [90,92,95,96,97]. This phenomenon was associated with increased carbohydrate intake.
The implications of higher TTG levels for clinical outcomes are not fully understood. However, it is known to correlate with the accumulation of visceral fat, which, in turn, is linked to a higher risk of developing type 2 diabetes and cardiovascular diseases [98,99]. This might also be relevant for PKU patients, as there are reports of obesity, dyslipidemia, hypertension, and oxidative stress in PKU patients [91,100,101,102,103]. Furthermore, in the present study, we found an increase in CerPE 38:2. CerPEs are produced in small amounts and located in the plasma membrane. To date, their function remains largely unknown [104]. However, higher levels of CerPE 38:2 have been suggested to increase the risk of type 2 diabetes [105]. Additionally, we report lower levels of LPC 17:0, which is also known to be inversely correlated with type 2 diabetes [106,107]. Higher TTG and CerPE 38:2 abundance and lower levels of LPC 17:0 might suggest that this cohort of PKU patients is at increased risk of CVD and type 2 diabetes. However, the question remains whether this risk is comparable to that of non-PKU individuals. Unfortunately, no anthropometric data were available for our study cohort. The link between CVD or diabetes and the plasma lipid composition may be of interest to future research.
Strengths of this study include a comprehensive metabolomic and lipidomic analysis of plasma/serum obtained in a well-defined population of PKU patients. Patients who participated in the COBESO study are well described with respect to their neuropsychological outcomes [24,25,26,27]. We are, therefore, able to correlate the biomarkers to clinically relevant outcomes. We chose to comprehensively describe the analytical methods and important biomarkers for PKU in this study, whereas we aim to describe the relation with neuropsychological outcomes in a separate paper, allowing both topics to receive the attention and depth they deserve. Another advantage of our combined untargeted-omics approach is that we found several new PKU biomarkers; in particular, the lipidome is rather unexplored. Further studies are necessary to explore the relevance of our findings.
A limitation of our study is that the control and patient samples were obtained from different studies, potentially introducing variability. We would like to emphasize that this study is exploratory and acknowledge that variables such as long-term sample storage may have affected the metabolome or lipid signature. It is well known that long-term storage significantly impacts the metabolome and the lipidome [108,109], including amino acids such as Phe [110,111], although the impact of storage of EDTA-plasma on the metabolome seems limited when samples are stored up to 7 years at −80 degrees Celcius [112]. Unfortunately, no long-term data on the stability of the PKU biomarkers found in our study are available in our laboratories. The impact of the difference in storage time between PKU patients and the non-PKU controls on the biomarkers is, therefore, uncertain. On the other hand, some of the biomarkers have been found by other groups as well, demonstrating the validity of these biomarkers. Another important aspect is that our results were obtained in a relatively small population. In combination with the other potential confounding variables, we realize that any conclusions based on this work are to be considered preliminary and should be further explored.
In the future, we intend to conduct correlation analyses between the biomarkers identified in this study and the previously obtained neurocognitive outcomes from the COBESO study. If these analyses reveal any interesting results, we will proceed to confirm the identity of the detected metabolomics and lipidomics features of interest. We intend to confirm the identity, preferably using standards, in addition to performing a DDA (data-dependent acquisition) experiment, allowing us to generate clean spectra and making feature annotation more reliable. Naturally, this study must be repeated with more recent data to verify the findings, preferably in a larger prospective study applying a quantitative analytical method, avoiding long-term storage effects. Additionally, future studies should incorporate anthropometrical and, most importantly, nutritional metrics to contextualize the metabolomic and lipidomic findings to diet.
To summarize, in this study, we have shown that the lipidomic and metabolic profiles of treated adult PKU patients differ from non-PKU individuals. Regarding lipidomics, total TG content was higher and alterations were found for lipid features belonging to the classes of glycerolipids, glycerophospholipids, and sphingolipids. Metabolomic results confirmed alterations in the lipidome. Furthermore, metabolomic results demonstrated that patients had higher levels of Phe and Phe catabolites, including N-lactoyl-Phe, Trp derivatives, pantothenic acid, and various dipeptides. Future studies should focus on the identification of unknown features and correlate significantly different lipid and metabolite features with outcome parameters. Altogether, this research provides insight into the metabolomic and lipid status of treated adult PKU patients. However, the clinical and/or therapeutic relevance and validity of the biomarkers need to be determined.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/metabo14090479/s1: Table S1: Overview lipid standards; Table S2: UPLC and MS settings in ESI(+) and ESI(−); Table S3: Qualitative fragments of lipids in negative ion mode; Table S4: Qualitative fragments of lipids in positive ion mode; Table S5: Information of the samples included in the lipidomics and metabolomics study; Table S6: Lipidomic features in ESI(+); Table S7: Lipidomic features in ESI(−); Table S8: Results from the Wilcoxon Rank Sum Test comparing control and PKU groups within the ESI(+) lipidomics dataset; Table S9: Results from the Wilcoxon Rank Sum Test comparing control and PKU groups within the ESI(−) lipidomics dataset; Table S10: Annotation metabolic features in ESI(+); Table S11: Annotation metabolic features in ESI(−); Table S12: Results from the Wilcoxon Rank Sum Test comparing control and PKU groups within the ESI(−) metabolomics dataset; Table S13: Results from the Wilcoxon Rank Sum Test comparing control and PKU groups within the ESI(+) metabolomics dataset; Table S14: Results Spearman correlations; Figure S1: Schematic representation of run order for lipidomics measurements; Figure S2: Example absence of isotopic pattern; Figure S3: Example of deconvolution failure; Figure S4: Example of chromatographic resolution failure; Figure S5: Examples of noise; Figure S6: Example of in-source fragmentation; Figure S7: Two-dimensional score plot of the PCA model based on the raw data; Figure S8: Time series plot of PC1 against the sample injection order in positive and negative ionization mode; Figure S9: Representative ion ESI-MS spectra for a selection of lipid classes detected; Figure S10: Statistically significant lipid features in ESI(+), at p < 0.05; Figure S11: Statistically significant lipid features in ESI(−), at p < 0.05; Figure S12: Metabolomics normalization and filtering effects on the data obtained in ESI(+); Figure S13: Metabolomics normalization and filtering effects on the data obtained in ESI(−); Figure S14: OPLS-DA results of metabolomic profiles in PKU patients compared to controls; Figure S15: Statistically significant metabolite features in ESI(+), at p < 0.05; Figure S16: Statistically significant metabolite features in ESI(−), at p < 0.05.

Author Contributions

Conceptualization, J.C.v.d.W., F.J.v.S. and M.R.H.-F.; Methodology, J.C.v.d.W., T.S.B., A.M.J.v.W. and U.F.H.E.; Software, J.C.v.d.W., P.d.B. and U.F.H.E.; Validation J.C.v.d.W., T.S.B., P.d.B. and U.F.H.E.; Formal Analysis, J.C.v.d.W. and A.M.J.v.W.; Investigation, J.C.v.d.W., A.M.J.v.W., U.F.H.E. and M.R.H.-F.; Resources, A.M.J.v.W., R.A.W., S.J.L.B., M.R.H.-F. and F.J.v.S.; Data Curation, J.C.v.d.W. and U.F.H.E.; Writing—Original Draft Preparation, J.C.v.d.W. and M.R.H.-F.; Writing—Review and Editing, all authors; Visualization, J.C.v.d.W. and A.M.J.v.W.; Supervision, F.J.v.S. and M.R.H.-F.; Project Administration, F.J.v.S. and M.R.H.-F.; Funding Acquisition, F.J.v.S. and M.R.H.-F. All authors have read and agreed to the published version of the manuscript.

Funding

A grant from the Dutch Association of Inherited Metabolic Diseases (ESN, 22 November 2018) and a grant from Nutricia (2019N0174) were received.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Medical Ethical Committee of the University Medical Center Groningen (METc2011/377 and METc2014/077, respectively).

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy and ethical restrictions.

Acknowledgments

We thank Ellis van Steenis for her assistance in collecting the blood sampling information on the samples that were part of the COBESO study. Figure 6 and the graphical abstract were drawn in part using images from Servier Medical Art. Servier Medical Art by Servier is licensed under a Creative Commons Attribution 3.0 Unported License (https://creativecommons.org/licenses/by/3.0/, accessed on 20 August 2024).

Conflicts of Interest

The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results. The UMCG, as an employer of A.M.J.v.W., received a research grant from Nutricia, honoraria as a speaker, and travel grants from Nutricia and Vitaflo. The UMCG as employer of F.J.v.S. has received research grants, advisory board fees, and/or speaker’s honoraria from Agios, Alexion, AlltRNA, Applied Pharma Research, Arla Food Int., Biomarin, Beatrix Research Fund, Codexis, Eurocept, ESPKU, Homology, Lucane Nestle-Codexis Alliance, Moderna, Nutricia, NPKUA, NPKUV, Orphan Europe, Pluvia Biotech, Rivium Medical BV, Vitaflo, Vivet, and ZONMW. J.C.v.d.W., T.S.B., P.d.B., J.G., U.F.H.E., K.L.M.C., S.J.L.B., R.A.W. and M.R.H.-F. declare no conflicts of interest.

References

  1. Hillert, A.; Anikster, Y.; Belanger-Quintana, A.; Burlina, A.; Burton, B.K.; Carducci, C.; Chiesa, A.E.; Christodoulou, J.; Đorđević, M.; Desviat, L.R.; et al. The Genetic Landscape and Epidemiology of Phenylketonuria. Am. J. Hum. Genet. 2020, 107, 234–250. [Google Scholar] [CrossRef]
  2. Blau, N.; van Spronsen, F.J.; Levy, H.L. Phenylketonuria. Lancet 2010, 376, 1417–1427. [Google Scholar] [CrossRef] [PubMed]
  3. Anikster, Y.; Haack, T.B.; Vilboux, T.; Pode-Shakked, B.; Thöny, B.; Shen, N.; Guarani, V.; Meissner, T.; Mayatepek, E.; Trefz, F.K.; et al. Biallelic Mutations in DNAJC12 Cause Hyperphenylalaninemia, Dystonia, and Intellectual Disability. Am. J. Hum. Genet. 2017, 100, 257–266. [Google Scholar] [CrossRef] [PubMed]
  4. Cleary, M.A.; Walter, J.H.; Wraith, J.E.; Jenkins, J.P.; Alani, S.M.; Tyler, K.; Whittle, D. Magnetic resonance imaging of the brain in phenylketonuria. Lancet 1994, 344, 87–90. [Google Scholar] [CrossRef]
  5. Fölling, A. Über Ausscheidung von Phenylbrenztraubensäure in den Harn als Stoffwechselanomalie in Verbindung mit Imbezillität. Biol. Chem. 1934, 227, 169–181. [Google Scholar] [CrossRef]
  6. Følling, I. The discovery of phenylketonuria. Acta Paediatr. Suppl. 1994, 407, 4–10. [Google Scholar] [CrossRef]
  7. Leuzzi, V.; Trasimeni, G.; Gualdi, G.F.; Antonozzi, I. Biochemical, clinical and neuroradiological (MRI) correlations in late-detected PKU patients. J. Inherit. Metab. Dis. 1995, 18, 624–634. [Google Scholar] [CrossRef]
  8. Penrose, L.; Quastel, J.H. Metabolic studies in phenylketonuria. Biochem. J. 1937, 31, 266–274. [Google Scholar] [CrossRef] [PubMed]
  9. Pietz, J.; Kreis, R.; Rupp, A.; Mayatepek, E.; Rating, D.; Boesch, C.; Bremer, H.J. Large neutral amino acids block phenylalanine transport into brain tissue in patients with phenylketonuria. J. Clin. Investig. 1999, 103, 1169–1178. [Google Scholar] [CrossRef]
  10. Smith, I.; Beasley, M.G.; Ades, A.E. Intelligence and quality of dietary treatment in phenylketonuria. Arch. Dis. Child. 1990, 65, 472–478. [Google Scholar] [CrossRef]
  11. Thompson, A.J.; Tillotson, S.; Smith, I.; Kendall, B.; Moore, S.G.; Brenton, D.P. Brain MRI changes in phenylketonuria. Associations with dietary status. Brain 1993, 116 Pt 4, 811–821. [Google Scholar] [CrossRef] [PubMed]
  12. VanZutphen, K.; Packman, W.; Sporri, L.; Needham, M.; Morgan, C.; Weisiger, K.; Packman, S. Executive functioning in children and adolescents with phenylketonuria. Clin. Genet. 2007, 72, 13–18. [Google Scholar] [CrossRef]
  13. van Spronsen, F.J.; van Wegberg, A.M.; Ahring, K.; Bélanger-Quintana, A.; Blau, N.; Bosch, A.M.; Burlina, A.; Campistol, J.; Feillet, F.; Giżewska, M.; et al. Key European guidelines for the diagnosis and management of patients with phenylketonuria. Lancet Diabetes Endocrinol. 2017, 5, 743–756. [Google Scholar] [CrossRef] [PubMed]
  14. Giovannini, M.; Verduci, E.; Salvatici, E.; Paci, S.; Riva, E. Phenylketonuria: Nutritional advances and challenges. Nutr. Metab. 2012, 9, 7. [Google Scholar] [CrossRef] [PubMed]
  15. van Wegberg, A.M.J.; MacDonald, A.; Ahring, K.; Bélanger-Quintana, A.; Blau, N.; Bosch, A.M.; Burlina, A.; Campistol, J.; Feillet, F.; Giżewska, M.; et al. The complete European guidelines on phenylketonuria: Diagnosis and treatment. Orphanet J. Rare Dis. 2017, 12, 162. [Google Scholar] [CrossRef]
  16. Burlina, A.; Biasucci, G.; Carbone, M.T.; Cazzorla, C.; Paci, S.; Pochiero, F.; Spada, M.; Tummolo, A.; Zuvadelli, J.; Leuzzi, V. Italian national consensus statement on management and pharmacological treatment of phenylketonuria. Orphanet J. Rare Dis. 2021, 16, 476. [Google Scholar] [CrossRef]
  17. Harding, C.O.; Longo, N.; Northrup, H.; Sacharow, S.; Singh, R.; Thomas, J.A.; Vockley, J.; Zori, R.T.; Bulloch Whitehall, K.; Lilienstein, J.; et al. Pegvaliase for the treatment of phenylketonuria: Final results of a long-term phase 3 clinical trial program. Mol. Genet. Metab. Rep. 2024, 39, 101084. [Google Scholar] [CrossRef]
  18. Krämer, J.; Baerwald, C.; Heimbold, C.; Kamrath, C.; Parhofer, K.G.; Reichert, A.; Rutsch, F.; Stolz, S.; Weinhold, N.; Muntau, A.C. Two years of pegvaliase in Germany: Experiences and best practice recommendations. Mol. Genet. Metab. 2023, 139, 107564. [Google Scholar] [CrossRef]
  19. Dijkstra, A.M.; van Vliet, N.; van Vliet, D.; Romani, C.; Huijbregts, S.C.J.; van der Goot, E.; Hovens, I.B.; van der Zee, E.A.; Kema, I.P.; Heiner-Fokkema, M.R.; et al. Correlations of blood and brain biochemistry in phenylketonuria: Results from the Pah-enu2 PKU mouse. Mol. Genet. Metab. 2021, 134, 250–256. [Google Scholar] [CrossRef]
  20. Miller, M.J.; Kennedy, A.D.; Eckhart, A.D.; Burrage, L.C.; Wulff, J.E.; Miller, L.A.; Milburn, M.V.; Ryals, J.A.; Beaudet, A.L.; Sun, Q.; et al. Untargeted metabolomic analysis for the clinical screening of inborn errors of metabolism. J. Inherit. Metab. Dis. 2015, 38, 1029–1039. [Google Scholar] [CrossRef]
  21. Coene, K.L.M.; Kluijtmans, L.A.J.; van der Heeft, E.; Engelke, U.F.H.; de Boer, S.; Hoegen, B.; Kwast, H.J.T.; van de Vorst, M.; Huigen, M.; Keularts, I.; et al. Next-generation metabolic screening: Targeted and untargeted metabolomics for the diagnosis of inborn errors of metabolism in individual patients. J. Inherit. Metab. Dis. 2018, 41, 337–353. [Google Scholar] [CrossRef] [PubMed]
  22. van Outersterp, R.E.; Moons, S.J.; Engelke, U.F.H.; Bentlage, H.; Peters, T.M.A.; van Rooij, A.; Huigen, M.; de Boer, S.; van der Heeft, E.; Kluijtmans, L.A.J.; et al. Amadori rearrangement products as potential biomarkers for inborn errors of amino-acid metabolism. Commun. Biol. 2021, 4, 367. [Google Scholar] [CrossRef]
  23. Guerra, I.M.S.; Ferreira, H.B.; Neves, B.; Melo, T.; Diogo, L.M.; Domingues, M.R.; Moreira, A.S.P. Lipids and phenylketonuria: Current evidences pointed the need for lipidomics studies. Arch. Biochem. Biophys. 2020, 688, 108431. [Google Scholar] [CrossRef] [PubMed]
  24. Jahja, R.; van Spronsen, F.J.; de Sonneville, L.M.J.; van der Meere, J.J.; Bosch, A.M.; Hollak, C.E.M.; Rubio-Gozalbo, M.E.; Brouwers, M.; Hofstede, F.C.; de Vries, M.C.; et al. Long-Term Follow-Up of Cognition and Mental Health in Adult Phenylketonuria: A PKU-COBESO Study. Behav. Genet. 2017, 47, 486–497. [Google Scholar] [CrossRef] [PubMed]
  25. Jahja, R.; van Spronsen, F.J.; de Sonneville, L.M.J.; van der Meere, J.J.; Bosch, A.M.; Hollak, C.E.M.; Rubio-Gozalbo, M.E.; Brouwers, M.; Hofstede, F.C.; de Vries, M.C.; et al. Social-cognitive functioning and social skills in patients with early treated phenylketonuria: A PKU-COBESO study. J. Inherit. Metab. Dis. 2016, 39, 355–362. [Google Scholar] [CrossRef]
  26. Jahja, R.; Huijbregts, S.C.J.; de Sonneville, L.M.J.; van der Meere, J.J.; Legemaat, A.M.; Bosch, A.M.; Hollak, C.E.M.; Rubio-Gozalbo, M.E.; Brouwers, M.; Hofstede, F.C.; et al. Cognitive profile and mental health in adult phenylketonuria: A PKU-COBESO study. Neuropsychology 2017, 31, 437–447. [Google Scholar] [CrossRef]
  27. Jahja, R.; Huijbregts, S.C.; de Sonneville, L.M.; van der Meere, J.J.; Bosch, A.M.; Hollak, C.E.; Rubio-Gozalbo, M.E.; Brouwers, M.C.; Hofstede, F.C.; de Vries, M.C.; et al. Mental health and social functioning in early treated Phenylketonuria: The PKU-COBESO study. Mol. Genet. Metab. 2013, 110, S57–S61. [Google Scholar] [CrossRef]
  28. Eisenga, M.F.; Gomes-Neto, A.W.; van Londen, M.; Ziengs, A.L.; Douwes, R.M.; Stam, S.P.; Osté, M.C.J.; Knobbe, T.J.; Hessels, N.R.; Buunk, A.M.; et al. Rationale and design of TransplantLines: A prospective cohort study and biobank of solid organ transplant recipients. BMJ Open 2018, 8, e024502. [Google Scholar] [CrossRef] [PubMed]
  29. Gil, A.; Zhang, W.; Wolters, J.C.; Permentier, H.; Boer, T.; Horvatovich, P.; Heiner-Fokkema, M.R.; Reijngoud, D.J.; Bischoff, R. One- vs. two-phase extraction: Re-evaluation of sample preparation procedures for untargeted lipidomics in plasma samples. Anal. Bioanal. Chem. 2018, 410, 5859–5870. [Google Scholar] [CrossRef]
  30. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef]
  31. Forest, A.; Ruiz, M.; Bouchard, B.; Boucher, G.; Gingras, O.; Daneault, C.; Robillard Frayne, I.; Rhainds, D.; Tardif, J.C.; Rioux, J.D.; et al. Comprehensive and Reproducible Untargeted Lipidomic Workflow Using LC-QTOF Validated for Human Plasma Analysis. J. Proteome Res. 2018, 17, 3657–3670. [Google Scholar] [CrossRef]
  32. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  33. Royal Society of Chemistry ChemSpider|Search and Share Chemistry. Available online: https://www.chemspider.com/ (accessed on 16 April 2024).
  34. Sud, M.; Fahy, E.; Cotter, D.; Brown, A.; Dennis, E.A.; Glass, C.K.; Merrill, A.H., Jr.; Murphy, R.C.; Raetz, C.R.; Russell, D.W.; et al. LMSD: LIPID MAPS structure database. Nucleic Acids Res. 2007, 35, D527–D532. [Google Scholar] [CrossRef] [PubMed]
  35. t’Kindt, R.; Telenga, E.D.; Jorge, L.; Van Oosterhout, A.J.; Sandra, P.; Ten Hacken, N.H.; Sandra, K. Profiling over 1500 lipids in induced lung sputum and the implications in studying lung diseases. Anal. Chem. 2015, 87, 4957–4964. [Google Scholar] [CrossRef] [PubMed]
  36. Fahy, E.; Subramaniam, S.; Brown, H.A.; Glass, C.K.; Merrill, A.H., Jr.; Murphy, R.C.; Raetz, C.R.; Russell, D.W.; Seyama, Y.; Shaw, W.; et al. A comprehensive classification system for lipids. J. Lipid Res. 2005, 46, 839–861. [Google Scholar] [CrossRef]
  37. Fahy, E.; Subramaniam, S.; Murphy, R.C.; Nishijima, M.; Raetz, C.R.; Shimizu, T.; Spener, F.; van Meer, G.; Wakelam, M.J.; Dennis, E.A. Update of the LIPID MAPS comprehensive classification system for lipids. J. Lipid Res. 2009, 50, S9–S14. [Google Scholar] [CrossRef]
  38. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef]
  39. Wishart, D.S.; Guo, A.; Oler, E.; Wang, F.; Anjum, A.; Peters, H.; Dizon, R.; Sayeeda, Z.; Tian, S.; Lee, B.L.; et al. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res. 2022, 50, D622–D631. [Google Scholar] [CrossRef]
  40. Wishart, D.S.; Jewison, T.; Guo, A.C.; Wilson, M.; Knox, C.; Liu, Y.; Djoumbou, Y.; Mandal, R.; Aziat, F.; Dong, E.; et al. HMDB 3.0—The Human Metabolome Database in 2013. Nucleic Acids Res. 2013, 41, D801–D807. [Google Scholar] [CrossRef]
  41. Wishart, D.S.; Knox, C.; Guo, A.C.; Eisner, R.; Young, N.; Gautam, B.; Hau, D.D.; Psychogios, N.; Dong, E.; Bouatra, S.; et al. HMDB: A knowledgebase for the human metabolome. Nucleic Acids Res. 2009, 37, D603–D610. [Google Scholar] [CrossRef]
  42. Wishart, D.S.; Tzur, D.; Knox, C.; Eisner, R.; Guo, A.C.; Young, N.; Cheng, D.; Jewell, K.; Arndt, D.; Sawhney, S.; et al. HMDB: The Human Metabolome Database. Nucleic Acids Res. 2007, 35, D521–D526. [Google Scholar] [CrossRef] [PubMed]
  43. Sumner, L.W.; Amberg, A.; Barrett, D.; Beale, M.H.; Beger, R.; Daykin, C.A.; Fan, T.W.; Fiehn, O.; Goodacre, R.; Griffin, J.L.; et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 2007, 3, 211–221. [Google Scholar] [CrossRef]
  44. Gao, K.; Mu, C.L.; Farzi, A.; Zhu, W.Y. Tryptophan Metabolism: A Link between the Gut Microbiota and Brain. Adv. Nutr. 2020, 11, 709–723. [Google Scholar] [CrossRef]
  45. Roth, W.; Zadeh, K.; Vekariya, R.; Ge, Y.; Mohamadzadeh, M. Tryptophan Metabolism and Gut-Brain Homeostasis. Int. J. Mol. Sci. 2021, 22, 2973. [Google Scholar] [CrossRef] [PubMed]
  46. Kałużna-Czaplińska, J.; Gątarek, P.; Chirumbolo, S.; Chartrand, M.S.; Bjørklund, G. How important is tryptophan in human health? Crit. Rev. Food Sci. Nutr. 2019, 59, 72–88. [Google Scholar] [CrossRef]
  47. O’Mahony, S.M.; Clarke, G.; Borre, Y.E.; Dinan, T.G.; Cryan, J.F. Serotonin, tryptophan metabolism and the brain-gut-microbiome axis. Behav. Brain Res. 2015, 277, 32–48. [Google Scholar] [CrossRef] [PubMed]
  48. Walther, D.J.; Bader, M. A unique central tryptophan hydroxylase isoform. Biochem. Pharmacol. 2003, 66, 1673–1680. [Google Scholar] [CrossRef]
  49. van Wegberg, A.M.J.; van der Weerd, J.C.; Engelke, U.F.H.; Coene, K.L.M.; Jahja, R.; Bakker, S.J.L.; Huijbregts, S.C.J.; Wevers, R.A.; Heiner-Fokkema, M.R.; van Spronsen, F.J. The clinical relevance of novel biomarkers as outcome parameter in adults with phenylketonuria. J. Inherit. Metab. Dis. 2024, 47, 524–635. [Google Scholar] [CrossRef]
  50. Andrade, F.; Cano, A.; Unceta Suarez, M.; Arza, A.; Vinuesa, A.; Ceberio, L.; López-Oslé, N.; de Frutos, G.; López-Oceja, R.; Aznal, E.; et al. Urine Phenylacetylglutamine Determination in Patients with Hyperphenylalaninemia. J. Clin. Med. 2021, 10, 3674. [Google Scholar] [CrossRef]
  51. Wild, J.; Shanmuganathan, M.; Hayashi, M.; Potter, M.; Britz-McKibbin, P. Metabolomics for improved treatment monitoring of phenylketonuria: Urinary biomarkers for non-invasive assessment of dietary adherence and nutritional deficiencies. Analyst 2019, 144, 6595–6608. [Google Scholar] [CrossRef]
  52. Xiong, X.; Sheng, X.; Liu, D.; Zeng, T.; Peng, Y.; Wang, Y. A GC/MS-based metabolomic approach for reliable diagnosis of phenylketonuria. Anal. Bioanal. Chem. 2015, 407, 8825–8833. [Google Scholar] [CrossRef]
  53. Berguig, G.Y.; Martin, N.T.; Creer, A.Y.; Xie, L.; Zhang, L.; Murphy, R.; Pacheco, G.; Bullens, S.; Olbertz, J.; Weng, H.H. Of mice and men: Plasma phenylalanine reduction in PKU corrects neurotransmitter pathways in the brain. Mol. Genet. Metab. 2019, 128, 422–430. [Google Scholar] [CrossRef] [PubMed]
  54. Stroup, B.M.; Ney, D.M.; Murali, S.G.; Rohr, F.; Gleason, S.T.; van Calcar, S.C.; Levy, H.L. Metabolomic Insights into the Nutritional Status of Adults and Adolescents with Phenylketonuria Consuming a Low-Phenylalanine Diet in Combination with Amino Acid and Glycomacropeptide Medical Foods. J. Nutr. Metab. 2017, 2017, 6859820. [Google Scholar] [CrossRef] [PubMed]
  55. Bokayeva, K.; Jamka, M.; Walkowiak, D.; Duś-Żuchowska, M.; Herzig, K.H.; Walkowiak, J. Vitamin Status in Patients with Phenylketonuria: A Systematic Review and Meta-Analysis. Int. J. Mol. Sci. 2024, 25, 5065. [Google Scholar] [CrossRef] [PubMed]
  56. Cannet, C.; Bayat, A.; Frauendienst-Egger, G.; Freisinger, P.; Spraul, M.; Himmelreich, N.; Kockaya, M.; Ahring, K.; Godejohann, M.; MacDonald, A.; et al. Phenylketonuria (PKU) Urinary Metabolomic Phenotype Is Defined by Genotype and Metabolite Imbalance: Results in 51 Early Treated Patients Using Ex Vivo (1)H-NMR Analysis. Molecules 2023, 28, 4916. [Google Scholar] [CrossRef]
  57. Rutkowski, B.; Rutkowski, P.; Słomińska, E.; Smolenski, R.T.; Swierczyński, J. Cellular toxicity of nicotinamide metabolites. J. Ren. Nutr. 2012, 22, 95–97. [Google Scholar] [CrossRef]
  58. Rutkowski, B.; Slominska, E.; Szolkiewicz, M.; Smolenski, R.T.; Striley, C.; Rutkowski, P.; Swierczynski, J. N-methyl-2-pyridone-5-carboxamide: A novel uremic toxin? Kidney Int. 2003, 63, S19–S21. [Google Scholar] [CrossRef]
  59. Slominska, E.M.; Kowalik, K.; Smolenski, R.T.; Szolkiewicz, M.; Rutkowski, P.; Rutkowski, B.; Swierczynski, J. Accumulation of poly(ADP-ribose) polymerase inhibitors in children with chronic renal failure. Pediatr. Nephrol. 2006, 21, 800–806. [Google Scholar] [CrossRef] [PubMed]
  60. Synesiou, E.; Fairbanks, L.D.; Simmonds, H.A.; Slominska, E.M.; Smolenski, R.T.; Carrey, E.A. 4-Pyridone-3-carboxamide-1-β-D-ribonucleoside triphosphate (4PyTP), a novel NAD metabolite accumulating in erythrocytes of uremic children: A biomarker for a toxic NAD analogue in other tissues? Toxins 2011, 3, 520–537. [Google Scholar] [CrossRef]
  61. Mancilla, V.J.; Mann, A.E.; Zhang, Y.; Allen, M.S. The Adult Phenylketonuria (PKU) Gut Microbiome. Microorganisms 2021, 9, 530. [Google Scholar] [CrossRef]
  62. Verduci, E.; Carbone, M.T.; Borghi, E.; Ottaviano, E.; Burlina, A.; Biasucci, G. Nutrition, Microbiota and Role of Gut-Brain Axis in Subjects with Phenylketonuria (PKU): A Review. Nutrients 2020, 12, 3319. [Google Scholar] [CrossRef]
  63. Chyan, Y.J.; Poeggeler, B.; Omar, R.A.; Chain, D.G.; Frangione, B.; Ghiso, J.; Pappolla, M.A. Potent neuroprotective properties against the Alzheimer beta-amyloid by an endogenous melatonin-related indole structure, indole-3-propionic acid. J. Biol. Chem. 1999, 274, 21937–21942. [Google Scholar] [CrossRef] [PubMed]
  64. Poeggeler, B.; Sambamurti, K.; Siedlak, S.L.; Perry, G.; Smith, M.A.; Pappolla, M.A. A novel endogenous indole protects rodent mitochondria and extends rotifer lifespan. PLoS ONE 2010, 5, e10206. [Google Scholar] [CrossRef] [PubMed]
  65. Bercik, P.; Verdu, E.F.; Foster, J.A.; Macri, J.; Potter, M.; Huang, X.; Malinowski, P.; Jackson, W.; Blennerhassett, P.; Neufeld, K.A.; et al. Chronic gastrointestinal inflammation induces anxiety-like behavior and alters central nervous system biochemistry in mice. Gastroenterology 2010, 139, 2102–2112.e1. [Google Scholar] [CrossRef] [PubMed]
  66. Braniste, V.; Al-Asmakh, M.; Kowal, C.; Anuar, F.; Abbaspour, A.; Tóth, M.; Korecka, A.; Bakocevic, N.; Ng, L.G.; Kundu, P.; et al. The gut microbiota influences blood-brain barrier permeability in mice. Sci. Transl. Med. 2014, 6, 263ra158. [Google Scholar] [CrossRef] [PubMed]
  67. Dickerson, F.; Severance, E.; Yolken, R. The microbiome, immunity, and schizophrenia and bipolar disorder. Brain Behav. Immun. 2017, 62, 46–52. [Google Scholar] [CrossRef]
  68. Xu, M.; Xu, X.; Li, J.; Li, F. Association between Gut Microbiota and Autism Spectrum Disorder: A Systematic Review and Meta-Analysis. Front. Psychiatry 2019, 10, 473. [Google Scholar] [CrossRef]
  69. Agostoni, C.; Marangoni, F.; Riva, E.; Giovannini, M.; Galli, C. Plasma arachidonic acid and serum thromboxane B2 concentrations in phenylketonuric children negatively correlate with dietary compliance. Prostaglandins Leukot. Essent. Fat. Acids 1997, 56, 219–222. [Google Scholar] [CrossRef]
  70. Galli, C.; Agostoni, C.; Mosconi, C.; Riva, E.; Salari, P.C.; Giovannini, M. Reduced plasma C-20 and C-22 polyunsaturated fatty acids in children with phenylketonuria during dietary intervention. J. Pediatr. 1991, 119, 562–567. [Google Scholar] [CrossRef]
  71. Giovannini, M.; Biasucci, G.; Agostoni, C.; Luotti, D.; Riva, E. Lipid status and fatty acid metabolism in phenylketonuria. J. Inherit. Metab. Dis. 1995, 18, 265–272. [Google Scholar] [CrossRef]
  72. Moseley, K.; Koch, R.; Moser, A.B. Lipid status and long-chain polyunsaturated fatty acid concentrations in adults and adolescents with phenylketonuria on phenylalanine-restricted diet. J. Inherit. Metab. Dis. 2002, 25, 56–64. [Google Scholar] [CrossRef] [PubMed]
  73. van Gool, C.J.; van Houwelingen, A.C.; Hornstra, G. The essential fatty acid status in phenylketonuria patients under treatment. J. Nutr. Biochem. 2000, 11, 543–547. [Google Scholar] [CrossRef]
  74. Pöge, A.P.; Bäumann, K.; Müller, E.; Leichsenring, M.; Schmidt, H.; Bremer, H.J. Long-chain polyunsaturated fatty acids in plasma and erythrocyte membrane lipids of children with phenylketonuria after controlled linoleic acid intake. J. Inherit. Metab. Dis. 1998, 21, 373–381. [Google Scholar] [CrossRef] [PubMed]
  75. Sanjurjo, P.; Perteagudo, L.; Rodríguez Soriano, J.; Vilaseca, A.; Campistol, J. Polyunsaturated fatty acid status in patients with phenylketonuria. J. Inherit. Metab. Dis. 1994, 17, 704–709. [Google Scholar] [CrossRef]
  76. Koletzko, B.; Sauerwald, T.; Demmelmair, H.; Herzog, M.; von Schenck, U.; Böhles, H.; Wendel, U.; Seidel, J. Dietary long-chain polyunsaturated fatty acid supplementation in infants with phenylketonuria: A randomized controlled trial. J. Inherit. Metab. Dis. 2007, 30, 326–332. [Google Scholar] [CrossRef] [PubMed]
  77. Campoy, C.; Escolano-Margarit, M.V.; Anjos, T.; Szajewska, H.; Uauy, R. Omega 3 fatty acids on child growth, visual acuity and neurodevelopment. Br. J. Nutr. 2012, 107 (Suppl. S2), S85–S106. [Google Scholar] [CrossRef]
  78. Decsi, T.; Koletzko, B. Role of Long-Chain Polyunsaturated Fatty Acids in Early Human Neurodevelopment. Nutr. Neurosci. 2000, 3, 293–306. [Google Scholar] [CrossRef] [PubMed]
  79. Koletzko, B.; Lien, E.; Agostoni, C.; Böhles, H.; Campoy, C.; Cetin, I.; Decsi, T.; Dudenhausen, J.W.; Dupont, C.; Forsyth, S.; et al. The roles of long-chain polyunsaturated fatty acids in pregnancy, lactation and infancy: Review of current knowledge and consensus recommendations. J. Perinat. Med. 2008, 36, 5–14. [Google Scholar] [CrossRef]
  80. Bazan, N.G. Omega-3 fatty acids, pro-inflammatory signaling and neuroprotection. Curr. Opin. Clin. Nutr. Metab. Care 2007, 10, 136–141. [Google Scholar] [CrossRef]
  81. Agostoni, C.; Massetto, N.; Biasucci, G.; Rottoli, A.; Bonvissuto, M.; Bruzzese, M.G.; Giovannini, M.; Riva, E. Effects of long-chain polyunsaturated fatty acid supplementation on fatty acid status and visual function in treated children with hyperphenylalaninemia. J. Pediatr. 2000, 137, 504–509. [Google Scholar] [CrossRef]
  82. Agostoni, C.; Riva, E.; Biasucci, G.; Luotti, D.; Bruzzese, M.G.; Marangoni, F.; Giovannini, M. The effects of n-3 and n-6 polyunsaturated fatty acids on plasma lipids and fatty acids of treated phenylketonuric children. Prostaglandins Leukot. Essent. Fat. Acids 1995, 53, 401–404. [Google Scholar] [CrossRef] [PubMed]
  83. Beblo, S.; Reinhardt, H.; Demmelmair, H.; Muntau, A.C.; Koletzko, B. Effect of fish oil supplementation on fatty acid status, coordination, and fine motor skills in children with phenylketonuria. J. Pediatr. 2007, 150, 479–484. [Google Scholar] [CrossRef]
  84. Beblo, S.; Reinhardt, H.; Muntau, A.C.; Mueller-Felber, W.; Roscher, A.A.; Koletzko, B. Fish oil supplementation improves visual evoked potentials in children with phenylketonuria. Neurology 2001, 57, 1488–1491. [Google Scholar] [CrossRef] [PubMed]
  85. Koletzko, B.; Beblo, S.; Demmelmair, H.; Müller-Felber, W.; Hanebutt, F.L. Does dietary DHA improve neural function in children? Observations in phenylketonuria. Prostaglandins Leukot. Essent. Fat. Acids 2009, 81, 159–164. [Google Scholar] [CrossRef]
  86. Agostoni, C.; Scaglioni, S.; Bonvissuto, M.; Bruzzese, M.G.; Giovannini, M.; Riva, E. Biochemical effects of supplemented long-chain polyunsaturated fatty acids in hyperphenylalaninemia. Prostaglandins Leukot. Essent. Fat. Acids 2001, 64, 111–115. [Google Scholar] [CrossRef] [PubMed]
  87. Agostoni, C.; Verduci, E.; Massetto, N.; Fiori, L.; Radaelli, G.; Riva, E.; Giovannini, M. Long term effects of long chain polyunsaturated fats in hyperphenylalaninemic children. Arch. Dis. Child. 2003, 88, 582–583. [Google Scholar] [CrossRef]
  88. Drzymała-Czyż, S.; Kałużny, Ł.; Krzyżanowska-Jankowska, P.; Walkowiak, D.; Mozrzymas, R.; Walkowiak, J. Deficiency of long-chain polyunsaturated fatty acids in phenylketonuria: A cross-sectional study. Acta Biochim. Pol. 2018, 65, 303–308. [Google Scholar] [CrossRef]
  89. Guerra, I.M.S.; Diogo, L.; Pinho, M.; Melo, T.; Domingues, P.; Domingues, M.R.; Moreira, A.S.P. Plasma Phospholipidomic Profile Differs between Children with Phenylketonuria and Healthy Children. J. Proteome Res. 2021, 20, 2651–2661. [Google Scholar] [CrossRef]
  90. Colomé, C.; Artuch, R.; Lambruschini, N.; Cambra, F.J.; Campistol, J.; Vilaseca, M. Is there a relationship between plasma phenylalanine and cholesterol in phenylketonuric patients under dietary treatment? Clin. Biochem. 2001, 34, 373–376. [Google Scholar] [CrossRef]
  91. Couce, M.L.; Vitoria, I.; Aldámiz-Echevarría, L.; Fernández-Marmiesse, A.; Roca, I.; Llarena, M.; Sánchez-Pintos, P.; Leis, R.; Hermida, A. Lipid profile status and other related factors in patients with Hyperphenylalaninaemia. Orphanet J. Rare Dis. 2016, 11, 123. [Google Scholar] [CrossRef]
  92. Schulpi, K.H.; Scarpalezou, A. Triglycerides, cholesterol, HDL, LDL, and VLDL cholesterol in serum of phenylketonuric children under dietary control. Clin. Pediatr. 1989, 28, 466–469. [Google Scholar] [CrossRef] [PubMed]
  93. Schulpis, K.H.; Karakonstantakis, T.; Bartzeliotou, A.; Karikas, G.A.; Papassotiriou, I. The association of serum lipids, lipoproteins and apolipoproteins with selected trace elements and minerals in phenylketonuric patients on diet. Clin. Nutr. 2004, 23, 401–407. [Google Scholar] [CrossRef] [PubMed]
  94. Azabdaftari, A.; van der Giet, M.; Schuchardt, M.; Hennermann, J.B.; Plöckinger, U.; Querfeld, U. The cardiovascular phenotype of adult patients with phenylketonuria. Orphanet J. Rare Dis. 2019, 14, 213. [Google Scholar] [CrossRef]
  95. Lavoie, S.M.; Harding, C.O.; Gillingham, M.B. Normal fatty acid concentrations in young children with phenylketonuria (PKU). Top. Clin. Nutr. 2009, 24, 333–340. [Google Scholar] [CrossRef]
  96. Fernandez-Crespo, S.; Vazquez-Agra, N.; Marques-Afonso, A.T.; Cruces-Sande, A.; Martinez-Olmos, M.A.; Araujo-Vilar, D.; Hermida-Ameijeiras, A. The value of waist circumference as a preditor of cardiovascular risk in adult patients with classic phenylketonuria. Med. Clin. 2023, 161, 470–475. [Google Scholar] [CrossRef]
  97. Rocha, J.C.; van Spronsen, F.J.; Almeida, M.F.; Soares, G.; Quelhas, D.; Ramos, E.; Guimarães, J.T.; Borges, N. Dietary treatment in phenylketonuria does not lead to increased risk of obesity or metabolic syndrome. Mol. Genet. Metab. 2012, 107, 659–663. [Google Scholar] [CrossRef]
  98. Després, J.P. Body fat distribution and risk of cardiovascular disease: An update. Circulation 2012, 126, 1301–1313. [Google Scholar] [CrossRef] [PubMed]
  99. Piché, M.E.; Tchernof, A.; Després, J.P. Obesity Phenotypes, Diabetes, and Cardiovascular Diseases. Circ. Res. 2020, 126, 1477–1500. [Google Scholar] [CrossRef]
  100. Hennermann, J.B.; Roloff, S.; Gellermann, J.; Vollmer, I.; Windt, E.; Vetter, B.; Plöckinger, U.; Mönch, E.; Querfeld, U. Chronic kidney disease in adolescent and adult patients with phenylketonuria. J. Inherit. Metab. Dis. 2013, 36, 747–756. [Google Scholar] [CrossRef]
  101. Rocha, J.C.; Martins, M.J. Oxidative stress in phenylketonuria: Future directions. J. Inherit. Metab. Dis. 2012, 35, 381–398. [Google Scholar] [CrossRef]
  102. Gündüz, M.; Çakar, S.; Kuyum, P.; Makay, B.; Arslan, N. Comparison of atherogenic risk factors among poorly controlled and well-controlled adolescent phenylketonuria patients. Cardiol. Young 2016, 26, 901–908. [Google Scholar] [CrossRef] [PubMed]
  103. Robertson, L.V.; McStravick, N.; Ripley, S.; Weetch, E.; Donald, S.; Adam, S.; Micciche, A.; Boocock, S.; MacDonald, A. Body mass index in adult patients with diet-treated phenylketonuria. J. Hum. Nutr. Diet. 2013, 26 (Suppl. S1), 1–6. [Google Scholar] [CrossRef] [PubMed]
  104. Bickert, A.; Ginkel, C.; Kol, M.; vom Dorp, K.; Jastrow, H.; Degen, J.; Jacobs, R.L.; Vance, D.E.; Winterhager, E.; Jiang, X.C.; et al. Functional characterization of enzymes catalyzing ceramide phosphoethanolamine biosynthesis in mice. J. Lipid Res. 2015, 56, 821–835. [Google Scholar] [CrossRef]
  105. Fall, T.; Salihovic, S.; Brandmaier, S.; Nowak, C.; Ganna, A.; Gustafsson, S.; Broeckling, C.D.; Prenni, J.E.; Kastenmüller, G.; Peters, A.; et al. Non-targeted metabolomics combined with genetic analyses identifies bile acid synthesis and phospholipid metabolism as being associated with incident type 2 diabetes. Diabetologia 2016, 59, 2114–2124. [Google Scholar] [CrossRef]
  106. Bao, L.; Zhang, Y.; Yan, S.; Yan, D.; Jiang, D. Lysophosphatidylcholine (17:0) Improves HFD-Induced Hyperglycemia & Insulin Resistance: A Mechanistic Mice Model Study. Diabetes Metab. Syndr. Obes. 2022, 15, 3511–3517. [Google Scholar] [CrossRef] [PubMed]
  107. Long, J.; Yang, H.; Yang, Z.; Jia, Q.; Liu, L.; Kong, L.; Cui, H.; Ding, S.; Qin, Q.; Zhang, N.; et al. Integrated biomarker profiling of the metabolome associated with impaired fasting glucose and type 2 diabetes mellitus in large-scale Chinese patients. Clin. Transl. Med. 2021, 11, e432. [Google Scholar] [CrossRef]
  108. Haid, M.; Muschet, C.; Wahl, S.; Römisch-Margl, W.; Prehn, C.; Möller, G.; Adamski, J. Long-Term Stability of Human Plasma Metabolites during Storage at −80 °C. J. Proteome Res. 2017, 17, 203–211. [Google Scholar] [CrossRef]
  109. Wagner-Golbs, A.; Neuber, S.; Kamlage, B.; Christiansen, N.; Bethan, B.; Rennefahrt, U.; Schatz, P.; Lind, L. Effects of Long-Term Storage at −80 °C on the Human Plasma Metabolome. Metabolites 2019, 9, 99. [Google Scholar] [CrossRef]
  110. Buchanan, J.L.; Vaquerano, J.T.; Taylor, E.B. Isolated Effects of Plasma Freezing versus Thawing on Metabolite Stability. Metabolites 2022, 12, 1098. [Google Scholar] [CrossRef]
  111. An, Z.; Shi, C.; Li, P.; Liu, L. Stability of Amino Acids and Related Amines in Human Serum under Different Preprocessing and Pre-Storage Conditions Based on iTRAQ®-LC-MS/MS. Biol. Open 2021, 10, bio055020. [Google Scholar] [CrossRef]
  112. Heiskanen, L.A.; Suoniemi, M.; Ta, H.X.; Tarasov, K.; Ekroos, K. Long-Term Performance and Stability of Molecular Shotgun Lipidomic Analysis of Human Plasma Samples. Anal. Chem. 2013, 85, 8757–8763. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flowchart of the data processing workflow. This figure illustrates the workflow applied to the datasets in this study to identify potential biomarkers using the Progenesis QI software Version 2.4. The process is divided into several steps, including data import, retention time alignment, peak picking, feature filtering and selection, and biomarker identification. Details of the workflow are extensively described in Figures S2–S6.
Figure 1. Flowchart of the data processing workflow. This figure illustrates the workflow applied to the datasets in this study to identify potential biomarkers using the Progenesis QI software Version 2.4. The process is divided into several steps, including data import, retention time alignment, peak picking, feature filtering and selection, and biomarker identification. Details of the workflow are extensively described in Figures S2–S6.
Metabolites 14 00479 g001
Figure 2. Three-dimensional PCA score plots of lipidomic profiles from ESI(+) and ESI(−) modes for control and PKU patients. Principal Component Analysis (PCA) score plots derived from lipidomic data obtained in both ESI(+) (A) and ESI(−) (B) ionization modes. The plots illustrate the spatial distribution of controls (n = 21) and PKU patients (n = 21), represented by blue and red circles, respectively. The axes of the plots correspond to the first three principal components (PC1 on the x-axis, PC2 on the y-axis, and PC3 on the z-axis), which capture most of the variances within the dataset. Panel (A), corresponding to the ESI(+) mode, reveals an R2X value of 0.50 and a cumulative Q2 of 0.19, indicating a moderate explanation of variance and predictability from the model. Meanwhile, Panel (B), illustrating ESI(−) mode, demonstrates a higher model efficacy with an R2X of 0.66 and a cumulative Q2 of 0.40.
Figure 2. Three-dimensional PCA score plots of lipidomic profiles from ESI(+) and ESI(−) modes for control and PKU patients. Principal Component Analysis (PCA) score plots derived from lipidomic data obtained in both ESI(+) (A) and ESI(−) (B) ionization modes. The plots illustrate the spatial distribution of controls (n = 21) and PKU patients (n = 21), represented by blue and red circles, respectively. The axes of the plots correspond to the first three principal components (PC1 on the x-axis, PC2 on the y-axis, and PC3 on the z-axis), which capture most of the variances within the dataset. Panel (A), corresponding to the ESI(+) mode, reveals an R2X value of 0.50 and a cumulative Q2 of 0.19, indicating a moderate explanation of variance and predictability from the model. Meanwhile, Panel (B), illustrating ESI(−) mode, demonstrates a higher model efficacy with an R2X of 0.66 and a cumulative Q2 of 0.40.
Metabolites 14 00479 g002
Figure 3. OPLS-DA results of lipidomic profiles in PKU patients compared to controls. The figure displays OPLSDA score plots for pairwise comparisons between PKU and control groups, utilizing lipidomics data. Panels (A,B) illustrate the results for ESI(+) and ESI(−), respectively. Each point represents the lipidomic profile of an individual, with PKU patients and controls plotted to show the clear distinction and grouping based on their lipid abundances. For ESI+, the model diagnostics are as follows: R2X = 0.249, R2Y(cum) = 0.943, and Q2 = 0.42, indicating a strong model with significant predictive capability and reliability in distinguishing between PKU patients and controls. The colored circles represent 95% confidence regions, demonstrating tight clustering and clear separation between the groups. In the negative ionization mode (B), the model shows a higher complexity and fit with R2X = 0.52, R2Y(cum) = 0.986, and Q2 = 0.673. Both models were statistically significant, with p-values below 0.05, underscoring the robustness of the lipidomic differences identified. Additionally, panels (C,D) display the top 50 features with the highest VIP scores for both models in ESI(+) and ESI(−), respectively.
Figure 3. OPLS-DA results of lipidomic profiles in PKU patients compared to controls. The figure displays OPLSDA score plots for pairwise comparisons between PKU and control groups, utilizing lipidomics data. Panels (A,B) illustrate the results for ESI(+) and ESI(−), respectively. Each point represents the lipidomic profile of an individual, with PKU patients and controls plotted to show the clear distinction and grouping based on their lipid abundances. For ESI+, the model diagnostics are as follows: R2X = 0.249, R2Y(cum) = 0.943, and Q2 = 0.42, indicating a strong model with significant predictive capability and reliability in distinguishing between PKU patients and controls. The colored circles represent 95% confidence regions, demonstrating tight clustering and clear separation between the groups. In the negative ionization mode (B), the model shows a higher complexity and fit with R2X = 0.52, R2Y(cum) = 0.986, and Q2 = 0.673. Both models were statistically significant, with p-values below 0.05, underscoring the robustness of the lipidomic differences identified. Additionally, panels (C,D) display the top 50 features with the highest VIP scores for both models in ESI(+) and ESI(−), respectively.
Metabolites 14 00479 g003aMetabolites 14 00479 g003b
Figure 4. Circular bar chart expressing the fold change of annotated lipid features. The fold change is calculated as the average ratio of lipid abundance in PKU patients to that in control subjects. Bars extending outward from the circle center represent the fold change ranging from 0 to 2, indicating higher (FC > 1) or lower (FC < 1) lipid abundances in PKU patients compared to controls. Significantly higher or lower (adjusted p-value < 0.05) lipids are marked with an asterisk (*). The ‘O-’ prefix is used for plasmanyl species to indicate the presence of an alkyl ether substituent, whereas the ‘P-’ prefix is used for plasmenyl species to indicate the alk-1-enyl ether substituent.
Figure 4. Circular bar chart expressing the fold change of annotated lipid features. The fold change is calculated as the average ratio of lipid abundance in PKU patients to that in control subjects. Bars extending outward from the circle center represent the fold change ranging from 0 to 2, indicating higher (FC > 1) or lower (FC < 1) lipid abundances in PKU patients compared to controls. Significantly higher or lower (adjusted p-value < 0.05) lipids are marked with an asterisk (*). The ‘O-’ prefix is used for plasmanyl species to indicate the presence of an alkyl ether substituent, whereas the ‘P-’ prefix is used for plasmenyl species to indicate the alk-1-enyl ether substituent.
Metabolites 14 00479 g004
Figure 5. Volcano plot displaying metabolomic data. Labeled points indicate metabolite or metabolite feature is significant. Metabolites starting with ‘N_’ were detected in negative ion mode electrospray ionization (ESI−), while those with ‘P_’ were identified in positive ion mode electrospray ionization (ESI+). The x-axis represents the log2 of the fold change, and the y-axis indicates the negative log10 of the adjusted p-values, with the dashed lines marking thresholds for statistical significance.
Figure 5. Volcano plot displaying metabolomic data. Labeled points indicate metabolite or metabolite feature is significant. Metabolites starting with ‘N_’ were detected in negative ion mode electrospray ionization (ESI−), while those with ‘P_’ were identified in positive ion mode electrospray ionization (ESI+). The x-axis represents the log2 of the fold change, and the y-axis indicates the negative log10 of the adjusted p-values, with the dashed lines marking thresholds for statistical significance.
Metabolites 14 00479 g005
Figure 6. Overview of detected metabolites and their location in the Trp metabolism. Simplified representation of the metabolic pathways of Trp in the human body. Trp is an essential amino acid that is catabolized by four key metabolic routes, including the serotonin, kynurenine, indole, and tryptamine pathway. The pathways have been thoroughly reviewed previously [44,45]. In short, after food consumption, most of the Trp is absorbed in the small intestine and enters the bloodstream bound to albumin or as free Trp. The remaining 4–6% of Trp reaches the large intestine and is the primary precursor of indole, tryptamine, and other molecules derived from the microbiome [46]. The indole and tryptamine pathway involves microbes, for example, bacteria belonging to Bacteroides, Clostridium, and Bifidobacterium. The kynurenine pathway is the primary route for Trp degradation in the liver. Trp is converted to kynurenine by the enzyme tryptophan 2,3-dioxygenase (TDO) in the liver or indoleamine 2,3-dioxygenase (IDO) in other tissues [47]. Kynurenine can then be further metabolized into several bioactive compounds, including kynurenic acid and quinolinic acid. This route is also the starting point for the synthesis of nicotinamide adenine dinucleotide (NAD) and other downstream metabolites, including 2PY and 4PY. In the serotonin pathway, the enzyme tryptophan hydroxylase (TPH) is predominantly expressed in the brain (TPH2) and enterochromaffin cells in the gut (TPH1) [48]. They convert Trp to 5-hydroxytryptophan (5-HTP), and thereafter, 5-HTP is converted to serotonin (5-hydroxytryptamine, 5-HT) that enters the circulation. Significantly elevated metabolites (2PY, I3LA, I3A) or almost reached significance (IPA, 4PY) are marked bold. 3-HAA = 3-hydroxyanthranilic acid, 5-HIAA = 5-hydroxyindoleacetic acid, 3-HK = 3-hydroxykynurenine, 5-HT = 5-hydroxytryptamine, 5-HTP = 5-hydroxytryptophan, IA = 3-indoleacrylic acid, IAA = indole-3-acetic acid, IAAld = indole-3-acetaldehyde, I3A = indole-3-carboxyaldehyde, IAM = indole-3-acetamide, I3LA = indole-3-lactic acid, IPA = indole-3-propionic acid, IPYA = indole-3-pyruvic acid, NAD = nicotinamide adenine dinucleotide, XA = xanthurenic acid.
Figure 6. Overview of detected metabolites and their location in the Trp metabolism. Simplified representation of the metabolic pathways of Trp in the human body. Trp is an essential amino acid that is catabolized by four key metabolic routes, including the serotonin, kynurenine, indole, and tryptamine pathway. The pathways have been thoroughly reviewed previously [44,45]. In short, after food consumption, most of the Trp is absorbed in the small intestine and enters the bloodstream bound to albumin or as free Trp. The remaining 4–6% of Trp reaches the large intestine and is the primary precursor of indole, tryptamine, and other molecules derived from the microbiome [46]. The indole and tryptamine pathway involves microbes, for example, bacteria belonging to Bacteroides, Clostridium, and Bifidobacterium. The kynurenine pathway is the primary route for Trp degradation in the liver. Trp is converted to kynurenine by the enzyme tryptophan 2,3-dioxygenase (TDO) in the liver or indoleamine 2,3-dioxygenase (IDO) in other tissues [47]. Kynurenine can then be further metabolized into several bioactive compounds, including kynurenic acid and quinolinic acid. This route is also the starting point for the synthesis of nicotinamide adenine dinucleotide (NAD) and other downstream metabolites, including 2PY and 4PY. In the serotonin pathway, the enzyme tryptophan hydroxylase (TPH) is predominantly expressed in the brain (TPH2) and enterochromaffin cells in the gut (TPH1) [48]. They convert Trp to 5-hydroxytryptophan (5-HTP), and thereafter, 5-HTP is converted to serotonin (5-hydroxytryptamine, 5-HT) that enters the circulation. Significantly elevated metabolites (2PY, I3LA, I3A) or almost reached significance (IPA, 4PY) are marked bold. 3-HAA = 3-hydroxyanthranilic acid, 5-HIAA = 5-hydroxyindoleacetic acid, 3-HK = 3-hydroxykynurenine, 5-HT = 5-hydroxytryptamine, 5-HTP = 5-hydroxytryptophan, IA = 3-indoleacrylic acid, IAA = indole-3-acetic acid, IAAld = indole-3-acetaldehyde, I3A = indole-3-carboxyaldehyde, IAM = indole-3-acetamide, I3LA = indole-3-lactic acid, IPA = indole-3-propionic acid, IPYA = indole-3-pyruvic acid, NAD = nicotinamide adenine dinucleotide, XA = xanthurenic acid.
Metabolites 14 00479 g006
Figure 7. Results of the pathway enrichment analysis. The enrichment ratio reflects the number of hits divided by the expected number of hits. The pathways are ordered from low to high raw p-value.
Figure 7. Results of the pathway enrichment analysis. The enrichment ratio reflects the number of hits divided by the expected number of hits. The pathways are ordered from low to high raw p-value.
Metabolites 14 00479 g007
Table 1. Characteristics of the study population 1.
Table 1. Characteristics of the study population 1.
MetabolomicsLipidomics
PKUControlPKUControl
n35202222
Age (years)29 ± 734 ± 533 ± 634 ± 5
Sex
male18 (51%)12 (60%)8 (36%)8 (36%)
female17 (49%)8 (40%)14 (64%)14 (64%)
Lifetime Phe (µmol/L)439 [223–1001], n = 33-446 [255–734], n = 20-
Concurrent Phe (µmol/L)616 [176–1250], n = 35 645 [176–1250], n = 22
BH4 treatment8/35-2/22-
Plasma total cholesterol (mmol/L) 4.39 ± 0.814.65 ± 0.66
Plasma total triglyceride (mmol/L) 1.52 ± 0.551.32 ± 0.74
Matrix for the current studySerumSerumEDTA-plasmaEDTA-plasma
1 Results are expressed as mean ± standard deviation (SD) or as median [range].
Table 2. The number of features retained after each filtering step.
Table 2. The number of features retained after each filtering step.
LipidomicsMetabolomics
ESI(+)ESI(−)ESI(+)ESI(−)
Batch ABatch BBatch ABatch B
Number of features detected15,220522336272258
Primary exclusion criteria
An isotopic abundance of 10035705271291571
Predominant adduct with z = 2 or 3328501570
Features with 80% missing values132012314430422314
CV > 30%26075689231145507562
Number of features passing primary criteria44384416725503
OPLS-DA
VIP > 0.532263243628281
Secondary exclusion criteria after OPLS-DA
The final number of unique features210645167
Table 3. Detailed results from the pathway enrichment analysis.
Table 3. Detailed results from the pathway enrichment analysis.
Total aExpectedHits bRaw p cHolm p dFDR e
Phenylalanine metabolism80.04737.64 × 10−60.0006120.000612
Phenylalanine, tyrosine and tryptophan biosynthesis40.02320.0001810.01430.00726
Nicotinate and nicotinamide metabolism150.08820.003070.2390.0819
Arginine and proline metabolism360.21120.017310.346
Caffeine metabolism100.05910.057110.914
Arginine biosynthesis140.08210.079211
Pantothenate and CoA biosynthesis200.11710.11111
Glutathione metabolism280.16410.15311
Glycine, serine and threonine metabolism330.19310.17811
a Total: the number of metabolites in the pathway b Hits is the matched number of metabolites from the data c Raw p: is the original p-value calculated from the enrichment analysis d Holm p: is the p-value adjusted by the Holm–Bonferroni method e FDR: p is the p-value adjusted using false discovery rate.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Weerd, J.C.v.d.; Wegberg, A.M.J.v.; Boer, T.S.; Engelke, U.F.H.; Coene, K.L.M.; Wevers, R.A.; Bakker, S.J.L.; Blaauw, P.d.; Groen, J.; Spronsen, F.J.v.; et al. Impact of Phenylketonuria on the Serum Metabolome and Plasma Lipidome: A Study in Early-Treated Patients. Metabolites 2024, 14, 479. https://doi.org/10.3390/metabo14090479

AMA Style

Weerd JCvd, Wegberg AMJv, Boer TS, Engelke UFH, Coene KLM, Wevers RA, Bakker SJL, Blaauw Pd, Groen J, Spronsen FJv, et al. Impact of Phenylketonuria on the Serum Metabolome and Plasma Lipidome: A Study in Early-Treated Patients. Metabolites. 2024; 14(9):479. https://doi.org/10.3390/metabo14090479

Chicago/Turabian Style

Weerd, Jorine C. van der, Annemiek M. J. van Wegberg, Theo S. Boer, Udo F. H. Engelke, Karlien L. M. Coene, Ron A. Wevers, Stephan J. L. Bakker, Pim de Blaauw, Joost Groen, Francjan J. van Spronsen, and et al. 2024. "Impact of Phenylketonuria on the Serum Metabolome and Plasma Lipidome: A Study in Early-Treated Patients" Metabolites 14, no. 9: 479. https://doi.org/10.3390/metabo14090479

APA Style

Weerd, J. C. v. d., Wegberg, A. M. J. v., Boer, T. S., Engelke, U. F. H., Coene, K. L. M., Wevers, R. A., Bakker, S. J. L., Blaauw, P. d., Groen, J., Spronsen, F. J. v., & Heiner-Fokkema, M. R. (2024). Impact of Phenylketonuria on the Serum Metabolome and Plasma Lipidome: A Study in Early-Treated Patients. Metabolites, 14(9), 479. https://doi.org/10.3390/metabo14090479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop